A Distance-Dependent Distribution of Presynaptic Boutons Tunes Frequency-Dependent Dendritic Integration.

How presynaptic inputs and neurotransmitter release dynamics are distributed along a dendritic tree is not well established. Here, we show that presynaptic boutons that form onto basal dendrites of CA1 pyramidal neurons display a decrease in active zone (AZ) size with distance from the soma, resulting in a distance-dependent increase in short-term facilitation. Our findings suggest that the spatial distribution of short-term facilitation serves to compensate for the electrotonic attenuation of subthreshold distal inputs during repeated stimulation and fine-tunes the preferred input frequency of dendritic domains.

Introduction

Pyramidal neurons receive thousands of excitatory inputs on their extensive dendritic arbors. As a consequence, neurons need strategies to balance the strength of their synaptic inputs so that signals arriving at distal synapses, far from the soma, have a meaningful contribution to neuronal output. Active dendritic integration can greatly boost synaptic signals by amplifying local depolarizations onto dendritic spikes that travel to the soma. Much has been done to describe this phenomenon, and its regulation, by focusing on postsynaptic compartments. In particular, previous studies have characterized the distribution of the size and strength of postsynaptic spines within dendritic domains and have established their contributions to dendritic integration (Branco and Häusser, 2010, Magee, 2000, Spruston, 2008). Overall, these findings showed that whereas apical dendrites show an increase in the strength of synaptic inputs with distance from the soma (Magee and Cook, 2000), the thinner basal and apical oblique dendrites show a decrease in spine size with distance (Katz et al., 2009, Menon et al., 2013, Walker et al., 2017). What emerges is a complex picture of synapse distribution along dendritic arbors that appear to favor either local (dendritic) or global (cell-wide) integration, depending on dendrite identity. On the other hand, we know much less about the structural and functional distribution of presynaptic boutons (de Jong et al., 2012) and the role they play in shaping the integration of synaptic inputs on dendrites and their subdomains (Chabrol et al., 2015). Critically, whereas postsynaptic strength establishes the amount of local dendritic depolarization, changes in presynaptic structure and function determine the dynamics of neurotransmitter release, through short-term forms of plasticity. This feature, in turn, governs the type of information that is transmitted across a synapse (Branco and Staras, 2009, Fioravante and Regehr, 2011). To understand how dendrites integrate synaptic inputs, it is therefore important to first uncover how presynaptic boutons and the dynamics of neurotransmitter release are distributed along dendrites. A non-random distribution that echoes that of postsynaptic spines could lead to specialized dendritic domains that code for specific streams of information. Thus, short-term forms of plasticity that govern the release of neurotransmitter across different presynaptic boutons may have profound effects on how information is processed in postsynaptic dendrites (Abrahamsson et al., 2012). More globally, short-term plasticity (STP) has been shown to play many important roles in vivo (Regehr, 2012) and, in excitatory synapses of the hippocampus, STP is thought to contribute to the transmission of information about an animal’s place field (Kandaswamy et al., 2010, Klyachko and Stevens, 2006). Understanding the dynamics of neurotransmitter release in hippocampal presynaptic boutons will therefore also have implications for how information is encoded by CA1 neurons during active behaviors. Here, we show that presynaptic boutons decrease in size along the basal dendrites of CA1 hippocampal neurons, resulting in a decrease in release probability and an increase in short-term facilitation with distance from the soma. We reveal that this spatial distribution in STP tunes dendritic domains to specific information frequencies, introducing a further level of specialization to dendritic computations.

Results

We set out to map the distribution of the structure and function of presynaptic terminals in CA1 stratum oriens, a region of the hippocampus that receives ordered axonal inputs that are mostly perpendicular to the basal dendrites of pyramidal neurons (Andersen et al., 1980). We first performed Serial Block-Face Scanning Electron Microscopy (SBFSEM) of three regions of CA1 stratum oriens from a P22 brain, ranging from deep (close to the stratum pyramidale) to superficial (close to alveus) areas and reconstructed basal dendrites together with their synaptic inputs (Figures 1A–1C). Although dendritic segments showed a wide distribution in the size of excitatory inputs, we found, as expected, a strong correlation between morphological measures of presynaptic and postsynaptic compartments (Figure S1) (Holderith et al., 2012, Schikorski and Stevens, 1997). In agreement with the notion that basal dendrites taper toward tip ends (Menon et al., 2013), we also saw a decrease in dendrite diameter along the stratum oriens with increasing distance away from stratum pyramidale (Figure 1D; n = 35, mean = prox 0.75 ± 0.03 μm, med 0.65 ± 0.02 μm, dist 0.61 ± 0.01 μm, p < 0.001 Kruskal-Wallis test, prox-dist adjusted p < 0.001 Dunn’s multiple comparison test). Importantly, of the many structural features measured at the synapse (Figure S1), AZ area, a good predictor of release probability (Holderith et al., 2012), showed a strong decrease in size with distance along stratum oriens (Figure 1E; n = 604, mean = prox 0.074 ± 0.003 μm2, med 0.058 ± 0.002 μm2, dist 0.061 ± 0.003 μm2, p < 0.0001 Kruskal-Wallis test; prox-med adjusted p < 0.001, prox-dist adjusted p < 0.001 Dunn’s multiple comparison test). Dendrite diameter, which can be taken as an indirect measure of distance from the soma along a tapering dendrite, also correlated well with AZ size (Figures 1F and 1G), providing further evidence that presynaptic terminals become smaller with distance along a basal dendrite. A reconstruction of proximal and distal dendrites performed on an adult brain (P100) showed a very similar distribution of presynaptic and postsynaptic properties with distance (Figures 1H–1K, Figure S1), indicating that the decrease in AZ size (and other synaptic parameters) was not unique to adolescent (P22) brains. Our dataset has some limitations, including the inability to follow the same dendrite from a single cell and to establish the identity of the inputs that arrive at different locations. We therefore turned to dissociated hippocampal neurons, where synaptic inputs can be labeled and followed along a single dendrite and where input identity is jumbled following the dissociation procedure. We find that measures of both presynaptic structure (vGlut labeling) and function (FM4-64 labeling of recycling vesicles) also decreased in a graded manner with distance along a dendrite (Figure S2). Together, these data suggest that the distance-dependent distribution of presynaptic inputs may be independent of the identity of the axon and is likely constrained by the properties of the postsynaptic dendrite.

