Mechanistic Mathematical Modeling Tests Hypotheses of the Neurovascular Coupling in fMRI

Functional magnetic resonance imaging (fMRI) measures brain activity by detecting the blood-oxygen-level dependent (BOLD) response to neural activity. The BOLD response depends on the neurovascular coupling, which connects cerebral blood flow, cerebral blood volume, and deoxyhemoglobin level to neuronal activity. The exact mechanisms behind this neurovascular coupling are not yet fully investigated. There are at least three different ways in which these mechanisms are being discussed. Firstly, mathematical models involving the so-called Balloon model describes the relation between oxygen metabolism, cerebral blood volume, and cerebral blood flow. However, the Balloon model does not describe cellular and biochemical mechanisms. Secondly, the metabolic feedback hypothesis, which is based on experimental findings on metabolism associated with brain activation, and thirdly, the neurotransmitter feed-forward hypothesis which describes intracellular pathways leading to vasoactive substance release. Both the metabolic feedback and the neurotransmitter feed-forward hypotheses have been extensively studied, but only experimentally. These two hypotheses have never been implemented as mathematical models. Here we investigate these two hypotheses by mechanistic mathematical modeling using a systems biology approach; these methods have been used in biological research for many years but never been applied to the BOLD response in fMRI. In the current work, model structures describing the metabolic feedback and the neurotransmitter feed-forward hypotheses were applied to measured BOLD responses in the visual cortex of 12 healthy volunteers. Evaluating each hypothesis separately shows that neither hypothesis alone can describe the data in a biologically plausible way. However, by adding metabolism to the neurotransmitter feed-forward model structure, we obtained a new model structure which is able to fit the estimation data and successfully predict new, independent validation data. These results open the door to a new type of fMRI analysis that more accurately reflects the true neuronal activity.     

Resting State Brain Function Analysis Using Concurrent BOLD in ASL Perfusion fMRI

The past decade has seen astounding discoveries about resting-state brain activity patterns in normal brain as well as their alterations in brain diseases. While the vast majority of resting-state studies are based on the blood-oxygen-level-dependent (BOLD) functional MRI (fMRI), arterial spin labeling (ASL) perfusion fMRI can simultaneously capture BOLD and cerebral blood flow (CBF) signals, providing a unique opportunity for assessing resting brain functions with concurrent BOLD (ccBOLD) and CBF signals. Before taking that benefit, it is necessary to validate the utility of ccBOLD signal for resting-state analysis using conventional BOLD (cvBOLD) signal acquired without ASL modulations. To address this technical issue, resting cvBOLD and ASL perfusion MRI were acquired from a large cohort (n = 89) of healthy subjects. Four widely used resting-state brain function analyses were conducted and compared between the two types of BOLD signal, including the posterior cingulate cortex (PCC) seed-based functional connectivity (FC) analysis, independent component analysis (ICA), analysis of amplitude of low frequency fluctuation (ALFF), and analysis of regional homogeneity (ReHo). Consistent default mode network (DMN) as well as other resting-state networks (RSNs) were observed from cvBOLD and ccBOLD using PCC-FC analysis and ICA. ALFF from both modalities were the same for most of brain regions but were different in peripheral regions suffering from the susceptibility gradients induced signal drop. ReHo showed difference in many brain regions, likely reflecting the SNR and resolution differences between the two BOLD modalities. The DMN and auditory networks showed highest CBF values among all RSNs. These results demonstrated the feasibility of ASL perfusion MRI for assessing resting brain functions using its concurrent BOLD in addition to CBF signal, which provides a potentially useful way to maximize the utility of ASL perfusion MRI.     

Assessment of single-vessel cerebral blood velocity by phase contrast fMRI

Current approaches to high-field functional MRI (fMRI) provide 2 means to map hemodynamics at the level of single vessels in the brain. One is through changes in deoxyhemoglobin in venules, i.e., blood oxygenation level–dependent (BOLD) fMRI, while the second is through changes in arteriole diameter, i.e., cerebral blood volume (CBV) fMRI. Here, we introduce cerebral blood flow–related velocity-based fMRI, denoted CBFv-fMRI, which uses high-resolution phase contrast (PC) MRI to form velocity measurements of flow. We use CBFv-fMRI in measure changes in blood velocity in single penetrating microvessels across rat parietal cortex. In contrast to the venule-dominated BOLD and arteriole-dominated CBV fMRI signals, CBFv-fMRI is comparable from both arterioles and venules. A single fMRI platform is used to map changes in blood pO₂ (BOLD), volume (CBV), and velocity (CBFv). This combined high-resolution single-vessel fMRI mapping scheme enables vessel-specific hemodynamic mapping in animal models of normal and diseased states and further has translational potential to map vascular dementia in diseased or injured human brains with ultra–high-field fMRI.

