{"id":"f51275b6-727d-43c4-a207-1e509a41e271","arxiv_id":"2606.22695","paper_version":2,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"SPIDER stitches overlapping power spectra and completes missing pairs to recover frequency-resolved directed connectivity from incomplete asynchronous brain recordings.","lead":"SPIDER reconstructs directed brain connectivity from incomplete, asynchronous multi-session recordings that never share a common clock. It could let labs pool partial datasets into whole-brain effective-connectivity maps that no single experiment can produce.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"Nuclear-norm completion of never-co-observed cross-spectra must recover phase structure accurate enough for spectral factorization to yield correct PDC; this is the least secure step for asynchronous data.","rationale":"The reader’s weakest assumption correctly isolates the stitching-plus-nuclear-norm-completion step as the single load-bearing point on which recovery of directed flow (and therefore the hierarchy claim) rests. No stronger internal contradiction is visible from the abstract alone; the concern is precisely the one flagged. Because proofs, rank conditions, overlap geometry, and simulation details remain unavailable, the abstract-only verdict stays UNVERDICTED. The proposed simulation directly tests whether completion preserves the phase information required by PDC for unobserved pairs—the precise condition that must hold for the strongest claim to be valid.","tokens_in":2186,"tokens_out":555,"duration_ms":17750,"concrete_test":"Generate a ground-truth linear Gaussian process on 20–50 nodes with known directed graph and realistic spectral rank; sample incomplete asynchronous sessions whose overlap geometry matches the IBL Neuropixels coverage (including never-co-observed pairs); apply the full SPIDER pipeline and compare Frobenius / rank-correlation error of recovered PDC on never-co-observed pairs versus co-observed pairs. If never-co-observed error exceeds co-observed error by >20 % or the recovered source of the theta hierarchy differs from ground truth, the completion step fails to preserve directed structure.","verdict_should_be":"UNCHANGED","load_bearing_attack":"SPIDER’s central claim requires that local spectral matrices (auto- and cross-spectra from co-recorded subsets) can be stitched and nuclear-norm-completed into a global Hermitian S(f) whose spectral factorization produces a transfer function whose PDC recovers true directed interactions, including for region pairs never observed together. Nuclear-norm completion recovers missing entries under a low-rank assumption, yet directed measures depend on the complex off-diagonal phase relationships that encode timing. Asynchronous sessions supply no shared clock, so absolute phases between non-overlapping blocks are never measured and must be invented solely from the geometry of the observed blocks. If the true neural spectral matrices are not sufficiently low-rank at the frequencies of interest, or if the overlap pattern fails to uniquely determine the missing blocks to the precision required by factorization, the completed S(f) can distort phases, the resulting PDC can be artifactual, and both the whole-brain maps and the theta-band hippocampal hierarchy become unreliable. The abstract asserts consistency guarantees but does not state the precise rank, stationarity, or overlap conditions under which those guarantees hold for complex spectral matrices.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript introduces SPIDER, a frequency-domain pipeline for estimating directed effective connectivity (via partial directed coherence) from multi-session neural recordings that only partially overlap in space and lack a shared clock. Local power-spectral density matrices from co-recorded channel subsets are stitched into a global spectral matrix; never-co-observed region pairs are filled by nuclear-norm (low-rank) completion; canonical spectral factorization then yields a transfer function from which frequency-resolved PDC is obtained. The abstract asserts consistency guarantees and reports validation on simulations, two-photon imaging, the IBL Neuropixels corpus (50 areas, 43 sessions, 12 labs), and human resting iEEG (43 patients). Beyond method validation, SPIDER is claimed to reveal largely recurrent spontaneous flow overall, but a significant theta-band feedforward hierarchy with the hippocampal formation as source, recovered in both mouse and human data.","tokens_in":2437,"tokens_out":1168,"duration_ms":20492,"significance":"If the method is sound, SPIDER would make whole-brain directed-connectivity analysis tractable for the multi-session, multi-animal, multi-lab datasets that currently cannot support Granger causality, DCM, or PDC because regions are never recorded together under a common clock. The cross-species, cross-modality recovery of a theta-band hippocampal feedforward hierarchy would be a substantive systems-neuroscience result. Strengths claimed in the abstract—consistency guarantees, multi-modal empirical validation, and application at the scale of 50 areas / 12 labs—would, if substantiated in the full methods and results, constitute a genuine methodological advance rather than a re-packaging of existing spectral tools.","major_comments":[{"comment":"The central load-bearing step is nuclear-norm completion of never-co-observed complex cross-spectra, followed by spectral factorization and PDC. Directed measures depend on off-diagonal phase structure that encodes relative timing; asynchronous sessions supply no shared clock, so absolute phases between non-overlapping blocks are never measured. The abstract asserts consistency guarantees but does not state the rank, stationarity, or overlap-graph conditions under which those guarantees hold for Hermitian spectral matrices at the frequencies used for PDC. Without those conditions (and without a proof or theorem statement visible here), it is unclear whether completed phases are uniquely determined to the precision factorization requires, or whether the resulting PDC for unobserved pairs can be artifactual.","section":"Abstract (method core: stitch + nuclear-norm completion + spectral factorization/PDC)"},{"comment":"Validation is described as recovering directed flow among 50 areas from 43 sessions never recorded together, plus a significant theta hierarchy. For the claim about never-co-observed pairs to be load-bearing, the manuscript must show recovery against ground truth or held-out simultaneous co-recordings (e.g., leave-out-overlap tests where a pair is treated as unobserved and then compared to its true co-recorded spectrum/PDC). Running the pipeline on IBL/iEEG and obtaining a hierarchy is not by itself evidence that completion preserved directed structure. The abstract does not indicate that such a held-out test was performed.","section":"Abstract (validation claims: simulations, IBL Neuropixels, human iEEG)"},{"comment":"The scientific headline—a significant theta-band feedforward hierarchy with hippocampal formation as source, recovered across mouse Neuropixels and human iEEG—depends on the statistical procedure used to declare hierarchy and significance (edge-wise tests, multiple-comparison control across frequencies and pairs, controls for session/lab/animal structure, and robustness to free parameters such as completion rank and band definitions). None of that procedure is stated in the abstract; if it is weak or confounded by the multi-lab design, the hierarchy claim does not hold even if the estimator is consistent under ideal conditions.","section":"Abstract (theta-band hierarchy result)"}],"minor_comments":[{"comment":"The abstract packs method, guarantees, four validation regimes, and a cross-species hierarchy into a single dense paragraph; a clearer separation of (i) estimator, (ii) theoretical conditions, (iii) validation design, and (iv) scientific findings would help readers assess each claim independently.","section":"Abstract"},{"comment":"Free parameters implied by the pipeline (nuclear-norm regularization or target rank, stitching weights / overlap thresholds, frequency-band definitions) are not named; even an abstract-level statement of what is fixed vs. tuned would reduce circularity concern.","section":"Abstract"},{"comment":"Terminology: “stitched power-spectra” and “global spectral matrix” should be clarified as complex Hermitian cross-spectral density matrices (not real power spectra alone), since phase is essential for PDC.","section":"Abstract"}],"recommendation":"uncertain","confidential_remarks":"This is an abstract-only review; the full text, proofs, figures, and tables were not available. Soundness cannot be scored above a provisional level. The skeptic concern about nuclear-norm completion of complex cross-spectra is the right load-bearing question and should be the first item checked if/when the full manuscript is supplied. Recommendation is uncertain pending full methods (consistency theorem statement, held-out co-recording recovery metrics, and hierarchy statistics). Fit for a methods-heavy q-bio.NC / computational neuroscience venue is high if those pieces land."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The one thing to know: this abstract pitches a genuine methods gap and a coherent pipeline for it. Existing Granger/DCM/PDC tools need simultaneous common-clock coverage; modern multi-session multi-lab data (IBL Neuropixels, multi-patient iEEG) do not. SPIDER stitches local power spectra from overlapping subsets, nuclear-norm-completes never-co-observed pairs, then does spectral factorization + PDC. If that works, it is useful for systems neuroscience.\n\nWhat looks new and solid on the page: the problem statement is clean, the pipeline is non-parametric and frequency-resolved, and they claim consistency guarantees plus validation on simulations, two-photon, IBL (50 areas / 43 sessions / 12 labs), and human iEEG. The empirical hook—largely recurrent spontaneous flow, but a significant theta feedforward hierarchy with hippocampal formation as source, recovered across mouse Neuropixels and human iEEG—is the kind of result people would actually use if the method holds. No invented entities; free parameters (completion regularization/rank, band definitions, stitching weights) are the usual ones for this class of estimator.\n\nSoft spots, in proportion: we only have the abstract. The load-bearing step is exactly what the stress-test flags—nuclear-norm completion of complex cross-spectra must preserve phase structure well enough that factorization yields correct PDC for pairs never observed together. Asynchronous sessions give no shared clock, so absolute phases between non-overlapping blocks are never measured. If the low-rank / stationarity / overlap conditions are too strong or unstated, the whole-brain maps and the hierarchy can be artifactual. Circular construction of the hierarchy is not forced by the abstract description (it is presented as an empirical output), but residual risk remains if completion hyperparameters or band choices are tuned to the finding. None of the proofs, overlap geometry, error bars, baselines, or code are inspectable here, so soundness is provisional.\n\nWho it is for: people who already work with multi-lab electrophysiology or imaging and want directed connectivity without requiring a single simultaneous recording. It deserves a serious referee if the full paper ships the guarantees, the overlap statistics, and honest failure cases. I would not desk-reject on the abstract alone; I would send it out and demand the math and the ablation on completion rank/phase recovery. Reading group: maybe, once the PDF is up. Cite in the next year: only after checking the full methods.","headline":"Abstract-only methods claim for directed connectivity on incomplete asynchronous multi-lab data; real gap, coherent pipeline, but soundness of stitching/completion for PDC is uncheckable here.","tokens_in":3055,"tokens_out":602,"would_cite":false,"duration_ms":6184,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"SPIDER recovers frequency-resolved directed brain connectivity from incomplete asynchronous multi-session recordings that never share a clock.","keywords":["effective connectivity","partial directed coherence","spectral factorization","power spectral density","asynchronous recordings","nuclear-norm completion","theta hierarchy","Neuropixels"],"falsifier":"On a held-out set of simultaneously recorded regions, compare SPIDER-recovered PDC (built only from incomplete asynchronous sessions) against ground-truth PDC computed from the simultaneous data; systematic mismatch in direction or frequency peak would falsify the claim.","tokens_in":3024,"feed_emoji":"🧠","tokens_out":770,"duration_ms":8118,"temperature":0.7,"pith_summary":"Mapping how information flows between brain regions usually requires every region to be recorded at the same time with a shared clock. Real experiments never deliver that: different sessions, animals, and labs cover different partially overlapping regions, so classical directed-connectivity tools cannot be applied. SPIDER stitches together the local power spectra estimated from those overlapping subsets into one global spectral matrix, uses nuclear-norm completion for region pairs that were never recorded together, and then recovers frequency-resolved directed interactions by spectral factorization and partial directed coherence. The method comes with consistency guarantees and is validated on simulations, calcium imaging, and large Neuropixels and human intracranial-EEG collections. If it works as claimed, multi-session atlases that were previously usable only for undirected or pairwise statistics can now yield whole-brain effective-connectivity maps, including directed hierarchies that no single recording can see.","feed_headline":"SPIDER maps directed brain flow from recordings that never share a clock","feed_subtitle":"Stitched spectra and nuclear-norm completion recover theta-band hierarchies across mouse and human datasets","key_machinery":"Stitched global spectral matrix: local power-spectral density estimates from overlapping channel subsets are assembled into one consistent matrix; nuclear-norm completion fills never-co-observed pairs; canonical spectral factorization then yields partial directed coherence (PDC) that quantifies directed frequency-resolved interactions.","core_discovery":"SPIDER recovers frequency-resolved directed information flow among brain regions that were never recorded together, by stitching local power-spectral estimates from partially overlapping multi-session recordings into a global spectral matrix, completing never-co-observed pairs by nuclear norm, and applying spectral factorization plus partial directed coherence, without any shared temporal reference.","pith_inferences":["If the low-rank structure of neural spectra is generic, SPIDER-style completion could be applied to other incomplete multi-modal datasets (e.g., combined fMRI and electrophysiology) that share only partial anatomical coverage.","Failure modes of nuclear-norm completion on spectra with strong high-rank noise would most likely appear first as inverted or frequency-smeared directed edges rather than as global scale errors.","The recovered hierarchy supplies a concrete prediction for closed-loop optogenetic or stimulation experiments: perturbing hippocampal theta should cascade more strongly into downstream areas than the reverse."],"forward_implications":["Whole-brain effective connectivity can be estimated from existing multi-session, multi-lab Neuropixels and iEEG collections that lack shared clocks.","Spontaneous activity is largely recurrent across frequencies, yet the theta band forms a significant feedforward hierarchy with hippocampal formation as source.","The same theta-band hippocampal-source hierarchy appears in both mouse Neuropixels and human resting iEEG, suggesting cross-species, cross-modality conservation.","Directed-flow maps become available for any multi-animal atlas whose regions have partial spatial overlap, without new simultaneous recordings."],"fun_headline_variants":["SPIDER stitches spectra to map directed flow without shared clocks","Directed brain connectivity from asynchronous multi-session recordings","SPIDER recovers theta hierarchies among regions never recorded together","Nuclear-norm completion yields whole-brain directed flow from partial data","SPIDER infers frequency-resolved flow across non-overlapping brain datasets"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"Local power spectra from partially overlapping subsets can be stitched and completed so that the resulting global spectral matrix still preserves the true directed phase and cross-spectral relationships needed for PDC.","fun_headline_variants_meta":{"raw":{"variants":["SPIDER stitches spectra to map directed flow without shared clocks","Directed brain connectivity from asynchronous multi-session recordings","SPIDER recovers theta hierarchies among regions never recorded together","Nuclear-norm completion yields whole-brain directed flow from partial data","SPIDER infers frequency-resolved flow across non-overlapping brain datasets"]},"model":"grok-4.5","effort":"low","cost_usd":0.004228,"raw_usage":{"total_tokens":1320,"prompt_tokens":827,"num_sources_used":0,"completion_tokens":66,"cost_in_usd_ticks":42280000,"prompt_tokens_details":{"text_tokens":827,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":427,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":827,"tokens_out":66,"duration_ms":4027,"temperature":1.0,"reasoning_tokens":427,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-12T12:50:12.699077+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"On a held-out set of simultaneously recorded regions, compare SPIDER-recovered PDC (built only from incomplete asynchronous sessions) against ground-truth PDC computed from the simultaneous data; systematic mismatch in direction or frequency peak would falsify the claim.","supporting_citations":[],"review_version":2}