Figure 1

Presynaptic Active Zone Size Scales with Dendrite Diameter and Distance from Cell Bodies

(A) SBFSEM low-magnification image showing experimental design: 3 areas (red, proximal; purple, medial; blue, distal) at increasing distances from the pyramidal cell layer in the stratum oriens region of the CA1; S.O., stratum oriens; S.P., stratum pyramidale. Scale bar, 50 μm.

(B) SBFSEM single section images (left panels), with corresponding 3D reconstructions (right panels), showing presynaptic (yellow) and postsynaptic (green) structures with PSDs (red) and AZs (light blue) highlighted. Top left panel shows a spine head (star) connected through the narrower spine neck to the main dendritic shaft (asterisk). Scale bars, 0.5 μm.

(C) Two dendrites (thin dendrite left, thick dendrite right) reconstructed in 3D with spine heads in purple and boutons in green. Scale bar, 1 μm.

(D–G) Data from a postnatal day 22 animal. (D) Cumulative fraction plot: dendrites reconstructed in the proximal area have larger diameters than distal and medial area dendrites; n = 35 dendrites, p < 0.001 ANOVA; proximal-distal adjusted p < 0.001; proximal-medial adjusted p < 0.01, Tukey’s multiple comparison test. (E) Cumulative fraction plot of AZ sizes, which are a larger in the proximal group; n = 604 AZs, p < 0.0001 Kruskal-Wallis test; proximal-medial adjusted p < 0.0001, proximal-distal adjusted p = 0.0001 Dunn’s multiple comparison test. (F) Smaller AZs tend to be found on thinner dendritic processes, n = 556 AZs, Spearman’s correlation; colors indicate area in which the dendrites were reconstructed. (G) Average AZ size (per dendrite) positively correlates with dendrite diameter, n = 35 dendrites, Spearman’s correlation.

(H–K) Data from a postnatal day 100 animal. (H) Proximal dendrites are thicker than distal dendrites, cumulative fraction plot n = 26 dendrites, p < 0.05 unpaired t test with Welch’s correction. (I) Cumulative fraction plot of AZ sizes, which are a larger in the proximal group; n = 505 AZs, p < 0.0001 Kolmogorov-Smirnov test. (J) Smaller AZs tend to be found on thinner dendritic processes, n = 505 AZs, Spearman’s correlation; colors indicate area in which the dendrites were reconstructed. (K) Average AZ size (per dendrite) positively correlates with dendrite diameter, n = 26 dendrites, Spearman’s correlation. See also Figure S1. Data are represented as mean ± SEM.

AZ size is thought to correlate well with release probability (Pr) (Holderith et al., 2012, Murthy et al., 2001, Schikorski and Stevens, 1997). In addition, Pr is tightly coupled to the dynamics of neurotransmitter release during a burst of action potentials, such that high Pr synapses undergo short-term depression, whereas low Pr synapses exhibit short-term facilitation (Dobrunz and Stevens, 1997). This correlation between Pr and STP prompted us to measure neurotransmitter release dynamics along basal dendrites. We therefore patch clamped and filled CA1 pyramidal neurons with a fluorescent dye (Alexa 594) to visualize the basal dendritic arbor and measure synaptic inputs along its dendrites. Two stimulating pipettes were positioned along the basal dendrite, one at a site proximal to the soma and one distal to it (Figure 2A). In agreement with our structural findings (Figure 1), synaptic AMPA receptor currents measured by stimulating either proximal or distal axons independently showed that distal inputs facilitated more than proximal ones (Figures 2B and 2C). The paired-pulse ratio (PPR) of the first two stimuli delivered at 20 Hz was larger for distal inputs (mean = 1.59 ± 0.05) compared to proximal ones (mean = 1.33 ± 0.05) and remained higher throughout the 5 pulse stimulus train (Figures 2C and 2D). To establish whether the observed differences in STP were affected by postsynaptic properties, such as the inactivation or saturation state of AMPA receptors, we used cyclothiazide (CTZ, a drug that prevents AMPAR desensitization [Partin et al., 1993, Patneau et al., 1993]) or γDGG (a competitive antagonist of glutamate receptors [Liu et al., 1999]) to directly assess any contributions to STP by postsynaptic receptors (Figure 2E). Neither drug had any effect on the difference in STP between distal and proximal inputs (CTZ PPR = prox 1.18 ± 0.09, dist 1.51 ± 0.09, n = 9, p = 0.01; γDGG PPR = prox 1.36 ± 0.07, dist 1.64 ± 0.09, n = 13, p < 0.01; two-way ANOVA for control, CTZ, and γDGG, interaction p = 0.84), although there was a positive trend of γDGG on facilitation (Figure S3) that could reflect multivesicular release during the train (Christie and Jahr, 2006, Oertner et al., 2002, Wadiche and Jahr, 2001). To further corroborate that these effects were indeed presynaptic in origin, we performed similar measures of PPR in transgenic mice where synaptotagmin 7 (Syt 7) was knocked out. Syt 7−/− mice have been shown to lack any facilitation (Jackman et al., 2016) and should therefore abolish the differences in PPR measured here. We find that whereas WT mice show the same distance-dependent increase in PPR observed above, Syt7−/− littermates show neither facilitation nor any obvious difference in PPR between proximal and distal domains (Figure 2F). Together, our data point to a presynaptic origin in the difference in PPR along basal dendrites.