This study presents a phase contrast-based, high field MRI-based approach for the functional mapping of cerebral blood velocity in individual cortical arterioles and venules in the rat cortex; this approach can be combined with previously established approaches to map BOLD, CBV, and blood velocity from penetrating microvessels.     

Relations between BOLD fMRI-Derived Resting Brain Activity and Cerebral Blood Flow

Consistent resting brain activity patterns have been repeatedly demonstrated using measures derived from resting BOLD fMRI data. While those metrics are presumed to reflect underlying spontaneous brain activity (SBA), it is challenging to prove that association because resting BOLD fMRI metrics are purely model-free and scale-free variables. Cerebral blood flow (CBF) is typically closely coupled to brain metabolism and is used as a surrogate marker for quantifying regional brain function, including resting function. Assessing the correlations between resting BOLD fMRI measures and CBF correlation should provide a means of linking of those measures to the underlying SBA, and a means to quantify those scale-free measures. The purpose of this paper was to examine the CBF correlations of 3 widely used neuroimaging-based SBA measures, including seed-region based functional connectivity (FC), regional homogeneity (ReHo), and amplitude of low frequency fluctuation (ALFF). Test-retest data were acquired to check the stability of potential correlations across time. Reproducible posterior cingulate cortex (PCC) FC vs regional CBF correlations were found in much of the default mode network and visual cortex. Dorsal anterior cingulate cortex (ACC) FC vs CBF correlations were consistently found in bilateral prefrontal cortex. Both ReHo and ALFF were found to be reliably correlated with CBF in most of brain cortex. None of the assessed SBA measures was correlated with whole brain mean CBF. These findings suggest that resting BOLD fMRI-derived measures are coupled with regional CBF and are therefore linked to regional SBA.     

Brain serotonin transporter binding of [123I]ADAM: within-subject variation between summer and winter data

The neurotransmitter serotonin (5-HT) controls several physiological functions, and a disturbance of the 5-HT system is implicated in many psychiatric conditions. Seasonal variation has been suggested in the 5-HT system. We investigated within-subject seasonal variation in brain serotonin transporter (SERT) binding with the SERT-ligand [(123)I]ADAM and single photon emission computed tomography (SPECT) in 12 healthy individuals. No systematic variation was found in the midbrain or thalamus areas between scans done in summer and winter. Our results suggest that factors other than season are more important in causing within-subject variation of brain SERT binding between summer and winter.     

Positron emission tomography ligand activation studies in the sports sciences: measuring neurochemistry in vivo

Functional neuroimaging with magnetic resonance imaging (fMRI) or positron emission tomography (PET) provides the methodology to unravel some of the fascinating, but hitherto largely unresolved interactions between physical exercise and brain function. Phenomena such as raised mood, pain modulation, and sport addiction associated with physical exercise are highly interesting psychophysical models that require further in depth understanding at the neurotransmitter level. PET ligand displacement studies allow in vivo monitoring of endogenous transmitter trafficking in the entire brain and, thereby, to identify the link between exercise-induced behavioral measures and the endogenous neurotransmitter release. This review focuses on the methodology of ligand displacement in the opioidergic system, which together with the dopaminergic system has been considered as a central neurotransmitter system underlying diverse sport-induced psychophysical effects. Understanding the basic principles of exercise-induced transmitter release in the brain will potentially aid clinical applications of endurance training, both as a preventative or therapeutic intervention.     

Functional mapping in the human brain using high magnetic fields.

An avidly pursued new dimension in magnetic resonance imaging (MRI) research is the acquisition of physiological and biochemical information non-invasively using the nuclear spins of the water molecules in the human body. In this trial, a recent and unique accomplishment was the introduction of the ability to map human brain function non-invasively. Today, functional images with subcentimetre resolution of the entire human brain can be generated in single subjects and in data acquisition times of several minutes using 1.5 tesla (T) MRI scanners that are often used in hospitals for clinical purposes. However, there have been accomplishments beyond this type of imaging using significantly higher magnetic fields such as 4 T. Efforts for developing high magnetic field human brain imaging and functional mapping using MRI (fMRI) were undertaken at about the same time. It has been demonstrated that high magnetic fields result in improved contrast and, more importantly, in elevated sensitivity to capillary level changes coupled to neuronal activity in the blood oxygenation level dependent (BOLD) contrast mechanism used in fMRI. These advantages have been used to generate, for example, high resolution functional maps of ocular dominance columns, retinotopy within the small lateral geniculate nucleus, true single-trial fMRI and early negative signal changes in the temporal evolution of the BOLD signal. So far these have not been duplicated or have been observed as significantly weaker effects at much lower field strengths. Some of these high-field advantages and accomplishments are reviewed in this paper.     

