Adaptive modulation of brain hemodynamics across stereotyped running episodes.

During locomotion, theta and gamma rhythms are essential to ensure timely communication between brain structures. However, their metabolic cost and contribution to neuroimaging signals remain elusive. To finely characterize neurovascular interactions during locomotion, we simultaneously recorded mesoscale brain hemodynamics using functional ultrasound (fUS) and local field potentials (LFP) in numerous brain structures of freely-running overtrained rats. Locomotion events were reliably followed by a surge in blood flow in a sequence involving the retrosplenial cortex, dorsal thalamus, dentate gyrus and CA regions successively, with delays ranging from 0.8 to 1.6 seconds after peak speed. Conversely, primary motor cortex was suppressed and subsequently recruited during reward uptake. Surprisingly, brain hemodynamics were strongly modulated across trials within the same recording session; cortical blood flow sharply decreased after 10-20 runs, while hippocampal responses strongly and linearly increased, particularly in the CA regions. This effect occurred while running speed and theta activity remained constant and was accompanied by an increase in the power of hippocampal, but not cortical, high-frequency oscillations (100-150 Hz). Our findings reveal distinct vascular subnetworks modulated across fast and slow timescales and suggest strong hemodynamic adaptation, despite the repetition of a stereotyped behavior.