Figure 2

Increases in Short-Term Facilitation with Distance along a Dendrite Boost Distal Synaptic Integration

(A) Whole-cell patch-clamp technique was used to record synaptic currents and fill CA1 pyramidal cells with a fluorescent dye to image its structure. Two stimulating pipettes were placed in the proximal (red) and distal (blue) extracellular domains of the basal dendritic tree of pyramidal cells to stimulate the local fibers.

(B) Single cell example of average EPSC responses (average of 20 individual sweeps) to trains of 5 pulses at 20 Hz delivered to the proximal region (red trace) and the distal region (blue trace). The distal response shows greater facilitation compared to the proximal one.

(C) Normalized average peak EPSC amplitudes for distal and proximal responses show greater sustained facilitation during a 5 pulse train (20 Hz) for distal synapses, n = 35 cells, multiple t tests with p values adjusted with the Holm-Sidak method, p < 0.05.

(D) Paired-pulse ratios (PPRs) for each individual cell recorded at distal and proximal synapses. The majority (29/35) of cells display greater facilitation in the distal domain, n = 35 cells, p = 0.0003 two-tailed paired t test.

(E) Distal increase in PPR is not ascribable to postsynaptic AMPA receptor desensitization (prevented by CTZ application) or to AMPA receptor saturation (avoided with γDGG application). Distal PPR is greater than proximal PPR with CTZ (n = 9 cells, p = 0.01), and γDGG (n = 13 cells, p < 0.01), two-tailed paired t test. Two-way ANOVA to test PPRs in control (D), CTZ, and γDGG conditions together shows no significant interaction, p = 0.84, indicating that the drugs have no effect on STP properties.

(F) Full Synaptotagmin7 KO eliminates facilitation and proximo-distal STP differences. Triangles in lighter colors are from Syt7KO mice, n = 9 cells, circles in darker colors are littermate wild-type mice, n = 12 cells. For WT mice, proximal facilitation is lower than distal, multiple t tests, p < 0.05. WT facilitation is greater than Syt7KO facilitation, p < 0.01 multiple t tests.

(G–I) Distal synaptic EPSCs take longer to reach the soma. (G) Left panel: normalized trace for a proximal and distal EPSC response showing the delayed kinetics of the distal compared to the proximal synaptic current. Right panel: the rise time constant of the EPSCs was significantly higher in distally triggered events, n = 49 cells, p < 0.001 Wilcoxon signed rank test. (H) Longer rise times correlate with the amount of facilitation, Spearman’s correlation. (I) PPR is higher when the stimulation electrode is placed further away from the soma, measured as distance along the dendrite, Spearman’s correlation.

(J) Proximal synapses display greater Pr than distal synapses. After MK-801 bath application, the normalized amplitude of EPSCs from proximally stimulated synapses decay faster, following successive stimulations, than distal ones; n = 8 cells. Data points were fit with a double exponential function (filled lines). Insets are example traces of 7 successive NMDA mediated EPSCs.

(K) Frequency tuning curve showing PPRs for all frequencies tested. Distal PPRs (second stimulus only) increase significantly in the 20 Hz (n = 35) same as (D), and 50 Hz range (n = 20), p < 0.05, multiple t tests with Holm-Sidak adjusted p values. 5 Hz (n = 13 cells), 10 Hz (n = 16), 80 Hz (n = 9).

(L–M) Short-term facilitation contributes to dendritic non-linear events in distal domains. (L) Current-clamp example traces (red proximal, blue distal stimulation) in response to a paired pulse, of increasing stimulus intensity (lighter color shades represent lower intensity). (M) The proportion of supra-linear events (at least 2 mV above the expected response) for distal synapses is greatly increased for the second pulse (P2) following facilitation at 20 Hz (0 events in P1, 8 events in P2, n = 10), while supra-linear events were detected to the first pulse (P1) for proximal stimulations (n = 10, 4 in P1, 5 in P2). At 5 Hz, where STP is absent, distal synapses had fewer supra-linear events (n = 5, 0 in P1, 2 in P2). See also Figures S2 and S3. Data are represented as mean ± SEM.

