Original Research ARTICLE
Astrocytic Coverage of Dendritic Spines, Dendritic Shafts, and Axonal Boutons in Hippocampal Neuropil
- 1UNN Institute of Neuroscience, N. I. Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod, Russia
- 2Department of Biophysics, Faculty of Biology, M. V. Lomonosov Moscow State University, Moscow, Russia
- 3Department of Biology, University at Albany, The State University of New York (SUNY), Albany, NY, United States
- 4Institute of Information Technologies, Mathematics and Mechanics, N. I. Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod, Russia
- 5Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, Moscow, Russia
- 6All-Russian Research Institute of Medicinal and Aromatic Plants, Moscow, Russia
Distal astrocytic processes have a complex morphology, reminiscent of branchlets and leaflets. Astrocytic branchlets are rod-like processes containing mitochondria and endoplasmic reticulum, capable of generating inositol-3-phosphate (IP3)-dependent Ca2+ signals. Leaflets are small and flat processes that protrude from branchlets and fill the space between synapses. Here we use three-dimensional (3D) reconstructions from serial section electron microscopy (EM) of rat CA1 hippocampal neuropil to determine the astrocytic coverage of dendritic spines, shafts and axonal boutons. The distance to the maximum of the astrocyte volume fraction (VF) correlated with the size of the spine when calculated from the center of mass of the postsynaptic density (PSD) or from the edge of the PSD, but not from the spine surface. This suggests that the astrocytic coverage of small and larger spines is similar in hippocampal neuropil. Diffusion simulations showed that such synaptic microenvironment favors glutamate spillover and extrasynaptic receptor activation at smaller spines. We used complexity and entropy measures to characterize astrocytic branchlets and leaflets. The 2D projections of astrocytic branchlets had smaller spatial complexity and entropy than leaflets, consistent with the higher structural complexity and less organized distribution of leaflets. The VF of astrocytic leaflets was highest around dendritic spines, lower around axonal boutons and lowest around dendritic shafts. In contrast, the VF of astrocytic branchlets was similarly low around these three neuronal compartments. Taken together, these results suggest that astrocytic leaflets preferentially contact synapses as opposed to the dendritic shaft, an arrangement that might favor neurotransmitter spillover and extrasynaptic receptor activation along dendritic shafts.
Synaptic contacts are highly specialized neuronal compartments. The presynaptic terminal is an axonal varicosity that is closely apposed to the postsynaptic terminal. The dendritic spine is the postsynaptic terminal of a glutamatergic synapse. Dendritic spines are highly dynamic structures that can be classified into several groups based on their morphology (e.g., mushroom spines, thin spines etc.) (Hering and Sheng, 2001). Spines have a neck and a head with an electron dense PSD. The space between synapses is occupied by astrocytic processes (Kettenmann et al., 1996; Chao et al., 2002; Witcher et al., 2007; Verkhratsky and Nedergaard, 2018).
The intimate anatomical and physiological interaction between spines, boutons and perisynapstic astrocytic process (PAP) has been referred to as the “synaptic triate” (Kettenmann et al., 1996), “tri-partite synapse” (Araque et al., 1999), or “astrocytic cradle” (Nedergaard and Verkhratsky, 2012; Verkhratsky and Nedergaard, 2014). PAPs have a high density of glutamate transporters and of inward rectifying K+ channels. They are thought to shape both synaptic transmission and extrasynaptic signaling via glutamate buffering/uptake and K+ clearance, respectively (Walz, 2000; Danbolt, 2001; Chever et al., 2010; Sibille et al., 2014). Astrocytic uptake limits both glutamate escape out of the synaptic cleft (Rusakov and Kullmann, 1998; Diamond, 2001; Huang and Bordey, 2004; Rose et al., 2018) and its diffusion into the synaptic cleft (Lozovaya et al., 2004; Wu et al., 2012). K+ clearance reduces local, activity-dependent extracellular K+ accumulation, which can depolarize the presynaptic terminal and facilitate the neurotransmitter release (Poolos et al., 1987; Shih et al., 2013; Sibille et al., 2014; Contini et al., 2017; Lebedeva et al., 2018). PAPs also contain glycogen granules and can metabolically interact with neighboring synapses (Phelps, 1972; Koizumi, 1974; Calı̀ et al., 2016). Astrocytic glycogen is converted into lactate and delivered to neurons through monocarboxylate transporters (MCTs) to be used as an energy substrate (Quistorff et al., 2008; Matsui et al., 2017). This exchange of metabolites, commonly referred to as the “astrocyte-neuron lactate shuttle” (ANLS), is upregulated during sleep deprivation and exhaustive exercise (Magistretti, 2011; Matsui et al., 2017). In addition to the energy substrate function, lactate released by astrocytes can serve as a signaling molecule/gliotransmitter to excite neurons and to modulate astrocytic Ca2+ activity (Tang et al., 2014; Lebedeva et al., 2015; Magistretti and Allaman, 2018). Astrocytes release several other gliotransmitters such as D-serine, glycine, glutamate, GABA, ATP, adenosine etc (Zorec et al., 2012; Araque et al., 2014). Astrocytic D-serine facilitates synaptic plasticity (Panatier et al., 2006; Henneberger et al., 2010; Papouin et al., 2017). Glutamate activates presynaptic metabotropic (Semyanov and Kullmann, 2000; Liu et al., 2004b) and axonal kainate receptors (Semyanov and Kullmann, 2001; Liu et al., 2004a), which in turn shape the profile of neurotransmitter release (Kullmann and Semyanov, 2002; Gordleeva et al., 2012). GABA release from astrocytes might contribute to increase of tonic GABAA receptor activation in neurons (Semyanov et al., 2004; Angulo et al., 2008; Héja et al., 2012). Despite the wealth of physiological data suggesting multiple types of interaction between the synapse and PAPs, there is no firm anatomical evidence for the existence of particular types of PAPs belonging to a particular type of “tri-partite” synapse. Our general understanding, based on 3D EM reconstructions in the hippocampal neuropil, is that astrocytes are sponge-like structures just filling all the space between synapses, dendrites, axons, etc. (Witcher et al., 2007; Patrushev et al., 2013; Medvedev et al., 2014). Nevertheless, in the rodent hippocampus the astrocyte VF is larger around dendritic spines than axonal boutons (Lehre and Rusakov, 2002). There is, however, an ongoing contentious. Although some ultrastructural studies indicate that there is a higher density of astrocytic processes around large mushroom spines (Jones and Greenough, 1996; Genoud et al., 2006; Lushnikova et al., 2009; Bernardinelli et al., 2014a,b), other works do not confirm these results, possibly due to the use of different analytical approaches (Medvedev et al., 2014; Heller and Rusakov, 2015). For this reason, a more detailed structural analysis of astrocytic coverage at hippocampal synapses is needed. Recent reports have indicated that astrocytic processes are heterogeneous and can be distinguished into thin organelle-free structures (leaflets) and thick processes that host endoplasmic reticulum (ER) and mitochondria (larger branches, smaller branchlets, and perivascular endfeet) (Patrushev et al., 2013; Khakh and Sofroniew, 2015).
Here we study how astrocytic leaflets and branchlets distributed around different types of dendritic spines, dendritic shafts and axonal boutons.
Materials and Methods
Sample Preparation and Electron Microscopy
Wistar rats (male, 250 ± 25 g) were transcardially perfused with a mixture of 3% paraformaldehyde and 0.5% glutaraldehyde in 0.1 M Na-cacodylate buffer (pH 7.2–7.4). The perfused brains were cut into 150 μm thick slices, which were fixed with 2.5% glutaraldehyde for 24 h and impregnated with 1% osmium tetroxide and 0.01% potassium dichromate for 1–2 h. The slices were then dehydrated in graded aqueous solutions of ethanol from 40 to 96% (10 min each) and finally in 100% acetone (three washes, 10 min each). Dehydrated slices were then incubated in a mixture of 50% epoxy resin and 50% pure acetone for 30 min. They were then embedded in a capsule of pure epoxy resin (Epon 812/AralditeM) at 60°C for 1 h and polymerized overnight at 80 °C.
Serial sections (60–70 nm) of CA1 str.radiatum were prepared using a Diatome diamond knife and were collected using Pioloform-coated slot copper grids. Sections were counterstained with saturated ethanolic uranyl acetate, followed by lead citrate, and then placed in a grid holder. Images were obtained with a JEOL 1010 EM at 6000x magnification.
EM negatives were scanned at 1200 dpi resolution. The scanned images were aligned using custom-made software. Complete reconstruction of astrocytic processes, dendrites, dendritic spines, PSDs and axons was performed using the software Reconstruct (by Dr. John Fiala, provided free of charge at SynapseWeb, Kristen M. Harris, PI) after manually tracing the contours of these structures (Supplementary Figure 1). The meshes generated by the Reconstruct software served as input data for our morphological analysis. The total volume of the reconstructed block was 393 μm3.
Equidistant Surface Analysis
PSDs, dendritic spines, axonal boutons and dendritic shaft chunks, were enclosed by a triangulated sphere for the equidistant surface analysis. The sphere was then transformed into the equidistant surface at a given distance (d) (Supplementary Figure 2): We calculated the shortest distance between the object and each vertex of the triangulated sphere. If this distance was smaller than d, then the vertex was moved away from the object to the distance equal to d. If this distance was larger than d, then the vertex was moved closer to the object and to the distance equal to d. Rearranging all the vertexes of the triangulated sphere gave a satisfactory approximation of the equidistant surface.
For astrocytic coverage analysis, equidistant surfaces were plotted (1) around the center of mass of the PSD, (2) around the edges of the PSD or (3) around the spine surface. In the first case, the equidistant surfaces were concentric spheres with d = 0 at the PSD center of mass. In the second case, d = 0 was at the edge of the PSD. In the third case, d = 0 was at the side spine surface: First, the equidistant surfaces were plotted around the entire reconstructed spine. Then parts of these surfaces which were inside presynaptic half-space and inside dendritic half-space were excluded (Supplementary Figure 2).
Dendritic shafts could be were recognized as long objects striding the whole reconstructed block. Equidistant surfaces were not generated around whole dendritic shafts because this would produce spatial averaging of their uneven microenvironment. Therefore, dendritic shafts were split into chunks of 1 μm length, which were treated as separate objects, and equidistant surfaces were plotted around their membrane-covered parts.
Axonal boutons radii were calculated as described in Supplementary Figure 2C.
Glutamate diffusion from synapses of different radii was modeled using Matlab 2016b (Mathworks; Natick, MA). We positioned one synapse at the center of a 10 μm3 world (Diamond, 2005). The pre- and post-synaptic terminals were represented as two hemispheres separated by a 20-nm thick synaptic cleft and surrounded by a 50 nm-wide extrasynaptic region. The radius of the PSD was calculated as 1/3 of the sum of the synapse radius and the extrasynaptic region. A perisynaptic annulus separated the PSD from the extrasynaptic region. The width of the perisynaptic region was calculated by subtracting the PSD area from the area covered by the synaptic and extrasynaptic regions. At the beginning of each simulation, we released 2,000 glutamate molecules from a point source located at the center of the synaptic cleft. Each glutamate molecule diffused with an apparent diffusion coefficient D* = 0.33 μm2/ms, according to the experimental measures and theoretical predictions at 36–37°C (Nielsen et al., 2004). Outside the cleft, D* was reduced to account for the tortuosity of the hippocampal neuropil, measured experimentally in P14-21 mice (λ = 1.45) (Scimemi et al., 2009). The space surrounding the extrasynaptic volume contained glutamate transporters at a concentration of 20 μM (Lehre and Danbolt, 1998), each of which was represented using a simplified Markov model of the glial glutamate transporter GLT1 (Bergles et al., 2002; Diamond, 2005; Scimemi et al., 2009). All reaction rates were corrected for Q10 = 3 to account for the temperature dependence of the glutamate transport process. The extracellular space was subdivided into 10-nm thick concentric shells and the simulations were run with a time step of 1 μs for a total of 10,000 steps (10 ms). At each time step, we calculated the number of glutamate molecules in the volume of the synaptic cleft above the PSD area and in the 50 nm-wide extrasynaptic volume. The obtained glutamate concentration waveforms were used to drive GluN (Lester and Jahr, 1992) and GluA Markov models (Jonas, 1993) in ChanneLab2 (Synaptosoft; Decatur, GA). The number of glutamate molecules bound by glutamate transporters were used to calculate the time course of the STCs.
The source code for these simulations is available at https://sites.google.com/site/scimemilab2013/software.
Complexity-Entropy Spectrum Analysis
The goal of the complexity-entropy analysis was to quantify the spatial properties of representative 2D images with respect to their balance between randomness and structural order, triviality and complexity. In mathematical terms, this is described as location in the complexity-entropy plane or their distribution along a complexity-entropy spectrum. Highly ordered structures (e.g., a regular grating) have near-zero entropy and near-zero complexity. In contrast, completely disordered structures (e.g., independent and identically distributed Gaussian samples) have maximal entropy and small statistical complexity. Intermediate values of entropy are associated with higher values of complexity if the underlying pattern contains features with preferred orientation. For example, signals generated by systems with deterministic chaos result in a complexity-entropy spectrum with near-0.5 peak complexity and near-0.7 entropy (Lamberti et al., 2004; Rosso et al., 2007).
Below we provide a concise description of the complexity-entropy analysis theory and implementation. More in-depth theoretical overview on the entropic complexity measures can be found in (Lamberti et al., 2004; Martin et al., 2006; Rosso et al., 2007) and the implementation is described in a method paper (Brazhe, 2018). The implementation code is available at DOI:10.5281/zenodo.1217636.
Complexity and entropy measures can be based on a feature distribution of the analyzed patterns, which can be compared to equiprobable or singular feature cases. Accordingly, the relative entropy (H) was defined as the Shannon entropy:
normalized by the entropy of an equally probable distribution (Smax), thus giving values in the range [0 … 1]. Here N is the number of features analyzed. The complexity measures were based on the notion of disequilibrium (Lamberti et al., 2004) and defined in terms of normalized Jensen-Shannon divergence from an equally probable distribution (Rosso et al., 2007):
In this expression, J[P,Pe] is Jensen-Shannon divergence of some distribution P from an equally probable distribution Pe:
and Jmax is obtained when the probability of just one feature is equal to 1, while the probability of all other features is zero.
Armed with the probability-based definitions of entropy and complexity, one needs a way to build descriptive feature probabilities from spatial patterns, preferably allowing for multiscale resolution and local analysis. A spatial pattern can be described in terms of local edge orientations and scales. Thus, statistics of the shearlet transform coefficients appears as a promising probability-inducing image decomposition (Kutyniok and Labate, 2012). Shearlets are a novel image analysis tool, an orientation-aware analog of wavelets for multivariate data. In short, the shearlet transform separates an input image into a spectrum of local spatial scales and orientations. Thereafter, we treated normalized power of shearlet coefficients as probability densities of given spatial feature orientation and scale at a given location.
The complexity-entropy analysis is illustrated on a toy example of a gradually corrupted stripe pattern in Supplementary Figure 3. Before introducing noise, the spatial pattern is ordered and corresponds to low entropy–low complexity region of the spectrum. Disrupting the image by swapping intensity values with between neighboring pixels increases entropy and is accompanied by a transient rise in complexity due to changes in the distribution of the power of shearlet coefficients.
We performed the complexity-entropy spectrum analysis of the spatial organization of leaflets and branchlets separately. Segmentation images (i.e., images with pixel values set to one within identified leaflets or branchlets or to zero otherwise) were sum-projected along the Z-axis and down-sampled 2 times in the X and Y directions to be treatable by the shearlet transform library (Gregory R. Lee, https://github.com/grlee77/PyShearlets). The resulting images represented the abundance of the given astrocytic structure in each point. Entropy and complexity calculations were only performed in regions containing non-zero pixels, while other areas (black in Figure 5A) were masked out.
The resulting pairs of H and complexity (C) values at each pixel were then binned and mapped onto the complexity-entropy plane (Figure 5B). To condense the apparent differences between the spectra of branchlets and leaflets spatial patterns, we plotted a “median curve,” which represents median complexity for all (branchlets + leaflets) points at each entropy value. This was achieved by binning all points into 51 intervals of entropy values, calculating the median complexity values for each interval, and then fitting the resulting set of points with a smooth cubic spline. This allowed us to use an “excess complexity” parameter to describe each (H,C) point by a single value of signed deviation from the “median curve.” The resulting distributions of excess complexity (kernel density estimates) for branchlets and leaflets spatial patterns are shown in Figure 5C.
The results of the presented complexity/entropy analysis can potentially depend on the orientation of the astrocyte relative to the imaged block. Astrocytes are 3D structures, and using a 3D measure of complexity and entropy might provide more information than analysis of projections along one axis. However, increasing dimensionality could deteriorate the shearlet-based definition of spatial feature probabilities because higher degrees of freedom would demand larger samples. Ideally, if resolution is equal along all three axes (X-Y-Z), one could draw a plane corresponding to the two first principal axes of the cloud of astrocytic voxels, and project all structures to this principal plane or to a set of planes, parallel to the principal one. However, in the case of serial section EM, Z-axis resolution is typically lower that X-Y resolution. Therefore, only the projection along Z-axis is possible.
The data are presented as mean ± standard error of mean (SEM). The statistical difference was tested with Mann–Whitney test using Origin 8 (OriginLab Corp.). The level of significance was set at P < 0.05.
Astrocytic Coverage of Thin and Mushroom Spines
Serial EM sections were used to generate a full 3D reconstruction of a 393 μm3 block from hippocampal str.radiatum. Astrocytic processes, dendritic shafts, dendritic spines, axons and axonal varicosities were visually identified, manually traced and 3D rendered. We analyzed the spatial distribution of astrocyte VF around dendritic spines using three different approaches (Figure 1A). (1) We calculated the astrocyte VF at increasing distances from the center of mass of the PSD (similar to the approach used by Medvedev et al., 2014). (2) We calculated the astrocyte VF at increasing distances from the edge of the PSD (similar to the approach used by Patrushev et al., 2013). (3) We calculated astrocyte VF at increasing distances from the surface of the dendritic spine (without including the volume occupied by the presynaptic half-space and dendritic shaft half-space). The first approach estimates the distribution of astrocyte around the presynaptic glutamate release site, which is closely apposed to the center of mass of the PSD. This method does not exclude the portion of the cleft around the PSD, where no astrocyte is typically present. Because synapses vary in size, the astrocyte-free area also varies with the size of the synapse. Therefore, astrocytic processes are located further away from the centers of larger synapses, suggesting these might have fewer astrocytic processes around them (Figures 1B,C). The second approach works around the fact that larger spines have larger PSDs, but the method is still prone to errors because the PSD does not cover the entire apposition zone, which is also bigger at bigger spines (Figures 1B,C). Thus, both approaches provide a measure of astrocyte VF without completely excluding the volume occupied by the synapse. In contrast, the third approach provides a genuine measure of astrocyte VF around dendritic spines. As expected, our estimates for the distance from the spine to the maximum astrocyte VF (VFmax) varied depending on the method of analysis. When using the first method, the distance to VFmax depended on the spine head radius (Figure 1D and Supplementary Figure 4). For the second method, this trend was less pronounced. For the third method, there was no correlation between the distance to VFmax and the spine radius. Specifically, 28.9% of spines had VFmax in immediate apposition to their membrane as if they formed as “astrocytic cradle” rather than a non-specific filling of the extracellular space (Verkhratsky and Nedergaard, 2014).
Figure 1. Astrocyte VF distribution around dendritic spines. (A) Astrocyte VF was measured between equidistant surfaces constructed around PSD center (top), around PSD surface (middle) and around side surfaces of the dendritic spine (bottom). (B) Changes of astrocyte VF with the distance from PSD center (top), from PSD surface (middle) and from spine surface (bottom) for 5 individual spines. Double-sided arrow indicates the distance to the VFmax. (C) Summary data as in (B) are presented as mean ± SEM for 207 spines. (D) Percentage of dendritic spines with VFmax at a different distance from PSD center (top), PSD surface (middle), and spine surface (bottom). The bin size of 30 nm was equal to the distance step in (A–C).
Based on morphological and physiological properties, dendritic spines can be subdivided into several classes (e.g., mushroom spines, thin spines) (Hering and Sheng, 2001). In our analysis, we considered spines with head radius < 0.2 μm as thin spines and spines with head radius > 0.2 μm as mushroom spines (Figures 2A,B). We then plotted astrocyte VF at increasing distances for spines of different sizes (Figure 2C).
Figure 2. Astrocytic coverage of mushroom and thin spines. Segmentation of mushroom and thin spines was based on spine head radius. (A) Examples of thin and mushroom spines. (B) Distribution of spines according to their head radii. Spine types were segmented by the radius of 0.2 μm: mushroom spines were defined as spines with radius ≥ 0.2 μm, thin spines were defined as spines with radius < 0.2 μm. (C) Changes of astrocyte VF with the distance from PSD center (left), from PSD surface (middle) and from spine surface (right) for thin (red circles) and mushroom (blue circles) spines. The data are presented as mean ± SEM.
Consistent with the above result, VF at a given distance (100 nm) from the center of mass of the PSD or from the edge of the PSD was larger at thin spines than at mushroom spines (from the center of mass: 5.5 ± 1.1% at thin spines, n = 128, and 2.2 ± 0.9% at mushroom spines, n = 79; p < 0.001 Mann–Whitney test; from the edge: 10.3 ± 0.9% at thin spines, n = 128, and 5.2 ± 0.8% at mushroom spines, n = 79; p < 0.001 Mann–Whitney test). In contrast, there was no significant difference in astrocyte VF at 100 nm from the spine surface between thin and mushroom spines (14.2 ± 1.2% at thin spines, n = 128, and 14.6 ± 1.4% at mushroom spines, n = 79; p = 0.44 Mann–Whitney test). This result suggests that the astrocytic coverage does not depend on the spine size. This finding, however, does not rule out that activity-dependent differences in the astrocytic coverage of spines might occur within each morphological group.
Functional Implications for Extrasynaptic Receptor Activation
If the level of astrocytic coverage does not scale with the spine size, what effect does this have on the activation of synaptic and extrasynaptic GluA and GluN receptors? To address this question, we performed computer simulations of glutamate diffusion. The pre- and postsynaptic terminals were represented as two hemispheres separated by a 20-nm thick synaptic cleft surrounded by a 50-nm wide extrasynaptic region. The synaptic radius (r) was varied from 50 nm to 450 nm and the radius of the PSD from 17 to 150 nm, to account for the presence of larger PSD areas at bigger spines (Ventura and Harris, 1999). We measured the time course of the glutamate concentration in the volume of the cleft above the PSD area and in the extrasynaptic volume, and used these measures to determine the open probability of GluA/N receptors and the time course of glutamate transporter current. The peak glutamate concentration became progressively smaller and prolonged in its rise and decay time at increasing values of r (Figures 3A,B), as the volume above the PSD and extrasynaptic regions increased (Figure 3A right panel, inset). The open probability of GluA and GluN receptors in the PSD did not vary appreciably among spines of different dimensions (Figures 3C,D, left and right panels). In contrast, the open probability of these receptors in the extrasynaptic region varied profoundly between small and larger spines (Figures 3C,D, middle and right panels). The open probability of extrasynaptic GluA receptors peaked at r = 0.15 μm and then progressively decreased (Figure 3C, right panel). GluN receptors have a higher steady-state affinity for glutamate with respect to GluA receptors (Tang et al., 1989; Trussell and Fischbach, 1989; Patneau and Mayer, 1990) (Lester et al., 1990). Their peak open probability was generally higher than that of GluA receptors (Figure 3D, middle panel), but decayed more abruptly in the extrasynaptic region of larger spines. Consequently, there was a progressive reduction in the relative activation of GluN versus GluA receptors in extrasynaptic space of bigger spines (Figure 3E). Glutamate clearance by glial transporters was delayed at larger spines, but its time course (measured as the centroid of synaptically-activated transporter currents, < t>) was similar to the one measured at smaller spines (Figure 3F). In the hippocampus, the activation of GluN receptors is crucial for the induction of long-term potentiation and depression (LTP and LTD, respectively). LTP/LTD expression is associated with spine enlargement/shrinkage, respectively (Lang et al., 2004; Nägerl et al., 2004; Zhou et al., 2004). Therefore, our findings suggest that the increased extrasynaptic GluN activation at the smaller spines might render these synapses more prone to express synaptic plasticity compared to larger spines.
Figure 3. Extrasynaptic GluN activation is favored at smaller spines. (A) Time courses of the glutamate concentration in the volume of the synaptic cleft above the PSD (red, left panel) and in the extrasynaptic region (green, middle panel), at synapses of increasing radius (r = 0.05–0.45 μm). Right panel: summary graph, showing the progressively smaller peak glutamate concentration in the cleft and in the extrasynaptic region of synapses of increasing r (inset). (B) Peak-normalized time courses of the glutamate concentration shown in (A). Right panel: summary graph showing the increased decay time of the glutamate concentration profile at larger synapses. (C) Open probability (Po) of GluA receptors in thePSD and extrasynaptic region. Right panel: summary graph of the GluA Po normalized by its value at r = 0.05 μm. The GluA Po in the PSD does not change at small and large synapses (red). In contrast, the Po of extrasynaptic GluA receptors (green) is larger at small than large synapses. (D) GluN Po in the PSD and extrasynaptic region. Right panel: summary graph of the GluN Po, normalized by its value at r = 0.05 μm. The GluN Po in the PSD region does not change at small and large synapses (red). In contrast, the Po of extrasynaptic GluN receptors (green) decreases abruptly as the synaptic radius increases. (E). Ratio of the GluN and GluA Po at synapses with different size. The GluN Po is larger than the GluA Po in and out of the PSD at small synapses. A reduction in the Po of GluN vs. GluA receptors occurs as the synaptic radius increases. Right panel: summary graph of the GluN/A Po normalized by its value at r = 0.05 μm. There is a progressive reduction in the relative Po of GluN and GluA receptors at larger synapses. (F) Peak normalized time course of synaptically-activated transporter currents (STCs) at astrocytes surrounding an active synapse with increasing radius. Right panel: summary graph of the STC latency and centroid (< t>) normalized by their values at synapses with r = 0.05 μm. A delay in the onset of the STC occurs at large synapses. This effect is not associated with changes in the time course of the STC.
Analysis of Astrocytic Branchlets and Leaflets
Since astrocytic leaflets are thin and flat, while branchlets are rod-like structures, they can be distinguished according to their surface-to-volume ratio (SVR) (Patrushev et al., 2013). To establish local SVRs we employed previously suggested secant spheres method (Patrushev et al., 2013). Following this method, a 300 nm-radius sphere is was centered along each point of the astrocyte surface. The we analyzed the volume and the surface area of astrocytic chunks confined inside of the sphere. We plotted the relation between the surface area and the volume of each chunk and used a K-means algorithm to identify two clusters of points (Figure 4A). This method allowed us to satisfactorily separate leaflets and branchlets according to their SVR distributions (Figures 4B,C and Supplementary Video).
Figure 4. Method for sorting of astrocytic branchlets and leaflets. (A) Secant sphere with a 600-nm radius centered at each point on the astrocyte surface was used. Surface and VF of the astrocyte process within this sphere was calculated. The relationship between surface and VF was plotted as circles for individual points on astrocyte surface. The circles were divided into two clusters with K-means algorithm. Circles to the left (71% of all points) were marked in red, circles to the right (29% of all points) were marked in green. (B) Distribution of “green” and “red” marked astrocyte processes according to their SVR. (C) 3D reconstruction of astrocyte processes marked in green and red according to the classification presented in (A,B). The green marked processes predominantly represent astrocytic branchlets, while the red marked processes represent leaflets.
The spatial organization of leaflets and branchlets was then analyzed using the complexity-entropy spectrum analysis. To estimate local complexity and entropy, we redefined the probability measure used by Zunino and Ribeiro (Ribeiro et al., 2012; Zunino and Ribeiro, 2016). Instead of using permutation entropy (Bandt and Pompe, 2002), we used local statistics of shearlet transform coefficients (Kutyniok and Labate, 2012). The analysis was performed on mean projections along the Z-axis of the segmentation stacks either of astrocytic leaflets or of astrocytic branchlets (Figure 5A). Well-structured astrocytic branchlets had low entropy and low complexity, while the space between the branchlets covered a wider spectrum of entropy-complexity pairs (H,C) (Figure 5A). The astrocytic leaflets formed a more chaotic image, which is reflected in higher entropy values. The presence of coincidental anisotropic spatial features across several spatial scales accounted for zones of higher spatial complexity (Figure 5A).
Figure 5. Spatial entropy and complexity of astrocytic processes. (A) Mean projection of segmented astrocytic processes (left), spatial entropy (middle), and spatial complexity (right) for astrocytic branchlets (top) and leaflets (bottom). Areas not containing any structures were assigned 0 entropy and complexity. (B) Middle, distributions of complexity and entropy for astrocytic branchlets (magenta) and leaflets (green). Color intensity is log-probability of the corresponding (H,C) pair to be found in the images, rescaled for better contrast. Black dots—overlap between the two color channels. Yellow curve: local median complexity for given entropy, taking all values (from branchlets and leaflets) together. Top, marginal distributions of entropy, kernel density estimate. Right, marginal distributions of complexity, kernel density estimate. (C) Distributions of “excess complexity,” i.e., difference between C values and local median curve, kernel density estimate.
Thus, branchlets and leaflets followed distinct complexity-entropy spectra (Figure 5B): (1) (H,C) pairs for branchlets extended to the area of low entropy and of low complexity, while (H,C) pairs for the leaflets reached higher maximal complexity and were clustered in the area of higher entropy and of higher complexity; (2) in the areas where the two distributions overlapped, the branchlets tended to have higher complexity for the same entropy values than leaflets. These differences become more evident if one compares deviations of complexity values from the median curve calculated for the pooled (H,C) values from the two groups (Figure 5C). Most of the (H,C) values for leaflets were below the median curve (negative deviations), whereas most values for branchlets were above the curve (positive deviations).
Next, we tested if astrocytic leaflets can be classified according to their morphological properties likewise dendritic spines in neurons (Hering and Sheng, 2001). We segmented whole leaflets as complete astrocytic processes with high SVR, devoid of organelles and connected to parent astrocytic branchlet (Figure 6A). Volume, surface and SVR of whole individual leaflets were analyzed. All these parameters had skewed distributions without obvious peaks which would indicate the existence of different classes of the leaflets (Figures 6B–D). Therefore, the leaflets were not subdivided further in this study. However, further studies based on other parameters may suggest a classification of astrocytic leaflets. For example, the distribution of the distances from leaflet to its nearest neighbor leaflet on the same parent branchlet could be fitted with two Gaussian distributions with means of 0.23 and 0.58 μm (Supplementary Figure 5), but no correlation of the distance to the nearest neighbor with the leaflet volume was observed (Supplementary Figure 5).
Figure 6. Leaflets parameters distributions. (A) Complete (uncut) structures identified as leaflets according to their morphological properties and having a connection to parent branchlet considered as individual leaflets. Whole leaflet volume, surface and SVR was measured. (B) Distribution of leaflet volumes. (C) Distribution of leaflet surfaces. (D) Distribution of leaflet SVRs.
Asymmetric Astrocytic Coverage of Axonal Boutons and Dendritic Shafts
Next, axonal boutons and dendritic shafts were reconstructed. Consistent with previous studies, the distribution of boutons radii pointed to the diversity of these structures (De Paola et al., 2006; Hara et al., 2011; Grillo et al., 2013). However, further classification of axonal boutons was beyond the scope of this study, and they were treated as a single pool (Figures 7A,B). Dendritic shafts were fragmented into 1 μm long chunks which were approximately in the size range of dendritic spines and axonal boutons (Figure 7C). We then generated equidistant surfaces around each bouton and chunk of dendritic shaft to estimate astrocyte VF distribution around these structures in comparison to dendritic spines. The highest VFmax of astrocyte processes was found around dendritic spines (VFmax: 19.2 ± 0.9%, n = 207, Figure 7D; Note, that VFmax was located at a different distance from each spine; therefore, mean VFmax is higher than the peak of the mean VF in Figure 7D). VFmax of astrocyte processes around axonal boutons was significantly smaller (VFmax: 15.6 ± 0.97%, n = 139, p = 0.007 for difference with spines, Mann-Whitney test, Figure 7D). Such asymmetry of glia near central synapses has been suggested as a mechanism favoring glutamate spillover onto presynaptic autoreceptors, which are responsible for activity-dependent regulation of glutamate release probability (Lehre and Rusakov, 2002). VFmax of astrocytic processes around dendritic shaft was smallest (VFmax: 11 ± 1%, n = 96, p < 0.001 for difference with spines, p < 0.001 for difference with boutons, Mann–Whitney test, Figure 7D). These results suggest that astrocytic processes preferentially approach synaptic structures, leaving the dendritic shaft less covered. Such arrangement of the processes may not only serve to increase the efficiency of glutamate uptake and K+ clearance but also to isolate synapse from the extrasynaptic glutamate (Lozovaya et al., 2004; Wu et al., 2012). On the other hand, lower coverage of dendritic shafts by astrocytic processes will make glutamate uptake less efficient around this compartment, and this will facilitate activation of extrasynaptic receptors by extracellular glutamate.
Figure 7. Astrocyte VF distribution around axonal boutons, dendritic shafts and dendritic spines. (A) 3D reconstruction of axonal boutons (yellow) and dendritic spines (gray). (B) Distribution of axonal boutons according to their radii. (C) Top, reconstruction of dendritic shafts. The shafts were fragmented to 1 μm long chunks. Bottom, equidistant surfaces were plotted around the side surface of each chunk. The VF of astrocytic processes was estimated between these surfaces. (D) Distribution of astrocyte VF with the distance from the spine (blue circles), axonal bouton (red circles) and dendritic shaft (green circles). (E) Distribution of leaflets VF with the distance from spine (blue circles), axonal bouton (red circles) and dendritic shaft (green circles). (F) Distribution of branchlets VF with the distance from spine (blue circles), axonal bouton (red circles) and dendritic shaft (green circles). The data on panels (D–F) are presented as mean ± SEM.
Astrocytic Leaflets But Not Branchlets Form PAP
Previously, we have already reported that hippocampal synapses are covered by organelle-free leaflets but not by astrocytic branchlets containing Ca2+ stores (Patrushev et al., 2013). Here we compared the distribution of astrocytic leaflets and of astrocytic branchlets around axonal boutons, dendritic spines and dendritic shafts. Interestingly, different VFmax of astrocyte around these neuronal compartments correlated with different VFmax of astrocytic leaflets (VFmax for spines: 15 ± 1%, n = 207; VFmax for boutons: 10.6 ± 0.7% p < 0.001 for difference with spines; VFmax for dendritic shafts: 7.4 ± 0.6%, n = 96, p < 0.001 for difference with spines, p < 0.001 for difference with boutons; Mann–Whitney test, Figure 7E). There were no compartment-specific differences in VFmax of astrocytic branchlets (VFmax for spines: 5.6 ± 0.7%, n = 207; VFmax for boutons: 6.3 ± 0.9% p = 0.96 for difference with spines; VFmax for dendritic shafts: 4.6 ± 0.8%, n = 96, p = 0.9 for difference with spines, p = 0.31 for difference with boutons; Mann–Whitney test, Figure 7F). This finding suggests that astrocytic branchlets host leaflets that extend toward synapses to cover them.
Because leaflets are devoid of organelles, including Ca2+ stores (ER and mitochondria), they cannot support Ca2+ release from endogenous Ca2+ stores as means of neuron-astrocyte signaling (Li and Rinzel, 1994; Agarwal et al., 2017; Sherwood et al., 2017). Nevertheless, there is a possibility that some synapses are located close to the astrocytic branchlets containing Ca2+ stores. We found that both thin and mushroom dendritic spines are predominantly located within the area of the leaflets, but not in proximity to astrocytic branchlets (Figure 8A). 72.7% of thin spines and 83.5% of mushroom spines had direct contact with leaflets, but only 22.8% of thin spines and 23.4% of mushroom spines had direct contact with the astrocytic branchlets (Figure 8B). This finding suggests that the majority of glutamate synapses cannot directly trigger Ca2+ release from endogenous Ca2+ stores; however, it does not rule out that they trigger Ca2+ entry through plasma membrane in leaflets (e.g., due to glutamate uptake-dependent activation of Na+/Ca2+ exchanger) (Reyes et al., 2012; Bazargani and Attwell, 2016).
Figure 8. Distribution of mushroom and thin spines around astrocytic branchlets and leaflets. (A) Top left, distribution of mushroom spines (gray) around astrocytic branchlets (green). Top right, distribution of thin spines (gray) around astrocytic branchlets (green). PSDs – red. Bottom left, distribution of mushroom spines around leaflets (red). Bottom right, distribution of thin spines around leaflets (red). (B) Percentage of thin (empty circles) and mushroom (filled circles) spines that have leaflets (red lines) or astrocytic branchlets (green traces) at different distances. Dashed line corresponds to 50 nm distance, which was considered a direct contact between spine and astrocytic process.
Astrocytic coverage of synapses plays an important role in shaping synaptic transmission and neurotransmitter spillover in the brain. We used three different methods to estimate the distribution of astrocyte VF around dendritic spines in CA1 hippocampal neuropil: from the center of mass of the PSD, from the edge of PSD and from the spine surface. When measured from the center of mass or from the edge of the PSD, the maximum astrocyte VF appeared closer to the PSD of smaller spines. However, astrocyte VF around the spine surface did not depend on the spine radius. This finding suggests that the spine size does not actually determine the astrocytic VF in it's microenvironment, it was previously suggested (Lushnikova et al., 2009; Bernardinelli et al., 2014a; Medvedev et al., 2014). When the spine size increases, it simply pushes the PAPs further away from the PSD and presynaptic active zone. This also leads to smaller and slower extrasynaptic glutamate concentration profiles. The net effect is a progressive reduction in the activation of extrasynaptic receptors around larger synapses, particularly pronounced for GluN receptors. This finding has multiple intriguing implications because it suggests that: (1) activation of extrasynaptic receptors because of glutamate spillover is more likely to occur around smaller synapses; (2) smaller spines are more likely to undergo long-term plasticity dependent on extrasynaptic GluN receptor activation; (3) larger spines are optimized for point-to-point communication, whereas smaller spines also mediate extrasynaptic signaling. Thus, the ratio between thin and mushroom spines shapes the propensity of a neuron to communicate through point-to-point or volume transmission in the hippocampus. In addition, larger spines have a lower SVR. A similar density of PAPs at small and larger spines may lead to reduced metabolic support per unit of cytoplasmic volume. This can potentially limit metabolic processes in larger spines relatively to thin spines, and thus serve as another plasticity limiting factor (Magistretti, 2011).
Astrocytic processes are morphologically highly plastic structures (Genoud et al., 2006; Theodosis et al., 2008; Tanaka et al., 2013; Bernardinelli et al., 2014a; Heller and Rusakov, 2015). Our findings suggest that the morphological plasticity of astrocytic processes may not be directly related to morphological plasticity of the associated synapses. Further studies are required to distinguish between different possibilities such as: (1) astrocytic morphological plasticity is not specific to particular synapses but is associated to groups of local synapses; or (2) it is synapse specific but is triggered independently of morphological plasticity of the spines. Alternatively, PAPs can regulate synaptic morphological plasticity themselves. For example, a contact with PAPs inhibits spine enlargement during memory consolidation in the amygdala (Ostroff et al., 2014). In fact, morphological remodeling is only one form of plasticity in PAPs. Two-fold increase in glutamate transporter density has been reported in mouse cortical astrocytes following whisker stimulation (Genoud et al., 2006).
The astrocytic processes can be categorized into primary branches, secondary branchlets (rod-like processes containing organelles, including ER and mitochondria) and leaflets (thin sheet-like organelle-free compartments) (Fernandez et al., 1984; Patrushev et al., 2013; Khakh and Sofroniew, 2015). This organization is reminiscent of that of spiny neuronal dendrites: astrocytic leaflets protrude from branchlets just like spines protrude from dendritic shafts. By using a K-means algorithm, we distinguished leaflets and branchlets based on their local SVR. Then we compared the spatial organization of astrocytic branchlets and leaflets using a spatial complexity-entropy spectra analysis. Because leaflets have different shapes and orientations, their 2D projection is characterized by a higher entropy than that of branchlets. In fact, branchlets introduced areas with a preferred characteristic scale and orientation, which ultimately lowered the entropy values. Interestingly, regions with similar entropy levels showed higher complexity values for branchlets than leaflets. This can be explained by the presence of a wider range of spatial scales in branchlets, whereas leaflets all tend to be small. The entropy-complexity spectrum analysis is a useful method not only to compare properties of astrocytic branchlets and leaflets, but also to characterize changes in astrocytic processes associated with different physiological and pathological states (Plata et al., 2018). In the future, this method of analysis can be extended to structures related to other cell types, such as dendrites and locations of synapses, and subcellular features, such as ER, mitochondria, vesicles.
Leaflets had skewed distributions of volume, surface area and SVR without obvious peaks. Thus, we did not classify the astrocytic leaflets into subclasses. However, we detected two distinguishable peaks when analyzing the distribution of distances to the nearest neighboring leaflets along an astrocytic branchlet. There was no correlation between these distance values and the size of the analyzed leaflet. The second peak roughly matched the diameter of dendritic spines. We speculate that leaflets grow toward spines to form astrocytic cradle around them; therefore, neighboring dendritic spines can space leaflets along the astrocytic branchlet in the hippocampal neuropil.
Dendritic spines are approached by axonal terminals to form a two-partite synapse. Astrocytic leaflets approach the synapses to form an astrocytic cradle (Verkhratsky and Nedergaard, 2014). Our findings show that the leaflet VF is larger around the spine surface than around axonal boutons. This result confirms previous report from 2D EM analysis, showing an asymmetric arrangement of PAPs around hippocampal pre- and post-synaptic terminals (Lehre and Rusakov, 2002). Interestingly, we could not detect a peak of astrocyte VF around the dendritic shaft. The reason why astrocytic leaflets do not extend toward the dendritic shaft is not entirely clear. Dendritic spines are located sparsely enough to form an “impenetrable wall” for the thin leaflets. However, the dendrite surface can be occupied by shaft synapses (e.g., GABAergic), limiting the space available for the leaflets (Megías et al., 2001). Nevertheless, two conclusions can be made from our results. First, in the hippocampal neuropil, astrocytic leaflets are specialized structures that form an astrocytic cradle around glutamatergic synapse. This cradle is asymmetric and has an opening on the presynaptic side. This structure might ensure glutamate and K+ clearance from the synapse and prevent the entry of extrasynaptic glutamate into the synaptic cleft (Lozovaya et al., 2004; Wu et al., 2012). Second, the low VF of astrocytic leaflets near dendritic shaft suggests less efficient neurotransmitter uptake around this compartment. The dendritic shaft contains extrasynaptic receptors and shaft synapses. This finding is consistent with lack of astrocytic cradle around shaft synapses, which may facilitate neurotransmitter spillover and activation of dendritic extrasynaptic receptors. Because many shaft synapses are GABAergic in CA1 pyramidal neurons, this may be a morphological mechanism promoting the activation of tonic GABAA conductances through lower astrocytic GABA uptake along the dendritic shaft (Kersante et al., 2013; Song et al., 2013). Unlike glutamate receptors, GABA receptors are not permeable to K+ (Wollmuth and Sobolevsky, 2004; Shih et al., 2013; Cheung et al., 2015). For this reason, GABAergic synapses may not require the presence of PAPs for K+ clearance.
Since astrocytic leaflets have no ER or mitochondria, they can only generate Ca2+ transients via Ca2+ entry through plasma membrane (Zorec et al., 2012; Bazargani and Attwell, 2016). Astrocytic branchlets have these intracellular organelles and can therefore generate/amplify Ca2+ signals via Ca2+ release from internal stores. Several reports have suggested that this is one the major signaling pathways to ensure neuron-astrocyte communication (Ullah et al., 2006; Agarwal et al., 2017; Sherwood et al., 2017). However, our morphological analysis indicates that less than a quarter of all synapses contribute to this form of signaling in hippocampal CA1 str.radiatum. This does not rule out that remote synapses can still modulate Ca2+ dynamics in astrocytic branchlets via glutamate spillover. Glutamate escaping synapses can reach high-affinity mGluRs on an astrocytic branchlet and trigger IP3 production. IP3 can also be generated in leaflets in response to local synaptic activity and then diffuse intracellularly to the astrocytic branchlet. Likewise, Ca2+ entering leaflets through plasma membrane can diffuse intracellularly to the parent branchlet. This scheme portrays the astrocytic branchlet with leaflets similarly to the dendritic shaft with dendritic spines. Dendritic spines receive synaptic signals (electric and Ca2+) which are integrated and amplified (dendritic spike) in the parent dendritic shaft (Larkum and Nevian, 2008; Spruston, 2008; Takahashi et al., 2012). Leaflets receive synaptic signals which are integrated and amplified in the parent astrocytic branchlet. This hypothesis, however, requires further experimental testing and mathematical modeling. Similar methods used for dendritic integration studies (e.g., local glutamate uncaging) can be employed to study the integration of Ca2+ signals in the astrocytic branchlet.
Most of modern cellular neuroscience methods have their caveats: Whole-cell recordings dialyze the cytoplasm preventing some forms of signaling and plasticity. Fluorescent Ca2+ sensors work as buffers and, therefore, strongly affect Ca2+ dynamics. Voltage-sensitive dyes affect the cell membrane properties, hence, the excitability of the cells. Super-resolution optical imaging requires intense laser illumination which increases local temperature and may affect metabolic processes. Mathematical models are based on numerous assumptions and simplifications. Overall, validity of animal research for understating human brain function is often questioned. Chemical fixation used in this study considered as a gold standard in the field for decades. Nevertheless, it affects the size and the structure of the extracellular space and distorts the morphological relationship between spines and astrocytic processes (Kinney et al., 2013; Korogod et al., 2015; Nicholson and Hrabetová, 2017). The use of electron microscopy to study extracellular space therefore has significant limitations. Thus, a biophysical approach based on diffusion should be used as a method of choice (Nicholson and Hrabetová, 2017). Although distortions in the neuropil happen because of chemical fixation, they likely affect the microenvironment around small and large spines to the same extent.
Finally, our conclusions are only applicable to the astrocytic coverage of the glutamatergic synapses in CA1 hippocampal neuropil. The scenario could differ in other brain regions: astrocytes cover 74% of the cerebellar Purkinje cell spines and 29% of the dendritic spines in the mouse visual cortex (Spacek, 1985; Ventura and Harris, 1999). In addition, astrocytic coverage of GABAergic synapses was not addressed in this work. However, recent reports highlighted important differences in the regulation of synaptic and extrasynaptic GABAergic signaling by neuronal and astrocytic GABA transporters (Kersante et al., 2013; Song et al., 2013). Therefore, the astrocytic coverage of GABAergic synapses also requires further morphological analysis.
All animal experiments have been performed with approval of Nizhny Novgorod State University ethics committee.
NG and IG performed the 3D reconstruction; NG, ASc, and AB analyzed the data and prepared the figures; ASc performed the diffusion simulations; VT discussed the results and consulted on the analysis; ASe, ASc, and AB wrote the paper; ASe supervised the project, planned the analysis, prepared the figures.
The research was supported by the Russian Science Foundation grant to ASe (project No. 16-14-00201, all parts of the paper except for diffusion simulations) and SUNY Albany and the NIH to ASc (NS102822, diffusion simulations part).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer, CC and the handling Editor declared their shared affiliation.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncel.2018.00248/full#supplementary-material
Agarwal, A., Wu, P. H., Hughes, E. G., Fukaya, M., Tischfield, M. A., Langseth, A. J., et al. (2017). Transient opening of the mitochondrial permeability transition pore induces microdomain calcium transients in astrocyte processes. Neuron 93, 587–605.e7. doi: 10.1016/j.neuron.2016.12.034
Bergles, D. E., Tzingounis, A. V., and Jahr, C. E. (2002). Comparison of coupled and uncoupled currents during glutamate uptake by GLT-1 transporters. J. Neurosci. 22, 10153–10162. doi: 10.1523/JNEUROSCI.22-23-10153.2002
Bernardinelli, Y., Randall, J., Janett, E., Nikonenko, I., Konig, S., Jones, E. V., et al. (2014b). Activity-dependent structural plasticity of perisynaptic astrocytic domains promotes excitatory synapse stability. Curr. Biol. 24, 1679–1688. doi: 10.1016/j.cub.2014.06.025
Calì, C., Baghabra, J., Boges, D. J., Holst, G. R., Kreshuk, A., Hamprecht, F. A., et al. (2016). Three-dimensional immersive virtual reality for studying cellular compartments in 3D models from EM preparations of neural tissues. J. Comp. Neurol. 524, 23–38. doi: 10.1002/cne.23852
Chao, T. I., Rickmann, M., and Wolff, J. R. (2002). “The synapse-astrocyte boundary: an anatomical basis for an integrative role of glia in synaptic transmission,” in The Tripartite Synapse: Glia in Synaptic Transmission, eds. A. Volterra, P. J. Magistretti, and P. G. Haydon (Oxford: Oxford Univ. Press), 4–23.
Chever, O., Djukic, B., McCarthy, K. D., and Amzica, F. (2010). Implication of Kir4.1 channel in excess potassium clearance: an in vivo study on anesthetized glial-conditional Kir4.1 knock-out mice. J. Neurosci. 30, 15769–15777. doi: 10.1523/JNEUROSCI.2078-10.2010.
Contini, D., Price, S. D., and Art, J. J. (2017). Accumulation of K+ in the synaptic cleft modulates activity by influencing both vestibular hair cell and calyx afferent in the turtle. J. Physiol. 595, 777–803. doi: 10.1113/JP273060
De Paola, V., Holtmaat, A., Knott, G., Song, S., Wilbrecht, L., Caroni, P., et al. (2006). Cell type-specific structural plasticity of axonal branches and boutons in the adult neocortex. Neuron 49, 861–875. doi: 10.1016/j.neuron.2006.02.017
Diamond, J. S. (2001). Neuronal glutamate transporters limit activation of NMDA receptors by neurotransmitter spillover on CA1 pyramidal cells. J. Neurosci. 21, 8328–8338. doi: 10.1523/JNEUROSCI.21-21-08328.2001
Diamond, J. S. (2005). Deriving the glutamate clearance time course from transporter currents in CA1 hippocampal astrocytes: transmitter uptake gets faster during development. J. Neurosci. 25, 2906–2916. doi: 10.1523/JNEUROSCI.5125-04.2005
Genoud, C., Quairiaux, C., Steiner, P., Hirling, H., Welker, E., and Knott, G. W. (2006). Plasticity of astrocytic coverage and glutamate transporter expression in adult mouse cortex. PLoS Biol. 4:e343. doi: 10.1371/journal.pbio.0040343
Gordleeva, S. Y., Stasenko, S. V., Semyanov, A. V., Dityatev, A. E., and Kazantsev, V. B. (2012). Bi-directional astrocytic regulation of neuronal activity within a network. Front. Comput. Neurosci. 6:92. doi: 10.3389/fncom.2012.00092
Grillo, F. W., Song, S., Teles-Grilo Ruivo, L. M., Huang, L., Gao, G., Knott, G. W., et al. (2013). Increased axonal bouton dynamics in the aging mouse cortex. Proc. Natl. Acad. Sci. U.S.A. 110, E1514–E1523. doi: 10.1073/pnas.1218731110
Hara, Y., Park, C. S., Janssen, W. G., Punsoni, M., Rapp, P. R., and Morrison, J. H. (2011). Synaptic characteristics of dentate gyrus axonal boutons and their relationships with aging, menopause, and memory in female rhesus monkeys. J. Neurosci. 31, 7737–7744. doi: 10.1523/JNEUROSCI.0822-11.2011
Héja, L., Nyitrai, G., Kekesi, O., Dobolyi, A., Szabo, P., Fiath, R., et al. (2012). Astrocytes convert network excitation to tonic inhibition of neurons. BMC Biol. 10:26. doi: 10.1186/1741-7007-10-26
Huang, H., and Bordey, A. (2004). Glial glutamate transporters limit spillover activation of presynaptic NMDA receptors and influence synaptic inhibition of Purkinje neurons. J. Neurosci. 24, 5659–5669. doi: 10.1523/JNEUROSCI.1338-04.2004
Jones, T. A., and Greenough, W. T. (1996). Ultrastructural evidence for increased contact between astrocytes and synapses in rats reared in a complex environment. Neurobiol. Learn. Mem. 65, 48–56. doi: 10.1006/nlme.1996.0005
Kersante, F., Rowley, S. C. S., Pavlov, I., Gutierrez-Mecinas, M., Semyanov, A., Reul, J. M., et al. (2013). A functional role for both gamma-aminobutyric acid (GABA) transporter-1 and GABA transporter-3 in the modulation of extracellular GABA and GABAergic tonic conductances in the rat hippocampus. J. Physiol. 591, 2429–2441. doi: 10.1113/jphysiol.2012.246298.
Kettenmann, H., Faissner, A., and Trotter, J. (1996). “Neuron-glia interactions in homeostasis and degeneration,” in Comprehensive Human Physiology: From Cellular Mechanisms to Integration, eds. R. Greger and U. Windhorst. (Berlin; Heidelberg: Springer Berlin Heidelberg), 533–543.
Kinney, J. P., Spacek, J., Bartol, T. M., Bajaj, C. L., Harris, K. M., and Sejnowski, T. J. (2013). Extracellular sheets and tunnels modulate glutamate diffusion in hippocampal neuropil. J. Comp. Neurol. 521, 448–464. doi: 10.1002/cne.23181
Kullmann, D., and Semyanov, A. (2002). Glutamatergic modulation of GABAergic signaling among hippocampal interneurons: novel mechanisms regulating hippocampal excitability. Epilepsia 43, 174–178. doi: 10.1046/j.1528-1157.43.s.5.12.x
Lang, C., Barco, A., Zablow, L., Kandel, E. R., Siegelbaum, S. A., and Zakharenko, S. S. (2004). Transient expansion of synaptically connected dendritic spines upon induction of hippocampal long-term potentiation. Proc. Natl. Acad. Sci. U.S.A. 101, 16665–16670. doi: 10.1073/pnas.0407581101
Lebedeva, A., Plata, A., Nosova, O., Tyurikova, O., and Semyanov, A. (2018). Activity-dependent changes in transporter and potassium currents in hippocampal astrocytes. Brain Res. Bull. 136, 37–43. doi: 10.1016/j.brainresbull.2017.08.015
Lebedeva, A. V., Dembitskaya, Y. V., Pimashkin, A. S., Zhuravleva, Z. D., Shishkova, E. A., and Semyanov, A. V. (2015). The role of energy substrates in astrocyte calcium activity of rat hippocampus in early postnatal ontogenesis. Sovrem. Tehnol. V Med. 7, 14–18. doi: 10.17691/stm2015.7.3.02
Lehre, K. P., and Danbolt, N. C. (1998). The number of glutamate transporter subtype molecules at glutamatergic synapses: chemical and stereological quantification in young adult rat brain. J. Neurosci. 18, 8751–8757. doi: 10.1523/JNEUROSCI.18-21-08751.1998
Lester, R. A., Clements, J. D., Westbrook, G. L., and Jahr, C. E. (1990). Channel kinetics determine the time course of NMDA receptor-mediated synaptic currents. Nature 346, 565–567. doi: 10.1038/346565a0.
Li, Y. X., and Rinzel, J. (1994). Equations for InsP3 receptor-mediated [Ca2+]i oscillations derived from a detailed kinetic model: a Hodgkin-Huxley like formalism. J. Theor. Biol. 166, 461–473. doi: 10.1006/jtbi.1994.1041
Liu, Q. S., Xu, Q., Arcuino, G., Kang, J., and Nedergaard, M. (2004a). Astrocyte-mediated activation of neuronal kainate receptors. Proc. Natl. Acad. Sci. U.S.A. 101, 3172–3177. doi: 10.1073/pnas.0306731101
Liu, Q. S., Xu, Q., Kang, J., and Nedergaard, M. (2004b). Astrocyte activation of presynaptic metabotropic glutamate receptors modulates hippocampal inhibitory synaptic transmission. Neuron Glia Biol. 1, 307–316. doi: 10.1017/S1740925X05000190
Lozovaya, N., Melnik, S., Tsintsadze, T., Grebenyuk, S., Kirichok, Y., and Krishtal, O. (2004). Protective cap over CA1 synapses: extrasynaptic glutamate does not reach the postsynaptic density. Brain Res. 1011, 195–205. doi: 10.1016/j.brainres.2004.03.023
Lushnikova, I., Skibo, G., Muller, D., and Nikonenko, I. (2009). Synaptic potentiation induces increased glial coverage of excitatory synapses in CA1 hippocampus. Hippocampus 19, 753–762. doi: 10.1002/hipo.20551
Martin, M. T., Plastino, A., and Rosso, O. A. (2006). Generalized statistical complexity measures: geometrical and analytical properties. Phys. A Stat. Mech. Appl. 369, 439–462. doi: 10.1016/j.physa.2005.11.053
Matsui, T., Omuro, H., Liu, Y. F., Soya, M., Shima, T., McEwen, B. S., et al. (2017). Astrocytic glycogen-derived lactate fuels the brain during exhaustive exercise to maintain endurance capacity. Proc. Natl. Acad. Sci. U.S.A. 114, 6358–6363. doi: 10.1073/pnas.1702739114
Medvedev, N., Popov, V., Henneberger, C., Kraev, I., Rusakov, D. A., and Stewart, M. G. (2014). Glia selectively approach synapses on thin dendritic spines. Philos. Trans. R. Soc. Lond. B Biol. Sci. 369:20140047. doi: 10.1098/rstb.2014.0047
Megías, M., Emri, Z., Freund, T. F., and Gulyas, A. I. (2001). Total number and distribution of inhibitory and excitatory synapses on hippocampal CA1 pyramidal cells. Neuroscience 102, 527–540. doi: 10.1016/S0306-4522(00)00496-6
Nägerl, U. V., Eberhorn, N., Cambridge, S. B., and Bonhoeffer, T. (2004). Bidirectional activity-dependent morphological plasticity in hippocampal neurons. Neuron 44, 759–767. doi: 10.1016/j.neuron.2004.11.016
Nielsen, T. A., DiGregorio, D. A., and Silver, R. A. (2004). Modulation of glutamate mobility reveals the mechanism underlying slow-rising AMPAR EPSCs and the diffusion coefficient in the synaptic cleft. Neuron 42, 757–771. doi: 10.1016/j.neuron.2004.04.003
Ostroff, L. E., Manzur, M. K., Cain, C. K., and Ledoux, J. E. (2014). Synapses lacking astrocyte appear in the amygdala during consolidation of Pavlovian threat conditioning. J. Comp. Neurol. 522, 2152–2163. doi: 10.1002/cne.23523
Panatier, A., Theodosis, D. T., Mothet, J. P., Touquet, B., Pollegioni, L., Poulain, D. A., et al. (2006). Glia-derived D-serine controls NMDA receptor activity and synaptic memory. Cell 125, 775–784. doi: 10.1016/j.cell.2006.02.051
Patneau, D. K., and Mayer, M. L. (1990). Structure-activity relationships for amino acid transmitter candidates acting at N-methyl-D-aspartate and quisqualate receptors. J. Neurosci. 10, 2385–2399. doi: 10.1523/JNEUROSCI.10-07-02385.1990
Patrushev, I., Gavrilov, N., Turlapov, V., and Semyanov, A. (2013). Subcellular location of astrocytic calcium stores favors extrasynaptic neuron-astrocyte communication. Cell Calcium 54, 343–349. doi: 10.1016/j.ceca.2013.08.003
Plata, A., Lebedeva, A., Denisov, P., Nosova, O., Postnikova, T. Y., Pimashkin, A., et al. (2018). Astrocytic atrophy following status epilepticus parallels reduced Ca2+ activity and impaired synaptic plasticity in the rat hippocampus. Front. Mol. Neurosci. 11:215. doi: 10.3389/fnmol.2018.00215
Poolos, N. P., Mauk, M. D., and Kocsis, J. D. (1987). Activity-evoked increases in extracellular potassium modulate presynaptic excitability in the CA1 region of the hippocampus. J. Neurophysiol. 58, 404–416. doi: 10.1152/jn.19126.96.36.1994
Reyes, R. C., Verkhratsky, A., and Parpura, V. (2012). Plasmalemmal Na+/Ca2+ exchanger modulates Ca2+-dependent exocytotic release of glutamate from rat cortical astrocytes. ASN Neuro 4:e00075. doi: 10.1042/AN20110059
Ribeiro, H. V., Zunino, L., Lenzi, E. K., Santoro, P. A., and Mendes, R. S. (2012). Complexity-entropy causality plane as a complexity measure for two-dimensional patterns. PLoS ONE 7:e40689. doi: 10.1371/journal.pone.0040689
Rose, C. R., Felix, L., Zeug, A., Dietrich, D., Reiner, A., and Henneberger, C. (2018). Astroglial glutamate signaling and uptake in the hippocampus. Front. Mol. Neurosci. 10:451. doi: 10.3389/fnmol.2017.00451
Rusakov, D. A., and Kullmann, D. M. (1998). Extrasynaptic glutamate diffusion in the hippocampus: ultrastructural constraints, uptake, and receptor activation. J. Neurosci. 18, 3158–3170. doi: 10.1523/JNEUROSCI.18-09-03158.1998
Scimemi, A., Tian, H., and Diamond, J. S. (2009). Neuronal transporters regulate glutamate clearance, NMDA receptor activation, and synaptic plasticity in the hippocampus. J. Neurosci. 29, 14581–14595. doi: 10.1523/JNEUROSCI.4845-09.2009
Semyanov, A., Walker, M., Kullmann, D., and Silver, R. (2004). Tonically active GABA(A) receptors: modulating gain and maintaining the tone. Trends Neurosci. 27, 262–269. doi: 10.1016/j.tins.2004.03.005
Sherwood, M. W., Arizono, M., Hisatsune, C., Bannai, H., Ebisui, E., Sherwood, J. L., et al. (2017). Astrocytic IP3 Rs: contribution to Ca2+ signalling and hippocampal LTP. Glia 65, 502–513. doi: 10.1002/glia.23107
Shih, P.-Y., Savtchenko, L. P., Kamasawa, N., Dembitskaya, Y., McHugh, T. J., Rusakov, D. A., et al. (2013). Retrograde synaptic signaling mediated by K+ efflux through postsynaptic NMDA receptors. Cell Rep. 5, 941–951. doi: 10.1016/j.celrep.2013.10.026
Sibille, J., Pannasch, U., and Rouach, N. (2014). Astroglial potassium clearance contributes to short-term plasticity of synaptically evoked currents at the tripartite synapse. J Physiol 592, 87–102. doi: 10.1113/jphysiol.2013.261735
Song, I., Volynski, K., Brenner, T., Ushkaryov, Y., Walker, M., and Semyanov, A. (2013). Different transporter systems regulate extracellular GABA from vesicular and non-vesicular sources. Front. Cell. Neurosci. 7 23. doi: 10.3389/fncel.2013.00023
Tanaka, M., Shih, P. Y., Gomi, H., Yoshida, T., Nakai, J., Ando, R., et al. (2013). Astrocytic Ca2+ signals are required for the functional integrity of tripartite synapses. Mol. Brain 6:6. doi: 10.1186/1756-6606-6-6
Tang, F., Lane, S., Korsak, A., Paton, J. F., Gourine, A. V., Kasparov, S., et al. (2014). Lactate-mediated glia-neuronal signalling in the mammalian brain. Nat. Commun. 5:3284. doi: 10.1038/ncomms4284
Theodosis, D. T., Poulain, D. A., and Oliet, S. H. (2008). Activity-dependent structural and functional plasticity of astrocyte-neuron interactions. Physiol. Rev. 88, 983–1008. doi: 10.1152/physrev.00036.2007
Ullah, G., Jung, P., and Cornell-Bell, A. H. (2006). Anti-phase calcium oscillations in astrocytes via inositol (1, 4, 5)-trisphosphate regeneration. Cell Calcium 39, 197–208. doi: 10.1016/j.ceca.2005.10.009
Wu, Y. W., Grebenyuk, S., McHugh, T. J., Rusakov, D. A., and Semyanov, A. (2012). Backpropagating action potentials enable detection of extrasynaptic glutamate by NMDA receptors. Cell Rep. 1, 495–505. doi: 10.1016/j.celrep.2012.03.007
Zorec, R., Araque, A., Carmignoto, G., Haydon, P. G., Verkhratsky, A., and Parpura, V. (2012). Astroglial excitability and gliotransmission: an appraisal of Ca2+ as a signalling route. ASN Neuro 4:61. doi: 10.1042/AN20110061
Keywords: astrocyte, branches, leaflets, dendritic shaft, axonal bouton, spatial entropy, spatial complexity
Citation: Gavrilov N, Golyagina I, Brazhe A, Scimemi A, Turlapov V and Semyanov A (2018) Astrocytic Coverage of Dendritic Spines, Dendritic Shafts, and Axonal Boutons in Hippocampal Neuropil. Front. Cell. Neurosci. 12:248. doi: 10.3389/fncel.2018.00248
Received: 28 February 2018; Accepted: 19 July 2018;
Published: 17 August 2018.
Edited by:Pierre J. Magistretti, King Abdullah University of Science and Technology (KAUST), Saudi Arabia
Reviewed by:Corrado Calì, King Abdullah University of Science and Technology, Saudi Arabia
Alexej Verkhratsky, University of Manchester, United Kingdom
Graham William Knott, École Polytechnique Fédérale de Lausanne, Switzerland
Copyright © 2018 Gavrilov, Golyagina, Brazhe, Scimemi, Turlapov and Semyanov. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alexey Semyanov, email@example.com