Dendritic branch structure compartmentalizes voltage-dependent calcium influx in cortical layer 2/3 pyramidal cells
Abstract
Back-propagating action potentials (bAPs) regulate synaptic plasticity by evoking voltage-dependent calcium influx throughout dendrites. Attenuation of bAP amplitude in distal dendritic compartments alters plasticity in a location-specific manner by reducing bAP-dependent calcium influx. However, it is not known if neurons exhibit branch-specific variability in bAP-dependent calcium signals, independent of distance-dependent attenuation. Here, we reveal that bAPs fail to evoke calcium influx through voltage-gated calcium channels (VGCCs) in a specific population of dendritic branches in mouse cortical layer 2/3 pyramidal cells, despite evoking substantial VGCC-mediated calcium influx in sister branches. These branches contain VGCCs and successfully propagate bAPs in the absence of synaptic input; nevertheless, they fail to exhibit bAP-evoked calcium influx due to a branch-specific reduction in bAP amplitude. We demonstrate that these branches have more elaborate branch structure compared to sister branches, which causes a local reduction in electrotonic impedance and bAP amplitude. Finally, we show that bAPs still amplify synaptically-mediated calcium influx in these branches because of differences in the voltage-dependence and kinetics of VGCCs and NMDA-type glutamate receptors. Branch-specific compartmentalization of bAP-dependent calcium signals may provide a mechanism for neurons to diversify synaptic tuning across the dendritic tree.
Editor's evaluation
Synaptic changes in forebrain neurons typically require the conjunction of dendritic action potentials and synaptic activation. Together these signals cause nonlinear changes in calcium influx that then drive plasticity. The strength of these interactions can vary in complex ways. This study uses state of the art methods to convincingly show how some of these complexities arise from known properties of neuronal dendrites and synaptic NMDA receptors.
https://doi.org/10.7554/eLife.76993.sa0Introduction
The synaptic tuning properties of neurons are regulated by calcium-dependent plasticity signals, including voltage-dependent calcium influx evoked by back-propagating action potentials (bAPs) (Magee and Johnston, 1997; Sjöström and Nelson, 2002; Häusser and Mel, 2003). Many investigations have focused on timing-dependent plasticity rules, such as spike-timing dependent plasticity (STDP), in which the delay between synaptic input and the bAP determines if an induction protocol results in long-term potentiation (LTP), depression (LTD), or maintenance of synapse strength (Feldman, 2012). However, other studies have shown that the outcome of plasticity induction protocols also depends on the synaptic location within the dendritic tree (Golding et al., 2002; Froemke et al., 2005; Sjöström and Häusser, 2006; Gordon et al., 2006). The same protocol can evoke a different sign, magnitude, or temporal profile of synaptic plasticity if applied to synapses in proximal or distal dendritic compartments. Such location-dependent plasticity rules may bolster the integrative capacity of individual neurons by supporting divergent tuning properties in basal and apical dendritic arbors (Larkum et al., 1999; Xu et al., 2012; Iacaruso et al., 2017).
Location-dependent plasticity rules vary along the proximodistal axis of dendrites because bAPs are attenuated as they propagate away from the soma (Regehr et al., 1989; Spruston et al., 1995; Golding et al., 2002; Waters et al., 2003; Froemke et al., 2005; Sjöström and Häusser, 2006). In distal dendritic compartments, smaller amplitude bAPs evoke less voltage-dependent calcium-influx, such that induction protocols that usually evoke LTP result in LTD (Nevian and Sakmann, 2006; Sjöström and Häusser, 2006). Experimentally boosting distal bAP amplitude can restore LTP by amplifying bAP-dependent calcium influx (Sjöström and Häusser, 2006), suggesting that location-dependent variation in synaptic plasticity rules is primarily determined by the local amplitude and timing of bAP-dependent calcium influx.
In principle, bAP-dependent calcium influx can also vary in a branch-specific manner, independent of distance from the soma. Dendritic bAP-dependent calcium influx could be shaped by local variability in ion channel expression, synaptic inputs, or electrotonic impedance due to dendritic branch structure (Rall and Rinzel, 1973; Vetter et al., 2001; Frick et al., 2004; Losonczy et al., 2008; Yaeger et al., 2019). Branch-specific variation in bAP-dependent calcium influx provides a potential mechanism for diversifying synaptic tuning across the dendritic tree, which is a core feature of computational models of neurons as hierarchical nonlinear integrators (Poirazi et al., 2003; Häusser and Mel, 2003; Francioni and Harnett, 2021; Bicknell and Häusser, 2021). The existence of branch-specific variation in bAP-dependent calcium influx would also affect the interpretation of dendritic calcium transients imaged in vivo, which can be evoked by local synaptic integration or global bAP-dependent influx (Wilson et al., 2016; Iacaruso et al., 2017; Yaeger et al., 2019; Beaulieu-Laroche et al., 2019; Kerlin et al., 2019; Francioni et al., 2019).
To discover potential dendrite branch-specific and distance-independent variation in bAP amplitude and bAP-dependent calcium influx, we measured bAP-dependent calcium influx in multiple dendrites of individual cortical L2/3 pyramidal cells in mice. We used somatic and dendritic electrical recordings, two-photon calcium imaging and glutamate uncaging, and dendritic voltage-imaging to experimentally identify dendrite branch-specific changes in bAP amplitude and its consequences for bAP-dependent calcium influx. In some branches, bAPs fail to evoke calcium influx through voltage-gated calcium channels (VGCCs), but still amplify synaptically-mediated calcium influx through NMDA receptors (NMDARs). We demonstrated that these branches successfully propagate bAPs in the absence of synaptic input and contain VGCCs, but exhibit more elaborate dendritic branching, which decreases their electrotonic impedance and reduces the amplitude of the bAP. Our results reveal branch-specific variation in bAP-dependent calcium signals, providing an additional mechanism by which calcium-dependent plasticity induction may vary throughout individual neurons.
Results
We measured calcium influx in the apical dendrites of cortical L2/3 pyramidal cells to investigate whether voltage-dependent calcium signals are regulated in a branch-specific manner. We acquired whole-cell current-clamp recordings from individual cells filled through the recording pipette with 10 µM Alexa Fluor 594 and 300 µM Fluo-5F to visualize neuronal morphology and monitor changes in intracellular calcium concentration using red and green fluorescence, respectively (Carter and Sabatini, 2004). Calcium-dependent fluorescence transients were measured as the change in green relative to red fluorescence (), which is linearly proportional to and comparable across recordings (Carter and Sabatini, 2004). We found that bAP-evoked calcium influx () varies across dendritic spines from within the same cells, even when matched for distance from the soma (Figure 1). Across many neurons and recording sites, we observed that was highly correlated in nearby spines ( < 6 µm) and between spines and their parent dendritic shafts (Figure 1—figure supplement 1), indicating that variations in are regulated across dendritic branches, rather than across individual dendritic spines. These data suggest that voltage-gated calcium channels (VGCCs), which mediate bAP-evoked calcium influx and are responsible for various forms of plasticity (Kapur et al., 1998; Dolmetsch et al., 2001; Yasuda et al., 2003; Nevian and Sakmann, 2006; Scheuss et al., 2006, Bloodgood and Sabatini, 2007; Brigidi et al., 2019), can be differentially activated by bAPs in a branch-specific manner within individual L2/3 pyramidal cells.
To investigate the mechanisms of branch-specific variation in voltage-dependent calcium signals, we measured calcium influx in response to glutamate uncaging-evoked excitatory postsynaptic potentials (uEPSPs) as a proxy for synaptically evoked calcium influx (). In addition, we measured bAP-dependent amplification of synaptically-evoked calcium influx due to transient relief of the Mg2+ block of NMDARs using uEPSP-bAP pairings (). In all uEPSP-bAP pairings, we delayed the bAP by 5ms relative to the glutamate uncaging laser pulse to maximize bAP-dependent amplification of postsynaptic calcium influx and mimic protocols that robustly induce NMDAR-dependent LTP (Froemke et al., 2005; Nevian and Sakmann, 2006; Holbro et al., 2010). Spines on branches with and without bAP-evoked calcium influx showed and supralinear bAP-dependent amplification in (Figure 1B–E; Holbro et al., 2010), even though one of the spines exhibited almost no bAP-evoked calcium influx when the bAP was evoked by itself (Figure 1D–E). These data suggest that bAP-evoked calcium influx through VGCCs and bAP-dependent amplification of synaptically evoked calcium influx through NMDARs can be regulated independently in a branch-specific manner (Bloodgood and Sabatini, 2007; Holbro et al., 2010).
bAP-evoked calcium influx and bAP-dependent amplification are decoupled in L2/3 pyramidal cells
To systematically investigate the regulation of bAP-evoked calcium influx and bAP-dependent amplification of synaptically evoked calcium influx, we recorded , , and in 97 dendritic spines from 33 cells and 15 mice (N = 97/33/15) spanning the full proximodistal extent of the dendrites. To focus our analysis on the supralinear component of , which is primarily carried by NMDARs (Holbro et al., 2010), we measured amplification as (Figure 2A). We found that and were uncorrelated across dendritic spines in L2/3 pyramidal cells (Figure 2B). To visualize the divergence between and , we grouped the data into high (0.1 < < 0.3, selected to avoid outliers), and low ( < 0.04) populations (Figure 2B) based on the amplitude of bAP-evoked calcium influx (Figure 2C–D). We plot data with this selection criterion throughout the figures to highlight differences across branches but focus our statistical analyses on the full distribution of data which is presented as scatter plots in each figure. The same criterion is used throughout all figures except in Figures 5B–D. Low branches exhibited almost no calcium influx evoked by bAPs but still exhibited large bAP-dependent amplification (Figure 2C–D). and showed similar patterns in dendritic spines and parent dendritic shafts (Figure 2C–D). We found that although attenuates with distance from the soma (Figure 2E and Figure 2—figure supplement 1B), it exhibits significant variability in distal branches, such that pairs of dendritic spines that are on different branches within the same neuron but at a similar distance from the soma sometimes had large differences in calcium influx (Figure 2E–F). varied over a similar range at all distances from the soma (Figure 2—figure supplement 1A). The divergence between and was not related to the degree of NMDAR activation (Figure 2—figure supplement 1C) nor the uEPSP amplitude or time-course (Figure 2—figure supplement 1D-H). Furthermore, we observed low branches in the presence of MNI-Glutamate, which is an antagonist of ionotropic GABA receptors, indicating that low branches are not caused by inhibition (Figure 2B). We confirmed this with the selective GABA-A receptor antagonist 10 µM SR 95531 (data not shown). Together, these data indicate that and are regulated independently in the dendrites of L2/3 pyramidal cells.
We considered three distinct mechanistic explanations for why low branches exhibit bAP-dependent amplification () but fail to exhibit bAP-evoked calcium influx ():
Low branches contain fewer (or no) voltage-gated calcium channels (VGCCs).
bAPs fail to backpropagate into low branches without the support of EPSPs.
bAP amplitude in low branches does not exceed the threshold for VGCC opening.
Although each of these mechanisms could be realized in several ways, they represent general explanations that we tested directly.
Low branches contain voltage-gated calcium channels
Low branches may contain functional VGCCs that are not opened by a single bAP but do open in response to stronger depolarizations. To examine this possibility, we imaged calcium influx evoked by a burst of 5 bAPs at 150 Hz in dendritic shafts located >75 µm from the soma in L2/3 pyramidal cells (Figure 3A; Larkum et al., 2007). In most low branches, calcium influx could be evoked by a burst of 5 bAPs (Figure 3A–B). However, the calcium influx evoked by a burst of 5 bAPs in high branches was higher than in low branches (Figure 3A–B). These data are consistent with the presence of VGCCs in low branches, although they do not distinguish whether there is a lower density of VGCCs (or a higher voltage threshold), or if bursts of bAPs are less effective at depolarizing low branches.
We attempted to make the depolarization evoked across the dendritic tree more uniform by triggering bAPs in the presence of the potassium channel antagonist 4-AP (Figure 3C–D; Hoffman et al., 1997). In all somatic recordings, application of 2 mM 4-AP increased the after-depolarization of the bAP (Figure 3C, inset), confirming the efficacy of the drug. Applying 4-AP significantly increased in every branch, regardless of how much calcium influx was evoked in control conditions (Figure 3C–D). Because evoked in the presence of 4-AP is similar in high and low branches (Figure 3C), we conclude that VGCCs are present in all dendritic compartments and mediate substantial calcium influx in response to sufficiently large depolarizations.
To determine if VGCCs can be effectively opened by local depolarization in low branches, we measured calcium influx evoked by direct application of 1 mM glutamate delivered via a patch pipette using a picospritzer (Figure 3E–H). We added 20 µM CPP and 50 µM MK-801 to the puff and bath solutions to block NMDARs and limit our measurements to calcium influx mediated by VGCCs. In every dendrite examined, direct application of glutamate led to significant calcium influx, regardless of (Figure 3F–H). The calcium signal was not an artifact of mechanical displacement due to the picospritzer (Figure 3—figure supplement 1A-B) and did not include an NMDAR component (Figure 3—figure supplement 1C-D), consistent with calcium influx mediated by VGCCs. We note that the large calcium signals evoked by glutamate demonstrate that the similarity in bAP-evoked calcium influx between high and low measured in 4-AP was not due to saturation of the calcium indicator. Collectively, our data indicate that all dendritic branches contain VGCCs, which can be opened by local depolarization. However, a specific population of low branches exhibit little to no calcium influx in response to back-propagating bAPs.
bAPs propagate to low branches in the absence of EPSPs
We only observe significant bAP-dependent calcium influx in low branches if the bAP is paired with synaptic input evoked by glutamate uncaging. Despite the small size of the uEPSP, it is possible that bAPs only back-propagate to low branches when assisted by local depolarization (Magee and Johnston, 1997; Gasparini et al., 2007). Since most of the depolarization evoked by EPSPs is carried by AMPARs, we measured , , and after applying the selective AMPAR antagonist, NBQX (Figure 4A–B). We observed bAP-dependent amplification of the remaining NMDAR-mediated calcium influx in all dendritic spines examined, despite fully blocking the depolarization evoked by glutamate uncaging (Figure 4A–B). As expected, glutamate uncaging evokes some calcium influx via NMDARs in the presence of NBQX, since NMDARs are only partially blocked by Mg2+ at resting potentials (Jahr and Stevens, 1990) and the driving force on calcium is extremely high. These data suggest that bAPs propagate to low branches without an uEPSP-mediated local depolarization.
To test whether bAPs propagate to low branches in the absence of glutamate receptor activation, we acquired loose-patch recordings in current-clamp from 19 dendrites in 15 neurons while simultaneously recording from the soma to measure the dendritic electrical signal evoked by a bAP (Figure 4C–I). Simultaneous two-photon imaging was used to monitor dendritic calcium influx. All recording sites were ≥75 µm from the soma. In the neuron shown in Figure 4C, we acquired consecutive loose-patch recordings from one high branch and one low branch. Only the high branch had detectable , but both dendritic recordings exhibited clear electrical signals evoked by the bAP (Figure 4D–F). We observed bAP-evoked electrical signals in all dendrites, even though only some exhibited (Figure 4G–I). These loose-patch recordings do not permit direct measurement of intracellular bAP amplitude (Figure 4—figure supplement 1); however, they demonstrate that bAPs propagate to low branches without the support of EPSPs.
bAP amplitude is attenuated in low branches
To measure the intracellular dendritic voltage evoked by bAPs, we performed interleaved dendritic voltage and calcium imaging while performing somatic whole-cell current-clamp recordings (Figure 5). Voltage imaging was achieved by expressing a highly-sensitive near-infrared genetically-encoded voltage indicator, QuasAr6a (Wang et al., 2021). Both voltage and calcium signals were recorded via structured illumination one-photon microscopy. All recording sites were ≥75 µm from the soma (except for the closest site in panel A). Consistent with our results from Figure 4, we observed bAP-dependent voltage signals in all branches, including those that did not exhibit bAP-evoked calcium influx (Figure 5). To visualize differences in high and low branches, we adjusted our selection criterion to account for differences in the range of measured with two-photon and one-photon calcium imaging (the selection range described above was scaled by 1.75, Figure 5B–D). Although low branches exhibited clear bAP-evoked voltage waveforms (Figure 5C), the amplitudes of the dendritic bAP-evoked voltage waveforms were attenuated in low branches (Figure 5C–D). The full-width half-max of the voltage waveforms was 6.2 ±0.9ms for high branches and 8.0 ±4.1ms for low branches. These data demonstrate that bAP amplitude is reduced in low branches.
Collectively, our data support the conclusion that low branches exhibit less bAP-evoked calcium influx because of branch-specific reductions in local bAP amplitude. First, our experiments with bursts of 5 bAPs, application of 4-AP, and glutamate application show that low branches contain VGCCs that open in response to large depolarizations (Figure 3). Second, our experiments with NBQX and dendritic loose-patch recordings show that bAPs successfully propagate to low branches in the absence of synaptic input (Figure 4). Third, our voltage imaging experiments show that the amplitude of bAPs is attenuated in low branches (Figure 5). These data are consistent with the mechanistic explanation that bAPs fail to evoke calcium influx in low branches because the local bAP amplitude does not exceed the threshold for VGCC opening. Therefore, we investigated why bAPs are selectively attenuated in low branches.
bAP-evoked calcium influx attenuates at branch points
bAP amplitude and bAP-evoked calcium influx can attenuate at branch points due to dendritic impedance mismatch and localized expression of potassium channels (Spruston et al., 1995; Williams and Stuart, 2000; Vetter et al., 2001; Frick et al., 2003; Cai et al., 2004; Gasparini et al., 2007; Harnett et al., 2013). Consistent with these observations, we found that often attenuates across branch points, as revealed by comparing differences in at pairs of sites within the same dendritic segment or across one dendritic branch point (Figure 6—figure supplement 1A-C). The ratio of (influx at distal site divided by proximal site) tends to be smaller for pairs of sites spanning a branch point than for those within the same segment (rank-sum: p = 0.0145, Figure 6—figure supplement 1C). This result was corroborated by the observation that low branches typically had higher dendritic branch order (defined as the number of on-path branch points between the soma and the recording site) than high branches (Figure 6—figure supplement 1D). While these data suggest a relationship between and the pattern of dendritic branching, they do not fully account for the variance in observed across branches in L2/3 pyramidal cells (Figure 6—figure supplement 1D, R2 = 0.178).
Low branches have a more elaborate branch structure than high branches
The electrotonic impedance (i.e., passive input resistance) of a dendrite scales inversely with the surface area of membrane. Therefore, if low branches are within sections of the dendritic tree that contain more extensive branching, then they have a smaller impedance, which might reduce bAP amplitude and . To examine this possibility, we developed a metric based on the dendritic morphology called ‘branch complexity’ that accounts for all nearby dendritic branches and is independent of branch order. This approach measures the total length of dendrite near a recording site with an exponential filter to discount branches that are further from the site (so ithas units of millimeters). We measured branch complexity over a length constant of λ = 145 µm because propagating bAPs simultaneously depolarize ~145 µm of dendrite in L2/3 pyramidal cells (conduction velocity is 154.9 µm/ms and AP duration is 0.94ms, Figure 6—figure supplement 2).
We compared the branch complexity of pairs of dendritic recording sites within the same cell that were located approximately the same distance from the soma but for which there was a large difference in (Figure 6A–E). The low site had a higher branch complexity than the high site in all pairs (Figure 6E), indicating that low sites are surrounded by more extensive branching than high sites. To confirm that branch complexity is inversely proportional to input resistance, we reconstructed all cells used in this analysis as biophysical compartment models in the NEURON simulation environment and directly measured the input resistance by injecting small hyperpolarizing currents into each recording site (Carnevale and Hines, 2006). In every pair, the low site had a smaller input resistance than the high site (Figure 6F). Additionally, branch complexity was inversely correlated with input resistance (R2 = 0.93 for λ = 145 µm), confirming its validity as a measure of electrotonic properties (Figure 6—figure supplement 3A-B). These results and conclusions hold for measurements of branch complexity over length constants ranging from 5 to 400 µm (Figure 6—figure supplement 3C).
To inspect why low branches had a higher branch complexity, we plotted the average number of dendritic branches distal and proximal to the target site (Figure 6G). Low branches are immediately proximal to dendritic branch elaboration (e.g. Figure 6A c.f. with Figure 6H, left), and closer to sister dendrites emerging proximal to the measurement site (c.f. with Figure 6H, middle and right). These results demonstrate that low branches are surrounded by more elaborate branching patterns than high branches, suggesting that dendritic branch structure may be sufficient to selectively reduce the amplitude of bAPs.
Dendritic branch structure is sufficient to explain reductions in bAP-evoked calcium influx
To determine if dendritic branch structure is sufficient for reducing , we reproduced our experimental results using a biophysical compartment model. We held all physiological parameters constant throughout each model except for the dendritic branch structure, allowing us to specifically measure the effects of dendritic morphology between high and low branches (Vetter et al., 2001). Using NEURON compartmental models of 8 reconstructed L2/3 pyramidal cells (all from the pairwise analysis described in Figure 6), we evoked a bAP with somatic current injection and measured the intracellular voltage waveform and the resulting calcium conductance at the same sites in which we had imaged calcium influx in acute brain slices (Figure 7A–E). For each cell, we fit the density of dendritic A-Type potassium channels such that the bAP peaked at –10 mV in the high site (Figure 7D), which is consistent with the amplitude of bAPs measured ∼100 µm from the soma in layer 2–3 pyramidal cells with dendritic whole-cell recordings (Larkum et al., 2007; Smith et al., 2013). We used this density of A-Type potassium channels throughout the entire dendritic tree of each cell. Under these conditions, for every within-cell pair, the low site had smaller bAP amplitude and peak calcium conductance, despite having the same channel densities as the high site (Figure 7B–E). As an additional measure of dendritic excitability, we plotted the A-Type potassium channel density required for bAPs to peak at –10 mV in each dendritic site (Figure 7F). Low sites required lower potassium channel densities for their bAPs to peak at –10 mV, indicating that the branch structure surrounding low sites is less excitable than for high sites (Figure 7F). To visualize the spread of bAP-evoked voltage and calcium signals throughout the dendritic tree, we plotted two example cells with the color of each segment indicating the peak voltage and calcium conductance (Figure 7—figure supplement 1). These data demonstrate that dendritic branch structure is sufficient to impact by altering the dendritic bAP amplitude in a branch-specific manner.
We performed a simulated ‘cutting’ experiment in which we removed dendritic branches from the reconstruction that were closer to the low site than the high site (Figure 7A), to determine if simplifying the dendritic branch structure can recover . After cutting dendritic branches near the low site, differences in the bAP amplitude, calcium channel conductance, and the A-Type potassium channel density required to generate bAPs that peak at –10 mV were all nullified (Figure 7B–F). We also computationally recapitulated our results with 4-AP (Figure 3C–D) by artificially removing A-Type potassium channels (Figure 7—figure supplement 2A-D). These data confirm that dendritic branch structure directly affects local dendritic excitability. Furthermore, these results support the model that increased dendritic branching around low recording sites is sufficient to reduce bAP amplitude and .
Branch structure-dependent changes in dendritic excitability may also affect the amplitude, kinetics, and efficacy of synaptic input arriving at each branch. Using our compartment models, we simulated synaptic input by injecting a conductance into each dendritic site (Figure 7G–H). As expected, the local EPSP amplitude was smaller in low branches because of their reduced input resistance (Figures 7G and 6F). Although EPSP size measured at the soma (but evoked in the dendrites) was lower for synapses on low branches, the difference was minimized due to frequency-dependent cable attenuation (Figure 7H, effect size dendritic EPSPs: –3.23, effect size somatic EPSPs: –1.4). The differences in local EPSP amplitude were abolished after cutting dendritic branches near the low site (Figure 7—figure supplement 2E-F). Therefore, although dendritic branch structure affects the local amplitude of synaptic potentials, it has a minimal impact on the efficacy of synaptic input at depolarizing the soma.
Reductions in bAP amplitude have a selective effect on VGCC-mediated calcium influx
Our data support a model in which elaborated dendritic branching reduces dendritic excitability, which lowers bAP amplitude and bAP-evoked calcium influx. However, it remains to be determined why bAP-dependent amplification of synaptic NMDAR-mediated calcium influx is spared by reductions in bAP amplitude that impact bAP-evoked calcium influx through VGCCs. Using published data (Jahr and Stevens, 1990; Reuveni et al., 1993), we plotted the voltage-dependent open-probabilities and time constants for VGCCs and NMDARs (Figure 8A–B). Using these parameters, we simulated VGCC and NMDAR conductances and the resulting voltage-dependent calcium influx in response to quadratic depolarizing voltage stimuli designed to resemble bAPs (Figure 8C). Due to the steep voltage-sensitivity and long time constant of VGCCs (Figure 8A–B), they exhibit a dramatic drop-off of calcium influx for voltage steps peaking below –10 mV (Figure 8C–E). Although the open-probability of NMDARs is highly sensitive to peak voltage (Figure 8C–D), the interaction between open-probability and driving force minimizes differences in calcium influx evoked across a wide range of voltages (Figure 8E). These data demonstrate how small reductions in bAP amplitude can eliminate mediated by VGCCs, while largely sparing bAP-dependent amplification of synaptically evoked, NMDAR-mediated calcium influx.
Discussion
Back-propagating APs regulate synaptic plasticity by evoking voltage-dependent calcium influx throughout dendrites. Here, we show that bAP-dependent calcium influx varies in a dendrite branch-specific manner in cortical L2/3 pyramidal cells due to branch-specific reductions in input resistance inherited from the local morphology of the dendritic tree. Sections of the dendritic tree with more elaborate branch patterns have lower input resistance, leading to a branch-specific reduction in bAP amplitude. These branches have a selective reduction in bAP-evoked calcium influx through VGCCs despite containing VGCCs and successfully propagating bAPs in the absence of synaptic input. However, these branches maintain bAP-dependent amplification of synaptically-evoked calcium influx through NMDARs due to the shallower voltage-dependence and faster kinetics of the Mg2+-block of NMDARs. Our results demonstrate that dendritic excitability, bAPs, and bAP-dependent calcium signals vary between branches, which may provide a mechanism for synaptic plasticity rules and the computational properties of dendrites to vary across subcompartments of individual neurons.
Branch-specific variation in dendritic morphology shapes voltage-dependent calcium signals
Many studies have demonstrated that bAP amplitude and bAP-evoked calcium influx are attenuated in the distal compartments of pyramidal cells (Regehr et al., 1989; Spruston et al., 1995; Magee and Johnston, 1997; Waters et al., 2003; Froemke et al., 2005; Sjöström and Häusser, 2006; Larkum et al., 2007). Our work additionally demonstrates that bAP amplitude and can vary across branches, independent of distance from the soma. This variance can be explained by the branch structure of the dendritic tree. Previous studies have focused on attenuation of as a function of dendritic branch order (Spruston et al., 1995; Williams and Stuart, 2000); however, we find that explaining the observed variance in requires consideration of the full dendritic morphology, including the number of branches distal to the recording site and on sister dendrites. Elaborate branching patterns reduce input resistance, local bAP amplitude, and bAP-evoked calcium influx because dendritic input resistance is proportional to the surface area of membrane. Although our experiments show that the effect of branch structure on impedance is sufficient to cause the branch-specific differences we observed, additional mechanisms, such as variation in potassium channel density, may also contribute. We note that because bAPs are time varying signals, it is necessary to consider the complex impedance of each dendrite to fully explain bAP attenuation. As an alternative, we use computational models that account for complex impedance to demonstrate that bAP attenuation is predicted by the dendritic branching structure.
Our data cannot precisely resolve how much bAPs are attenuated in low branches because loose-patch recordings do not directly measure intracellular voltage (Figure 4) and voltage imaging signals are not precisely calibrated to absolute voltage (Figure 5). However, due to the differential impact on VGCC- and NMDAR-mediated calcium influx, we hypothesize that bAPs peak between –40 mV and –15 mV in low branches, based on our analysis in Figure 8. Testing this prediction with voltage indicators would require precise calibration of dendritic voltage-dependent fluorescence, which requires knowledge of the dendritic resting potential and a perturbation to a known voltage (such as 0 mV). These results provide a key insight into the mechanisms underlying the intracellular heterogeneity of dendritic biophysical properties and bolster our understanding of dendritic physiology.
Activation of VGCCs in low branches
Our data indicate that bAP-dependent activation of VGCCs is branch-specific in cortical L2/3 pyramidal cells. Low branches exhibit little to no calcium influx evoked by single bAPs or even bursts of 5 bAPs at 150 Hz. In vivo measurements of L2/3 pyramidal cells have not reported high frequency bursts of bAPs (Smith et al., 2013; Wei et al., 2020), suggesting that in general, somatic depolarization does not lead to activation of VGCCs in low branches without additional depolarization from local synaptic input. We demonstrated that dendritic glutamate application can evoke VGCC-mediated calcium influx in low dendrites; however, this approach was designed to maximize local depolarization, rather than to resemble a physiological stimulus.
The lack of bAP-dependent activation of VGCCs in low branches may have important implications for neuronal function. First, because VGCC-mediated influx is responsible for multiple forms of synaptic and cellular plasticity (Kapur et al., 1998; Dolmetsch et al., 2001; Yasuda et al., 2003; Nevian and Sakmann, 2006; Scheuss et al., 2006, Brigidi et al., 2019), these data suggest that each of these forms of plasticity may be branch-specific in L2/3 pyramidal cells, at least with respect to their activation by bAPs. Second, these results have clear significance for the interpretation of dendritic calcium signals in vivo. A fundamental challenge in measuring dendritic activity with calcium imaging in vivo is distinguishing between calcium influx evoked by local synaptic input and global calcium signals, such as those evoked by bAPs (Xu et al., 2012; Beaulieu-Laroche et al., 2019; Francioni et al., 2019; Kerlin et al., 2019). Our results indicate that in low branches, bAPs do not evoke calcium signals without the addition of local synaptic input, which could make them a useful site for measuring local dendritic processing.
Implications for synaptic plasticity
Previous work has demonstrated that bAP-evoked calcium influx selectively drives LTD while bAP-dependent amplification of synaptic calcium influx selectively drives LTP (Feldman, 2000; Matsuzaki et al., 2004; Nevian and Sakmann, 2006; Lee et al., 2009). Therefore, identifying branches with a selective reduction in bAP-evoked calcium influx while maintaining bAP-dependent amplification of synaptically mediated calcium influx suggests that the ratio of LTD to LTP varies across branches in L2/3 pyramidal cells. Since the ratio of STDP-dependent depression to potentiation determines the cooperativity, strength, and sparsity of synaptic inputs (Song et al., 2000; Froemke et al., 2005; Sjöström and Häusser, 2006), our work suggests that synaptic tuning properties may vary in a branch-specific manner based on the local dendritic branch structure. These findings complement a series of papers showing that calcium-dependent plasticity rules vary as a function of distance from the soma (Golding et al., 2002; Froemke et al., 2005; Sjöström and Häusser, 2006; Gordon et al., 2006).
Our experiments were conducted in juvenile mice (between P21 and P28). Although the dendritic branch structure is static after this age (Koleske, 2013; Richards et al., 2020), it is possible that changes in channel composition dampen the effects of dendritic branch structure on excitability. We show that such effects are possible with compartment models (Figure 7F), and it is known that changes in dendritic excitability can occur as a result of synaptic plasticity induction (Frick et al., 2003; Frick et al., 2004; Losonczy et al., 2008). Therefore, our findings may suggest a developmental function of low branches in shaping synaptic connectivity, consistent with the lower incidence of calcium influx in specific branches of L2/3 pyramidal cells during development in visual cortex (Yaeger et al., 2019).
Implications for dendritic computation
Our work suggests that the dendritic branch structure of low branches can expand the representational capacity of neurons (Häusser and Mel, 2003). First, the lower input resistance of low branches makes them more compartmentalized from the rest of the cell, such that they may be more biased toward performing local computations (Polsky et al., 2004). Interestingly, our compartment models indicate that compartmentalization is asymmetric (Williams, 2004), such that forward propagation of synaptic input from low branches to the soma may be equivalent to high branches, despite their different biophysical properties. Second, branch-specific variation in calcium-dependent plasticity signals may diversify synaptic tuning, which would increase dendritic computational capacity (Poirazi et al., 2003; Häusser and Mel, 2003; Francioni and Harnett, 2021; Bicknell and Häusser, 2021). Together, these considerations indicate that low branches may be hotspots for dendritic computation due to their enhanced dependence on cooperative synaptic processing and higher likelihood of containing heterogeneous synaptic tuning.
Materials and methods
Genetic constructs
Request a detailed protocolFor voltage imaging, we used QuasAr6a, an improved Archaerhodopsin-derived near-infrared genetically encoded voltage indicator (GEVI) (Wang et al., 2021). To improve expression and membrane trafficking we created the construct CAG::QuasAr6a-TS-dmCitrine-TSx3-ER2-p2a-jRGECO1a-CAAX, following the design in Adam et al., 2019, where TS is the trafficking sequence from Kir2.1 (Gradinaru et al., 2010), dmCitrine is the non-fluorescent Y66G mutant of mCitrine (Adam et al., 2019), and ER2 is the endoplasmic reticulum export signal FCYENEV (Gradinaru et al., 2010). The self-cleaving p2a linker enabled bicistronic co-expression of a membrane-targeted Ca2+ indicator, jRGECO1a-CAAX. In the present studies, this indicator was only used for locating expressing neurons; for consistency with other experiments in this manuscript, we used the blue-shifted dye, Fluo-5F for Ca2+ imaging.
The genes were cloned into a second-generation lentiviral backbone with a CAG promoter (Addgene: #124775) using standard Gibson Assembly. Briefly, the vector was linearized by double digestion using BamHl and EcoRl (New England Biolabs, Ipswich, MA) and purified by the GeneJET gel extraction kit (ThermoFisher, Waltham, MA). DNA fragments were generated by PCR amplification and then fused with the backbones using NEBuilder HiFi DNA assembly kit (New England Biolabs). All plasmids were verified by sequencing (GeneWiz, Cambridge, MA).
In utero electroporation (IUE)
Request a detailed protocolAll procedures involving animals were in accordance with the National Institutes of Health Guide for the care and use of laboratory animals and were approved by the Harvard University Institutional Animal Care and Use Committee (IACUC, Protocol #: IS00000571-3). The IUE surgery was made as described previously (Kwon and Sabatini, 2011). Embryonic day 15.5 (E15.5) timed-pregnant female CD1 mice (Charles River, Wilmington, MA) were deeply anesthetized and maintained with 2% isoflurane. The animal body temperature was maintained at 37 °C. Uterine horns were carefully exposed, and periodically rinsed with warm PBS. The plasmid DNA was diluted with PBS (2 μg/μL; 0.005% fast green), and 1 µL of the mixture was injected into the left lateral ventricle of pups. Electrical pulses (40 V, 50ms duration) were delivered five times at 1 Hz using a tweezers electroporation electrode (CUY650P5; Nepa Gene, Ichikawa, Japan). Injected embryos were placed back into the abdominal cavity, and the surgical wound was sutured with PGCL25 absorbable sutures (Patterson, Saint Paul, MN).
Slice preparation
Request a detailed protocolAcute coronal slices were prepared from the somatosensory cortex of young adult, C57Bl/6 j wild-type mice (Jackson Labs, Bar Harbor, ME) between postnatal days P21-P28. Animals were anesthetized via inhalation of isoflurane then immediately decapitated. The brain was rapidly removed and placed in ice-cold cutting artificial cerebrospinal fluid (ACSF) containing (in mM): 125 NaCl, 2.5 KCl, 25 NaHCO3, 1.25 NaH2PO4, 1 CaCl2, 10 MgCl2, and 25 glucose, saturated with 95% O2 and 5% CO2. We cut coronal slices with a VT1200S Vibratome (Leica, Buffalo Grove, IL) while maintained in the cutting solution. Slices were then incubated at 34 °C for 30 min in a recovery ACSF containing (in mM): 92 NaCl, 28.5 NaHCO3, 2.5 KCl, 1.25 HaH2PO4, 2 CaCl2, 4 MgCl2, 25 glucose, 20 HEPES, 3 Sodium Pyruvate, and 5 Sodium Ascorbate, saturated with 95% O2 and 5% CO2. Slices were next transferred to a 20 °C room in recovery ACSF until use ( > 30 min after slicing) and not held for longer than 8 hr. Experiments were conducted at 34 °C in recording ACSF containing (in mM): 125 NaCl, 2.5 KCl, 25 NaHCO3, 1.25 NaH2PO4, 1.5 CaCl2, 1 MgCl2, and 25 glucose, saturated with 95% O2 and 5% CO2. All solutions were prepared with osmolality between 300–310 mOsm/kg, adjusted with either water or glucose. For voltage imaging experiments in Figure 5, coronal slices were prepared from CD1 mice between P21-P28. The slicing solution contained (in mM): 210 sucrose, 3 KCl, 26 NaHCO3, 1.25 NaH2PO4, 5 MgCl2, 10 D-glucose, 3 sodium ascorbate and 0.5 CaCl2, and was saturated with 95% O2 and 5% CO2. The slices were transferred to an incubation chamber containing the recording solution (in mM): 124 NaCl, 3 KCl, 26 NaHCO3, 1.25 NaH2PO4, 2 MgCl2, 15 D-glucose and 2 CaCl2, and was saturated with 95% O2 and 5% CO2.
Electrophysiology
Request a detailed protocolSomatic whole-cell recordings were acquired from cortical excitatory L2/3 pyramidal cells using IR-DIC. Patch pipettes (2–4.5 MΩ) were filled with an internal solution containing (in mM): 130 K-Gluconate, 10 KCl, 10 HEPES, 4 MgATP, 0.5 Na2GTP, 10 Phosphocreatine-disodium salt, 0.3 Fluo-5F and 0.01 Alexa 594. For voltage-imaging experiments in Figure 5, the internal solution contained (in mM): 8 NaCl, 130 KMeSO3, 10 HEPES, 5 KCl, 0.5 EGTA, 4 Mg-ATP, 0.3 Na3-GTP, 0.3 Fluo-5F, and 0.01 JF549i. The pH was adjusted to 7.3 using KOH and osmolality was adjusted to 285–295 mOsm/kg with water. In a subset of experiments shown in Figure 2E, 300 µM Fluo-4 was used instead of Fluo-5F, but amplitudes of calcium signals were comparable, so we merged the data. All recordings were performed with a Multiclamp 700B amplifier (Molecular Devices, San Jose, CA). Pipette capacitance was neutralized prior to break-in and the series resistance was fully balanced for all recordings. Series resistance ranged from 7 to 25 MΩ. We elicited APs by injecting a 1–2ms current of 0.5–3.5 nA. The average resting membrane potential was –79.5 mV with a standard deviation of 6.13 mV (N = 420). Recordings were aborted if the cell’s membrane potential exceeded –60 mV, if the dendrites begun filling with calcium (as indicated by the baseline Fluo-5F signal), or if the series resistance became too high for us to reliably evoke action potentials. For experiments in which we performed glutamate uncaging, we also aborted recordings if the input resistance of the cell exceeded 200 MΩ. When we performed glutamate uncaging, we added the following drugs via bath application: 3.75 mM MNI-Glutamate (for uncaging), 10 µM 2-CA (to reduce presynaptic release probability and maintain quiescent recording conditions), 1 unit/mL glutamate pyruvate transaminase (which catalyzes free glutamate into -ketoglutarate), and 3 mM Sodium Pyruvate (a necessary cofactor for GPT), to the recording ACSF. In experiments where we blocked uncaging-evoked depolarization in Figure 4, we added 10 µM NBQX to the recording ACSF.
Dendritic loose-patch recordings were acquired from cortical excitatory L2/3 pyramidal cell dendrites after acquiring somatic whole-cell recordings and filling the cells with 10 µM Alexa 594. The open-tip resistance of dendritic pipettes was 8–12 MΩ, and the pipettes were filled with recording ACSF. Dendritic recordings were made with a Multiclamp 700B amplifier (Molecular Devices) in current-clamp mode. We used scanning DIC and two-photon imaging to target pipettes to dendritic branches. Before acquiring a loose-patch seal, we pushed the dendritic pipette into the dendrite until we observed a visible kink in the dendrite. For each dendritic recording, we recorded the electrical signal evoked by an AP before and after applying a brief pulse of negative pressure to achieve a loose-patch seal (Figure 4—figure supplement 1). All dendritic analyses in Figure 4 are derived from the difference between these two signals, which removes small electrical artifacts of somatic current injection and focuses our analysis on the signal arriving directly from the dendritic membrane.
The glutamate puff solution used in Figure 3 was applied with a large patch-pipette ( < 1 MΩ) connected to a picospritzer. The puff solution was composed of recording ACSF in addition to 1 mM glutamate, 20 µM CPP and 50 µM MK-801. For glutamate puff experiments, we also added 20 µM CPP and 50 µM MK-801 to the recording ACSF to block NMDARs. In some experiments, we puffed recording ACSF onto the dendrite without glutamate or NMDAR blockers to control for displacement artifacts evoked by the positive pressure (Figure 3—figure supplement 1). We used scanning DIC and two-photon imaging to bring puff pipettes near dendrites and used a short 5–15ms puff to apply glutamate to dendritic branches. The pressure of the puff was increased gradually from 5 to 20 PSI until a clear somatic depolarization and dendritic calcium signal were observed.
Two-photon imaging and glutamate uncaging
Request a detailed protocolTwo-photon imaging was performed with a custom-built microscope described previously (Carter and Sabatini, 2004). We tuned a mode-locked femtosecond laser to 840 nm to excite both Fluo-5F and Alexa 594. All calcium imaging was performed >10 min after break-in to allow fluorophores to passively diffuse into the dendritic tree. Stimulus-evoked changes in fluorescence were quantified as the change in green (Fluo-5F) fluorescence relative to a baseline period, divided by the baseline red (Alexa 594) fluorescence (ΔG/R) (Bloodgood and Sabatini, 2007). We used a custom algorithm to center the field of view on each trial (Carter and Sabatini, 2004). We measured fluorescence with a 2ms line-scan through the recording sites (at least 1 dendritic shaft and 0–3 dendritic spines depending on the experiment) and averaged over the spatial coordinates of each compartment to extract a time-varying trace (Figure 1). We waited at least 10 s in between trials to minimize photodamage and allow the calcium signal to fully return to baseline. In the experiments in Figure 6—figure supplement 2, we performed point-scans on dendritic spines to maximize the sampling rate (recorded at 125 kHz and downsampled to 8 kHz to minimize shot noise).
To perform laser photolysis of MNI-glutamate, we used a short, 500 µs femtosecond laser pulse tuned to 725 nm. For each experiment, we gradually increased laser intensity until either a 0.5 mV EPSP was generated or a clear calcium signal in the dendritic spine was evoked. We moved the uncaging spot around the dendritic spine to maximize the EPSP size as an attempt to uncage glutamate directly apposed to the postsynaptic density (Figure 1B and D).
Simultaneous voltage and calcium imaging
Request a detailed protocolInterleaved voltage- and calcium-imaging experiments in Figure 5 were conducted on modified home-built structured illumination epifluorescence microscope previously reported (Adam et al., 2019). Briefly, blue illumination was patterned by a digital micromirror device (DMD) and used for structured illumination calcium imaging via the calcium-sensitive dye Fluo-5F. Yellow illumination was also patterned by a DMD and used for structured illumination mapping of dendritic morphology via a dye JF549i. Red illumination was patterned by a holographic spatial light modulator (SLM) and used for structured illumination voltage imaging.
Laser lines from a blue laser (488 nm, 150 mW, Obis LS) and green laser (561 nm, 150 mW, Obis LS) were combined, sent through an acousto-optic modulator for amplitude control, and expanded onto a DMD (V-7000 VIS, ViALUX) for spatial modulation. The DMD was then re-imaged onto the sample via a 25 x water-immersion objective, NA 1.05, Olympus XLPLN-MP. A red laser (MLL-FN-639, 1 W, CNI lasers) was expanded onto an SLM (P1920-0635-HSP8, Meadowlark, Frederick, CO), combined with the blue and yellow lasers via a polarizing beamsplitter, and re-imaged onto the back aperture of the objective. Fluorescence was collected by the objective and imaged onto a sCMOS camera (Hamamatsu Orca Flash 4.0) with the appropriate emission filter for each indicator. Voltage-imaging recordings were acquired at a 1 kHz frame rate. Calcium imaging was recorded at 20 Hz.Two-photon imaging and reconstruction was performed with a custom-built microscope adapted to be combined with the 1 P illumination. Maximum intensity projections of z-stacks were used to form images of the dendritic arbor.
Physiology and imaging analysis
Request a detailed protocolAll analysis was performed using custom software written in MATLAB and Python. All plots with error bars show mean +/- standard error. To compute the amplitude and peak time of uEPSPs, we first averaged across trials, subtracted the baseline voltage, and applied a median filter with a 1ms window. Then, we found the maximum voltage between 0 and 25ms after uncaging and defined this as the amplitude and peak time. The amplitude of stimulus-evoked calcium signals (ΔG/R) was computed using a 10ms average of the fluorescence signal centered around the average peak time for each recording site. To compute the amplitude of dendritic electrical signals in Figure 4, we first computed the average electrical signal in a 0.2ms window surrounding the peak voltage for each trial. The amplitude of each recording was defined as the average peak voltage from the trials in the top 50th percentile signal-to-noise ratio (peak divided by baseline standard deviation). While this method is unsatisfactory for proper comparisons of amplitude across recordings, we only use these data to determine if somatic APs evoke a measurable electrical signal in the dendrites. For measurements of the ratio of from within segment or across branch point pairs (Figure 6—figure supplement 1A-C), we only considered pairs if the proximal site had a because the ratio of is not well-defined if the proximal site failed to exhibit AP-evoked calcium influx. never dropped to ∼0 and then recovered in more distal sites along the same branch. Effect sizes in Figure 7G–H were computed using Cohen’s d.
For analysis of epifluorescence voltage- and calcium-imaging data in Figure 5, we created a spike-triggered average (STA) movie comprising the average bAP-evoked signal. An initial estimate of the STA waveform for each dendritic branch was calculated by averaging the pixel values in a region of interest comprising the 5% of pixels with greatest mean brightness. Maps of the baseline (F) and spike-dependent () signals were then calculated pixel-by-pixel by linear regression of the STA movie against a normalized copy of the initial STA waveform estimate, after accounting for a pixel-dependent offset. To account for spatially variable background, watershed segmentation was applied to the map to identify regions with high spike-dependent signal. For the region surrounding pixels with the highest , a pixel-by-pixel plot of vs. was fit to a line via linear regression. The inverse of the slope was the calculated .
Morphological analysis
Request a detailed protocolAfter all experiments, we collected two-photon z-stacks of each cell to map the dendritic morphology between the soma and the recording sites. In the neurons used in Figure 6, we imaged the entire apical dendritic tree after finishing experiments to reconstruct the cells for morphological analyses and compartmental modeling in NEURON. For every recording site, we measured the on-path distance from the soma using maximum z-stack projections of each cell. The distance for each recording site is a slight underestimate of the true on-path distance from the soma because we ignored changes along the z-axis. All imaging was conducted close to the surface of the slice and the maximal angle between the soma and any recording site was 25°, so errors in measured distance were <10%.
We developed a metric for branch complexity that fully accounts for the extent of nearby dendritic branching, independent of the structure of the dendritic tree. First, we defined the term as the number of dendritic branches at a given on-path distance from a recording site (Figure 6G). always starts as 2 (one on each side of the recording site) and increases by 1 after every branch point. This is distinct from Scholl analysis as it measures distance based on the on-path distance rather than absolute distance. Next, we performed point-wise multiplication between and an exponential filter to discount for dendritic branches that are further away from the recording site. Finally, we computed branch complexity as , which is largest at sites that are closest to large numbers of dendritic branches and is inversely proportional to electrotonic input resistance.
The true electrotonic input resistance is analytically determined by the cable length constant, which is not feasible to measure in thin caliber L2/3 pyramidal cell dendrites. As an alternative, we estimated how much membrane is simultaneously depolarized by a bAP as an estimate of how much of the dendritic tree can affect the depolarization at a particular site. We used a rapid point-scan imaging approach that allowed us to measure in dendritic spines at >8 kHz (Figure 6—figure supplement 2A). We plotted the latency to peak influx (measured from the derivative of fluorescence) as a function of the distance from the soma and used linear regression to estimate the conduction velocity (Figure 6—figure supplement 2A-B). We found that the conduction velocity of bAPs in L2/3 dendrites is 154.9 µm/ms. In the same cells used to measure conduction velocity, we plotted the somatic AP waveforms and measured their full-width half-max (Figure 6—figure supplement 2C-D). The somatic AP has a duration of about 0.94ms (Figure 6—figure supplement 2C-D). To compute the spatial width of a propagating bAP, we multiplied the conduction velocity by the duration (154.9 µm/ms * 0.94ms = 145.6 µm). This estimate of bAP width led to a choice of length constant that maximized the correlation between branch complexity and input resistance (Figure 6, Figure 6—figure supplement 3).
Biophysical compartmental modeling
Request a detailed protocolSimulations used in Figures 6F and 7, and Figure 7—figure supplement 1 were performed in the NEURON simulation environment (Carnevale and Hines, 2006) using Python (Version 3.7). We reconstructed all 8 neurons used in Figure 6 with maximum z-projections of the Alexa 594 fluorescence signal. Each compartmental model was composed of a cylindrical soma with length and diameter equal to 12 µm, a 30 µm axon initial segment (AIS) with a 2 µm diameter connected to the soma, a 300 µm axon with a 1 µm diameter connected to the AIS, and a dendritic tree connected to the soma that with identical branch structure to the 8 neurons used in Figure 6 and a uniform diameter of 1 µm. Our two-photon imaging data did not permit precise estimation of dendritic diameter but from inspection all dendrites had similar diameters. In some experiments, we performed simulated cuts of the dendritic tree in NEURON. For each cell, we removed branches close to the low recording site until the input resistance of the high and low sites were comparable.
We used standard passive parameters: Cm = 1 µF/cm2, Rm = 7000 Ω-cm2, Ri = 100 Ω-cm. Active conductances were based on a previous model of L2/3 pyramidal cells that was fit to in vivo somatic and dendritic recordings (Smith et al., 2013). Dendritic active conductances were completely uniform, such that all differences across dendrites were determined by the distance from the soma and the dendritic branch structure. Each model contained the following channels with given conductance density in pS/µm2: voltage-gated sodium channels (axon: 170, AIS: 2,550 soma: 85, dendrites: 85), voltage-gated potassium channels (axon: 33, AIS: 100, soma: 100, dendrites: 3), M-Type potassium channels (soma: 2.2, dendrites: 1), calcium-dependent potassium channels (soma: 3, dendrites: 3), high voltage-activated calcium channels (soma: 0.5, dendrites: 0.5), and low voltage-activated calcium channels (soma: 3, dendrites: 1.5). Additionally, we added an A-Type potassium channel (Migliore et al., 1999) to the soma and dendrites with a variable conductance density that was set for each cell independently. For experiments in Figures 6F, 7A–E and G–H, and Figure 7—figure supplement 1, we used a Nelder-Mead simplex algorithm (from the SciPy toolbox in Python) to set the A-Type potassium channel conductance density such that backpropagating APs peaked at –10 mV in the high recording site within each pair (values for each site plotted in left-most column of points in Figure 7F). The computed density was used for the soma and entire dendritic tree. We chose to set the AP peak at –10 mV because it is consistent with the AP amplitude measured with whole-cell recordings from L2/3 pyramidal cell dendrites > 100 µm from the soma (Waters et al., 2003; Larkum et al., 2007; Smith et al., 2013). The same values for A-Type potassium channel density were used in experiments with simulated cut dendrites. The values shown in Figure 7F use the same algorithm to select A-Type potassium channel density for high and low recording sites independently, for both the intact and cut dendritic branch structure.
Conductance simulations
Request a detailed protocolWe computed the voltage-dependent open probability and time constant of VGCCs and NMDARs using published data (Jahr and Stevens, 1990; Reuveni et al., 1993). For VGCCs, we used the following voltage-dependent forward and backward rate constants ( and , respectively) of the activation and inactivation gates (in mV):
To compute open probability and time constant for the activation (denoted ‘’) and inactivation (denoted ‘’) gates of VGCCs, we used the following equations: , . Because the time constant of the inactivation gate is so slow (the minimum time constant is 170ms between –80 and +40 mV), it acts as a constant, so we did not plot it in Figure 8B, although it was accounted for in the simulations in Figure 8C–E.
For NMDARs (denoted ‘’), we calculated the voltage-dependent open-probability and time constant as and based on the measured on and off rates of the Mg2+ block (Jahr and Stevens, 1990), with [Mg2+] = 1 µM:
For numerical integration in Figure 8C, we used Euler’s method to compute changes in the state of each gate with . We used the following differential equations:
The open probability of NMDARs is equal to , and the open probability of VGCCs is equal to . To convert the open probability into a voltage-dependent calcium current, we used a modified version of the Goldman-Hodgkin-Katz current equation with and :
We used quadratic voltage depolarizations to mimic the shape of APs. For each stimulus, we used a preset amplitude () and duration ). Then, we solved the following equation for : , such that the stimulus would have a maximum height of along the domain . Then, we added the equation: to the baseline voltage of –70 mV.
To determine the integrated calcium influx in Figure 8E, we integrated the calcium current evoked by the voltage step after subtracting the baseline current (NMDARs are partially open at rest).
Data availability
All data and code is posted on the Harvard Dataverse (https://doi.org/10.7910/DVN/ZHNKGE).
-
Harvard DataverseDendritic branch structure compartmentalizes voltage-dependent calcium influx in cortical layer 2/3 pyramidal cells.https://doi.org/10.7910/DVN/ZHNKGE
References
-
Signaling to the nucleus by an L-type calcium channel-calmodulin complex through the MAP kinase pathwayScience (New York, N.Y.) 294:333–339.https://doi.org/10.1126/science.1063395
-
Normalization of Ca2+ signals by small oblique dendrites of CA1 pyramidal neuronsThe Journal of Neuroscience 23:3243–3250.
-
LTP is accompanied by an enhanced local excitability of pyramidal neuron dendritesNature Neuroscience 7:126–135.https://doi.org/10.1038/nn1178
-
Plasticity compartments in basal dendrites of neocortical pyramidal neuronsThe Journal of Neuroscience 26:12717–12726.https://doi.org/10.1523/JNEUROSCI.3502-06.2006
-
Dendrites: bug or feature?Current Opinion in Neurobiology 13:372–383.https://doi.org/10.1016/s0959-4388(03)00075-8
-
A quantitative description of NMDA receptor-channel kinetic behaviorThe Journal of Neuroscience 10:1830–1837.
-
L-Type calcium channels are required for one form of hippocampal mossy fiber LTPJournal of Neurophysiology 79:2181–2190.https://doi.org/10.1152/jn.1998.79.4.2181
-
Molecular mechanisms of dendrite stabilityNature Reviews. Neuroscience 14:536–550.https://doi.org/10.1038/nrn3486
-
Dendritic spikes in apical dendrites of neocortical layer 2/3 pyramidal neuronsThe Journal of Neuroscience 27:8999–9008.https://doi.org/10.1523/JNEUROSCI.1717-07.2007
-
A synaptically controlled, associative signal for Hebbian plasticity in hippocampal neuronsScience (New York, N.Y.) 275:209–213.https://doi.org/10.1126/science.275.5297.209
-
Role of an A-type K+ conductance in the back-propagation of action potentials in the dendrites of hippocampal pyramidal neuronsJournal of Computational Neuroscience 7:5–15.https://doi.org/10.1023/a:1008906225285
-
Spine Ca2+ signaling in spike-timing-dependent plasticityThe Journal of Neuroscience 26:11001–11013.https://doi.org/10.1523/JNEUROSCI.1749-06.2006
-
Computational subunits in thin dendrites of pyramidal cellsNature Neuroscience 7:621–627.https://doi.org/10.1038/nn1253
-
Stepwise repolarization from Ca2+ plateaus in neocortical pyramidal cells: evidence for nonhomogeneous distribution of HVA Ca2+ channels in dendritesThe Journal of Neuroscience 13:4609–4621.
-
Experience-Dependent Development of Dendritic Arbors in Mouse Visual CortexThe Journal of Neuroscience 40:6536–6556.https://doi.org/10.1523/JNEUROSCI.2910-19.2020
-
Nonlinear [Ca2+] signaling in dendrites and spines caused by activity-dependent depression of Ca2+ extrusionThe Journal of Neuroscience 26:8183–8194.https://doi.org/10.1523/JNEUROSCI.1962-06.2006
-
Spike timing, calcium signals and synaptic plasticityCurrent Opinion in Neurobiology 12:305–314.https://doi.org/10.1016/s0959-4388(02)00325-2
-
Competitive Hebbian learning through spike-timing-dependent synaptic plasticityNature Neuroscience 3:919–926.https://doi.org/10.1038/78829
-
Activity-dependent action potential invasion and calcium influx into hippocampal CA1 dendritesScience (New York, N.Y.) 268:297–300.https://doi.org/10.1126/science.7716524
-
Propagation of action potentials in dendrites depends on dendritic morphologyJournal of Neurophysiology 85:926–937.https://doi.org/10.1152/jn.2001.85.2.926
-
Supralinear Ca2+ influx into dendritic tufts of layer 2/3 neocortical pyramidal neurons in vitro and in vivoThe Journal of Neuroscience 23:8558–8567.
-
A comparison of neuronal population dynamics measured with calcium imaging and electrophysiologyPLOS Computational Biology 16:e1008198.https://doi.org/10.1371/journal.pcbi.1008198
-
Backpropagation of physiological spike trains in neocortical pyramidal neurons: implications for temporal coding in dendritesThe Journal of Neuroscience 20:8238–8246.
-
Spatial compartmentalization and functional impact of conductance in pyramidal neuronsNature Neuroscience 7:961–967.https://doi.org/10.1038/nn1305
-
Orientation selectivity and the functional clustering of synaptic inputs in primary visual cortexNature Neuroscience 19:1003–1009.https://doi.org/10.1038/nn.4323
-
Plasticity of calcium channels in dendritic spinesNature Neuroscience 6:948–955.https://doi.org/10.1038/nn1112
Article and author information
Author details
Funding
National Institute of Neurological Disorders and Stroke (F31NS113353)
- Andrew T Landau
National Institute of Neurological Disorders and Stroke (R37NS046579)
- Bernardo L Sabatini
National Institute of Mental Health (1RF1MH117042-01)
- Adam E Cohen
Defense Advanced Research Projects Agency (Vannevar Bush Faculty Fellowship)
- Adam E Cohen
Brain Research Foundation (Scientific Innovation)
- Adam E Cohen
Harvard Medical School (Harvard Brain Initiative)
- Adam E Cohen
- Bernardo L Sabatini
Life Sciences Research Foundation (Merck Awardee)
- J David Wong-Campos
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank all members of the Sabatini Lab for comments and advice, in particular R Hakim, E Sayed, P Capelli, S Melzer, K Reinhold, T Kula, A Granger, and S Kim for suggestions on the manuscript. We also thank J Reggiani, N Hahn, W Regehr, C Harvey, and B Bean for helpful suggestions and critical insights. This work was supported by 5F31NS113353-03 from NINDS to ATL, R37NS046579 from NINDS to BLS, NIH R01 1RF1MH117042-01 to AEC, a Vannevar Bush Faculty Fellowship and a Brain Research Foundation Scientific Innovation Award to AEC, and a grant from the Harvard Brain Initiative to BLS and AEC. JDW is a Merck Awardee of the Life Sciences Research Foundation.
Ethics
All procedures involving animals were in accordance with the National Institutes of Health Guide for the care and use of laboratory animals and were approved by the Harvard University Institutional Animal Care and Use Committee (IACUC) (Protocol #: IS00000571-3).
Copyright
© 2022, Landau et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,703
- views
-
- 420
- downloads
-
- 22
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Neuroscience
Williams syndrome (WS; OMIM#194050) is a rare disorder, which is caused by the microdeletion of one copy of 25–27 genes, and WS patients display diverse neuronal deficits. Although remarkable progresses have been achieved, the mechanisms for these distinct deficits are still largely unknown. Here, we have shown that neural progenitor cells (NPCs) in WS forebrain organoids display abnormal proliferation and differentiation capabilities, and synapse formation. Genes with altered expression are related to neuronal development and neurogenesis. Single cell RNA-seq (scRNA-seq) data analysis revealed 13 clusters in healthy control and WS organoids. WS organoids show an aberrant generation of excitatory neurons. Mechanistically, the expression of transthyretin (TTR) are remarkably decreased in WS forebrain organoids. We have found that GTF2IRD1 encoded by one WS associated gene GTF2IRD1 binds to TTR promoter regions and regulates the expression of TTR. In addition, exogenous TTR can activate ERK signaling and rescue neurogenic deficits of WS forebrain organoids. Gtf2ird1-deficient mice display similar neurodevelopmental deficits as observed in WS organoids. Collectively, our study reveals critical function of GTF2IRD1 in regulating neurodevelopment of WS forebrain organoids and mice through regulating TTR-ERK pathway.
-
- Computational and Systems Biology
- Neuroscience
Hypothalamic kisspeptin (Kiss1) neurons are vital for pubertal development and reproduction. Arcuate nucleus Kiss1 (Kiss1ARH) neurons are responsible for the pulsatile release of gonadotropin-releasing hormone (GnRH). In females, the behavior of Kiss1ARH neurons, expressing Kiss1, neurokinin B (NKB), and dynorphin (Dyn), varies throughout the ovarian cycle. Studies indicate that 17β-estradiol (E2) reduces peptide expression but increases Slc17a6 (Vglut2) mRNA and glutamate neurotransmission in these neurons, suggesting a shift from peptidergic to glutamatergic signaling. To investigate this shift, we combined transcriptomics, electrophysiology, and mathematical modeling. Our results demonstrate that E2 treatment upregulates the mRNA expression of voltage-activated calcium channels, elevating the whole-cell calcium current that contributes to high-frequency burst firing. Additionally, E2 treatment decreased the mRNA levels of canonical transient receptor potential (TPRC) 5 and G protein-coupled K+ (GIRK) channels. When Trpc5 channels in Kiss1ARH neurons were deleted using CRISPR/SaCas9, the slow excitatory postsynaptic potential was eliminated. Our data enabled us to formulate a biophysically realistic mathematical model of Kiss1ARH neurons, suggesting that E2 modifies ionic conductances in these neurons, enabling the transition from high-frequency synchronous firing through NKB-driven activation of TRPC5 channels to a short bursting mode facilitating glutamate release. In a low E2 milieu, synchronous firing of Kiss1ARH neurons drives pulsatile release of GnRH, while the transition to burst firing with high, preovulatory levels of E2 would facilitate the GnRH surge through its glutamatergic synaptic connection to preoptic Kiss1 neurons.