To confirm the location specificity of our stimulus pipettes, we measured the kinetics of synaptic transmission at the soma. In agreement with the electrotonic decay of signals along a dendrite (Mainen et al., 1996, Rall, 1962, Spruston et al., 1994), we confirmed that distal stimuli (mean = 1.66 ± 0.07 ms) elicited slower events than proximal stimuli (mean = 1.23 ± 0.09 ms; p < 0.001, Wilcoxon sign rank test; Figure 2G). Since both the rise time and time to peak of events are directly related to the distance from the recording pipette at the soma, we looked at their relationship with PPR across different cells. A positive correlation emerged between synaptic kinetics and PPR (Figure 2H and Figures S3F and S3G, n = 49), suggesting a gradual change in PPR along a dendrite. Indeed, a positive correlation was also observed when PPR was plotted against the absolute distance of the stimulating electrode from the soma (Figure 2I). Our structural observations (Figure 1) predict that distal synapses will have a lower release probability (Pr) than proximal ones (Holderith et al., 2012), which may, in turn, help account for the increased levels of distal facilitation. To test for this, we used the irreversible open-channel blocker of NMDA receptors (NMDARs), MK-801, to measure Pr at proximal and distal compartments (Figure 2J). We recorded NMDAR currents in the presence of MK-801 in response to successive stimuli to measure the gradual block of NMDAR channels. The decay was best fit by a double exponential function, with a fast and a slow phase, indicating at least two groups of synapses with high and low release probabilities, respectively (Hessler et al., 1993). The two time constants were very similar for proximal and distal dendrites (proximal, τ1: 1.44 and τ2: 14.64; distal, τ1: 1.37 and τ2: 16.14), but in distal dendrites a larger fraction of the decay was explained by the slower time constant (proximal: 49%; distal: 69%). The overall slower decay observed in distal dendrites (mean half-life, proximal: 7.28 ± 1.2, distal: 13.99 ± 3.03, n = 8, p = 0.027, paired t test) indicates that synaptic NMDARs experienced less neurotransmitter in response to successive stimuli and therefore belong to synapses with, on average, a lower Pr than those found in proximal dendrites (Figure 2J). These findings mirror the structural correlation of presynaptic AZ size with dendrite diameter (Figures 1G and 1K) and strengthen our hypothesis that both the structure and function of presynaptic boutons display a distance-dependent distribution along basal dendrites. Furthermore, the lack of cross-inhibition by MK-801 between proximal and distal inputs is further support for the specificity of our stimulating electrodes in recruiting axons that form connections locally (Figure S3). Finally, we also established the frequency range over which inputs showed facilitation. We found that distal inputs showed significantly larger facilitation than proximal ones over a limited frequency range, between 10 and 50 Hz (Figures 2K and S3), with a peak at 20 Hz (PPR means = 5 Hz, prox 1.13 ± 0.04, dist 1.05 ± 0.04; 10 Hz, prox 1.18 ± 0.06, dist 1.35 ± 0.05; 20 Hz, see above; 50 Hz, prox 1.24 ± 0.07, p < 0.01 dist 1.53 ± 0.07, p < 0.01; 80 Hz, prox 1.17 ± 0.12, dist 1.22 ± 0.08). So far, our data show that distal inputs are tuned to band-pass frequencies in the near gamma frequency range, which curiously match the high-frequency discharges measured in vivo in hippocampal place cells when an animal passes through a place field (Huxter et al., 2003).

To establish whether the biased distribution in short-term facilitation can influence dendritic integration and neuronal output, we performed current-clamp experiments and locally stimulated either proximal or distal afferents, as above. We found that stimulation of distal inputs at frequencies that showed increased facilitation (20 Hz), also showed a tendency to produce non-linear summations to multiple stimuli across the range of stimulation intensities (Figures 2L and 2M). This effect was less pronounced in proximal dendrites, where non-linear events were only observed for high stimulus intensities. The difference in integration properties between proximal and distal events has been shown to depend on the opening of NMDARs (Ariav et al., 2003, Branco and Häusser, 2011, Major et al., 2008, Schiller et al., 2000), which are more likely to be activated at the high impedance distal dendrites, where levels of depolarization to a given input are higher. We found that distal inputs are better suited to respond to multiple stimuli than to a single stimulus. This was apparent when comparing responses along a paired pulse. The second stimulus of a 20 Hz pair was consistently more likely to summate in a supralinear manner, when compared either to the first stimulus or to stimuli delivered at 5 Hz (Figures 2L and 2M). Together, our data show that the integration properties of distal domains are markedly different from proximal ones and depend on input frequency. Although there are many factors that likely play a role in driving these differences in dendritic integration, STP is an obvious candidate. To explore this further, we turned to a computational model where we could directly assess the role played by STP on the integration properties of different dendritic domains.

We built a model consisting of a cluster of proximal and distal synapses (15 each; Figure 3A), where STP was adjusted to match our data (Figure 3B). This simple model aims to first determine the impact that a distance-dependent STP gradient can have on dendritic integration, when all other variables are the same. Activating either proximal or distal synapses separately with Poisson input trains delivered at different mean frequencies elicited non-linear responses (Figure 3C). Distal inputs, however, showed larger levels of membrane depolarization compared to proximal ones, over a large frequency range (continuous lines in Figure 3D). This distal amplification of synaptic inputs was partially lost when distal synapses were switched to proximal STP properties, suggesting that increased levels of facilitation enhance supralinear integration in distal dendrites. Interestingly, the difference in STP along a dendrite not only played a role in modulating dendritic non-linearities, it also contributed to distance-dependent input normalization; removal of active conductances from the model showed similar levels of membrane depolarization for distal and proximal inputs, which was lost when distal synapses were tuned to proximal STP properties (dashed lines in Figure 3D; note that the red and blue dashed lines overlap). This result suggests that the gradient in STP can boost distal inputs to counteract the effects of passive decay. Finally, increased levels of facilitation also resulted in an increased variance of membrane responses over a limited frequency domain, given by the short-term dynamics of the synapse (Figure 3E). This model shows that the STP gradient can have a significant impact on synaptic input integration, mainly exploiting differences in the recruitment of dendritic non-linearities. As such, it will act alongside other mechanisms that impact on non-linear integration, such as synaptic conductances, impedance, or release probability, all of which can also display gradients. However, models that include additional gradients of these other variables showed that the input-output curve remains biased toward larger depolarizations at higher frequencies in distal dendrites (Figure S4), indicating that the magnitude of the STP gradient that we observe is sufficient to have a significant impact on dendrite integration across a variety of conditions. Together, these simulations support the idea that distal inputs have an increased capacity for information transfer that is boosted by the short-term dynamics of neurotransmitter release from presynaptic inputs.

Figure 3

PPR Gradient Enhances Distal Supra-linear Integration

(A) Compartmental model illustration with distal and proximal synapses along one dendritic branch (left) and respective single synapse AMPA conductance traces for 20 Hz stimulation (right).

(B) Model data showing PPR dynamics for each pulse in a train of 20 Hz (left) and for PPR between 5 and 100 Hz for the first two pulses in the train.

