Frequency Domain Connectivity Identification: An Application of Partial Directed Coherence in fMRI

Functional magnetic resonance imaging (fMRI) has become an important tool in Neuroscience due to its noninvasive and high spatial resolution properties compared to other methods like PET or EEG. Characterization of the neural connectivity has been the aim of several cognitive researches, as the interactions among cortical areas lie at the heart of many brain dysfunctions and mental disorders. Several methods like correlation analysis, structural equation modeling, and dynamic causal models have been proposed to quantify connectivity strength. An important concept related to connectivity modeling is Granger causality, which is one of the most popular definitions for the measure of directional dependence between time series. In this article, we propose the application of the partial directed coherence (PDC) for the connectivity analysis of multisubject fMRI data using multivariate bootstrap. PDC is a frequency domain counterpart of Granger causality and has become a very prominent tool in EEG studies. The achieved frequency decomposition of connectivity is useful in separating interactions from neural modules from those originating in scanner noise, breath, and heart beating. Real fMRI dataset of six subjects executing a language processing protocol was used for the analysis of connectivity. Hum Brain Mapp 30:452-461, 2009. (C) 2007 Wiley-Liss, Inc.

Increased coherence among striatal regions in the theta range during attentive wakefulness

The striatum, the largest component of the basal ganglia, is usually subdivided into associative, motor and limbic components. However, the electrophysiological interactions between these three subsystems during behavior remain largely unknown. We hypothesized that the striatum might be particularly active during exploratory behavior, which is presumably associated with increased attention. We investigated the modulation of local field potentials (LFPs) in the striatum during attentive wakefulness in freely moving rats. To this end, we implanted microelectrodes into different parts of the striatum of Wistar rats, as well as into the motor, associative and limbic cortices. We then used electromyograms to identify motor activity and analyzed the instantaneous frequency, power spectra and partial directed coherence during exploratory behavior. We observed fine modulation in the theta frequency range of striatal LFPs in 92.5 ± 2.5% of all epochs of exploratory behavior. Concomitantly, the theta power spectrum increased in all striatal channels (P < 0.001), and coherence analysis revealed strong connectivity (coefficients >0.7) between the primary motor cortex and the rostral part of the caudatoputamen nucleus, as well as among all striatal channels (P < 0.001). Conclusively, we observed a pattern of strong theta band activation in the entire striatum during attentive wakefulness, as well as a strong coherence between the motor cortex and the entire striatum. We suggest that this activation reflects the integration of motor, cognitive and limbic systems during attentive wakefulness.

Increased coherence among striatal regions in the theta range during attentive wakefulness

The striatum, the largest component of the basal ganglia, is usually subdivided into associative, motor and limbic components. However, the electrophysiological interactions between these three subsystems during behavior remain largely unknown. We hypothesized that the striatum might be particularly active during exploratory behavior, which is presumably associated with increased attention. We investigated the modulation of local field potentials (LFPs) in the striatum during attentive wakefulness in freely moving rats. To this end, we implanted microelectrodes into different parts of the striatum of Wistar rats, as well as into the motor, associative and limbic cortices. We then used electromyograms to identify motor activity and analyzed the instantaneous frequency, power spectra and partial directed coherence during exploratory behavior. We observed fine modulation in the theta frequency range of striatal LFPs in 92.5 ± 2.5% of all epochs of exploratory behavior. Concomitantly, the theta power spectrum increased in all striatal channels (P < 0.001), and coherence analysis revealed strong connectivity (coefficients >0.7) between the primary motor cortex and the rostral part of the caudatoputamen nucleus, as well as among all striatal channels (P < 0.001). Conclusively, we observed a pattern of strong theta band activation in the entire striatum during attentive wakefulness, as well as a strong coherence between the motor cortex and the entire striatum. We suggest that this activation reflects the integration of motor, cognitive and limbic systems during attentive wakefulness.

Dynamic directed interictal connectivity in left and right temporal lobe epilepsy.

