Recovering directional brain networks under temporal undersampling: an application to schizophrenia

Two decades of functional connectivity research have established schizophrenia as a disorder of distributed dysconnectivity, with a robust thalamocortical signature: reduced prefrontal and increased sensory coupling. A fundamental issue is that functional connectivity is undirected, operates at a single slow timescale, and cannot reveal causal direction. Moreover, the mismatch between BOLD sampling speed and neural dynamics can hide edges, fabricate spurious ones, and reverse the apparent orientation of causal relationships. To overcome these limitations, we introduce a general framework for estimating directed causal graphs from fMRI that explicitly accounts for temporal undersampling. We apply RnR, a causal discovery method built on the rate-agnostic RASL framework, to resting-state fMRI from the multi-site FBIRN cohort. Rather than returning a single directed graph, RnR recovers an equivalence class of directed graphs consistent with the observed data, each annotated with the sampling rate that would produce it and classifies each estimated orientation by its stability across inferred rates. This provides a principled basis for distinguishing directed interpretations that are safe to trust from those that are timescale-contingent. Benchmarking against five single-timescale estimators on schizophrenia data, RnR recovered substantially more group-differentiating directed edges, reproducing the fields most replicated finding in directed form: a sensory-to-visual hyperconnectivity hub oriented from the post-central gyrus component to primary visual cortex. Through simulation, we show that coarse spatial resolution has shielded single-timescale methods from the full effects of undersampling, whereas finer parcellations will require explicit undersampling modeling. This work reframes fMRI effective connectivity estimation by treating undersampling as a fundamental property of the measurement, enabling directed interpretations that are grounded in the data-generating process.