Compare runs Export
Dataset Figures Sensitivity Report
Recordings / M13Rec001 / Signature

M13Rec001, Mouse heart, spiral wave and optogenetic termination

64 × 55 px, 9,267 frames, 1 kHz, 133 µm/px, deposit checks: negative values, nan mask all frames, temporal smoothing, spatial smoothing

Run 0b805159, 2026-10-27 22:13 UTC, protocol parameters v2, code 752b95b (uncommitted changes), None workers, signature 080dd569, run view 0b805159, 13 signatures of this recording

Not reportable: range check passed, code uncommitted changes, parameters unsigned, reviewer sign-off pending; 19 rules open

Parameters

% Mouse whole heart optical mapping (Luther lab), v2, read-only
% set by developer (not a domain scientist); provisional, REVIEW_NEEDED 16
baseline_frames  = 641;     % no baseline rule recorded for this run
pad_frames       = 400;     % each side, outside the window
spatial_sigma    = 1.0;     % px, Gaussian, phase only
spatial_truncate = 4.0;     % kernel radius, sigma
lowpass          = [];      % none
polarity         = -1;      % measured, skew of dV/dt
rotor_threshold  = 0.9;     % turns, |winding| per plaquette
rotor_erode_px   = 3;       % px from tissue edge ignored
link_px          = 4.0;     % tracking, max step per frame
max_gap          = 5;       % frames a track may vanish
min_cycles       = 1.0;     % own rotations, min_frames = 64
pair_dist_px     = 3.0;     % opposite-charge pairing
paired_frac      = 0.5;     % max paired fraction of a track
render_step      = 4;       % movie only

Segments

SegmentDefinedAnalysedCycle, framesExpected RotorsBounded By
Illumination with arrhythmia0–3,170400–3,170[not measured]presentLED on before frame 0 (DAQ) → termination step (image, ±17 frames)
Illumination after termination3,170–4,9773,170–4,977[not measured]absenttermination step (image, ±17 frames) → illumination off (DAQ, ±11 frames)
Light off4,977–9,2674,977–8,867[not measured]absentillumination off (DAQ, ±11 frames) → recording end