The neural basis of the blood-oxygen-level-dependent functional magnetic resonance imaging signal

Magnetic resonance imaging (MRI) has rapidly become an important tool in clinical medicine and biological research. Its functional variant (functional magnetic resonance imaging; fMRI) is currently the most widely used method for brain mapping and studying the neural basis of human cognition. While the method is widespread, there is insufficient knowledge of the physiological basis of the fMRI signal to interpret the data confidently with respect to neural activity. This paper reviews the basic principles of MRI and fMRI, and subsequently discusses in some detail the relationship between the blood-oxygen-level-dependent (BOLD) fMRI signal and the neural activity elicited during sensory stimulation. To examine this relationship, we conducted the first simultaneous intracortical recordings of neural signals and BOLD responses. Depending on the temporal characteristics of the stimulus, a moderate to strong correlation was found between the neural activity measured with microelectrodes and the BOLD signal averaged over a small area around the microelectrode tips. However, the BOLD signal had significantly higher variability than the neural activity, indicating that human fMRI combined with traditional statistical methods underestimates the reliability of the neuronal activity. To understand the relative contribution of several types of neuronal signals to the haemodynamic response, we compared local field potentials (LFPs), single- and multi-unit activity (MUA) with high spatio-temporal fMRI responses recorded simultaneously in monkey visual cortex. At recording sites characterized by transient responses, only the LFP signal was significantly correlated with the haemodynamic response. Furthermore, the LFPs had the largest magnitude signal and linear systems analysis showed that the LFPs were better than the MUAs at predicting the fMRI responses. These findings, together with an analysis of the neural signals, indicate that the BOLD signal primarily measures the input and processing of neuronal information within a region and not the output signal transmitted to other brain regions.     

Functional magnetic resonance imaging: imaging techniques and contrast mechanisms.

Functional magnetic resonance imaging (fMRI) is a widely used technique for generating images or maps of human brain activity. The applications of the technique are widespread in cognitive neuroscience and it is hoped they will eventually extend into clinical practice. The activation signal measured with fMRI is predicated on indirectly measuring changes in the concentration of deoxyhaemoglobin which arise from an increase in blood oxygenation in the vicinity of neuronal firing. The exact mechanisms of this blood oxygenation level dependent (BOLD) contrast are highly complex. The signal measured is dependent on both the underlying physiological events and the imaging physics. BOLD contrast, although sensitive, is not a quantifiable measure of neuronal activity. A number of different imaging techniques and parameters can be used for fMRI, the choice of which depends on the particular requirements of each functional imaging experiment. The high-speed MRI technique, echo-planar imaging provides the basis for most fMRI experiments. The problems inherent to this method and the ways in which these may be overcome are particularly important in the move towards performing functional studies on higher field MRI systems. Future developments in techniques and hardware are also likely to enhance the measurement of brain activity using MRI.     

La stimulation électrique pallidale dans le traitement des mouvements anormaux : l’expérience montpelliéraine dans l’utilisation de l’imagerie fonctionnelle comme critères de sélection des patients

Magnetic resonance anatomical imaging gives to the practician an important set of criteria for selecting patients with dystonic and dyskinetic syndromes who can get benefice from a treatment by electrical deep brain stimulation. However, anatomical anomalies seen on MRI, whatever their number and morphological extension, do not reflect the effective (residual) functional activity of the corresponding structures. We describe, in this short report, the way how, in our experience, functional imaging techniques used for measuring cerebral perfusion and neurotransmission in Single Photon Emission Computerized Tomography (SPECT), cerebral metabolism in Positron Emission Tomography (PET) and motor neuroactivations in functional Magnetic Resonance Imaging (fMRI) are capable of progressively fulfil this deficiency.     

Simultaneous fMRI and Electrophysiology in the Rodent Brain

