The brain is active even in the absence of explicit input or output as demonstrated from electrophysiological as well as imaging studies. Using a combined approach we measured spontaneous fluctuations in the blood oxygen level dependent (BOLD) signal along with electroencephalography (EEG) in eleven healthy subjects during relaxed wakefulness (eyes closed). In contrast to other studies which used the EEG frequency information to guide the functional MRI (fMRI) analysis, we opted for transient EEG events, which identify and quantify brain electric microstates as time epochs with quasi-stable field topography. We then used this microstate information as regressors for the BOLD fluctuations. Single trial EEGs were segmented with a specific module of the LORETA (low resolution electromagnetic tomography) software package in which microstates are represented as normalized vectors constituted by scalp electric potentials, i.e., the related 3-dimensional distribution of cortical current density in the brain. Using the occurrence and the duration of each microstate, we modeled the hemodynamic response function (HRF) which revealed BOLD activation in all subjects. The BOLD activation patterns resembled well known resting-state networks (RSNs) such as the default mode network. Furthermore we "cross validated" the data performing a BOLD independent component analysis (ICA) and computing the correlation between each ICs and the EEG microstates across all subjects. This study shows for the first time that the information contained within EEG microstates on a millisecond timescale is able to elicit BOLD activation patterns consistent with well known RSNs, opening new avenues for multimodal imaging data processing.