The superior temporal sulcus (STS) encompasses a complex set of regions involved in a wide range of cognitive functions. To understand its functional properties, neuromodulation approaches such brain stimulation or neurofeedback can be used. We investigated whether the posterior STS (pSTS), a core region in the face perception and imagery network, could be specifically identified based on the presence of dynamic facial expressions (and not just on simple motion or static face signals), and probed with neurofeedback. Recognition of facial expressions is critically impaired in autism spectrum disorder, making this region a relevant target for future clinical neurofeedback studies. We used a stringent localizer approach based on the contrast of dynamic facial expressions against static neutral faces plus moving dots. The target region had to be specifically responsive to dynamic facial expressions instead of mere motion and/or the presence of a static face. The localizer was successful in selecting this region across subjects. Neurofeedback was then performed, using this region as a target, with two novel feedback rules (mean or derivative-based, using visual or auditory interfaces). Our results provide evidence that a facial expression-selective cluster in pSTS can be identified and may represent a suitable target for neurofeedback approaches, aiming at social and emotional cognition. These findings highlight the presence of a highly selective region in STS encoding dynamic aspects of facial expressions. Future studies should elucidate its role as a mechanistic target for neurofeedback strategies in clinical disorders of social cognition such as autism.

FMRI-based neurofeedback (NF) interventions represent the method of choice for the neuromodulation of localized brain areas. Although we have already validated an fMRI-NF protocol targeting the facial expressions processing network (FEPN), its dissemination is hampered by the economical and logistical constraints of fMRI-NF interventions, which may be however surpassed by transferring it to EEG setups, due to their low cost and portability. One of the major challenges of this procedure is then to reconstruct the BOLD-fMRI signal measured at the FEPN using only EEG signals. Because these types of approaches have been poorly explored so far, here we systematically investigated the extent at which the BOLD-fMRI signal recorded from the FEPN during a fMRI-NF protocol could be reconstructed from the simultaneously recorded EEG signal. Several features from both scalp and source spaces (the latter estimated using continuous EEG source imaging) were extracted and used as predictors in a regression problem using random forests. Furthermore, three different approaches to deal with the hemodynamic delay of the BOLD signal where tested. The resulting models where compared with the only approach already proposed in the literature that uses spectral features and considers different time delays. Our results show that the combination of different kinds of features extracted from the scalp, and convolving them with multiple HRF functions peaking at different latencies, increases significantly the reconstruction accuracy (defined as the correlation between the measured and reconstructed BOLD signal) from 20% (the current literature) to 52%. With this pipeline, a more accurate reconstruction of the BOLD signal can be obtained, which will positively impact the transfer of fMRI-based neurofeedback interventions to EEG setups, and more importantly, their dissemination and efficacy in modulating the activity of the desired brain areas.