To examine the neural basis of the blood oxygenation level dependent (BOLD) magnetic resonance imaging (MRI) signal, we have developed a rodent model in which functional MRI data and in vivo intracortical recording can be performed simultaneously. The combination of MRI and electrical recording is technically challenging because the electrodes used for recording distort the MRI images and the MRI acquisition induces noise in the electrical recording. To minimize the mutual interference of the two modalities, glass microelectrodes were used rather than metal and a noise removal algorithm was implemented for the electrophysiology data. In our studies, two microelectrodes were separately implanted in bilateral primary somatosensory cortices (SI) of the rat and fixed in place. One coronal slice covering the electrode tips was selected for functional MRI. Electrode shafts and fixation positions were not included in the image slice to avoid imaging artifacts. The removed scalp was replaced with toothpaste to reduce susceptibility mismatch and prevent Gibbs ringing artifacts in the images. The artifact structure induced in the electrical recordings by the rapidly-switching magnetic fields during image acquisition was characterized by averaging all cycles of scans for each run. The noise structure during imaging was then subtracted from original recordings. The denoised time courses were then used for further analysis in combination with the fMRI data. As an example, the simultaneous acquisition was used to determine the relationship between spontaneous fMRI BOLD signals and band-limited intracortical electrical activity. Simultaneous fMRI and electrophysiological recording in the rodent will provide a platform for many exciting applications in neuroscience in addition to elucidating the relationship between the fMRI BOLD signal and neuronal activity.Download video file.^{(95M, mp4)}     

Reproducibility of swallow-induced cortical BOLD positive and negative fMRI activity

Functional MRI (fMRI) studies have demonstrated that a number of brain regions (cingulate, insula, prefrontal, and sensory/motor cortices) display blood oxygen level-dependent (BOLD) positive activity during swallow. Negative BOLD activations and reproducibility of these activations have not been systematically studied. The aim of our study was to investigate the reproducibility of swallow-related cortical positive and negative BOLD activity across different fMRI sessions. We studied 16 healthy volunteers utilizing an fMRI event-related analysis. Individual analysis using a general linear model was used to remove undesirable signal changes correlated with motion, white matter, and cerebrospinal fluid. The group analysis used a mixed-effects multilevel model to identify active cortical regions. The volume and magnitude of a BOLD signal within each cluster was compared between the two study sessions. All subjects showed significant clustered BOLD activity within the known areas of cortical swallowing network across both sessions. The cross-correlation coefficient of percent fMRI signal change and the number of activated voxels across both positive and negative BOLD networks were similar between the two studies (r ≥ 0.87, P < 0.0001). Swallow-associated negative BOLD activity was comparable to the well-defined "default-mode" network, and positive BOLD activity had noticeable overlap with the previously described "task-positive" network. Swallow activates two parallel cortical networks. These include a positive and a negative BOLD network, respectively, correlated and anticorrelated with swallow stimulus. Group cortical activity maps, as well as extent and amplitude of activity induced by volitional swallowing in the cortical swallowing network, are reproducible between study sessions.     

Imaging neurodegeneration in Parkinson's disease

Neuroimaging techniques have evolved over the past several years giving us unprecedented information about the degenerative process in Parkinson's disease (PD) and other movement disorders. Functional imaging approaches such as positron emission tomography (PET) and single photon emission computerised tomography (SPECT) have been successfully employed to detect dopaminergic dysfunction in PD, even while at a preclinical stage, and to demonstrate the effects of therapies on function of intact dopaminergic neurons within the affected striatum. PET and SPECT can also monitor PD progression as reflected by changes in brain levodopa and glucose metabolism and dopamine transporter binding. Structural imaging approaches include magnetic resonance imaging (MRI) and transcranial sonography (TCS). Recent advances in voxel-based morphometry and diffusion-weighted MRI have provided exciting potential applications for the differential diagnosis of parkinsonian syndromes. Substantia nigra hyperechogenicity, detected with TCS, may provide a marker of susceptibility to PD, probably reflecting disturbances of iron metabolism, but does not appear to correlate well with disease severity or change with disease progression. In the future novel radiotracers may help us assess the involvement of non-dopaminergic brain pathways in the pathology of both motor and non-motor complications in PD.     

利用功能磁共振研究记忆功能的研究进展