(C) Example simulation traces for Poisson input trains delivered independently at each synapse at two different rates. Top raster shows input times for all synapses, and traces below show responses to the same input for proximal (red) and distal (blue) synapses, and for distal synapses equipped with PPR identical to proximal synapses (purple).

(D and E) (D) Mean depolarization during stimulation for different input frequencies and (E) respective standard deviation of the membrane potential, showing that increased PPF at distal synapses enhances supra-linear integration. Dashed lines in (D) and (E) are for passive models for each position (red proximal, blue distal, and purple distal with proximal PPR properties); note that in (D) the red and blue dashed lines overlap. See also Figure S4. Data are represented as mean ± SEM.

Discussion

Our findings uncovered a gradient in the distribution of presynaptic terminals along dendrites that dictates the short-term dynamics of synaptic transmission. This spatial gradient in short-term facilitation serves to both normalize the distance-dependent decay in the amplitude of subthreshold inputs to repeated stimuli, as well as tune the preferred input frequency of different dendritic domains through supralinear integration. As a result, the type of information transmitted across a synapse will depend on its location along a dendrite, an arrangement that could be exploited by different input streams to achieve input-specific differential integration on the same dendrite.

Previous studies have shown that pre- and post-synaptic compartments are well matched, both structurally and functionally (Kay et al., 2011, Murthy et al., 2001). Our findings of a spatial bias in the distribution of presynaptic boutons would predict that postsynaptic structures should follow suit. Indeed, work focusing on postsynaptic spines has shown that along the basal dendrites of CA1 pyramidal neurons, non-perforated spines, which represent the great majority of excitatory inputs, show a decrease in size (including PSD size) with dendritic distance (Menon et al., 2013, Walker et al., 2017). This reduction in size is thought to locally compensate for the increased impedance of thinner, distal dendrites, normalizing responses locally, rather than cell-wide (Katz et al., 2009). By including presynaptic boutons and neurotransmitter release dynamics, our findings uncover spatially segregated domains for the transfer of information within dendrites. However, our experiments provide no information on the identity of the axons that innervate different dendritic compartments along the stratum oriens. In general, it is thought that axons in this region of the hippocampus arrive mainly from pyramidal cells in areas CA3 and CA2, although other brain areas may be involved. There is, in fact, a biased topographic projection of CA3 pyramidal cells to the stratum oriens, such that CA3 cells whose soma lie closer to the dentate gyrus (DG) project preferentially to distal dendritic domains of CA1 neurons, whereas those CA3 cells further away from the DG project to more proximal dendritic domains (Ishizuka et al., 1990). It is therefore possible that the properties of distal and proximal boutons are dictated by the identity of the presynaptic neuron itself. However, our in vitro findings argue against this view. Using dissociated hippocampal neurons, we show that a similar distance-dependent decrease in presynaptic bouton properties also occurs in vitro (Figure S2) (de Jong et al., 2012) and suggests that the distribution of synapses observed in vivo may well be independent of the circuit or the identity of the synaptic input. It remains a likely possibility that the biased synaptic distribution along dendrites is therefore specified by the postsynaptic neuron in a cell-autonomous manner. Indeed, work in dissociated hippocampal neurons has shown that local dendritic depolarization is a major determinant of presynaptic release probability (Pr), such that increases in depolarization induce a homeostatic decrease in Pr (Branco et al., 2008). It is therefore tempting to speculate that the increase in impedance in distal dendrites, which would result in larger synaptic amplitudes (Spruston, 2008), would in turn act to reduce Pr locally. In this way, dendritic impedance could act as a readout of dendritic distance. Possible molecular mechanisms may include retrograde messengers, some of which have been shown to modulate presynaptic function in response to postsynaptic membrane depolarization (Regehr et al., 2009).

Our data show that synapses are distributed in a distance-dependent manner along dendrites and we provide multiple lines of evidence to suggest that this distribution is also graded with distance. First, we find a graded correlation between synapse morphology and dendrite diameter (Figures 1G and 1K), a measure that has been previously shown to be a good proxy for distance along a given dendrite (Walker et al., 2017). Second, measures of PPR are correlated to the time course of the EPSC (Figure 2H), reflecting the distance-dependent filtering suffered by an EPSC as it travels to the soma. However, these measures of distance are indirect and incur a certain amount of noise that likely arises from the fact that dendrites are heterogenous entities, especially across different cells, with variable degrees of tapering and where the passive electrotonic decay of synaptic events can be influenced by other factors. Importantly, direct measures of dendritic distance in either acute slices (Figure 2I) or in primary neuronal cells (Figure S2) both show a correlation with presynaptic properties. Together, our data support the notion of a graded distribution of synaptic structure and function along dendrites.

Our study finds a strong correlation between presynaptic structure and function, where boutons with smaller AZs have a lower Pr and an increased PPR. However, although PPR has been shown to correlate with the overall Pr of a bouton (Dobrunz and Stevens, 1997), this correlation may well be driven by the release probability of individual vesicles (Pvr) along an active zone. Differences in Pvr can arise from a number of different scenarios, ranging from structural features (e.g., the relative distribution of vesicles and calcium channels) to molecular heterogeneity (e.g., the expression of specific proteins that influence exocytosis). Although we have not explored these possibilities in our study, they could also underlie some of the differences in PPR observed here.

