OSCR

Robust activity-dependent mitochondrial calcium dynamics at the AIS is dispensable for action potential generation.

Code ↔ Paper

The paper beside its authors' code: matches between them have not been computed for this paper yet.

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

  1. function options = NoRMCorreSetParms(varargin)
  2. % Struct for setting the NoRMCorre algorithm parameters. Any parameter that is
  3. % not set gets a default value
  4. % Author: Eftychios A. Pnevmatikakis
  5. % Simons Foundation, 2016
  6. Names = [
  7. % dataset info
  8. 'd1 ' % number of rows
  9. 'd2 ' % number of cols
  10. 'd3 ' % number of planes (for 3d imaging, default: 1)
  11. % patches
  12. 'grid_size ' % size of non-overlapping regions (default: [d1,d2,d3])
  13. 'overlap_pre ' % size of overlapping region (default: [32,32,16])
  14. 'min_patch_size ' % minimum size of patch (default: [32,32,16])
  15. 'us_fac ' % upsampling factor for subpixel registration (default: 20)
  16. 'mot_uf ' % degree of patches upsampling (default: [4,4,1])
  17. 'max_dev ' % maximum deviation of patch shift from rigid shift (default: [3,3,1])
  18. 'overlap_post ' % size of overlapping region after upsampling (default: [32,32,16])
  19. 'max_shift ' % maximum rigid shift in each direction (default: [15,15,5])
  20. 'phase_flag ' % flag for using phase correlation (default: false)
  21. 'shifts_method ' % method to apply shifts ('FFT','cubic','linear')
  22. % template updating
  23. 'upd_template ' % flag for online template updating (default: true)
  24. 'init_batch ' % length of initial batch (default: 100)
  25. 'bin_width ' % width of each bin (default: 10)
  26. 'buffer_width ' % number of local means to keep in memory (default: 50)
  27. 'method ' % method for averaging the template (default: {'median';'mean})
  28. 'iter ' % number of data passes (default: 1)
  29. 'boundary ' % method of boundary treatment 'NaN','copy','zero','template' (default: 'copy')
  30. % misc
  31. 'add_value ' % add dc value to data (default: 0)
  32. 'use_parallel ' % for each frame, update patches in parallel (default: false)
  33. 'memmap ' % flag for saving memory mapped motion corrected file (default: false)
  34. 'mem_filename ' % name for memory mapped file (default: 'motion_corrected.mat')
  35. 'mem_batch_size ' % batch size during memory mapping for speed (default: 5000)
  36. % plotting
  37. 'plot_flag ' % flag for plotting results in real time (default: false)
  38. 'make_avi ' % flag for making movie (default: false)
  39. 'name ' % name for movie (default: 'motion_corrected.avi')
  40. 'fr ' % frame rate for movie (default: 30)
  41. % output type
  42. 'output_type ' % 'mat' (load in memory), 'memmap', 'tiff', 'hdf5', 'bin' (default:mat)
  43. 'h5_groupname ' % name for hdf5 dataset (default: 'mov')
  44. 'h5_filename ' % name for hdf5 saved file (default: 'motion_corrected.h5')
  45. 'tiff_filename ' % name for saved tiff stack (default: 'motion_corrected.tif')
  46. % use windowing
  47. 'use_windowing ' % flag for windowing data before fft (default: false)
  48. 'window_length ' % length of window on each side of the signal as a fraction of signal length
  49. % total length = length(signal)(1 + 2*window_length). (default: 0.5)
  50. % bitsize for reading .raw files
  51. 'bitsize ' % (default: 2 (uint16). other choices 1 (uint8), 4 (single), 8 (double))
  52. % offset from bidirectional sampling
  53. 'correct_bidir ' % check for offset due to bidirectional scanning (default: true)
  54. 'nFrames ' % number of frames to average (default: 50)
  55. 'bidir_us ' % upsampling factor for bidirectional sampling (default: 10)
  56. 'col_shift ' % known bi-directional offset provided by the user (default: [])
  57. ];
  58. [m,n] = size(Names);
  59. names = lower(Names);
  60. % Combine all leading options structures o1, o2, ... in l1Set(o1,o2,...).
  61. options = [];
  62. for j = 1:m
  63. eval(['options.' Names(j,:) '= [];']);
  64. end
  65. i = 1;
  66. while i <= nargin
  67. arg = varargin{i};
  68. if ischar(arg), break; end
  69. if ~isempty(arg) % [] is a valid options argument
  70. if ~isa(arg,'struct')
  71. error(sprintf(['Expected argument %d to be a string parameter name ' ...
  72. 'or an options structure\ncreated with OPTIMSET.'], i));
  73. end
  74. for j = 1:m
  75. if any(strcmp(fieldnames(arg),deblank(Names(j,:))))
  76. eval(['val = arg.' Names(j,:) ';']);
  77. else
  78. val = [];
  79. end
  80. if ~isempty(val)
  81. eval(['options.' Names(j,:) '= val;']);
  82. end
  83. end
  84. end
  85. i = i + 1;
  86. end
  87. % A finite state machine to parse name-value pairs.
  88. if rem(nargin-i+1,2) ~= 0
  89. error('Arguments must occur in name-value pairs.');
  90. end
  91. expectval = 0; % start expecting a name, not a value
  92. while i <= nargin
  93. arg = varargin{i};
  94. if ~expectval
  95. if ~ischar(arg)
  96. error(sprintf('Expected argument %d to be a string parameter name.', i));
  97. end
  98. lowArg = lower(arg);
  99. j = strmatch(lowArg,names);
  100. if isempty(j) % if no matches
  101. error(sprintf('Unrecognized parameter name ''%s''.', arg));
  102. elseif length(j) > 1 % if more than one match
  103. % Check for any exact matches (in case any names are subsets of others)
  104. k = strmatch(lowArg,names,'exact');
  105. if length(k) == 1
  106. j = k;
  107. else
  108. msg = sprintf('Ambiguous parameter name ''%s'' ', arg);
  109. msg = [msg '(' deblank(Names(j(1),:))];
  110. for k = j(2:length(j))'
  111. msg = [msg ', ' deblank(Names(k,:))];
  112. end
  113. msg = sprintf('%s).', msg);
  114. error(msg);
  115. end
  116. end
  117. expectval = 1; % we expect a value next
  118. else
  119. eval(['options.' Names(j,:) '= arg;']);
  120. expectval = 0;
  121. end
  122. i = i + 1;
  123. end
  124. if expectval
  125. error(sprintf('Expected value for parameter ''%s''.', arg));
  126. end
  127. Values = [
  128. % dataset info
  129. {[]}
  130. {[]}
  131. {1}
  132. % patches
  133. {[]} % size of non-overlapping regions (default: [d1,d2,d3])
  134. {[32,32,16]} % size of overlapping region (default: [32,32,16])
  135. {[32,32,16]} % minimum size of patch (default: [32,32,16])
  136. {50} % upsampling factor for subpixel registration (default: 50)
  137. {[4,4,1]} % degree of patches upsampling (default: [4,4,1])
  138. {[3,3,1]} % maximum deviation of patch shift from rigid shift (default: [3,3,1])
  139. {[32,32,16]} % size of overlapping region after upsampling (default: [32,32,16])
  140. {[15,15,5]} % maximum rigid shift in each direction
  141. {false} % use phase correlation (good for high SNR)
  142. {'FFT'} % method for applying shifts ('FFT', 'linear', 'cubic')
  143. % template updating
  144. {true} % flag for online template updating (default: true)
  145. {100} % length of initial batch (default: 100)
  146. {50} % width of each bin (default: 10)
  147. {50} % number of local means to keep in memory (default: 50)
  148. {{'median';'mean'}} % method for averaging the template (default: {'median';'mean'}
  149. {1} % number of data passes (default: 1)
  150. {'copy'} % method of boundary treatment (default: 'copy')
  151. % misc
  152. {0} % add dc value to data (default: 0)
  153. {false} % for each frame, update patches in parallel (default: false)
  154. {false} % flag for saving memory mapped motion corrected file (default: false)
  155. {'motion_corrected.mat'} % name for memory mapped file (default: 'motion_corrected.mat')
  156. {1000} % batch size used during memory mapping for faster mapping
  157. % plotting
  158. {false} % flag for plotting results in real time (default: false)
  159. {false} % flag for making movie (default: false)
  160. {'motion_corrected.avi'} % name for movie (default: 'motion_corrected.avi')
  161. {30} % frame rate for movie (default: 30)
  162. % output_type
  163. {'mat'}
  164. {'mov'}
  165. {'motion_corrected.h5'}
  166. {'motion_corrected.tif'}
  167. % use_windowing
  168. {false}
  169. {0.5}
  170. % bitsize for reading .raw files
  171. {2}
  172. % offset from bidirectional sampling
  173. {true}
  174. {50}
  175. {10}
  176. {[]}
  177. ];
  178. for j = 1:m
  179. if eval(['isempty(options.' Names(j,:) ')'])
  180. eval(['options.' Names(j,:) '= Values{j};']);
  181. end
  182. end
  183. if isempty(options.d1); options.d1 = input('What is the total number of rows? \n'); end
  184. if isempty(options.d2); options.d2 = input('What is the total number of columns? \n'); end
  185. %if options.d3 == 1; nd = 2; else nd = 3; end
  186. if isempty(options.grid_size); options.grid_size = [options.d1,options.d2,options.d3]; end
  187. if length(options.grid_size) == 1; options.grid_size = options.grid_size*ones(1,3); end
  188. if length(options.grid_size) == 2; options.grid_size(3) = 1; end
  189. if length(options.min_patch_size) == 1; options.min_patch_size = options.min_patch_size*ones(1,3); end
  190. if length(options.min_patch_size) == 2; options.min_patch_size(3) = 1; end
  191. if length(options.overlap_pre) == 1; options.overlap_pre = options.overlap_pre*ones(1,3); end
  192. if length(options.overlap_pre) == 2; options.overlap_pre(3) = 1; end
  193. if length(options.overlap_post) == 1; options.overlap_post = options.overlap_post*ones(1,3); end
  194. if length(options.overlap_post) == 2; options.overlap_post(3) = 1; end
  195. if length(options.max_shift) == 1; options.max_shift = options.max_shift*ones(1,3); end
  196. if length(options.max_shift) == 2; options.max_shift(3) = 1; end
  197. if length(options.max_dev) == 1; options.max_dev = options.max_dev*ones(1,3); end
  198. if length(options.max_dev) == 2; options.max_dev(3) = 1; end
  199. if length(options.mot_uf) == 1; options.mot_uf = options.mot_uf*ones(1,3); end
  200. if length(options.mot_uf) == 2; options.mot_uf(3) = 1; end
  201. options.mot_uf(options.grid_size >= [options.d1,options.d2,options.d3]) = 1;

NoRMCorreSetParms.m at commit c329c3b, no license · at the source

Overview

Authors: Koen Kole1,2, Maarten HP Kole1,3
  1. Department of Axonal Signaling, Netherlands Institute for Neuroscience, Royal Netherlands Academy of Arts and Science, Amsterdam, The Netherlands
  2. Donders Center for Neuroscience, Donders Institute for Brain, Cognition and Behaviour, Nijmegen, The Netherlands
  3. Cell Biology, Neurobiology and Biophysics, Department of Biology, Faculty of Science, University of Utrecht, Utrecht, The Netherlands
Journal: The Journal of physiology, volume 604, issue 7, pages 3032-3052
Dates: received 24 May 2025; accepted 11 February 2026; published online 10 March 2026; in print 1 April 2026
Type: Research article · Language: English
License: CC BY
Identifiers: DOI 10.1113/jp289290 · PMID 41806299 · PMCID PMC13039255 · OpenAlex W7134897345
Open access: hybrid, a free copy (OpenAlex)
Status: code verified
Categories: optical imaging (calcium, voltage, 2-photon) (modality), mouse (organism), cellular / molecular (subfield)
Methods: Statistics, Spectral & time-frequency, Evoked potentials, Connectivity, Single-unit activity, calcium imaging
Keywords: action potential, axon initial segment, calcium buffer, mitochondria, pyramidal neuron
MeSH: Action Potentials*, Axon Initial Segment*, Calcium*, Mitochondria*, Pyramidal Cells*, Animals, Calcium Channels, Calcium Signaling, Mice (* major topic)
Topic: Mitochondrial Function and Pathology (Molecular Biology, Biochemistry, Genetics and Molecular Biology), according to OpenAlex
Funding: ZonMw (0 451 001 201 0066); Exacte en Natuurwetenschappen (865.17.003); Netherlands Organisation for Health Research and Development (0 451 001 201 0066)
Citations: not cited yet (Europe PMC); 67 references in the paper
Research resources: we used wild‐type mice RRID:IMSR_RJ:C57BL, RRID:MMRRC_03, RRID:SCR_0, RRID:SCR_01

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

License: none: the authors keep all their rights
State: the link answers, verified on 30 September 2026
Evidence: files inventoried
Commit: c329c3b7e46844ef53edc10c8fc2af124ad32aff, 14 April 2022
Languages: MATLAB (18)
Size: 54 files, 18 scripts
Software Heritage: not archived
Found in: “Data availability statement”
Holds: README
Not found: license file, CITATION.cff, environment file, tests, continuous integration, documentation
Availability: 1 check, the latest on 30 September 2026: the link answers
  • 30 September 2026: the link answers
18 files

microns-explorer.org

License: none: the authors keep all their rights
State: the link answers, verified on 30 September 2026
Evidence: the link answers
Software Heritage: not checked
Found in: “Data availability statement”
Not found: README, license file, CITATION.cff, environment file, tests, continuous integration, documentation
Availability: 1 check, the latest on 30 September 2026: the link answers (HTTP 200)
  • 30 September 2026: the link answers (HTTP 200)
At the source: microns-explorer.org

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 (http://www.microns-explorer.org). The script used to analyse Ca2+ imaging data is available via https://github.com/Kolelab/mtGCaMP.

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://doi.org/10.1113/jp289290

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/jp289290},
url = {https://doi.org/10.1113/jp289290},
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/03/10
VL - 604
IS - 7
SP - 3032
EP - 3052
SN - 0022-3751
PB - Wiley
DO - 10.1113/jp289290
UR - https://doi.org/10.1113/jp289290
LA - en
ER -

CSL-JSON

{
"id": "10.1113/jp289290",
"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": "J Physiol",
"volume": "604",
"issue": "7",
"page": "3032-3052",
"DOI": "10.1113/jp289290",
"PMID": "41806299",
"PMCID": "PMC13039255",
"ISSN": "0022-3751",
"publisher": "Wiley",
"URL": "https://doi.org/10.1113/jp289290",
"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: Biomedicines
In 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: iScience
In 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 communications
In 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 reports
In 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 communications
In 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: Neuron
In 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 advances
In 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 neuroscience
In 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 biology
In 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 communications
In 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.

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.