功能磁共振成像(fMRI)作为一种无创伤,可反复实验的成像技术,已被广泛地应用于各项脑功能的研究。用fMRI进行脑功能研究的主要依据是血流敏感性和BOLD对比增强原理。记忆是人脑的高级功能,其过程分为编码加工、固化、存储和提取几个阶段。大脑皮质、海马、乳头体、丘脑是参与记忆的主要解剖结构。记忆的刺激方式对各个脑区的激活是具有差异性的,记忆功能的测量和分析方法也在不断的改进。年龄和性别的不同都会对记忆能力产生影响,同时激活的脑区也会有相应的改变。此外,情感和记忆的关系正越来越受到人们的关注。本文阐述利用功能磁共振成像研究记忆功能的最新进展。     

A New Functional MRI Approach for Investigating Modulations of Brain Oxygen Metabolism

Functional MRI (fMRI) using the blood oxygenation level dependent (BOLD) signal is a common technique in the study of brain function. The BOLD signal is sensitive to the complex interaction of physiological changes including cerebral blood flow (CBF), cerebral blood volume (CBV), and cerebral oxygen metabolism (CMRO₂). A primary goal of quantitative fMRI methods is to combine BOLD imaging with other measurements (such as CBF measured with arterial spin labeling) to derive information about CMRO₂. This requires an accurate mathematical model to relate the BOLD signal to the physiological and hemodynamic changes; the most commonly used of these is the Davis model. Here, we propose a new nonlinear model that is straightforward and shows heuristic value in clearly relating the BOLD signal to blood flow, blood volume and the blood flow-oxygen metabolism coupling ratio. The model was tested for accuracy against a more detailed model adapted for magnetic fields of 1.5, 3 and 7T. The mathematical form of the heuristic model suggests a new ratio method for comparing combined BOLD and CBF data from two different stimulus responses to determine whether CBF and CMRO₂ coupling differs. The method does not require a calibration experiment or knowledge of parameter values as long as the exponential parameter describing the CBF-CBV relationship remains constant between stimuli. The method was found to work well for 1.5 and 3T but is prone to systematic error at 7T. If more specific information regarding changes in CMRO₂ is required, then with accuracy similar to that of the Davis model, the heuristic model can be applied to calibrated BOLD data at 1.5T, 3T and 7T. Both models work well over a reasonable range of blood flow and oxygen metabolism changes but are less accurate when applied to a simulated caffeine experiment in which CBF decreases and CMRO₂ increases.     

Hemodynamic traveling waves in human visual cortex

Functional MRI (fMRI) experiments rely on precise characterization of the blood oxygen level dependent (BOLD) signal. As the spatial resolution of fMRI reaches the sub-millimeter range, the need for quantitative modelling of spatiotemporal properties of this hemodynamic signal has become pressing. Here, we find that a detailed physiologically-based model of spatiotemporal BOLD responses predicts traveling waves with velocities and spatial ranges in empirically observable ranges. Two measurable parameters, related to physiology, characterize these waves: wave velocity and damping rate. To test these predictions, high-resolution fMRI data are acquired from subjects viewing discrete visual stimuli. Predictions and experiment show strong agreement, in particular confirming BOLD waves propagating for at least 5-10 mm across the cortical surface at speeds of 2-12 mm s-1. These observations enable fundamentally new approaches to fMRI analysis, crucial for fMRI data acquired at high spatial resolution.     

A Dual Tracer PET-MRI Protocol for the Quantitative Measure of Regional Brain Energy Substrates Uptake in the Rat

We present a method for comparing the uptake of the brain''s two key energy substrates: glucose and ketones (acetoacetate [AcAc] in this case) in the rat. The developed method is a small-animal positron emission tomography (PET) protocol, in which ¹¹C-AcAc and ¹⁸F-fluorodeoxyglucose (¹⁸F-FDG) are injected sequentially in each animal. This dual tracer PET acquisition is possible because of the short half-life of ¹¹C (20.4 min). The rats also undergo a magnetic resonance imaging (MRI) acquisition seven days before the PET protocol. Prior to image analysis, PET and MRI images are coregistered to allow the measurement of regional cerebral uptake (cortex, hippocampus, striatum, and cerebellum). A quantitative measure of ¹¹C-AcAc and ¹⁸F-FDG brain uptake (cerebral metabolic rate; μmol/100 g/min) is determined by kinetic modeling using the image-derived input function (IDIF) method. Our new dual tracer PET protocol is robust and flexible; the two tracers used can be replaced by different radiotracers to evaluate other processes in the brain. Moreover, our protocol is applicable to the study of brain fuel supply in multiple conditions such as normal aging and neurodegenerative pathologies such as Alzheimer''s and Parkinson''s diseases.     

