IntroductionSmall animal fMRI is an essential part of translational research in the cognitive neurosciences. Due to small dimensions and animal physiology preclinical fMRI is prone to artifacts that may lead to misinterpretation of the data. To reach unbiased translational conclusions, it is, therefore, crucial to identify potential sources of experimental noise and to develop correction methods for contributions that cannot be avoided such as physiological noise. Aim of this study was to assess origin and prevalence of hemodynamic oscillations (HDO) in preclinical fMRI in rat, as well as their impact on data analysis.MethodsFollowing the development of algorithms for HDO detection and suppression, HDO prevalence in fMRI measurements was investigated for different anesthetic regimens, comprising isoflurane and medetomidine, and for both gradient echo and spin echo fMRI sequences. In addition to assessing the effect of vasodilation on HDO, it was studied if HDO have a direct neuronal correlate using local field potential (LFP) recordings. Finally, the impact of HDO on analysis of fMRI data was assessed, studying both the impact on calculation of activation maps as well as the impact on brain network analysis. Overall, 303 fMRI measurements and 32 LFP recordings were performed in 71 rats.ResultsIn total, 62% of the fMRI measurements showed HDO with a frequency of (0.20 ± 0.02) Hz. This frequent occurrence indicated that HDO cannot be generally neglected in fMRI experiments. Using the developed algorithms, HDO were detected with a specificity of 95%, and removed efficiently from the signal time courses. HDO occurred brain-wide under vasoconstrictive conditions in both small and large blood vessels. Vasodilation immediately interrupted HDO, which, however, returned within 1 h under vasoconstrictive conditions. No direct neuronal correlate of HDO was observed in LFP recordings. HDO significantly impacted analysis of fMRI data, leading to altered cluster sizes and F-values for activated voxels, as well as altered brain networks, when comparing data with and without HDO.DiscussionWe therefore conclude that HDO are caused by vasomotion under certain anesthetic conditions and should be corrected during fMRI data analysis to avoid bias.