One other recent study has also shown a graded distribution of short-term plasticity along dendrites (Abrahamsson et al., 2012), although both the mechanism and direction of the gradient were opposed to that shown here. Stellate cells, a type of interneuron in the cerebellum, show a decrease in short-term facilitation along their thin dendrites (Abrahamsson et al., 2012). However, the effect is purely postsynaptic and arises from the large local depolarization in distal dendritic domains that reduces the driving force of postsynaptic receptors, resulting in the sub-linear integration of clustered inputs. As a result, stellate cells become ideal integrators of decorrelated inputs in both space and time, preferring the arrival of sparsely distributed, asynchronous events. This is clearly distinct form the mechanisms preferentially used by basal dendrites in the hippocampus and cortex, where clustered inputs drive local supralinear dendritic events (Branco and Häusser, 2011, Losonczy and Magee, 2006, Makara and Magee, 2013). Indeed, our data show that in CA1 basal dendrites, facilitation is distributed in the opposite direction to stellate cells, increasing with distance from the soma. Furthermore, according to our experimental data and computational model, increased distal facilitation is needed to trigger supralinear dendritic integration as this requires fast concomitant activation of multiple clustered synapses (Branco and Häusser, 2011), the likelihood of which is boosted by presynaptic facilitation. Collectively, these studies underscore the importance of short-term forms of plasticity on dendritic integration.

Finally, the difference in STP in proximal versus distal dendrites will also have important consequences on the way each domain encodes information. Proximal synapses that facilitate less will be better suited to respond to more isolated or lower-frequency inputs, suggesting that they respond better to temporally decorrelated events. Distal dendrites, on the other hand, with their higher levels of facilitation will act as spatio-temporal filters that favor high-frequency, clustered inputs. These highly selective distal inputs may, therefore, carry salient information for working memory or place field location. Facilitation of synaptic transmission in excitatory hippocampal synapses has been proposed to act as an adaptive high-pass filter that transmits and amplifies the signals encoding place-field information (Kandaswamy et al., 2010, Klyachko and Stevens, 2006). Our data suggest that the distal domains of CA1 basal dendrites are better suited to perform this computation. Future work will need to establish the identity of the inputs that arrive along different dendritic locations and understand the type of information they encode.

STAR★Methods

Key Resources Table

Contact for Reagent and Resource Sharing

As Lead Contact, Juan Burrone is responsible for all reagent and resource requests. Please contact Juan Burrone at juan.burrone@kcl.ac.uk with requests and inquiries.

Method Details

Animals

All animal procedures were approved by the local ethics committee and licensed under the UK Animals (Scientific Procedures) Act of 1986. Male and Female SV-129 mice were housed grouped in standard cages and provided with ad libitum food and water. The Syt 7 knockout mice were obtained from The Jackson Laboratory (Chakrabarti et al., 2003). Sprague-Dawley rats were obtained from Charles River Laboratory.

Dissociated hippocampal cultures

Primary hippocampal cultures were prepared from embryonic day 18 Sprague-Dawley rats (Charles River Laboratory). Dissociated cells were plated onto 18 mm diameter coverslips (Menzel Gläser, Germany) pre-treated with poly-D-lysine (50 μg/mL) and laminin (20 μg/mL) at a density of 350 cells/mm2 in Neurobasal media containing 1% fetal calf serum, 1% B-27 supplement, 0.5% glutamax and 0.5% penicillin/streptomycin. Neurons were kept for 17-21 days in vitro.

Electron Microscopy

Two mice (post-natal day 22 and 100) were transcardially perfused with 20 mL of ice-cold saline solution followed by 200 mL of ice-cold fixative (2% PFA and 0.2% glutaraldehyde mixture in 0.1 M phosphate buffer), followed by incubation overnight in fresh fixative at 4°C. Coronal vibratome sections (60 μm) were cut using a Leica VT1000S vibratome and further fixed in 1.5% potassium ferrocyanide: 2% osmium tetroxide in cacodylate buffer for 30 min at 4°C. Tissue was then thoroughly rinsed in distilled water and incubated in 1% aqueous thiocarbohydrazide for 4 min. After further rinsing, the samples were treated with 2% aqueous osmium tetroxide for 30 min, rinsed and en-bloc stained in 1% uranyl acetate for 2 hr. To further enhance contrasts in the samples, one last treatment with Walton’s Lead was carried out for 30 min at 60°C, before proceeding to dehydration in an ethanol series and infiltration with Durcupan ACM resin (Sigma). After embedding and curing, tissue blocks were mounted on Gatan 3View aluminum pins using conductive glue (CircuitWorks Conductive Epoxy) and trimmed accordingly. Before imaging, samples were gold coated to increase electron conductivity. The specimens were then placed inside a Jeol field emission scanning electron microscope (JSM-7100F) equipped with a 3View 2XP system (Gatan). Section thickness was set at 40 nm (Z resolution). Samples were imaged at 2.5kV under high vacuum using a 2048x2048 scan rate, which gave a final pixel size of 4.4 nm.

Electron microscope images were registered and manually segmented using the ImageJ plugin TrakEM2 (Cardona et al., 2012). Extracted 3D structures were exported to the Blender software with the Neuromorph toolset (Jorstad et al., 2015), which was used to compute surface, volume, and length measurements and render 3D reconstructions shown in Figure 1.

Electrophysiology

