Frequency- and spike-timing-dependent mitochondrial Ca2+ signaling regulates the metabolic rate and synaptic efficacy in cortical neurons
Abstract
Mitochondrial activity is crucial for the plasticity of central synapses, but how the firing pattern of pre- and postsynaptic neurons affects the mitochondria remains elusive. We recorded changes in the fluorescence of cytosolic and mitochondrial Ca2+ indicators in cell bodies, axons, and dendrites of cortical pyramidal neurons in mouse brain slices while evoking pre- and postsynaptic spikes. Postsynaptic spike firing elicited fast mitochondrial Ca2+ responses that were about threefold larger in the somas and apical dendrites than in basal dendrites and axons. The amplitude of these responses and metabolic activity were extremely sensitive to the firing frequency. Furthermore, while an EPSP alone caused no detectable Ca2+ elevation in the dendritic mitochondria, the coincidence of EPSP with a backpropagating spike produced prominent, highly localized mitochondrial Ca2+ hotspots. Our results indicate that mitochondria decode the spike firing frequency and the Hebbian temporal coincidences into the Ca2+ signals, which are further translated into the metabolic output and most probably lead to long-term changes in synaptic efficacy.
Introduction
For the neuronal circuit to function properly, the energy demand in all compartments of the individual neurons needs to be precisely matched by local, primarily mitochondrial (Celsi et al., 2009; Rangaraju et al., 2019), ATP production. During periods of enhanced neuronal activity, mitochondria accelerate ATP production by allowing the cytosolic Ca2+ elevations to propagate into the mitochondrial matrix (Ashrafi et al., 2020; Díaz-García et al., 2021). When cytosolic [Ca2+]i rises, Ca2+ ions, powered by the steep mitochondrial membrane potential, flow into the mitochondrial matrix via the Ca2+ uniporter MCU (Baughman et al., 2011; De Stefani et al., 2011), and are then extruded back into the cytosol by the mitochondrial Na+/Ca2+ antiporter NCLX (Palty et al., 2010). The mitochondrial Ca2+ elevation increases ATP production by activating at least three Krebs cycle enzymes (De Stefani et al., 2016; Wescott et al., 2019). When Ca2+ signaling is disrupted, mitochondria in presynaptic terminals fail to maintain the stable ATP concentration during enhanced activity periods (Ashrafi et al., 2020; Giorgio et al., 2013; Glancy and Balaban, 2012). The link between the distinctive types of the neuronal electrical activity, mitochondrial Ca2+ signaling and metabolism in other neuronal compartments is poorly understood, however.
Recent evidence using photo-uncaging of glutamate on spines of cultured hippocampal neurons indicates that dendritic mitochondria play a critical role in long-term regulation of synaptic strength (Rangaraju et al., 2019). Thus, local inhibition of the mitochondria by a phototoxic protein abolishes the synaptic plasticity in the affected dendritic segment. Under physiological conditions, the plastic changes in synaptic efficacy are believed to be controlled by the relative timing of pre- and postsynaptic action potentials (APs) (Bi and Poo, 1998; Holtmaat and Svoboda, 2009; Markram et al., 1997b). Spike-timing–dependent plasticity (STDP) relies on activation of the postsynaptic NMDA receptors, which are Ca2+ permeable and require both glutamate and depolarization to open (Nowak et al., 1984). Mitochondrial handling of Ca2+ ions is known to have a significant effect on cytosolic Ca2+ dynamics (Rizzuto et al., 2004; Szabadkai and Duchen, 2008). The mitochondrial depletion, however, had no effect on glutamate evoked cytosolic Ca2+ dynamics (Rangaraju et al., 2019), indicating that these organelles play a yet unexplored downstream role in the cascade of reactions leading to STDP.
Here, using whole-cell electrical recordings from Layer 5 pyramidal neurons in cortical slices and fluorescence imaging of cytosolic and mitochondrial Ca2+ indicators, we show that single or few spikes trigger rapidly rising and decaying mitochondrial Ca2+ elevations in all neuronal compartments, with kinetics similar to cytosolic Ca2+ transients. Our evidence indicates that the mitochondria’s Ca2+ signaling and metabolic rate depend critically on spike firing frequency. We further report that, in dendrites, the coincidence of unitary EPSP and a backpropagating action potential produces a localized [Ca2+]m transient, which, in addition to enhancing local ATP synthesis, could play a role in STDP.
Results
Spike-elicited mitochondrial Ca2+ transients
To determine how neuronal electrical activity affects mitochondrial Ca2+ dynamics, we performed somatic whole-cell recordings from L5 pyramidal neurons expressing the mitochondria-targeted, neuron-specific Ca2+ indicator, mitoGCaMP6m. The intracellular solution was supplemented with the cytosolic Ca2+ indicator, Fura-2, that could be excited separately from mitoGCaMP6m. Series of thin, high-resolution optical sections over the vertical extent of the neuron, representing Fura-2 and mitoGCaMP6m fluorescence elicited by two-photon excitation at 760 and 960 nm, respectively, were used to reconstruct the morphology of the neuron and to reveal the position of the individual mitochondria within its soma and processes. To improve the temporal resolution of the optical signals, dynamic fluorescence measurements were obtained from the smaller regions of interest in soma, axon, and dendrites. In a typical experiment, two APs elicited by the injection of two brief current pulses via the patch pipette caused cytosolic and mitochondrial Ca2+ elevations in the soma of an L5 cell (Figure 1a). Although the cytosolic Ca2+ signals had relatively even intensity, the change in mitoGCaMP6m fluorescence occurred at ‘hotspots’, each representing an individual mitochondrion. The mitochondrial Ca2+ elevations began with a short delay after the beginning of the spike train, and they were observed only in the electrically active neurons. At the same time, the fluorescence of the nearby mitoGCaMP6m expressing mitochondria belonging to the non-active cells did not change (Figure 1—figure supplement 1).
The AP-elicited mitochondrial Ca2+ elevations required Ca2+ influx from the extracellular space. It is likely that Ca2+ ions enter via the voltage-gated Ca2+ channels, as bath application of Cd2+ (200 µM) inhibited both cytosolic and mitochondrial Ca2+ transients (Figure 1—figure supplement 2). Previous studies have shown that Ca2+ depletion of the endoplasmic reticulum (ER) in cultured neurons does not influence mitochondrial Ca2+ elevations (Ashrafi et al., 2020). Consistent with these results, we found that blockade of ER Ca2+ channels by Ryanodine (100 µM) and Dantrolene (100 µM) (Figure 1—figure supplement 3) had no significant effect on mitochondrial Ca2+ transients.
An increase in the number of spikes produced progressively larger mitoGCaMP6m transients in the soma (Figure 1b) and the apical dendrite (Figure 1—figure supplement 4) with a progressively slower rising phase. Thus, the mito-Ca2+ signals elicited by 50 APs grew throughout the spike train duration, reaching a peak ΔF/F value of 59% ± 11% (n = 29) and 54% ± 11% (n = 22) for soma and apical dendrite, respectively, at 120 ± 18ms (n = 10) after the train end. In contrast with previous reports (Ashrafi et al., 2020; Lin et al., 2019), mito-Ca2+ transients elicited by 2–20 spikes decayed as rapidly or even faster than cytosolic Ca2+ transients. However, following large elevations, [Ca2+]m remained high for tens of seconds, reaching the resting level long after the cytosolic Ca2+ concentration completely recovered. We systematically examined the relationship between the peak amplitude and the decay rate of the mito-Ca2+ transients in 42 pyramidal neurons (Figure 1c, Figure 1—figure supplement 5). Decay of smaller transients (peak ΔF/F < 40%) was always monoexponential, with τ = 0.86 ± 0.11 s (n = 25) as well as decay of some middle-sized (ΔF/F 40–75%) transients (τ = 0.82 ± 0.10 s, n = 10). The decay of other middle size and large (ΔF/F > 75%) transients was double exponential, with τ1 = 0.86 ± 0.13 s and τ2 = 22.0 ± 2.85 s (n = 23). The decay of cytosolic Ca2+ transients always followed a single exponential time course which was characterized by τ = 1.03 ± 0.14 s (n = 22) for small and τ = 1.19 ± 0.4 s (n = 20, p = 0.43) for large transients.
To elucidate the differences in mitochondrial Ca2+ signaling between distinct neuronal compartments, we monitored mitoGCaMP6m fluorescence in ~10 µm long regions of interest in basal, apical dendrites, and axon initial segments (AISs) during trains of five APs (Figure 1d). While the amplitude and time course of the cytosolic Ca2+ responses in all these compartments were similar, the magnitude of mitochondrial signals was remarkably polar. The amplitude of the mitochondrial Ca2+ elevations in apical dendrite was as prominent as in the soma. In contrast, in the AIS and thin basal dendrites, the mito-Ca2+ responses were dramatically smaller.
We next sought to evaluate the relationship between the number of APs and the mean peak amplitude of the mitoGCaMP6m transients in the soma, AIS, apical and basal dendrites of 29 neurons (Figure 1e). The steepness of this relationship was greater in soma and apical dendrites (1.3% and 1.1% ΔF/F per spike, respectively) compared with AIS and basal dendrites (0.4% and 0.3% ΔF/F per spike, respectively). The compartmental differences in mitochondrial signals were not due to a different magnitude of cytosolic Ca2+ elevations. The pattern of [Ca2+]i during the neuronal activity is known to be complex. However, the differences in peak cytosolic Ca2+ levels were subtle, and they poorly correlated with the magnitude of mitochondrial Ca2+ elevations. For example, the peak ΔF/F amplitude of Fura-2 transients elicited by twenty spikes was 26% ± 3.6% (n = 22) for soma, 37% ± 4% (n = 13) for apical dendrite, 38% ± 4% (n = 13) for basal dendrites, and 25% ± 4% in the AIS (n = 7). The compartmental differences in magnitude of the mitoGCaMP6m transients could, at least partially, be explained by the lower expression level of fluorescence probe in the mitochondria localized within the thinner neuronal processes. This seems to be unlikely, however, since, in all neuronal compartments, maximal mitoGCaMP6m fluorescence of the individual mitochondria following prolonged depolarization of the neuronal membrane was similar (Figure 1—figure supplement 6). Our evidence, therefore, points to the existence of an as-yet-unidentified mechanism that differentially regulates the mitochondrial Ca2+ entry in distinct neuronal compartments.
Frequency-dependent amplification of mitochondrial Ca2+ elevations
Next, we examined whether mitochondrial Ca2+ elevations are sensitive to firing frequency. Figure 2a shows an optical recording from the soma of a representative L5 neuron in which cytosolic and mitochondrial Ca2+ transients were elicited by trains of twenty APs at 20, 50, and 100 Hz. As spike frequency increased, the rise of cytosolic Ca2+ transients became progressively faster and their peak amplitude modestly increased (Figure 2—figure supplement 1). In contrast, the mitochondrial Ca2+ signals showed a very different frequency dependence. Firing at a frequency of 20 Hz or lower elicited a minimal elevation in [Ca2+]m, whereas the response to spikes at a frequency of 50 Hz or higher was dramatically greater. The steep frequency dependence of the mitochondrial signals was observed in all neuronal compartments, including somas, apical, and basal dendrites, making it unlikely that it reflects the frequency-dependent failure of AP backpropagation (Spruston et al., 1995). The frequency-dependent amplification of mitochondrial Ca2+ transients was observed in all 21 neurons tested with either 50, 20, or 5 APs (Figure 2b). The mean ratio of peak amplitudes of the mitochondrial Ca2+ transients elicited by trains of spikes at 50 and 20 Hz was larger when neurons were subjected to longer (3.36 ± 0.46 times, n = 22 for 50 APs) than to shorter (1.74 ± 0.18 times, n = 5 for 5 APs) trains (Figure 2c). Systematic varying of AP frequency in a range from 10 to 100 Hz revealed that the peak amplitude of mito-Ca2+ transients behaved as Bolzmannian function of the frequency, with mean half-amplitude of 38 Hz and steepness of 22 Hz–1 (n = 12).
The decay time course of the mitochondrial Ca2+ transients elicited by high-frequency spike trains was significantly slower than those produced by the same number of spikes at a lower frequency (Figure 2—figure supplement 2). Thus, trains of fifty APs at 50 Hz elicited the extremely slowly decaying mito-Ca2+ transients (τ > 10 s) in 11/15 somas and 7/10 apical dendrites. In contrast, all but one mito-Ca2+ transient produced in the same neurons by fifty APs at 20 Hz decayed rapidly.
Frequency-dependent acceleration of the mitochondrial NAD(P)H metabolism
To determine whether the frequency-dependent amplification of Ca2+ signals in the mitochondria affects their metabolic activity, we monitored the changes in NAD(P)H autofluorescence elicited by trains of brief, just suprathreshold antidromic stimuli at a different frequency (Figure 3a). Whole-cell recording from a single pyramidal neuron within the region of interest was obtained to tune the stimuli intensity such that each stimulus would elicit only one AP. In cortical neurons, NAD(P)H signals primarily reflect changes in mitochondrial NAD(P)H pool (Díaz-García et al., 2021). In response to electrical stimulation, we observed a negative deflection in the NAD(P)H autofluorescence (‘dip’) which indicates an increased rate of electron transfer reflected in NAD(P)H oxidation, followed by a positive transient (‘overshoot’) which indicates Krebs-cycle-dependent replenishment of NAD(P)H pool. At higher stimulation frequency, the magnitude of both dip and overshoot of the NAD(P)H signals were enhanced, consistent with the previous reports that the rates of NAD(P)H oxidation and synthesis are dependent on the Ca2+ level in the mitochondrial matrix (Díaz-García et al., 2021). Spatio-temporal analysis of the NAD(P)H autofluorescence dynamics at two stimulation frequencies (Figure 3b) revealed that the changes in the fluorescence were spatially restricted to Layer 5 of a single cortical column and that the frequency-dependent amplification of both dip and overshoot amplitude was relatively uniform within the stimulated region. A comparison of NAD(P)H autofluorescence changes in response to a train of 50 stimuli at 20 Hz and 50 Hz in ten cortical slices (Figure 3c) revealed the same frequency dependence as with [Ca2+]m. The dip amplitude increased from –15 ± 2 a.u. at 20 Hz to –21 ± 2 a.u. at 50 Hz (n = 15 ROIs, p < 0.001), the overshoot’s peak amplitude increased from 8 ± 2 a.u. to 19 ± 3 a.u. (p < 0.001) and the overshoot area increased from 156 ± 46 a.u.∙s to 441 ± 71 a.u.∙s (p < 0.001), respectively. Because the glial responses might partially contaminate NAD(P)H signals obtained in the cortical slices, we tested the frequency dependence of NAD(P)H responses in the in Stratum pyramidale of CA1 area of the hippocampus (Figure 3—figure supplement 1) which predominantly contains neuronal cell bodies. As in the neocortex, NAD(P)H signals in response to trains of stimuli delivered to the Stratum oriens were significantly enhanced at the higher stimulation frequency.
Localized dendritic [Ca2+]m elevations elicited by the coincidence of postsynaptic AP and EPSP
We next sought to elucidate how synaptic activity affects mitochondrial Ca2+ dynamics. After filling the cell for ~20 min to allow the diffusion of Fura-2 into the dendrites, we positioned the bipolar electrode close to an apical or basal dendrite. Delivery of a single brief stimulus (0.1ms), with its amplitude adjusted to keep the subsequent EPSP below the threshold for postsynaptic spike generation produced no detectable cytosolic or mitochondrial Ca2+ signals (n = 7). The single AP elicited by brief somatic current pulse injection, however, produced a cytosolic but no mitochondrial Ca2+ response in the dendrites. Remarkably, the coincidence of the EPSP and backpropagating AP synergized to elicit a robust cytosolic and mitochondrial Ca2+ response (Figure 4a). While most dendritic mitochondria were silent, the EPSP and AP coincidence created single localized mitochondrial Ca2+ “hotspots” with peak ΔF/F amplitude of 6.4% ± 0.8% (n = 8) in the dendritic regions of interest. We interpreted the appearance of these hotspots as evidence for a highly restricted, probably single mitochondrial Ca2+ elevation, in the vicinity of the active spine. The unitary character of the mitochondrial signals under this experimental paradigm can be explained by a spatial sparseness of the synapses formed by the presynaptic fiber on the dendrites of the postsynaptic cortical cell (Markram et al., 1997a) such that only one out of a few currently active spines could be found in the relatively short segment of a dendritic branch that we examined. The cytosolic Ca2+ elevation, measured in the same hotspot at which the mitoGCaMP6m signal was detected, was not significantly larger than in the nearby dendrite. The failure to see the cytosolic Ca2+ “hotspots” is, most probably, due to temporal and amplitude resolution of our optical recording that was insufficient to reveal Ca2+ elevation in the tiny volume between the spine neck and the mitochondrion during the fast single spine Ca2+ transient (Miyazaki and Ross, 2017; Svoboda et al., 1996). We extended our analysis by measuring the cytosolic and mitochondrial Ca2+ signals elicited by 20 unpaired APs and APs paired with EPSP. In both apical and basal dendrites, the paired APs produced a significantly larger mitochondrial signal. In contrast, the cytosolic Ca2+ elevation amplitude was not significantly different for unpaired and paired stimulation (Figure 4b and c). Hence, the peak amplitude of the mitoGCaMP6m transients elicited by paired APs was 2.73 ± 0.32 (n = 8) and 2.11 ± 0.32 (n = 9) times higher than of those evoked by the unpaired APs in the basal and apical dendrites, respectively (Figure 4d). The ability of postsynaptic neurons to generate an AP was crucial for triggering the dendritic mito-Ca2+ transients. Thus, intracellular dialysis with a solution containing a blocker of voltage-gated Na+ channels, QX-314, which prevented the firing of the postsynaptic cells (Connors and Prince, 1982) while producing only a minor effect on the EPSP generation, dramatically reduced the amplitude of synaptically evoked mito-Ca2+ transients (Figure 4—figure supplement 1). As with single subthreshold EPSPs, even the significantly larger EPSPs evoked in QX-314 dialyzed neurons by a train of strong synaptic stimuli, in the absence of AP, produced no significant Ca2+ elevation in dendritic mitochondria.
It is well established that the main route of Ca2+ entry into the spines is via the NMDARs (Mainen et al., 1999; Miyazaki and Ross, 2017), which require glutamate and depolarization to relieve the Mg2+ block (Nowak et al., 1984). Blockade of NMDARs by bath applied APV (50 µM) almost completely and reversibly abolished the mitoGCaMP6m transients evoked by the paired APs in the dendrites (Figure 4—figure supplement 2).
Discussion
Using simultaneous electrical recordings, cytosolic and mitochondrial Ca2+ imaging in L5 pyramidal neurons, we found that relatively rare, singular spike firing, the firing mode associated with cortical circuit-based information processing in vivo (Brecht and Sakmann, 2002), produces fast-rising, rapidly decaying mitochondrial Ca2+ transients in all neuronal compartments. In contrast to a recent report of ‘loose’ coupling of mitochondrial Ca2+ transients to neuronal activity and cytosolic Ca2+ transients in in vivo cortical neurons (Lin et al., 2019), our data suggest a tight, causative relationship between neuronal electrical activity, cytosolic and mitochondrial Ca2+ levels.
Several studies (Ashrafi et al., 2020; Devaraju et al., 2017; Kwon et al., 2016; Lewis et al., 2018) have revealed the tens-of-seconds long mitochondrial Ca2+ transients in presynaptic terminals of cultured central neurons. While this extremely slow rate of Ca2+ clearance may reflect compartmental specialization of mitochondrial Ca2+ handling at the presynaptic sites, it seems more likely to be due to excessively intense stimulation. Indeed, small mitochondrial Ca2+ signals in axons of cultured hippocampal neurons elicited by 1–5 antidromic stimuli decayed rapidly (Gazit et al., 2016). Our results indicate that the prolonged mitochondrial Ca2+ elevations occur only following the long period of robust high-frequency firing. Although such intense firing is not typically observed in cortical pyramidal cells under physiological conditions, it may occur during various neurological diseases (Makinson et al., 2017; Zott et al., 2019) contributing to the mitochondrial Ca2+ overload and disease progression.
Our evidence indicates that the summation of the unitary AP-evoked mitochondrial Ca2+ transients steeply depends on firing frequency. Thus, in contrast to cytosolic Ca2+ signaling, which is less sensitive to firing rate, the mitochondrial Ca2+ elevations were strongly amplified at firing frequencies of >40 Hz. As expected, the enhanced mitochondrial Ca2+ uptake accelerated the rate of NAD(P)H production and consumption via the tricarboxylic acid cycle and the electron transport chain, respectively (Díaz-García et al., 2021), thereby upregulating the ATP synthesis.
The finding that action potential induced mitochondrial Ca2+ transients are about threefold greater in the soma and apical dendrites than in proximal axon and basal dendrites, most likely, indicates the compartmental difference in expression of mitochondrial Ca2+ channel, MCU. Alternatively, the subcellular differences in mitochondrial Ca2+ handling may reflect region-specific MCU molecular tuning (Ashrafi et al., 2020; Patron et al., 2019) or mitochondrial morphology (Lewis et al., 2018). From the functional viewpoint, the greater mitochondrial Ca2+ transients in thick neuronal processes might be necessary to compensate for the higher energetic cost of AP generation (Attwell and Laughlin, 2001).
A recent study suggests that the metabolic activity of mitochondria plays a pivotal role in synaptic plasticity (Rangaraju et al., 2019). How neuronal activity is linked to this process is poorly understood, however. Our results suggest that mitochondria can detect the Hebbian time coincidences between the pre- and postsynaptic spikes (Markram et al., 1997b). The resultant Ca2+ elevations in the mitochondrial matrix could be an essential part of the cascade of events underlying spike-time-dependent synaptic plasticity. This cascade might involve the initiation of fission of the dendritic mitochondria, as has been proposed for chemically induced NMDAR-dependent LTP in hippocampal neuronal culture (Divakaruni et al., 2018), although it remains unclear whether the LTP associated burst of fission events occurs at a physiologically relevant time scale. The involvement of the frequency-dependent mitochondrial Ca2+ signaling, most probably, explains the observed repetition rate requirement for LTP induction (Inglebert et al., 2020; Lisman and Spruston, 2005; Sjöström et al., 2001). Indeed, the spike-timing-dependent potentiation of cortical synapses, in addition to precise spike timing, requires a sufficiently high repetition frequency (Sjöström et al., 2001). Our data indicate that the unique ability of the mitochondria to decode firing frequency and Hebbian timing code of neuronal activity make this organelle a long-thought link between the firing pattern, metabolism, and plasticity.
Materials and methods
Experimental animals
Request a detailed protocolAll experiments were approved by the Animal Care and Use Committee of Ben Gurion University of the Negev. C57BL/6 mice were obtained from Envigo (Israel).
Viral constructs production and purification
Request a detailed protocolcDNA of 2MT-GCaMP6m (mitoGCaMP6m) was subcloned by restriction/ligation (restriction enzymes and T4-ligase were from Fermentas/Thermo Scientific Life Science Research) into a plasmid containing adeno-associated virus 2 (AAV2) inverted terminal repeats flanking a cassette consisting of the neuronal-specific human synapsin one promoter (hSyn), the woodchuck post-transcriptional regulatory element (WPRE) and the bovine growth hormone polyA signal.
Viral particles were produced in HEK293T cells (ATCC) as previously described (Tevet and Gitler, 2016), using pAdDelta5 helper plasmids (a kind gift from Dr. Adi Mizrahi) and the pAAV2/9 n plasmid (Addgene #112865) which encode the rep/cap proteins of AAV2 and AAV9, respectively. Viral particles were then purified over iodixanol (Sigma-Aldrich) step gradients and concentrated using Amicon filters (EMD). Virus titers were measured by determining the number of DNase I–resistant vector genomes (vg) using qPCR with a linearized genome plasmid as a standard (Challis et al., 2019).
Stereotaxic injections
Request a detailed protocolMice at the age of P21-25 were deeply anesthetized with Ketamine/Xylazine and then stereotactic bilateral injections were performed into the Layer 5 of the somatosensory cortex using a microliter syringe (Hamilton, Israel) at a rate of 0.25 μl/minute, with 500 nl of AAV9-hSyn-Mito-GCaMP6m containing 1 × 1010 vg. After the injection, the needle was left in place for additional 3 min before being slowly removed from the brain. Coordinates for injections were (in mm): 4.1 rostral to lambda, ± 1.8 left/right of midline, –0.5 ventral to the pial surface.
Acute coronal brain slices preparation
Request a detailed protocolCoronal slices were prepared from mice three weeks post-injection (at the age of 6–7 weeks). The 300-µm-thick coronal cortical or horizontal hippocampal slices were prepared using standard techniques, as previously described (Baranauskas et al., 2013; Fleidervish et al., 2010). Mice were anesthetized with isoflurane (5%) and decapitated. The slices were cut on a vibratome (VT1200; Leica) and placed in a holding chamber containing oxygenated artificial cerebrospinal fluid (ACSF) at room temperature; they were transferred to a recording chamber after more than 1 hr of incubation. The composition of the ACSF (in mM): 124 NaCl, 3 KCl, 2 CaCl2, 2 MgSO4, 1.25 NaH2PO4, 26 NaHCO3, and 10 glucose; pH 7.4 when saturated with 95% O2/CO2.
Fluorescence imaging
Request a detailed protocolMost experiments were performed on L5 pyramidal neurons in somatosensory neocortical slices. The cells were viewed with a 40 or 60× Olympus water-immersion lens of Ultima IV two-photon microscope (Bruker) equipped with a Mai Tai Deep See pulsed laser (Spectra-Physics). MitoGCaMP6m was excited at 940–950 nm. L5 pyramidal cells with low resting fluorescence that responded to electrical stimulation delivered via the nearby placed bipolar electrode were selected for whole cell recording (see below). In order to measure the cytosolic Ca2+ transients, the intracellular solution was supplemented by Ca2+ indicator, Fura-2 (100 µmol/l). The indicator was selected to minimize the interference with the mitoGCaMP6m fluorescence measurements. The Fura-2 fluorescence was elicited by two-photon excitation at 780 nm, and it declined as a function of the cytosolic Ca2+ concentration. The neuronal morphology and the labelled mitochondria localization was obtained by scanning a Z-series of 30–40 high-resolution images at interval of 0.5 µm. The dynamic mitoGCaMP6m and Fura-2 imaging were performed from small regions of interest at frame rate of 10–50 Hz.
Electrophysiology
Request a detailed protocolSomatic whole-cell recordings were obtained using patch pipettes pulled from thick-walled borosilicate glass capillaries (1.5 mm outer diameter; Science Products, Germany). All recordings were at 30°C ± 0.5°C maintained with a temperature control unit (Luigs & Neumann, Rattingen). For current-clamp experiments the pipette solution contained (in mM): 130 K–gluconate, 6 KCl, 2 MgCl2, 4 NaCl, and 10 Hepes, with pH adjusted to 7.25 with KOH. Pipettes had resistances of 5–7 MΩ when filled with this solution supplemented with Fura-2 (Molecular Probes). Recordings were made using a Multiclamp 700B amplifier (Molecular Devices) equipped with CV-7B headstage (Molecular Devices). Data were low–pass–filtered at 10 kHz (−3 dB), single-pole Bessel filter and digitized at 20 kHz using Digidata 1,322 A digitizer driven by PClamp 10 software (Molecular Devices). Care was taken to maintain the access resistance below 10 MΩ.
Wide-field fluorescence imaging
Request a detailed protocolThe NAD(P)H auto-fluorescence signals were obtained using a 40× water-immersion lens (Olympus) in a BX51WI microscope (Olympus). The fluorescence was excited by using a high-intensity LED device (385 ± 4 nm, Prizmatix), and the emission was collected by using a modified Olympus U-MNU2 filter set (DC = 400 nm; EM = 420 nm). Images were collected with the Orca Flash 4.0 CMOS camera (Hamamatsu), using a pixel binning of 512 × 512, at a rate of 300ms per frame. A bipolar stimulating electrode (WPI, 0.01 MΩ) was placed ~100 µM below the region of interest, at the L-5/6 border in cortical slices or in CA1 stratum oriens in the hippocampal slices. The 0.1ms long extracellular stimuli were delivered using an optically coupled stimulus isolation unit (A.M.P.I) driven via the pClamp 10 software. Somatic whole-cell recordings (see above) were made from a pyramidal neuron in the middle of the region of interest. The stimulation intensity was carefully controlled so that each stimulus triggered only a single antidromic spike with a latency of <1ms post-stimulus. The baseline fluorescence was kept around 1500 a.u. throughout the experiments by regulating the intensity of LED emitted light.
Data analysis
Request a detailed protocolElectrophysiological data analysis was accomplished using pCLAMP10 software (Molecular Devices) and Origin 6.0 (OriginLab). The figures were created using CorelDraw X7 suite (Corel Corporation).
Statistical analysis
Request a detailed protocolIf not otherwise noted, data are expressed as mean ± SE. Student t-test for paired or unpaired data was used for statistical analysis.
Data availability
A representative subset of the raw electrical recording and imaging data has been deposited to Dryad (https://doi.org/10.5061/dryad.sxksn0348). The dataset contains the Microcal Origin opj files of the electrical and optical recordings and quantitative analysis of the data. We are unable to make all raw electrical recording and imaging data publicly available as due to the large size of our raw dataset (>10TB). Interested researchers should contact the corresponding author to gain access to the raw data.
-
Dryad Digital RepositoryFrequency- and spike-timing-dependent mitochondrial Ca2+ signaling regulates the metabolic rate and synaptic efficacy in cortical neurons.https://doi.org/10.5061/dryad.sxksn0348
References
-
An energy budget for signaling in the grey matter of the brainJournal of Cerebral Blood Flow and Metabolism 21:1133–1145.https://doi.org/10.1097/00004647-200110000-00001
-
Mitochondria, calcium and cell death: A deadly triad in neurodegenerationBiochimica et Biophysica Acta (BBA) - Bioenergetics 1787:335–344.https://doi.org/10.1016/j.bbabio.2009.02.021
-
Effects of local anesthetic QX-314 on the membrane properties of hippocampal pyramidal neuronsThe Journal of Pharmacology and Experimental Therapeutics 220:476–481.
-
Enjoy the Trip: Calcium in Mitochondria Back and ForthAnnual Review of Biochemistry 85:161–192.https://doi.org/10.1146/annurev-biochem-060614-034216
-
Role of mitochondrial Ca2+ in the regulation of cellular energeticsBiochemistry 51:2959–2973.https://doi.org/10.1021/bi2018909
-
Experience-dependent structural synaptic plasticity in the mammalian brainNature Reviews. Neuroscience 10:647–658.https://doi.org/10.1038/nrn2699
-
Physiology and anatomy of synaptic connections between thick tufted pyramidal neurones in the developing rat neocortexThe Journal of Physiology 500 (Pt 2):409–440.https://doi.org/10.1113/jphysiol.1997.sp022031
-
MICU3 is a tissue-specific enhancer of mitochondrial calcium uptakeCell Death and Differentiation 26:179–195.https://doi.org/10.1038/s41418-018-0113-8
-
Direct measurement of coupling between dendritic spines and shaftsScience (New York, N.Y.) 272:716–719.https://doi.org/10.1126/science.272.5262.716
-
Mitochondria: the hub of cellular Ca2+ signalingPhysiology (Bethesda, Md.) 23:84–94.https://doi.org/10.1152/physiol.00046.2007
-
Using FRAP or FRAPA to Visualize the Movement of Fluorescently Labeled Proteins or Cellular Organelles in Live Cultured Neurons Transformed with Adeno-Associated VirusesMethods in Molecular Biology (Clifton, N.J.) 1474:125–151.https://doi.org/10.1007/978-1-4939-6352-2_8
-
Voltage-energized Calcium-sensitive ATP Production by MitochondriaNature Metabolism 1:975–984.https://doi.org/10.1038/s42255-019-0126-8
-
A vicious cycle of β amyloid-dependent neuronal hyperactivationScience (New York, N.Y.) 365:559–565.https://doi.org/10.1126/science.aay0198
Article and author information
Author details
Funding
Israel Science Foundation (1384/19)
- Ilya Fleidervish
Israel Science Foundation (1310/19)
- Daniel Gitler
National Institutes of Health (RF1 AG065628)
- Israel Sekler
- Grace Stutzmann
Israel Science Foundation (1763/21)
- Israel Sekler
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This research was supported by The Israel Science Foundation (grant No. 1384/19 for IF, grant No.1763/21 for IS, and grant No. 1310/19 for DG) and by US National Institutes of Health grant (RF1 AG065628).
Ethics
All experiments were approved by the Animal Care and Use Committee of Ben Gurion University of the Negev (protocol #IL-68-09-2019).
Copyright
© 2022, Stoler 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
-
- 1,737
- views
-
- 365
- downloads
-
- 15
- 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
-
- Computational and Systems Biology
- Neuroscience
Diffusional kurtosis imaging (DKI) is a methodology for measuring the extent of non-Gaussian diffusion in biological tissue, which has shown great promise in clinical diagnosis, treatment planning, and monitoring of many neurological diseases and disorders. However, robust, fast, and accurate estimation of kurtosis from clinically feasible data acquisitions remains a challenge. In this study, we first outline a new accurate approach of estimating mean kurtosis via the sub-diffusion mathematical framework. Crucially, this extension of the conventional DKI overcomes the limitation on the maximum b-value of the latter. Kurtosis and diffusivity can now be simply computed as functions of the sub-diffusion model parameters. Second, we propose a new fast and robust fitting procedure to estimate the sub-diffusion model parameters using two diffusion times without increasing acquisition time as for the conventional DKI. Third, our sub-diffusion-based kurtosis mapping method is evaluated using both simulations and the Connectome 1.0 human brain data. Exquisite tissue contrast is achieved even when the diffusion encoded data is collected in only minutes. In summary, our findings suggest robust, fast, and accurate estimation of mean kurtosis can be realised within a clinically feasible diffusion-weighted magnetic resonance imaging data acquisition time.
-
- Neuroscience
The retina transforms patterns of light into visual feature representations supporting behaviour. These representations are distributed across various types of retinal ganglion cells (RGCs), whose spatial and temporal tuning properties have been studied extensively in many model organisms, including the mouse. However, it has been difficult to link the potentially nonlinear retinal transformations of natural visual inputs to specific ethological purposes. Here, we discover a nonlinear selectivity to chromatic contrast in an RGC type that allows the detection of changes in visual context. We trained a convolutional neural network (CNN) model on large-scale functional recordings of RGC responses to natural mouse movies, and then used this model to search in silico for stimuli that maximally excite distinct types of RGCs. This procedure predicted centre colour opponency in transient suppressed-by-contrast (tSbC) RGCs, a cell type whose function is being debated. We confirmed experimentally that these cells indeed responded very selectively to Green-OFF, UV-ON contrasts. This type of chromatic contrast was characteristic of transitions from ground to sky in the visual scene, as might be elicited by head or eye movements across the horizon. Because tSbC cells performed best among all RGC types at reliably detecting these transitions, we suggest a role for this RGC type in providing contextual information (i.e. sky or ground) necessary for the selection of appropriate behavioural responses to other stimuli, such as looming objects. Our work showcases how a combination of experiments with natural stimuli and computational modelling allows discovering novel types of stimulus selectivity and identifying their potential ethological relevance.