Frames at 9,267 total, 1.0 ms per frame. Cycle from the MAP channel of the DAQ log (median interval) when one exists, else from per-pixel dominant frequency. Boundaries, expectations and the segment design are provisional (project lead (not a domain scientist), following the paper's design, 10.1371/journal.pcbi.1011660).

Metrics

MetricUnitIllumination with arrhythmiaIllumination after terminationLight offPhysiological Range
Rhythm
Dominant frequencyHz 15.5232.2142.314–4.884 18.75–20 Hz, one recording, 10.1371/journal.pcbi.1011660, ±25% provisional
Rotors
Persistent rotorscount 4–1700 none cited
Max concurrent persistent rotorscount 1–200 none cited
Frames with a singularityfraction 0.4890.001–0.0070.003–0.031 none cited
Beats, where a global rhythm exists
Beatscount 5449 none cited
Global amplitude ratioratio 0.700.960.96 gate ≥ 0.30 provisional
APD50ms 141617 none cited
APD90ms 267678 none cited
Conduction velocitycm/s 29.758.459.4 none cited

Values at the run setting. Red: range across the sensitivity grid (12 baseline × sigma settings, 3 trackers) where the spread exceeds 25 % of the median; not reportable as a value. Beat metrics are not varied over the grid. Dominant frequency over the whole analysed window: 2.362 Hz.

CriterionRuleResult
Singularities present where the design expects them, absent elsewhereeach 'present' segment: frac >= 0.1 and >= 2.0 x max frac over 'absent' segments, at every settingholds at 12 filter settings provisional
Persistent rotors present where the design expects them, absent elsewhereeach 'present' segment: count >= 1; each 'absent' segment: count <= 0, at every settingholds at 12 filter settings and 3 trackers provisional
Persistence threshold, paper contextlifetime ≥ 64 frames (measured cycle) and ≥ 50 frames (cited 20 Hz)agree provisional

Report

The PDF carries the parameters as run, segments, metrics, criteria, reportable values, every provisional rule, the review items, the opening frame as camera, voltage and phase with the parameter footer, and the provenance block. The JSON carries the same plus the full stored metrics, sensitivity and critic records.

Methods

Using the provisional parameter set version 2 for protocol f65d786f-f194-4be2-88b8-10d3dd25d674, the fixed parameters were baseline_frames=641, spatial_sigma=1.0, spatial_truncate=4.0, pad_frames=400, rotor_threshold=0.9, rotor_erode_px=3, link_px=4.0, max_gap=5, min_cycles=1.0, pair_dist_px=3.0, paired_frac=0.5, and render_step=4. The frame rate was 1000 Hz from the publication source recorded for this movie, and the pixel size was 133 um from the publication source. The analysed window was frames 400 to 8867, with provisional segment definitions and provisional criteria: illuminated_arrhythmia from frames 0 to 3170, illuminated_after_termination from frames 3170 to 4977, and light_off from frames 4977 to 9267; within the analysed window these correspond to analysed frames 400 to 3170, 3170 to 4977, and 4977 to 8867.

Reportable values

MetricSegmentValueUnitStatus
persistent rotor tracksall4–17countunstable, not reportable
dominant freq (Hz)light off2.31–4.88Hzunstable, not reportable
persistent rotorslight off0countno cited range
singularity frames fractionlight off0.00283–0.0311fraction of framesunstable, not reportable
persistent rotor count maxlight off0countno cited range
dominant freq (Hz)illuminated arrhythmia15.52Hzwithin range
persistent rotorsilluminated arrhythmia4–17countunstable, not reportable
dominant freq (Hz)illuminated after termination2.214Hzno cited range
singularity frames fractionilluminated arrhythmia0.4888fraction of framesno cited range
persistent rotorsilluminated after termination0countno cited range
persistent rotor count maxilluminated arrhythmia1–2countunstable, not reportable
singularity frames fractionilluminated after termination0.000553–0.00719fraction of framesunstable, not reportable
persistent rotor count maxilluminated after termination0countno cited range

physiology_ok: True; send back: False. Verdict: report only stable/descriptive findings and the stable in-range illuminated-arrhythmia dominant frequency; physiology_ok is true because the only ranged check, illuminated_arrhythmia.dominant_freq_hz, is stable at 15.523 Hz and within the checked 14.0625–25.0 Hz contextual range derived from the cited mouse spiral-wave passages reporting approximately 20.0 Hz and 18.75 Hz (DOI 10.1371/journal.pcbi.1011660). Do not report single values for unstable metrics: light_off.singularity_frames_frac 0.002827763496143959–0.031105398457583547, light_off.dominant_freq_hz 2.314–4.884 Hz, illuminated_arrhythmia.persistent_rotors 4.0–17.0, illuminated_arrhythmia.persistent_rotor_count_max 1.0–2.0, illuminated_after_termination.singularity_frames_frac 0.0005534034311012728–0.007194244604316547, and persistent_rotor_tracks 4.0–17.0. Persistent-rotor presence holds provisionally with the 64-frame measured-cycle criterion and the 50-frame paper-rotation criterion, and strict/default/loose trackers agree on the rule verdict despite varying counts. Review the pre-deposit processing/raw-data label, the frame-3170 event interpretation, and baseline/illumination artifact before using this for critical analysis; the ranges, the stability threshold, and the criteria are provisional.

Provisional rules

Rule, LocationDescriptionExpertStatus
protocol_parameter_set
worker/voltage/protocol_params.py (developer_defaults, set_params)
Each protocol's fixed, versioned analysis parameters (version 1 = developer defaults); no run or agent may deviateoptical mapping analystprovisional
filter_defaults
worker/voltage/preprocess.py (Params)
Default baseline window (frames) and spatial Gaussian sigma/truncate for phaseoptical mapping analyst (cardiac electrophysiology)provisional
tissue_mask
worker/voltage/io.py (tissue_mask)
Tissue = not NaN in frame 0 and nonzero temporal stdoptical mapping analystprovisional
polarity
worker/voltage/preprocess.py (polarity)
Dye polarity from sign of median per-pixel skew of dV/dtoptical mapping analyst / voltage dye chemistprovisional
singularity_detection
worker/voltage/preprocess.py (Params.rotor_threshold, rotor_erode_px)
Phase singularity = |winding| > threshold turns, >= erode px from tissue edgecomputational cardiac electrophysiologist (phase mapping)provisional
rotor_persistence
worker/voltage/rotors.py (TrackParams)
Track linking distance/gap and minimum lifetime for a persistent rotorcomputational cardiac electrophysiologist (phase mapping)provisional
chirality_pairing
worker/voltage/rotors.py (TrackParams.pair_dist_px, paired_frac)
Opposite-charge singularities within pair distance are a pair; tracks paired most of their life are excluded from rotorscomputational cardiac electrophysiologist (phase mapping)provisional
sensitivity_grid
worker/voltage/sensitivity.py (SIGMAS, BASELINE_MULTIPLES, baseline_levels)
Filter levels tested: spatial sigma and baseline window valuesoptical mapping analystprovisional
analysis_window
worker/voltage/runs.py (default_window)
Frames analysed: whole recording minus one filter pad at each end (unless the plan narrows it)cardiac electrophysiologist familiar with the datasetprovisional
prior_processing_flags
worker/voltage/inspect.py (prior_processing)
Thresholds for flagging depositor processing (normalization, masking, temporal/spatial smoothing)imaging scientist (camera acquisition and noise)provisional
baseline_step_detection
worker/voltage/events.py (detect_baseline_step)
Baseline step = consensus upward step in tissue-mean raw baseline (cause not inferred; first misread as light onset, contradicted by the DAQ)optical mapping experimentalistprovisional
daq_interpretation
worker/voltage/daq.py
DAQ channel names from depositor notebook; frames = first N triggers of closest burst; midpoint thresholdsexperimentalist who recorded the data (or the depositor)provisional
segment_definitions
worker/voltage/annotations.py (MOUSE_SEGMENTS)
Mouse segments bounded by DAQ light events and the baseline step, interpreted as spiral termination under illumination (project lead, from Hussaini et al. 2023)cardiac electrophysiologist (ideally the depositing authors)provisional
rotor_finding_criterion
worker/voltage/sensitivity.py (MIN_FRAC, MIN_RATIO)
Singularity before/after finding holds if frac(before) >= min_frac and >= min_ratio x frac(after) at every setting (also kept, superseded, for persistent rotors)cardiac electrophysiologist + biostatisticianprovisional
persistent_rotor_criterion
worker/voltage/sensitivity.py (MIN_ROTORS_PRESENT, MAX_ROTORS_ABSENT)
Persistent-rotor finding holds if every 'present' segment has >= 1 persistent rotor and every 'absent' segment has none, at every settingcardiac electrophysiologistprovisional
pipeline_check
worker/agents/pipeline.py (INSTRUCTIONS, validator)
The pipeline agent's hold/proceed decision from protocol mismatches and inspector anomalies; it chooses no parameteroptical mapping analystprovisional
physiology_reference_rates
worker/voltage/physiology.py (RANGES, TOLERANCE_FRAC)
Cited physiological ranges (single-study values widened by a 25% tolerance), applied to segments by their expected rhythm; used by the critic and for frame-rate sanity checkscardiac electrophysiologist per speciesprovisional
stability_rule
worker/agents/critic.py (STABLE_FRAC, stability)
A metric is reportable only if its spread across all sensitivity settings is <= 25% of its median (zero spread when all values are 0); otherwise 'unstable, not reportable'biostatistician + optical mapping analystprovisional
persistence_threshold_variants
worker/agents/critic.py (persistence_variants)
The persistent-rotor rule is run at two lifetimes, one measured cycle of the fastest segment and the paper's cited rotation period; both verdicts stored, disagreement flaggedcomputational cardiac electrophysiologist (phase mapping)provisional

Review items

TopicDetailSeverityReviewer
unstable:light_off.singularity_frames_fracFilter-grid sensitivity gives 0.002827763496143959–0.031105398457583547; do not report a single light-off singularity-frames fraction.warnanalysis
unstable:light_off.dominant_freq_hzBaseline-window sensitivity gives 2.314–4.884 Hz; do not report a single light-off dominant frequency.warnanalysis
unstable:illuminated_arrhythmia.persistent_rotorsFilter and tracker sensitivity gives 4.0–17.0 persistent rotors; report presence only, not a single count.warnanalysis
unstable:illuminated_arrhythmia.persistent_rotor_count_maxFilter and tracker sensitivity gives 1.0–2.0 maximum simultaneous persistent rotors; do not report a single maximum count.warnanalysis
unstable:illuminated_after_termination.singularity_frames_fracFilter-grid sensitivity gives 0.0005534034311012728–0.007194244604316547; do not report a single post-termination singularity-frames fraction.warnanalysis
unstable:persistent_rotor_tracksAggregate persistent-rotor tracks vary from 4.0 to 17.0 across filters and trackers; do not report a single track count.warnanalysis
Pre-deposit processing / raw-data labelThe recording is labeled as raw experimental data in the README, but documentation and measurements indicate masking, spatial smoothing, temporal smoothing, and normalization; analyses that require raw photon counts or independent pixel noise should not treat this as raw camera data.warndata_inspector
Event interpretation at frame 3170Frame 3170 is a baseline_step only; illumination began before frame 0 and ended at frame 4977, so frame 3170 should not be treated as light onset. The termination label is provisional.warndata_inspector
Baseline and illumination artifactThe paper reports stimulation-light crosstalk and re-normalization; this recording also has a large slow baseline range and a baseline step, so absolute baseline or amplitude analyses require artifact review.warndata_inspector
Provisional tool thresholdsEvent detection, processing flags, SNR, rhythm, and reference-rate thresholds have not yet been reviewed by a domain scientist.infodata_inspector

Data and code availability

Recording10.5281/zenodo.10228296, CC-BY-4.0, Hussaini, Sayedeh; Luther, Stefan; Lädke, Sarah L.
Publication10.1371/journal.pcbi.1011660
Analysis code752b95b01dde60f86eb71b47343448c7af4551c1 (uncommitted changes)
Stored arrays[not stored]
Rundone
Review items30 open