Robust activity-dependent mitochondrial calcium dynamics at the AIS is dispensable for action potential generation.
Paper
Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
MATLAB · 215 lines · 9.9 KB · no license
- function options = NoRMCorreSetParms(varargin)
- % Struct for setting the NoRMCorre algorithm parameters. Any parameter that is
- % not set gets a default value
- % Author: Eftychios A. Pnevmatikakis
- % Simons Foundation, 2016
- Names = [
- % dataset info
- 'd1 ' % number of rows
- 'd2 ' % number of cols
- 'd3 ' % number of planes (for 3d imaging, default: 1)
- % patches
- 'grid_size ' % size of non-overlapping regions (default: [d1,d2,d3])
- 'overlap_pre ' % size of overlapping region (default: [32,32,16])
- 'min_patch_size ' % minimum size of patch (default: [32,32,16])
- 'us_fac ' % upsampling factor for subpixel registration (default: 20)
- 'mot_uf ' % degree of patches upsampling (default: [4,4,1])
- 'max_dev ' % maximum deviation of patch shift from rigid shift (default: [3,3,1])
- 'overlap_post ' % size of overlapping region after upsampling (default: [32,32,16])
- 'max_shift ' % maximum rigid shift in each direction (default: [15,15,5])
- 'phase_flag ' % flag for using phase correlation (default: false)
- 'shifts_method ' % method to apply shifts ('FFT','cubic','linear')
- % template updating
- 'upd_template ' % flag for online template updating (default: true)
- 'init_batch ' % length of initial batch (default: 100)
- 'bin_width ' % width of each bin (default: 10)
- 'buffer_width ' % number of local means to keep in memory (default: 50)
- 'method ' % method for averaging the template (default: {'median';'mean})
- 'iter ' % number of data passes (default: 1)
- 'boundary ' % method of boundary treatment 'NaN','copy','zero','template' (default: 'copy')
- % misc
- 'add_value ' % add dc value to data (default: 0)
- 'use_parallel ' % for each frame, update patches in parallel (default: false)
- 'memmap ' % flag for saving memory mapped motion corrected file (default: false)
- 'mem_filename ' % name for memory mapped file (default: 'motion_corrected.mat')
- 'mem_batch_size ' % batch size during memory mapping for speed (default: 5000)
- % plotting
- 'plot_flag ' % flag for plotting results in real time (default: false)
- 'make_avi ' % flag for making movie (default: false)
- 'name ' % name for movie (default: 'motion_corrected.avi')
- 'fr ' % frame rate for movie (default: 30)
- % output type
- 'output_type ' % 'mat' (load in memory), 'memmap', 'tiff', 'hdf5', 'bin' (default:mat)
- 'h5_groupname ' % name for hdf5 dataset (default: 'mov')
- 'h5_filename ' % name for hdf5 saved file (default: 'motion_corrected.h5')
- 'tiff_filename ' % name for saved tiff stack (default: 'motion_corrected.tif')
- % use windowing
- 'use_windowing ' % flag for windowing data before fft (default: false)
- 'window_length ' % length of window on each side of the signal as a fraction of signal length
- % total length = length(signal)(1 + 2*window_length). (default: 0.5)
- % bitsize for reading .raw files
- 'bitsize ' % (default: 2 (uint16). other choices 1 (uint8), 4 (single), 8 (double))
- % offset from bidirectional sampling
- 'correct_bidir ' % check for offset due to bidirectional scanning (default: true)
- 'nFrames ' % number of frames to average (default: 50)
- 'bidir_us ' % upsampling factor for bidirectional sampling (default: 10)
- 'col_shift ' % known bi-directional offset provided by the user (default: [])
- ];
- [m,n] = size(Names);
- names = lower(Names);
- % Combine all leading options structures o1, o2, ... in l1Set(o1,o2,...).
- options = [];
- for j = 1:m
- eval(['options.' Names(j,:) '= [];']);
- end
- i = 1;
- while i <= nargin
- arg = varargin{i};
- if ischar(arg), break; end
- if ~isempty(arg) % [] is a valid options argument
- if ~isa(arg,'struct')
- error(sprintf(['Expected argument %d to be a string parameter name ' ...
- 'or an options structure\ncreated with OPTIMSET.'], i));
- end
- for j = 1:m
- if any(strcmp(fieldnames(arg),deblank(Names(j,:))))
- eval(['val = arg.' Names(j,:) ';']);
- else
- val = [];
- end
- if ~isempty(val)
- eval(['options.' Names(j,:) '= val;']);
- end
- end
- end
- i = i + 1;
- end
- % A finite state machine to parse name-value pairs.
- if rem(nargin-i+1,2) ~= 0
- error('Arguments must occur in name-value pairs.');
- end
- expectval = 0; % start expecting a name, not a value
- while i <= nargin
- arg = varargin{i};
- if ~expectval
- if ~ischar(arg)
- error(sprintf('Expected argument %d to be a string parameter name.', i));
- end
- lowArg = lower(arg);
- j = strmatch(lowArg,names);
- if isempty(j) % if no matches
- error(sprintf('Unrecognized parameter name ''%s''.', arg));
- elseif length(j) > 1 % if more than one match
- % Check for any exact matches (in case any names are subsets of others)
- k = strmatch(lowArg,names,'exact');
- if length(k) == 1
- j = k;
- else
- msg = sprintf('Ambiguous parameter name ''%s'' ', arg);
- msg = [msg '(' deblank(Names(j(1),:))];
- for k = j(2:length(j))'
- msg = [msg ', ' deblank(Names(k,:))];
- end
- msg = sprintf('%s).', msg);
- error(msg);
- end
- end
- expectval = 1; % we expect a value next
- else
- eval(['options.' Names(j,:) '= arg;']);
- expectval = 0;
- end
- i = i + 1;
- end
- if expectval
- error(sprintf('Expected value for parameter ''%s''.', arg));
- end
- Values = [
- % dataset info
- {[]}
- {[]}
- {1}
- % patches
- {[]} % size of non-overlapping regions (default: [d1,d2,d3])
- {[32,32,16]} % size of overlapping region (default: [32,32,16])
- {[32,32,16]} % minimum size of patch (default: [32,32,16])
- {50} % upsampling factor for subpixel registration (default: 50)
- {[4,4,1]} % degree of patches upsampling (default: [4,4,1])
- {[3,3,1]} % maximum deviation of patch shift from rigid shift (default: [3,3,1])
- {[32,32,16]} % size of overlapping region after upsampling (default: [32,32,16])
- {[15,15,5]} % maximum rigid shift in each direction
- {false} % use phase correlation (good for high SNR)
- {'FFT'} % method for applying shifts ('FFT', 'linear', 'cubic')
- % template updating
- {true} % flag for online template updating (default: true)
- {100} % length of initial batch (default: 100)
- {50} % width of each bin (default: 10)
- {50} % number of local means to keep in memory (default: 50)
- {{'median';'mean'}} % method for averaging the template (default: {'median';'mean'}
- {1} % number of data passes (default: 1)
- {'copy'} % method of boundary treatment (default: 'copy')
- % misc
- {0} % add dc value to data (default: 0)
- {false} % for each frame, update patches in parallel (default: false)
- {false} % flag for saving memory mapped motion corrected file (default: false)
- {'motion_corrected.mat'} % name for memory mapped file (default: 'motion_corrected.mat')
- {1000} % batch size used during memory mapping for faster mapping
- % plotting
- {false} % flag for plotting results in real time (default: false)
- {false} % flag for making movie (default: false)
- {'motion_corrected.avi'} % name for movie (default: 'motion_corrected.avi')
- {30} % frame rate for movie (default: 30)
- % output_type
- {'mat'}
- {'mov'}
- {'motion_corrected.h5'}
- {'motion_corrected.tif'}
- % use_windowing
- {false}
- {0.5}
- % bitsize for reading .raw files
- {2}
- % offset from bidirectional sampling
- {true}
- {50}
- {10}
- {[]}
- ];
- for j = 1:m
- if eval(['isempty(options.' Names(j,:) ')'])
- eval(['options.' Names(j,:) '= Values{j};']);
- end
- end
- if isempty(options.d1); options.d1 = input('What is the total number of rows? \n'); end
- if isempty(options.d2); options.d2 = input('What is the total number of columns? \n'); end
- %if options.d3 == 1; nd = 2; else nd = 3; end
- if isempty(options.grid_size); options.grid_size = [options.d1,options.d2,options.d3]; end
- if length(options.grid_size) == 1; options.grid_size = options.grid_size*ones(1,3); end
- if length(options.grid_size) == 2; options.grid_size(3) = 1; end
- if length(options.min_patch_size) == 1; options.min_patch_size = options.min_patch_size*ones(1,3); end
- if length(options.min_patch_size) == 2; options.min_patch_size(3) = 1; end
- if length(options.overlap_pre) == 1; options.overlap_pre = options.overlap_pre*ones(1,3); end
- if length(options.overlap_pre) == 2; options.overlap_pre(3) = 1; end
- if length(options.overlap_post) == 1; options.overlap_post = options.overlap_post*ones(1,3); end
- if length(options.overlap_post) == 2; options.overlap_post(3) = 1; end
- if length(options.max_shift) == 1; options.max_shift = options.max_shift*ones(1,3); end
- if length(options.max_shift) == 2; options.max_shift(3) = 1; end
- if length(options.max_dev) == 1; options.max_dev = options.max_dev*ones(1,3); end
- if length(options.max_dev) == 2; options.max_dev(3) = 1; end
- if length(options.mot_uf) == 1; options.mot_uf = options.mot_uf*ones(1,3); end
- if length(options.mot_uf) == 2; options.mot_uf(3) = 1; end
- options.mot_uf(options.grid_size >= [options.d1,options.d2,options.d3]) = 1;
NoRMCorreSetParms.m at commit c329c3b, no license · at the source
Overview
- Department of Axonal Signaling, Netherlands Institute for Neuroscience, Royal Netherlands Academy of Arts and Science, Amsterdam, The Netherlands
- Donders Center for Neuroscience, Donders Institute for Brain, Cognition and Behaviour, Nijmegen, The Netherlands
- Cell Biology, Neurobiology and Biophysics, Department of Biology, Faculty of Science, University of Utrecht, Utrecht, The Netherlands
Abstract
Abstract: Mitochondria are diverse and multifaceted intracellular organelles regulating oxidative energy supply, lipid metabolism and calcium (Ca2+) signalling. In neurons the spatial sequestration of cytoplasmic Ca2+ by mitochondria plays a critical role in determining activity‐dependent spine plasticity, shaping the presynaptic transmitter release characteristics and contributing to sustained action potential firing. Here, we tested the hypothesis that mitochondria at the axon initial segment (AIS) affect the microdomain cytoplasmic Ca2+ transients, thereby regulating Ca2+‐dependent voltage‐gated ion channels at the plasma membrane and initiation of action potentials. Using 3D electron microscopy reconstructions and virally injecting genetically encoded fluorescence indicators we visualized the ultrastructure and distribution of mitochondria selectively in thick‐tufted layer 5 pyramidal neurons. We found that most mitochondria were stably clustered to the proximal AIS, while few were observed at distal sites. Simultaneous two‐photon imaging of action potential‐dependent cytoplasmic and mitochondrial Ca2+, combined with electrophysiological recordings showed that AIS mitochondria exhibit powerful activity‐dependent cytosolic Ca2+ uptake. However, while intracellular application of the mitochondrial Ca2+ uniporter inhibitor Ru360 fully blocked mitochondrial Ca2+ import and increased the slow afterhyperpolarization duration, it did not affect action potential input–output function, action potential dynamics nor the ability to produce high‐frequency burst output. Together, the results indicate that AIS mitochondria are dispensable for temporal and rate encoding, suggesting that mt‐Ca2+ buffering at the AIS may be involved in non‐electrical roles.
Key points: Mitochondrial Ca2+ buffering controls multiple Ca2+‐dependent intracellular processes and their subcellular location of the organelles defines local physiological properties in neurons.
Recent studies implicate mitochondrial Ca2+ uptake in the slow afterhyperpolarization and maintenance of action potential firing.
Using electron microscopy and virally delivered genetically encoded tools we examined mitochondria in the layer 5 pyramidal neuron axon initial segment (AIS), the site where action potentials initiate, and found that cytoplasmic Ca2+ influx is powerfully buffered by proximally clustered mitochondria.
Electrophysiological recordings during the block of the mitochondrial calcium uniporter reveal a role in the slow afterhyperpolarization, while AIS action potential initiation and action potential waveforms are independent from mitochondria.
These findings indicate AIS mitochondria under physiological conditions exert non‐electrical roles.
Reproduced under the paper's license (CC BY), from the paper cited above.
Repositories
Its files are read in the Code ↔ Paper reader above.
Kolelab/mtGCaMP
c329c3b7e46844ef53edc10c8fc2af124ad32aff, 14 April 2022Availability: 1 check, the latest on 30 September 2026: the link answers
- 30 September 2026: the link answers
18 files
- Functions/
NoRMCorreSetParms.m , MATLAB, 215 lines - Functions/
cell2mat_ov.m , MATLAB, 27 lines - Functions/
cell2mat_ov_sum.m , MATLAB, 42 lines - Functions/
construct_grid.m , MATLAB, 33 lines - Functions/
construct_weights.m , MATLAB, 8 lines - Functions/
correct_bidirectional_of , MATLAB, 77 linesfset.m - Functions/
dftregistration_min_max. , MATLAB, 267 linesm - Functions/
dftregistration_min_max_ , MATLAB, 315 lines3d.m - Functions/
mat2cell_ov.m , MATLAB, 29 lines - Functions/
natsort.m , MATLAB, 344 lines - Functions/
natsortfiles.m , MATLAB, 199 lines - Functions/
normcorre.m , MATLAB, 524 lines - Functions/
read_file.m , MATLAB, 89 lines - Functions/
refreshdisp.m , MATLAB, 12 lines - Functions/
remove_boundaries.m , MATLAB, 40 lines - Functions/
saveastiff.m , MATLAB, 304 lines - Functions/
shift_reconstruct.m , MATLAB, 57 lines - mtGCaMP_analysis.m, MATLAB, 1,413 lines
microns-explorer.org
Availability: 1 check, the latest on 30 September 2026: the link answers (HTTP 200)
- 30 September 2026: the link answers (HTTP 200)
The paper's code and data availability statement is in the Data section.
Tracing map
Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.
What the map holds:
- 2 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 18 scripts, each with its path and the digest of its content;
- no match between paragraphs and code yet;
- neither the text of the paper nor the code itself.
Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.
Data
No dataset and no data link were found in the paper.
Data availability statement
All data are available upon request to the corresponding authors. The Microns dataset used for 3D EM analysis is available at www.microns‐explorer.org
Reproduced under the paper's license (CC BY), from the paper cited above.
Versions
The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.
Version 1, 30 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 2 authors, 5 keywords, 9 MeSH terms, 3 funders, 67 references, 4 RRIDs.
Cite
This paper
Kole, K., & Kole, M. H. (2026). Robust activity-dependent mitochondrial calcium dynamics at the AIS is dispensable for action potential generation. The Journal of physiology, 604(7), 3032-3052. https://
BibTeX
@article{kole2026robust,
author = {Kole, Koen and Kole, Maarten HP},
title = {{Robust activity-dependent mitochondrial calcium dynamics at the AIS is dispensable for action potential generation}},
journal = {The Journal of physiology},
year = {2026},
month = mar,
volume = {604},
number = {7},
pages = {3032--3052},
publisher = {Wiley},
issn = {0022-3751},
doi = {10.1113/
url = {https://
pmid = {41806299},
pmcid = {PMC13039255}
}
RIS
TY - JOUR
AU - Kole, Koen
AU - Kole, Maarten HP
TI - Robust activity-dependent mitochondrial calcium dynamics at the AIS is dispensable for action potential generation
T2 - The Journal of physiology
J2 - J Physiol
PY - 2026
DA - 2026/
VL - 604
IS - 7
SP - 3032
EP - 3052
SN - 0022-3751
PB - Wiley
DO - 10.1113/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1113/
"type": "article-journal",
"title": "Robust activity-dependent mitochondrial calcium dynamics at the AIS is dispensable for action potential generation",
"container-title": "The Journal of physiology",
"author": [
{
"family": "Kole",
"given": "Koen"
},
{
"family": "Kole",
"given": "Maarten HP"
}
],
"container-title-short":
"volume": "604",
"issue": "7",
"page": "3032-3052",
"DOI": "10.1113/
"PMID": "41806299",
"PMCID": "PMC13039255",
"ISSN": "0022-3751",
"publisher": "Wiley",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
3,
10
]
]
}
}
The tracing map gets a citation of its own once an author has validated it and it has a DOI.
Similar papers
The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.
- [1] doi:10.3390/biomedicines14081670 [code]
- Limited Detectability of Network Functional Alterations in a Tauopathy Model Using Mouse Primary Cortical Cultures.Journal: BiomedicinesIn common: Curve Fitting Toolbox, Parallel Computing Toolbox, Image Processing Toolbox, 2 other tools, optical imaging (calcium, voltage, 2-photon), mouse
- [2] doi:10.1016/j.isci.2026.117187 [code]
- Functional and structural characterization of dendritic spine pathology in a mouse model of tauopathy.Journal: iScienceIn common: Curve Fitting Toolbox, Parallel Computing Toolbox, Image Processing Toolbox, 2 other tools, mouse
- [3] doi:10.1038/s41467-026-75347-4 [code]
- Sleep reveals dynamics integrating and segregating movement and stimulus representations in V1.Journal: Nature communicationsIn common: Curve Fitting Toolbox, Parallel Computing Toolbox, Image Processing Toolbox, 2 other tools, mouse
- [4] doi:10.1016/j.celrep.2026.117646 [code]
- Medial entorhinal-hippocampal desynchronization parallels the emergence of memory impairment in a mouse model of Alzheimer's disease pathology.Journal: Cell reportsIn common: Curve Fitting Toolbox, Parallel Computing Toolbox, Image Processing Toolbox, 2 other tools, mouse
- [5] doi:10.1038/s41467-026-73106-z [code]
- Respiratory pauses highlight sleep architecture in mice.Journal: Nature communicationsIn common: Curve Fitting Toolbox, Parallel Computing Toolbox, Image Processing Toolbox, 2 other tools, mouse
- [6] doi:10.1016/j.neuron.2026.03.034 [code]
- Dentate gyrus interneurons modulate winner-take-all network dynamics in freely behaving mice.Journal: NeuronIn common: Curve Fitting Toolbox, Parallel Computing Toolbox, Image Processing Toolbox, 2 other tools, mouse
- [7] doi:10.1126/sciadv.aeh7220 [code]
- Central complex representations of self-movement are sufficient to compute wind direction in flight.Journal: Science advancesIn common: Curve Fitting Toolbox, Parallel Computing Toolbox, Image Processing Toolbox, 2 other tools
- [8] doi:10.1038/s41593-026-02357-2 [code]
- Experience reorganizes content-specific memory traces in macaques.Journal: Nature neuroscienceIn common: Curve Fitting Toolbox, Parallel Computing Toolbox, Image Processing Toolbox, 2 other tools
- [9] doi:10.1371/journal.pcbi.1014563 [code]
- Single pulse electrical stimulation in white matter modulates iEEG visual responses in human early visual cortex.Journal: PLoS computational biologyIn common: Curve Fitting Toolbox, Parallel Computing Toolbox, Image Processing Toolbox, 2 other tools
- [10] doi:10.1038/s41467-026-75490-y [code]
- Topographically organized dorsal raphe activity modulates forebrain sensory-motor representations and contributes to defensive behaviors.Journal: Nature communicationsIn common: Curve Fitting Toolbox, Parallel Computing Toolbox, Image Processing Toolbox, 2 other tools
Contribute
The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.
Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 2 repositories of the authors' code, each at its verified commit and with its license, 18 scripts, and 0 matches between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:04e7bffdc04e1cb3…
Add the badge to its README
The badge links the code to this page. Copy one of these into the README of the paper's code: only you decide where it goes, and nothing is changed for you.
Markdown
[, paste the snippet at the top, then “Commit changes…” and, to review it first, “Create a new branch and start a pull request”. You open the pull request; OSCR asks for no permission.
Request its removal
To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).
Discussion, reproductions, activity
Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.
Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