Mice (21–33 days old) were sacrificed by decapitation following Isoflurane anesthesia, the brain was immediately extracted in ice cold high sucrose solution (in mM: 240 Sucrose, 5 KCl, 1.25 Na2PO4, 2 MgSO4, 1 CaCl2, 26 NaHCO3, 10 D-glucose, Saturated with 95% O2 and 5% CO2). In the same solution acute 300 μm thick coronal hippocampal slices were cut using a Leica vibratome (VT1000 S, Leica Microsystems). Slices were then transferred to a holding chamber containing room temperature ACSF (in mM: 125 NaCl, 5 KCl, 1.25 Na2PO4, 1 MgSO4, 2 CaCl2, 26 NaHCO3, 20 D-glucose) incubated for 1 hr and kept for up to 6 hr of experiments. Cells were visualized with a Scientifica two-photon microscope equipped with a water immersion 40X 0.8 numerical aperture Olympus lens. Dodt Gradient Contrast was used to approach and patch the neurons while a Chameleon femtosecond pulsed laser (Coherent) was used for two-photon imaging of the dendritic arbors. Whole-cell recordings were performed using a Multiclamp 700B amplifier (Molecular Devices), traces were filtered at 3 KHz and digitized at 50KHz. Series resistance was < 20 MΩ. Patch pipettes were pulled (Sutter Puller P-97; Sutter Instruments) from thick-walled borosilicate glass capillaries with an inner filament (1.5 mm outer diameter, 0.86 mm inner diameter; Sutter Instruments). Pipette resistance was 3–4 MΩ after fire polishing. Voltage clamp experiments were performed in the same extracellular ACSF with the addition of 25 μM AP-5 (Cambridge Bioscience) and 20 μM SR95531 (Cambridge Bioscience). The intracellular solution contained in mM: 135 CsMeSO3, 10 HEPES, 10 Na2-Phosphocreatine, 5 Glutathione, 4 MgCl2, 4 Na2ATP, 0.4 NaGTP, 5 QX-314 (Cambridge Bioscience) and 20 μM Alexa Fluor 594 (Molecular Probes). Cells were held at −65 mV and visualized with the two-photon laser tuned at 840 nm, ≈10 min after membrane rupture to allow the dye to spread throughout. Two stimulating unipolar glass electrodes were placed in the distal and proximal dendritic region. Stimulus intensity ranged between 0.1 and 0.8 mA using an Iso-Flex stimulator (Intracel). Fibers were stimulated 20 times with 5 pulse trains every 30 s to calculate the average EPSC. For the NMDAR depletion experiments the extracellular solution contained the AMPAR channel blocker NBQX 10 μM (Santa Cruz Biotechnology) instead of AP-5. Baseline recordings of synaptic inputs, stimulated at proximal and distal sites, were obtained as above but with a single stimulation pulse repeated 20 times every 10 s. MK-801(40 μM, Cayman Chemical Company) was bath applied, and allowed to equilibrate for 5 min, after which the stimulation of either the proximal or distal site was repeated a further 40 times to obtain the curves in Figure 2J. For cross-depletion experiments in Figure S3, we continued the experiment and delivered another 40 stimuli to the site that had not been depleted. Care was taken to vary the site (proximal or distal) that was depleted first. In current clamp experiments the extracellular ACSF contained 20 μM SR95531 while the intracellular solution contained in mM: 115 K-MeSO4, 20 KCl, 10 Na2-Phosphocreatine, 10 HEPES, 2 MgATP, 2 Na2ATP, 0.4 mM Na2GTP, and 20 μM Alexa Fluor 594. The stimulus was delivered with a single bipolar glass theta electrode either distally or proximally with intensity ranging between 1 and 6.2 V. All recordings were performed at 30°C. Fibers were stimulated every 1 min with increasing voltage intensity consisting of 0.2 V steps. All electrophysiology experiments were analyzed in IGOR Pro software (Wavemetrics) with the NeuroMatic 2.7 package and MATLAB (MathWorks). To detect non-linear events, a linear fit was calculated for the gradual increase in EPSP amplitude in response to increasing current steps of the first pulse. The expected amplitude of the second pulse was calculated from this fit. If the recorded EPSP was 2 V above the expected amplitude it was then classified as a non-linearity. All chemicals were from Sigma unless otherwise stated.

Modeling

Simulations were performed with the NEURON simulation environment (Hines and Carnevale, 1997). The model consisted of a soma connected to one dendrite with length of one length constant distributed over 50 segments (diameter = 1 μm). Passive parameters were Cm = 1 μF/cm2, Rm = 10,000 Ω.cm2, Ra = 80 Ω·cm and a leak conductance with a reversal of −65 mV. Unless otherwise noted, active conductances in the dendrite were (in mS/cm2): voltage-activated sodium channels (4), voltage-activated potassium channels (0.8), M-type potassium channels (0.005), high-threshold voltage-activated calcium channels (0.05), low-threshold voltage-activated calcium channels (0.15 × 10-3). Fifteen synapses containing AMPA and NMDA receptor conductances were distributed uniformly over the first or last 50% of the dendrite for the proximal and distal scenarios, respectively. Synaptic conductances were modeled as double exponential functions with gmax = 0.5 nS, τ1 = 0.1 ms, τ2 = 1 ms for AMPA and τ1 = 1 ms, τ2 = 20 ms for NMDA. For fitting short-term plasticity dynamics we initially used the Tsodyks and Markram model (Tsodyks and Markram, 1997), and fitted it by systematically varying the model parameters and comparing the model fits using the peak postsynaptic conductance for each presynaptic pulse on a train of 5 pulses according to , where modelG and dataG are the peak postsynaptic conductance for model and data respectively, for pulse i. While this model produced good fits for the first half of the pulse train, it failed to satisfactorily capture the profile of the remaining pulses. We then fitted the Varela et al. model (Varela et al., 1997) following the same procedure, which produced a good fit to all 5 pulses, mainly because of separately accounting for fast and slow depression kinetics. The Varela et al. model (Varela et al., 1997) was then used in all simulations. Short-term plasticity was implemented with time constants adjusted to match the experimental data and applied to both AMPA and NMDA conductances (facilitation: 2.5 distal, 1.6 proximal; fast depression: 0.4; slow depression: 0.94).