OBJECTIVE There is increasing evidence that epileptic activity involves widespread brain networks rather than single sources and that these networks contribute to interictal brain dysfunction. We investigated the fast-varying behavior of epileptic networks during interictal spikes in right and left temporal lobe epilepsy (RTLE and LTLE) at a whole-brain scale using directed connectivity. METHODS In 16 patients, 8 with LTLE and 8 with RTLE, we estimated the electrical source activity in 82 cortical regions of interest (ROIs) using high-density electroencephalography (EEG), individual head models, and a distributed linear inverse solution. A multivariate, time-varying, and frequency-resolved Granger-causal modeling (weighted Partial Directed Coherence) was applied to the source signal of all ROIs. A nonparametric statistical test assessed differences between spike and baseline epochs. Connectivity results between RTLE and LTLE were compared between RTLE and LTLE and with neuropsychological impairments. RESULTS Ipsilateral anterior temporal structures were identified as key drivers for both groups, concordant with the epileptogenic zone estimated invasively. We observed an increase in outflow from the key driver already before the spike. There were also important temporal and extratemporal ipsilateral drivers in both conditions, and contralateral only in RTLE. A different network pattern between LTLE and RTLE was found: in RTLE there was a much more prominent ipsilateral to contralateral pattern than in LTLE. Half of the RTLE patients but none of the LTLE patients had neuropsychological deficits consistent with contralateral temporal lobe dysfunction, suggesting a relationship between connectivity changes and cognitive deficits. SIGNIFICANCE The different patterns of time-varying connectivity in LTLE and RTLE suggest that they are not symmetrical entities, in line with our neuropsychological results. The highest outflow region was concordant with invasive validation of the epileptogenic zone. This enhanced characterization of dynamic connectivity patterns could better explain cognitive deficits and help the management of epilepsy surgery candidates.

Anticipation of pain enhances the nociceptive transmission and functional connectivity within pain network in rats

Background: Expectation is a very potent pain modulator in both humans and animals. There is evidence that pain transmission neurons are modulated by expectation preceding painful stimuli. Nonetheless, few studies have examined the influence of pain expectation on the pain-related neuronal activity and the functional connectivity within the central nociceptive network. Results: This study used a tone-laser conditioning paradigm to establish the pain expectation in rats, and simultaneously recorded the anterior cingulate cortex (ACC), the medial dorsal thalamus (MD), and the primary somatosensory cortex (SI) to investigate the effect of pain expectation on laser-induced neuronal responses. Cross-correlation and partial directed coherence analysis were used to determine the functional interactions within and between the recorded areas during nociceptive transmission. The results showed that under anticipation condition, the neuronal activity to the auditory cue was significantly increased in the ACC area, whereas those to actual noxious stimuli were enhanced in all the recorded areas. Furthermore, neuronal correlations within and between these areas were significantly increased under conditions of expectation compared to those under non-expectation conditions, indicating an enhanced synchronization of neural activity within the pain network. In addition, information flow from the medial (ACC and MD) to the lateral (SI cortex) pain pathway increased, suggesting that the emotion-related neural circuits may modulate the neuronal activity in the somatosensory pathway during nociceptive transmission. Conclusion: These results demonstrate that the nociceptive processing in both medial and lateral pain systems is modulated by the expectation of pain.

Information theoretic interpretation of frequency domain connectivity measures

In order to provide adequate multivariate measures of information flow between neural structures, modified expressions of partial directed coherence (PDC) and directed transfer function (DTF), two popular multivariate connectivity measures employed in neuroscience, are introduced and their formal relationship to mutual information rates are proved.

Analyzing the connectivity between regions of interest: An approach based on cluster Granger causality for NRI data analysis

