Identification of Negative BOLD Responses in Epilepsy Using Windkessel Models

Alongside positive blood oxygenation level–dependent (BOLD) responses associated with interictal epileptic discharges, a variety of negative BOLD responses (NBRs) are typically found in epileptic patients. Previous studies suggest that, in general, up to four mechanisms might underlie the genesis of NBRs in the brain: (i) neuronal disruption of network activity, (ii) altered balance of neurometabolic/vascular couplings, (iii) arterial blood stealing, and (iv) enhanced cortical inhibition. Detecting and classifying these mechanisms from BOLD signals are pivotal for the improvement of the specificity of the electroencephalography–functional magnetic resonance imaging (EEG-fMRI) image modality to identify the seizure-onset zones in refractory local epilepsy. This requires models with physiological interpretation that furnish the understanding of how these mechanisms are fingerprinted by their BOLD responses. Here, we used a Windkessel model with viscoelastic compliance/inductance in combination with dynamic models of both neuronal population activity and tissue/blood O2 to classify the hemodynamic response functions (HRFs) linked to the above mechanisms in the irritative zones of epileptic patients. First, we evaluated the most relevant imprints on the BOLD response caused by variations of key model parameters. Second, we demonstrated that a general linear model is enough to accurately represent the four different types of NBRs. Third, we tested the ability of a machine learning classifier, built from a simulated ensemble of HRFs, to predict the mechanism underlying the BOLD signal from irritative zones. Cross-validation indicates that these four mechanisms can be classified from realistic fMRI BOLD signals. To demonstrate proof of concept, we applied our methodology to EEG-fMRI data from five epileptic patients undergoing neurosurgery, suggesting the presence of some of these mechanisms. We concluded that a proper identification and interpretation of NBR mechanisms in epilepsy can be performed by combining general linear models and biophysically inspired models.

The cholinergic system modulates negative BOLD responses in the prefrontal cortex once electrical perforant pathway stimulation triggers neuronal afterdischarges in the hippocampus

Repeated high-frequency pulse-burst stimulations of the rat perforant pathway elicited positive BOLD responses in the right hippocampus, septum and prefrontal cortex. However, when the first stimulation period also triggered neuronal afterdischarges in the hippocampus, then a delayed negative BOLD response in the prefrontal cortex was generated. While neuronal activity and cerebral blood volume (CBV) increased in the hippocampus during the period of hippocampal neuronal afterdischarges (h-nAD), CBV decreased in the prefrontal cortex, although neuronal activity did not decrease. Only after termination of h-nAD did CBV in the prefrontal cortex increase again. Thus, h-nAD triggered neuronal activity in the prefrontal cortex that counteracted the usual neuronal activity-related functional hyperemia. This process was significantly enhanced by pilocarpine, a mACh receptor agonist, and completely blocked when pilocarpine was co-administered with scopolamine, a mACh receptor antagonist. Scopolamine did not prevent the formation of the negative BOLD response, thus mACh receptors modulate the strength of the negative BOLD response.

Ageing and the Ipsilateral M1 BOLD Response: A Connectivity Study

Young people exhibit a negative BOLD response in ipsilateral primary motor cortex (M1) when making unilateral movements, such as button presses. This negative BOLD response becomes more positive as people age. In this study, we investigated why this occurs, in terms of the underlying effective connectivity and haemodynamics. We applied dynamic causal modeling (DCM) to task fMRI data from 635 participants aged 18–88 from the Cam-CAN dataset, who performed a cued button pressing task with their right hand. We found that connectivity from contralateral supplementary motor area (SMA) and dorsal premotor cortex (PMd) to ipsilateral M1 became more positive with age, explaining 44% of the variability across people in ipsilateral M1 responses. In contrast, connectivity from contralateral M1 to ipsilateral M1 was weaker and did not correlate with individual differences in rM1 BOLD. Neurovascular and haemodynamic parameters in the model were not able to explain the age-related shift to positive BOLD. Our results add to a body of evidence implicating neural, rather than vascular factors as the predominant cause of negative BOLD—while emphasising the importance of inter-hemispheric connectivity. This study provides a foundation for investigating the clinical and lifestyle factors that determine the sign and amplitude of the M1 BOLD response in ageing, which could serve as a proxy for neural and vascular health, via the underlying neurovascular mechanisms.

FMRI and intra-cranial electrocorticography recordings in the same human subjects reveals negative BOLD signal coupled with silenced neuronal activity

Positive blood oxygenation level-dependent (BOLD) responses (PBR), as measured by functional Magnetic Resonance Imaging (fMRI), are the most utilized measurements to non-invasively map activity in the brain. Recent studies have consistently shown that BOLD responses are not exclusively positive. Negative BOLD responses (NBR) have been reported in response to specific sensory stimulations and tasks. However, the exact relationship between NBR and the underlying metabolic and neuronal demand is still under debate. In this study, we investigated the neurophysiological basis of negative BOLD using fMRI and intra-cranial electrophysiology (electrocorticography, ECoG) measurements from the same human participants. We show that, for those electrodes that responded to visual stimulation, PBR are correlated with high-frequency band (HFB) responses. Crucially, NBR were associated with an absence of HFB power responses and an unpredicted decrease in the alpha power responses.

Ageing and the ipsilateral M1 BOLD response: a connectivity study