Simultaneous BOLD fMRI and fiber-optic calcium recording in rat neocortex

Functional magnetic resonance imaging (fMRI) based on blood oxygen level-dependent (BOLD) contrast is widely used for probing brain activity, but its relationship to underlying neural activity remains elusive. Here, we combined fMRI with fiber-optic recordings of fluorescent calcium indicator signals to investigate this relationship in rat somatosensory cortex. Electrical forepaw stimulation (1-10 Hz) evoked fast calcium signals of neuronal origin that showed frequency-dependent adaptation. Additionally, slower calcium signals occurred in astrocyte networks, as verified by astrocyte-specific staining and two-photon microscopy. Without apparent glia activation, we could predict BOLD responses well from simultaneously recorded fiber-optic signals, assuming an impulse response function and taking into account neuronal adaptation. In cases with glia activation, we uncovered additional prolonged BOLD signal components. Our findings highlight the complexity of fMRI BOLD signals, involving both neuronal and glial activity. Combined fMRI and fiber-optic recordings should help to clarify cellular mechanisms underlying BOLD signals.     

Effect of induced high myopia on functional MRI signal changes

PurposeThe current study evaluated the effect of lens-induced high myopia (IHM) on the activity of the occipital visual cortex during two visual stimuli presentations to the subjects. This was done by measuring the Blood Oxygenation Level Dependent (BOLD) signal using functional MRI (fMRI).MethodsBOLD contrast fMRI was performed with a 1.5T MRI scanner on 12 emmetropic subjects (refractive error <±0.25Diopter) with no history of neurologic disorder. IHM conditions were applied to subjects by three convex lenses of +5D, +7D and +10D. Visual stimuli with 0.34 cpd and 1.84 cpd spatial frequencies (SF) were presented as a block paradigm to the participants in three IHM states and normal vision state during fMRI data acquisition. Resultant fMRI data were compared among different refractive states.ResultsData analysis showed that IHM did not cause a significant change in the visual cortex activity throughout the presentation of 0.34 cpd SF visual stimulus and BOLD signal intensity remained approximately constant (p = 0.17). Although, fMRI responses to visual stimuli with spatial frequency of 1.84 cpd demonstrated that visual cortex activity was significantly reduced in IHM states compared to normal vision (p = 0.01), the results showed no significant differences between three different values of IHM.ConclusionsThis study shows severe blurring caused by lens induced high myopia can decrease BOLD signal intensity depending on the visual stimulus pattern details. However in the low and moderate range of spatial frequencies, blur increment from +5D up to +10D is not associated with further reduction in the BOLD signal of the occipital visual cortex.     

Interictal functional connectivity of human epileptic networks assessed by intracerebral EEG and BOLD signal fluctuations

In this study, we aimed to demonstrate whether spontaneous fluctuations in the blood oxygen level dependent (BOLD) signal derived from resting state functional magnetic resonance imaging (fMRI) reflect spontaneous neuronal activity in pathological brain regions as well as in regions spared by epileptiform discharges. This is a crucial issue as coherent fluctuations of fMRI signals between remote brain areas are now widely used to define functional connectivity in physiology and in pathophysiology. We quantified functional connectivity using non-linear measures of cross-correlation between signals obtained from intracerebral EEG (iEEG) and resting-state functional MRI (fMRI) in 5 patients suffering from intractable temporal lobe epilepsy (TLE). Functional connectivity was quantified with both modalities in areas exhibiting different electrophysiological states (epileptic and non affected regions) during the interictal period. Functional connectivity as measured from the iEEG signal was higher in regions affected by electrical epileptiform abnormalities relative to non-affected areas, whereas an opposite pattern was found for functional connectivity measured from the BOLD signal. Significant negative correlations were found between the functional connectivities of iEEG and BOLD signal when considering all pairs of signals (theta, alpha, beta and broadband) and when considering pairs of signals in regions spared by epileptiform discharges (in broadband signal). This suggests differential effects of epileptic phenomena on electrophysiological and hemodynamic signals and/or an alteration of the neurovascular coupling secondary to pathological plasticity in TLE even in regions spared by epileptiform discharges. In addition, indices of directionality calculated from both modalities were consistent showing that the epileptogenic regions exert a significant influence onto the non epileptic areas during the interictal period. This study shows that functional connectivity measured by iEEG and BOLD signals give complementary but sometimes inconsistent information in TLE.