The identification, modeling, and analysis of interactions between nodes of neural systems in the human brain have become the aim of interest of many studies in neuroscience. The complex neural network structure and its correlations with brain functions have played a role in all areas of neuroscience, including the comprehension of cognitive and emotional processing. Indeed, understanding how information is stored, retrieved, processed, and transmitted is one of the ultimate challenges in brain research. In this context, in functional neuroimaging, connectivity analysis is a major tool for the exploration and characterization of the information flow between specialized brain regions. In most functional magnetic resonance imaging (fMRI) studies, connectivity analysis is carried out by first selecting regions of interest (ROI) and then calculating an average BOLD time series (across the voxels in each cluster). Some studies have shown that the average may not be a good choice and have suggested, as an alternative, the use of principal component analysis (PCA) to extract the principal eigen-time series from the ROI(s). In this paper, we introduce a novel approach called cluster Granger analysis (CGA) to study connectivity between ROIs. The main aim of this method was to employ multiple eigen-time series in each ROI to avoid temporal information loss during identification of Granger causality. Such information loss is inherent in averaging (e.g., to yield a single ""representative"" time series per ROI). This, in turn, may lead to a lack of power in detecting connections. The proposed approach is based on multivariate statistical analysis and integrates PCA and partial canonical correlation in a framework of Granger causality for clusters (sets) of time series. We also describe an algorithm for statistical significance testing based on bootstrapping. By using Monte Carlo simulations, we show that the proposed approach outperforms conventional Granger causality analysis (i.e., using representative time series extracted by signal averaging or first principal components estimation from ROIs). The usefulness of the CGA approach in real fMRI data is illustrated in an experiment using human faces expressing emotions. With this data set, the proposed approach suggested the presence of significantly more connections between the ROIs than were detected using a single representative time series in each ROI. (c) 2010 Elsevier Inc. All rights reserved.

Hyperdopaminergia and NMDA Receptor Hypofunction Disrupt Neural Phase Signaling

Neural phase signaling has gained attention as a putative coding mechanism through which the brain binds the activity of neurons across distributed brain areas to generate thoughts, percepts, and behaviors. Neural phase signaling has been shown to play a role in various cognitive processes, and it has been suggested that altered phase signaling may play a role in mediating the cognitive deficits observed across neuropsychiatric illness. Here, we investigated neural phase signaling in two mouse models of cognitive dysfunction: mice with genetically induced hyperdopaminergia [dopamine transporter knock-out (DAT-KO) mice] and mice with genetically induced NMDA receptor hypofunction [NMDA receptor subunit-1 knockdown (NR1-KD) mice]. Cognitive function in these mice was assessed using a radial-arm maze task, and local field potentials were recorded from dorsal hippocampus and prefrontal cortex as DAT-KO mice, NR1-KD mice, and their littermate controls engaged in behavioral exploration. Our results demonstrate that both DAT-KO and NR1-KD mice display deficits in spatial cognitive performance. Moreover, we show that persistent hyperdopaminergia alters interstructural phase signaling, whereas NMDA receptor hypofunction alters interstructural and intrastructural phase signaling. These results demonstrate that dopamine and NMDA receptor dependent glutamate signaling play a critical role in coordinating neural phase signaling, and encourage further studies to investigate the role that deficits in phase signaling play in mediating cognitive dysfunction.

Multimodal (EEG-fMRI) functional connectivity study of levodopa effect in Parkinson’s disease

Aim: To assess if the intake of levodopa in patients with Parkinson’s Disease (PD) changes cerebral connectivity, as revealed by simultaneous recording of hemodynamic (functional MRI, or fMRI) and electric (electroencephalogram, EEG) signals. Particularly, we hypothesize that the strongest changes in FC will involve the motor network, which is the most impaired in PD. Methods: Eight patients with diagnosis of PD “probable”, therapy with levodopa exclusively, normal cognitive and affective status, were included. Exclusion criteria were: moderate-severe rest tremor, levodopa induced dyskinesia, evidence of gray or white matter abnormalities on structural MRI. Scalp EEG (64 channels) were acquired inside the scanner (1.5 Tesla) before and after the intake of levodopa. fMRI functional connectivity was computed from four regions of interest: right and left supplementary motor area (SMA) and right and left precentral gyrus (primary motor cortex). Weighted partial directed coherence (w-PDC) was computed in the inverse space after the removal of EEG gradient and cardioballistic artifacts. Results and discussion: fMRI group analysis shows that the intake of levodopa increases hemodynamic functional connectivity among the SMAs / primary motor cortex and: sensory-motor network itself, attention network and default mode network. w-PDC analysis shows that EEG connectivity among regions of the motor network has the tendency to decrease after the intake the levodopa; furthermore, regions belonging to the DMN have the tendency to increase their outflow toward the rest of the brain. These findings, even if in a small sample of patients, suggest that other resting state physiological functional networks, beyond the motor one, are affected in patients with PD. The behavioral and cognitive tasks corresponding to the affected networks could benefit from the intake of levodopa.