Young people exhibit a negative BOLD response in ipsilateral primary motor cortex (M1) when making unilateral movements, such as button presses. This negative BOLD response becomes more positive as people age. Here we investigated why this occurs, in terms of the underlying effective connectivity and haemodynamics. We applied dynamic causal modelling (DCM) to task fMRI data from 635 participants aged 18-88 from the Cam-CAN dataset, who performed a cued button pressing task with their right hand. We found that connectivity from contralateral supplementary motor area (SMA) and dorsal premotor cortex (PMd) to ipsilateral M1 became more positive with age, explaining 44% of the variability across people in ipsilateral M1 responses. Neurovascular and haemodynamic parameters in the model were not able to explain the age-related shift to positive BOLD. Our results add to a body of evidence implicating neural, rather than vascular factors as the predominant cause of negative BOLD - while emphasising the importance of inter-hemispheric connectivity. This study provides a foundation for investigating the clinical and lifestyle factors that determine the sign and amplitude of the M1 BOLD response, which could serve as a proxy for neural and vascular health, via the underlying neurovascular mechanisms.

Dissociations between glucose metabolism and blood oxygenation in the human default mode network revealed by simultaneous PET-fMRI

The finding of reduced functional MRI (fMRI) activity in the default mode network (DMN) during externally focused cognitive control has been highly influential to our understanding of human brain function. However, these negative fMRI responses, measured as relative decreases in the blood-oxygenation-level–dependent (BOLD) response between rest and task, have also prompted major questions of interpretation. Using hybrid functional positron emission tomography (PET)-MRI, this study shows that task-positive and -negative BOLD responses do not reflect antagonistic patterns of synaptic metabolism. Task-positive BOLD responses in attention and control networks were accompanied by concomitant increases in glucose metabolism during cognitive control, but metabolism in widespread DMN remained high during rest and task despite negative BOLD responses. Dissociations between glucose metabolism and the BOLD response specific to the DMN reveal functional heterogeneity in this network and demonstrate that negative BOLD responses during cognitive control should not be interpreted to reflect relative increases in metabolic activity during rest. Rather, neurovascular coupling underlying BOLD response patterns during rest and task in DMN appears fundamentally different from BOLD responses in other association networks during cognitive control.

Do we really understand the role of the prefrontal cortex in placebo analgesia?

Several reviews have strongly implicated prefrontal cortical engagement in expectation-based placebo analgesia. We recently found a robust placebo analgesic response and associated decreases in pain-related cortical activations, without observable prefrontal engagement. We hypothesized our substantial conditioning and weak verbal instructions diminished expectation-related prefrontal activation. To test this, we examined the same subjects during a conditioning procedure, in which expectancy of pain relief was high. In two conditioning sessions, noxious heat was applied to a leg region treated with an ″analgesic″ cream and another treated with a ″moisturizing″ cream. In reality, both creams were inert, but the temperature applied to the moisturizing-cream area was 2°C higher than that applied to the analgesic-cream area. Functional MRI was acquired during the second conditioning session. Pain ratings were lower for the low heat than the high heat, with corresponding reduced activations in pain-related regions. Similar to previous studies with strong expectation for pain relief, we observed more prefrontal activations during the ″analgesic″ than the control condition. Nevertheless, contrary to the idea of active prefrontal engagement, the relative activation was based on differences in negative BOLD signals. A literature review revealed that only a few studies conclusively showed active engagement of prefrontal cortex, i.e. increased positive BOLD signal during high expectation compared to a control, with variable timing and spatial-specificity. We suggest that this variability is due to the heterogeneous influence of cognitive, emotional and motivational factors. Future studies should attempt to unravel the multiple contributions to placebo responsiveness in the prefrontal cortex.

Laminar-specific functional connectivity mapping with multi-slice line-scanning fMRI

ABSTRACTLaminar BOLD-fMRI has been applied to better depict the neuronal input and output circuitry and functional connectivity across cortical layers by measuring local hemodynamic changes. Despite extensive studies detecting laminar fMRI signals to illustrate the canonical microcircuit, the spatiotemporal characteristics of laminar-specific information flow across different cortical regions remain to be fully investigated in both evoked and resting states. Here, we developed a multi-slice line-scanning fMRI (MS-LS) method to detect laminar fMRI signals in adjacent cortical regions with high spatial (50 µm) and temporal resolution (100 ms) in anesthetized rats. Across different scanning trials, we detected both laminar-specific positive and negative BOLD responses in the surrounding cortical region adjacent to the most activated cortex under evoked condition. Specifically, in contrast to the typical Layer (L) 4 correlation across different regions due to the thalamocortical projections for trials with positive BOLD, a strong correlation pattern specific in L2/3 was detected for the trials with negative BOLD in adjacent regions, which indicate a brain state-dependent laminar-fMRI responses based on cortiocotical interaction from different trials. Also, we acquired the laminar-specific rs-fMRI signals across different cortical regions, of which the high spatiotemporal resolution allows us to estimate lag times based on the maximal cross-correlation of laminar-specific rs-fMRI signals. In contrast to the larger variability of lag times in L1 and 6, robust lag time differences in L2/3, 4, and 5 across multiple cortices represented the low-frequency rs-fMRI signal propagation from the caudal to the rostral slice. In summary, our work provides a unique laminar fMRI mapping scheme to better characterize trial-specific intra- and inter-laminar functional connectivity with MS-LS, presenting layer-specific spatiotemporal variation of both evoked and rs-fMRI signal.