For the passive scenario, active conductances were turned off and synapses only contained AMPA conductances. For the data shown in Figures 3A and 3B a single synapse was activated with a train of ten pulses at 5-100 Hz, and in Figures 3C–3E independent Poisson trains 500 ms long were delivered at each synapse with frequencies varying from 5-80 Hz. Mean depolarization and standard deviation in Figures 3D and 3E were measured during the stimulation period. All simulations were performed at 35°C. For the model variations shown in Figure S4, the following variables were changed for distal synapses: release probability was reduced by 30%, AMPA conductance was reduced by 19%, and dendrite diameter was reduced by 15% (all of the values were estimated from experimental data). The decrease in mean Pr was taken from the double exponential fit in Figure 2J, where the time-constant is inversely proportional to Pr (Hessler et al., 1993). We took the mean time-constant for a double exponential fit, taking into account the fraction of the curve each one accounted for (τmean = τfast∗frfast + τslow∗frslow, where τfast and τslow are the fast and slow time-constants and frfast and frslow are the fractions of the curve they account for), and calculated the change in τmean (proximal: 8.09 ms; distal: 11.49 ms), which is inversely proportional to mean Pr. We find that distal synapses show a ∼30% decrease in mean Pr when compared to proximal compartments. The code for the mathematical model is provided as supplementary information (see supplementary text).

Immuno-histochemistry and FM4-64 staining of primary neurons

Dissociated hippocampal neurons expressing GCaMP3 were immunostained using the following primary antibodies: rabbit αVGlut1 (1:500 Synaptic Systems) and chicken αGFP (1:1000, Abcam). Neurons were fixed in 4% PFA for 20 min and permeabilized using 0.25% Triton X-100 (Sigma, UK) in PBS for 5 min. Cells were incubated with 10% goat serum (Sigma, UK) for 1 hr at room temperature, then incubated with primary antibodies in 2% goat serum overnight at 4°C, and finally with Alexa-conjugated secondary antibodies (1:1000, Molecular Probes) for 1 hr at room temperature. Coverslips were mounted onto glass microscope slides using mowiol. Imaging of immunostained neurons was performed using an Olympus FV1000 confocal microscope equipped with a 40X/0.8 NA water-immersion objective (Olympus).

FM4-64 staining was performed by incubating neurons in high K+ HBS (78.5mM NaCl, 60mM KCl, 10mM HEPES, 10mM Glucose, 2mM CaCl, 1.3mM MgCl) supplemented with 10 μM FM4-64 (Molecular Probes), 1 μM TTX (Alamone Labs), 2.5nM AP-5 and 2nM CNQX (Tocris) for 90 s to load the entire releasable pool of vesicles. Cells were then washed twice for 4 min in HBS (139mM NaCl, 2.5mM KCl, 10mM HEPES, 10mM Glucose, 2mM CaCl, 1.3mM MgCl, 1 μM TTX) to remove FM4-64 from all external membranes. Neurons were then imaged using an Olympus IX71 inverted microscope with a CCD camera (Coolsnap HQ) controlled by Slidebook software (Intelligent Imaging Innovations), equipped with a 40X/1.0 NA oil-immersion objective (Olympus). The excitation light source was a xenon-arc lamp (Lambda LS; Sutter Instruments), in which light exposure was regulated by a rapid shutter (smartShutter; Sutter Instruments) controlled by a Sutter Instruments lambda 10-3 controller. Filtering was provided by a 470 ± 20 nm band pass excitation and 515 ± 20 nm band pass emission (Chroma Technology Corporation) filter set for GCaMP3 and a 565 ± 22 nm band pass excitation and 590-nm long pass dichroic plus 650 ± 36 nm band pass emission (Chroma Technology Corporation) filters for FM4-64.

Statistics

Statistical analysis was performed in MATLAB (MathWorks) and Prism (Graphpad), all the data analyzed with parametric tests was first tested for normality with the D’Agostino and Pearson normality test.

Further information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact, Juan Burrone (juan.burrone@kcl.ac.uk).

REAGENT or RESOURCE	SOURCE	IDENTIFIER
Antibodies
Rabbit polyclonal anti-VGlut1	Synaptic Systems	Cat#153303; RRID: AB_887875
Chicken polyclonal anti-GFP	Abcam	Cat# ab13970; RRID: AB_300798
Chemicals, Peptides, and Recombinant Proteins
AP-5	Cambridge Bioscience	CAY14539; CAS: 79055-68-8
SR 95531	Cambridge Bioscience	CAY14585; CAS: 104104-50-9
QX-314	Cambridge Bioscience	CAY10011032; CAS: 24003-58-5
Alexa Fluor 594	Molecular Probes	A10438
NBQX	Santa Cruz Biotechnology	sc-222048; CAS: 479347-86-9
MK-801	Cayman Chemical Company	10009019; CAS: 77086-22-7
CTZ	Cambridge Bioscience	CAY16335; CAS: 2259-96-3
γDGG	Santa Cruz Biotechnology	sc-203728; CAS: 6729-55-1
FM4-64	Molecular Probes	T13320
CNQX	Tocris	1045; CAS: 479347-85-8
TTX	Alomone Labs	T-550; CAS: 18660-81-6
Experimental Models: Organisms/Strains
Mouse: Syt7 KO: B6.129S1-Syt7tm1Nan/J	The Jackson Laboratory	JAX:004950
Mouse: wt: 129S2/SvPasCrl	Charles River Laboratories	129 Mice
Oligonucleotides
Primers: Syt7mut Forward: CTT GGG TGG AGA GGC TAT TC, Reverse: AGG TGA GAT GAC AGG AGA TCSyt7wt Forward: CAT CCT CCA CTG GCC ATG AAT G Reverse: GCT TCA CCT TGG TCT CCA G	N/A	N/A
Software and Algorithms
TrakEM2	Cardona et al., 2012	https://imagej.net/TrakEM2
Neuromorph	Jorstad et al., 2015	https://neuromorph.epfl.ch/index.html
NEURON	Hines and Carnevale, 1997	https://www.neuron.yale.edu/neuron/download