Multiple frequency functional connectivity in the hand somatosensory network:an EEG study

Objective: To investigate the dynamics of communication within the primary somatosensory neuronal network. Methods: Multichannel EEG responses evoked by median nerve stimulation were recorded from six healthy participants. We investigated the directional connectivity of the evoked responses by assessing the Partial Directed Coherence (PDC) among five neuronal nodes (brainstem, thalamus and three in the primary sensorimotor cortex), which had been identified by using the Functional Source Separation (FSS) algorithm. We analyzed directional connectivity separately in the low (1-200. Hz, LF) and high (450-750. Hz, HF) frequency ranges. Results: LF forward connectivity showed peaks at 16, 20, 30 and 50. ms post-stimulus. An estimate of the strength of connectivity was modulated by feedback involving cortical and subcortical nodes. In HF, forward connectivity showed peaks at 20, 30 and 50. ms, with no apparent feedback-related strength changes. Conclusions: In this first non-invasive study in humans, we documented directional connectivity across subcortical and cortical somatosensory pathway, discriminating transmission properties within LF and HF ranges. Significance: The combined use of FSS and PDC in a simple protocol such as median nerve stimulation sheds light on how high and low frequency components of the somatosensory evoked response are functionally interrelated in sustaining somatosensory perception in healthy individuals. Thus, these components may potentially be explored as biomarkers of pathological conditions. © 2012 International Federation of Clinical Neurophysiology.

Graphical modelling for brain connectivity via partial coherence

Spectral and coherence methodologies are ubiquitous for the analysis of multiple time series. Partial coherence analysis may be used to try to determine graphical models for brain functional connectivity. The outcome of such an analysis may be considerably influenced by factors such as the degree of spectral smoothing, line and interference removal, matrix inversion stabilization and the suppression of effects caused by side-lobe leakage, the combination of results from different epochs and people, and multiple hypothesis testing. This paper examines each of these steps in turn and provides a possible path which produces relatively ‘clean’ connectivity plots. In particular we show how spectral matrix diagonal up-weighting can simultaneously stabilize spectral matrix inversion and reduce effects caused by side-lobe leakage, and use the stepdown multiple hypothesis test procedure to help formulate an interaction strength.

Effects Of Misalignments In The Retinal Nerve Fiber Layer Thickness Measurements With Spectral Domain Optical Coherence Tomography.

Purpose. To investigate misalignments (MAs) on retinal nerve fiber layer thickness (RNFLT) measurements obtained with Cirrus(©) SD-OCT. Methods. This was a retrospective, observational, cross-sectional study. Twenty-seven healthy and 29 glaucomatous eyes of 56 individuals with one normal exam and another showing MA were included. MAs were defined as an improper alignment of vertical vessels in the en face image. MAs were classified in complete MA (CMA) and partial MA (PMA), according to their site: 1 (superior, outside the measurement ring (MR)), 2 (superior, within MR), 3 (inferior, within MR), and 4 (inferior, outside MR). We compared RNFLT measurements of aligned versus misaligned exams in all 4 sectors, in the superior area (sectors 1 + 2), inferior area (sectors 3 + 4), and within the measurement ring (sectors 2 + 3). Results. RNFLT measurements at 12 clock-hour of eyes with MAs in the superior area (sectors 1 + 2) were significantly lower than those obtained in the same eyes without MAs (P = 0.043). No significant difference was found in other areas (sectors 1 + 2 + 3 + 4, sectors 3 + 4, and sectors 2 + 3). Conclusion. SD-OCT scans with superior MAs may present lower superior RNFLT measurements compared to aligned exams.

Prospective ultramorphological characterization of human hair by optical coherence tomography

