TE-aware state analysis and evolution in dynamic network connectivity

Multi-echo functional magnetic resonance imaging (fMRI) acquires signals with distinct contrast profiles across echo times (TEs). Empirical evidence suggests that different TEs capture distinct signal contributions with varying sensitivity to BOLD and non-BOLD processes from long to short TEs (Stroman, 2002; Krishnamurthy, 2023; Dong, 2024). Previous literature has not yet examined how brain states extracted from temporal dynamic functional connectivity (dFNC) evolve across TEs. Here we present a methodological and empirical study of TE-dependent evolution of dynamic connectivity states, with a phenomenological two-component model for interpretation. Multiple echoes were acquired with an echo planar time-resolved imaging (EPTI) sequence. TE-specific components were extracted via atlas-based group ICA (Neuromark 1.0), and brain connectivity states were generated via windowed dFNC and spatial ICA clustering of windows. Within-subject and group-level states evolved across TEs in network organization and variance. States exhibited TE-dependent variance profiles, with distinct peak echo times; and group-level points of state transition were independently mapped to TE-dependent transitions in variance. We modified a bi-exponential model to simulate and model changes in the dFNC time-course across TEs. The bi-exponential model revealed high consistency with subject states (R2 fit=0.86, res. error=0.55, 86% model convergence rate), suggesting a two-component signal system contributing to dynamic TE-dependent connectivity changes in state. A significant correlation between relative model terms was estimated on a cell-wise basis across states (r(24) = -0.43, p=0.03, 95% CI [-0.70, -0.04]). This suggests that specific bi-exponential signal parameters (i.e., short FNC1,TE=0, long FNC2,TE=0, and long FNC2,decay) were associated with the observed TE-dependent variations in connectivity states across TEs. In sum, our analyses provide downstream implications for TE-aware state analysis dependent on connectivity magnitude (FNCTE=0) and rate of connectivity decay (FNCdecay), in a two-component signal model, with future implications for TE-dependent state biomarkers.