Background/purpose The continuous advancement in cosmetic science has led to an increasing demand for the development of non-invasive, reliable scientific techniques directed toward claim substantiation, which is of utmost relevance, to obtain data regarding the efficacy and safety of cosmetic products. Methods In this work, we used the optical coherence tomography (OCT) technique to produce in vitro transversal section-images of human hair. We also compared the OCT signal before and after chemical treatment with an 18% w/w ammonium thioglycolate solution. Results The mean diameter of the medulla was 29 +/- 7 mu m and the hair diameter was 122 +/- 16 mu m in our samples of standard Afro-ethnic hair. A three-dimensional (3D) image was constructed starting from 601 cross-sectional images (slices). Each slice was taken in steps of 6.0 mu m at eight frames per second, and the entire 3D image was constructed in 60 s. Conclusion It was possible to identify, using the A-scan protocol, the principal structures: the cuticle, cortex and medulla. After chemical treatment, it was not possible to identify the main structures of hair fiber due to index matching promoted by deleterious action of the chemical agent.

Inhibition of glial cell proinflammatory activities by peroxisome proliferator-activated receptor gamma agonist confers partial protection during antimyelin oligodendrocyte glycoprotein demyelination in vitro.

Peroxisome proliferator-activated receptor gamma (PPAR-gamma) is a member of the nuclear hormone superfamily originally characterized as a regulator of adipocyte differentiation and lipid metabolism. In addition, PPAR-gamma has important immunomodulatory functions. If the effect of PPAR-gamma's activation in T-cell-mediated demyelination has been recently demonstrated, nothing is known about the role of PPAR-gamma in antibody-induced demyelination in the absence of T-cell interactions and monocyte/macrophage activation. Therefore, we investigated PPAR-gamma's involvement by using an in vitro model of inflammatory demyelination in three-dimensional aggregating rat brain cell cultures. We found that PPAR-gamma was not constitutively expressed in these cultures but was strongly up-regulated following demyelination mediated by antibodies directed against myelin oligodendrocyte glycoprotein (MOG) in the presence of complement. Pioglitazone, a selective PPAR-gamma agonist, partially protected aggregates from anti-MOG demyelination. Heat shock responses and the expression of the proinflammatory cytokine tumor necrosis factor-alpha were diminished by pioglitazone treatment. Therefore, pioglitazone protection seems to be linked to an inhibition of glial cell proinflammatory activities following anti-MOG induced demyelination. We show that PPAR-gamma agonists act not only on T cells but also on antibody-mediated demyelination. This may represent a significant benefit in treating multiple sclerosis patients.

Characterization of HbpR binding by site-directed mutagenesis of its DNA-binding site and by deletion of the effector domain.

In the presence of 2-hydroxybiphenyl, the enhancer binding protein, HbpR, activates the sigma54-dependent P(hbpC) promoter and controls the initial steps of 2-hydroxybiphenyl degradation in Pseudomonas azelaica. In the activation process, an oligomeric HbpR complex of unknown subunit composition binds to an operator region containing two imperfect palindromic sequences. Here, the HbpR-DNA binding interactions were investigated by site-directed mutagenesis of the operator region and by DNA-binding assays using purified HbpR. Mutations that disrupted the twofold symmetry in the palindromes did not affect the binding affinity of HbpR, but various mutations along a 60 bp region, and also outside the direct palindromic sequences, decreased the binding affinity. Footprints of HbpR on mutant operator fragments showed that a partial loss of binding contacts occurs, suggesting that the binding of one HbpR 'protomer' in the oligomeric complex is impaired whilst leaving the other contacts intact. An HbpR variant, devoid of its N-terminal sensing A-domain, was unable to activate transcription from the hbpC promoter while maintaining protection of the operator DNA in footprints. Wild-type HbpR was unable to activate transcription from the hbpC promoter when delta A-HbpR was expressed in the same cell, suggesting the formation of (repressing) hetero-oligomers. This model implies that HbpR can self-associate on its operator DNA without effector recognition or ATP binding. Furthermore, our findings suggest that the N-terminal sensing domain of HbpR is needed to activate the central ATPase domain rather than to repress a constitutively active C domain, as is the case for the related regulatory protein XylR.