Estimation of functional connectivity from electromagnetic signals and the amount of empirical data required

An increasing number of neuroimaging studies are concerned with the identification of interactions or statistical dependencies between brain areas. Dependencies between the activities of different brain regions can be quantified with functional connectivity measures such as the cross-correlation coefficient. An important factor limiting the accuracy of such measures is the amount of empirical data available. For event-related protocols, the amount of data also affects the temporal resolution of the analysis. We use analytical expressions to calculate the amount of empirical data needed to establish whether a certain level of dependency is significant when the time series are autocorrelated, as is the case for biological signals. These analytical results are then contrasted with estimates from simulations based on real data recorded with magnetoencephalography during a resting-state paradigm and during the presentation of visual stimuli. Results indicate that, for broadband signals, 50-100 s of data is required to detect a true underlying cross-correlations coefficient of 0.05. This corresponds to a resolution of a few hundred milliseconds for typical event-related recordings. The required time window increases for narrow band signals as frequency decreases. For instance, approximately 3 times as much data is necessary for signals in the alpha band. Important implications can be derived for the design and interpretation of experiments to characterize weak interactions, which are potentially important for brain processing.

Abnormal left and right amygdala-orbitofrontal cortical functional connectivity to emotional faces:state versus trait vulnerability markers of depression in bipolar disorder

Background - Amygdala-orbitofrontal cortical (OFC) functional connectivity (FC) to emotional stimuli and relationships with white matter remain little examined in bipolar disorder individuals (BD). Methods - Thirty-one BD (type I; n = 17 remitted; n = 14 depressed) and 24 age- and gender-ratio-matched healthy individuals (HC) viewed neutral, mild, and intense happy or sad emotional faces in two experiments. The FC was computed as linear and nonlinear dependence measures between amygdala and OFC time series. Effects of group, laterality, and emotion intensity upon amygdala-OFC FC and amygdala-OFC FC white matter fractional anisotropy (FA) relationships were examined. Results - The BD versus HC showed significantly greater right amygdala-OFC FC (p = .001) in the sad experiment and significantly reduced bilateral amygdala-OFC FC (p = .007) in the happy experiment. Depressed but not remitted female BD versus female HC showed significantly greater left amygdala-OFC FC (p = .001) to all faces in the sad experiment and reduced bilateral amygdala-OFC FC to intense happy faces (p = .01). There was a significant nonlinear relationship (p = .001) between left amygdala-OFC FC to sad faces and FA in HC. In BD, antidepressants were associated with significantly reduced left amygdala-OFC FC to mild sad faces (p = .001). Conclusions - In BD, abnormally elevated right amygdala-OFC FC to sad stimuli might represent a trait vulnerability for depression, whereas abnormally elevated left amygdala-OFC FC to sad stimuli and abnormally reduced amygdala-OFC FC to intense happy stimuli might represent a depression state marker. Abnormal FC measures might normalize with antidepressant medications in BD. Nonlinear amygdala-OFC FC–FA relationships in BD and HC require further study.

Differential patterns of activity and functional connectivity in emotion processing neural circuitry to angry and happy faces in adolescents with and without suicide attempt

Background - Neural substrates of emotion dysregulation in adolescent suicide attempters remain unexamined. Method - We used functional magnetic resonance imaging to measure neural activity to neutral, mild or intense (i.e. 0%, 50% or 100% intensity) emotion face morphs in two separate emotion-processing runs (angry and happy) in three adolescent groups: (1) history of suicide attempt and depression (ATT, n = 14); (2) history of depression alone (NAT, n = 15); and (3) healthy controls (HC, n = 15). Post-hoc analyses were conducted on interactions from 3 group × 3 condition (intensities) whole-brain analyses (p < 0.05, corrected) for each emotion run. Results - To 50% intensity angry faces, ATT showed significantly greater activity than NAT in anterior cingulate gyral–dorsolateral prefrontal cortical attentional control circuitry, primary sensory and temporal cortices; and significantly greater activity than HC in the primary sensory cortex, while NAT had significantly lower activity than HC in the anterior cingulate gyrus and ventromedial prefrontal cortex. To neutral faces during the angry emotion-processing run, ATT had significantly lower activity than NAT in the fusiform gyrus. ATT also showed significantly lower activity than HC to 100% intensity happy faces in the primary sensory cortex, and to neutral faces in the happy run in the anterior cingulate and left medial frontal gyri (all p < 0.006,corrected). Psychophysiological interaction analyses revealed significantly reduced anterior cingulate gyral–insula functional connectivity to 50% intensity angry faces in ATT v. NAT or HC. Conclusions - Elevated activity in attention control circuitry, and reduced anterior cingulate gyral–insula functional connectivity, to 50% intensity angry faces in ATT than other groups suggest that ATT may show inefficient recruitment of attentional control neural circuitry when regulating attention to mild intensity angry faces, which may represent a potential biological marker for suicide risk.

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.

Hypnotic induction is followed by state-like changes in the organization of EEG functional connectivity in the theta and beta frequency bands in high-hypnotically susceptible individuals

Altered state theories of hypnosis posit that a qualitatively distinct state of mental processing, which emerges in those with high hypnotic susceptibility following a hypnotic induction, enables the generation of anomalous experiences in response to specific hypnotic suggestions. If so then such a state should be observable as a discrete pattern of changes to functional connectivity (shared information) between brain regions following a hypnotic induction in high but not low hypnotically susceptible participants. Twenty-eight channel EEG was recorded from 12 high susceptible (highs) and 11 low susceptible (lows) participants with their eyes closed prior to and following a standard hypnotic induction. The EEG was used to provide a measure of functional connectivity using both coherence (COH) and the imaginary component of coherence (iCOH), which is insensitive to the effects of volume conduction. COH and iCOH were calculated between all electrode pairs for the frequency bands: delta (0.1-3.9 Hz), theta (4-7.9 Hz) alpha (8-12.9 Hz), beta1 (13-19.9 Hz), beta2 (20-29.9 Hz) and gamma (30-45 Hz). The results showed that there was an increase in theta iCOH from the pre-hypnosis to hypnosis condition in highs but not lows with a large proportion of significant links being focused on a central-parietal hub. There was also a decrease in beta1 iCOH from the pre-hypnosis to hypnosis condition with a focus on a fronto-central and an occipital hub that was greater in high compared to low susceptibles. There were no significant differences for COH or for spectral band amplitude in any frequency band. The results are interpreted as indicating that the hypnotic induction elicited a qualitative change in the organization of specific control systems within the brain for high as compared to low susceptible participants. This change in the functional organization of neural networks is a plausible indicator of the much theorized "hypnotic-state". © 2014 Jamieson and Burgess.

Functional connectivity of gamma EEG activity is modulated at low frequency during conscious recollection

We examined two subjectively distinct memory states that are elicited during recognition memory in humans and compared them in terms of the gamma oscillations (20–60 Hz) in the electroencepahalogram (EEG) that they induced. These subjective states, ‘recollection’ and ‘familiarity’ both entail correct recognition but one involves a clear and conscious recollection of the event including memory for contextual detail whilst the other involves a sense of familiarity without clear recollection. Here we show that during a verbal recognition memory test, the subjective experience of ‘recollection’ induced higher amplitude gamma oscillations than the subjective experience of ‘familiarity’ in the time period 300–500 ms after stimulus presentation. Recollection, but not familiarity, was also associated with greater functional connectivity in the gamma frequency range between frontal and parietal sites. Furthermore, the magnitude of the gamma functional connectivity varied over time and was modulated at 3 Hz. Previous studies in animals have shown local theta frequency modulation (3–7 Hz) of gamma-oscillations but this is the first time that a similar effect has been reported in the human EEG.

Altered resting-state EEG source functional connectivity in schizophrenia:the effect of illness duration

Despite the increasing body of evidence supporting the hypothesis of schizophrenia as a disconnection syndrome, studies of resting-state EEG Source Functional Connectivity (EEG-SFC) in people affected by schizophrenia are sparse. The aim of the present study was to investigate resting-state EEG-SFC in 77 stable, medicated patients with schizophrenia (SCZ) compared to 78 healthy volunteers (HV). In order to study the effect of illness duration, SCZ were divided in those with a short duration of disease (SDD; n = 25) and those with a long duration of disease (LDD; n = 52). Resting-state EEG recordings in eyes closed condition were analyzed and lagged phase synchronization (LPS) indices were calculated for each ROI pair in the source-space EEG data. In delta and theta bands, SCZ had greater EEG-SFC than HV; a higher theta band connectivity in frontal regions was observed in LDD compared with SDD. In the alpha band, SCZ showed lower frontal EEG-SFC compared with HV whereas no differences were found between LDD and SDD. In the beta1 band, SCZ had greater EEG-SFC compared with HVs and in the beta2 band, LDD presented lower frontal and parieto-temporal EEG-SFC compared with HV. In the gamma band, SDD had greater connectivity values compared with LDD and HV. This study suggests that resting state brain network connectivity is abnormally organized in schizophrenia, with different patterns for the different EEG frequency components and that EEG can be a powerful tool to further elucidate the complexity of such disordered connectivity.

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.

The lateral and ventromedial prefrontal cortex work as a dynamic integrated system:evidence from FMRI connectivity analysis

Recent functional magnetic resonance imaging (fMRI) investigations of the interaction between cognition and reward processing have found that the lateral prefrontal cortex (PFC) areas are preferentially activated to both increasing cognitive demand and reward level. Conversely, ventromedial PFC (VMPFC) areas show decreased activation to the same conditions, indicating a possible reciprocal relationship between cognitive and emotional processing regions. We report an fMRI study of a rewarded working memory task, in which we further explore how the relationship between reward and cognitive processing is mediated. We not only assess the integrity of reciprocal neural connections between the lateral PFC and VMPFC brain regions in different experimental contexts but also test whether additional cortical and subcortical regions influence this relationship. Psychophysiological interaction analyses were used as a measure of functional connectivity in order to characterize the influence of both cognitive and motivational variables on connectivity between the lateral PFC and the VMPFC. Psychophysiological interactions revealed negative functional connectivity between the lateral PFC and the VMPFC in the context of high memory load, and high memory load in tandem with a highly motivating context, but not in the context of reward alone. Physiophysiological interactions further indicated that the dorsal anterior cingulate and the caudate nucleus modulate this pathway. These findings provide evidence for a dynamic interplay between lateral PFC and VMPFC regions and are consistent with an emotional gating role for the VMPFC during cognitively demanding tasks. Our findings also support neuropsychological theories of mood disorders, which have long emphasized a dysfunctional relationship between emotion/motivational and cognitive processes in depression.

Functional interconnectivity between the globus pallidus and the subthalamic nucleus in the mouse brain slice

In accordance with its central role in basal ganglia circuitry, changes in the rate of action potential firing and pattern of activity in the globus pallidus (GP)-subthalamic nucleus (STN) network are apparent in movement disorders. In this study we have developed a mouse brain slice preparation that maintains the functional connectivity between the GP and STN in order to assess its role in shaping and modulating bursting activity promoted by pharmacological manipulations. Fibre-tract tracing studies indicated that a parasagittal slice cut 20 deg to the midline best preserved connectivity between the GP and the STN. IPSCs and EPSCs elicited by electrical stimulation confirmed connectivity from GP to STN in 44/59 slices and from STN to GP in 22/33 slices, respectively. In control slices, 74/76 (97%) of STN cells fired tonically at a rate of 10.3 ± 1.3 Hz. This rate and pattern of single spiking activity was unaffected by bath application of the GABAA antagonist picrotoxin (50 μM, n = 9) or the glutamate receptor antagonist (6-cyano-7-nitroquinoxaline-2, 3-dione (CNQX) 10 μM, n = 8). Bursting activity in STN neurones could be induced pharmacologically by application of NMDA alone (20 μM, 3/18 cells, 17%) but was more robust if NMDA was applied in conjunction with apamin (20-100 nM, 34/77 cells, 44%). Once again, neither picrotoxin (50 μM, n = 5) nor CNQX (10 μM, n = 5) had any effect on the frequency or pattern of the STN neurone activity while paired STN and GP recordings of tonic and bursting activity show no evidence of coherent activity. Thus, in a mouse brain slice preparation where functional GP-STN connectivity is preserved, no regenerative synaptically mediated activity indicative of a dynamic network is evident, either in the resting state or when neuronal bursting in both the GP and STN is generated by application of NMDA/apamin. This difference from the brain in Parkinson's disease may be attributed either to insufficient preservation of cortico-striato-pallidal or cortico-subthalamic circuitry, and/or an essential requirement for adaptive changes resulting from dopamine depletion for the expression of network activity within this tissue complex. © The Physiological Society 2005.

Disconnected brains:what is the role of fMRI in connectivity research?

In this paper we consider how functional Magnetic Resonance Imaging (fMRI) has been used to study cortical connectivity in autism and autistic spectrum disorders (ASD). We discuss those studies that have contributed to the evidence supporting a model of disordered cortical connectivity in autism and (ASD), with a focusing emphasis on the application to research into the underconnectivity model. We note that the analytical techniques employed are limited and do not allow interpretation in terms of effective, or directional connectivity, nor do they provide information about the temporal or spectral characteristics of the functional networks being studied. We highlight how currently the features of neural generators that are being assessed by functional connectivity in fMRI are unclear. In addition, we note the importance in clinical studies of considering the consequences of task choice for the nature of the imaging data that can be collected and also of individual differences in participant state and trait characteristics for the accurate interpretation of imaging data. We discuss how alternative techniques such as EEG/MEG may address the limitations of fMRI in assessing brain connectivity, and additionally consider the potential of multimodal approaches. We conclude that fMRI has made significant contributions towards our understanding of the brain in terms of neural systems but that the conclusions drawn from its application in the sphere of autism research need to be approached with caution. It is important in research of this kind that we are aware of the need to examine the methodological and analytical techniques closely when applying findings in clinical populations, not only when they are used to support the development of theoretical models but also to inform diagnostic or treatment decisions.

Brain connectivity in positive and negative syndrome schizophrenia

If, as is widely believed, schizophrenia is characterized by abnormalities of brain functional connectivity, then it seems reasonable to expect that different subtypes of schizophrenia could be discriminated in the same way. However, evidence for differences in functional connectivity between the subtypes of schizophrenia is largely lacking and, where it exists, it could be accounted for by clinical differences between the patients (e.g. medication) or by the limitations of the measures used. In this study, we measured EEG functional connectivity in unmedicated male patients diagnosed with either positive or negative syndrome schizophrenia and compared them with age and sex matched healthy controls. Using new methodology (Medkour et al., 2009) based on partial coherence, brain connectivity plots were constructed for positive and negative syndrome patients and controls. Reliable differences in the pattern of functional connectivity were found with both syndromes showing not only an absence of some of the connections that were seen in controls but also the presence of connections that the controls did not show. Comparing connectivity graphs using the Hamming distance, the negative-syndrome patients were found to be more distant from the controls than were the positive syndrome patients. Bootstrap distributions of these distances were created which showed a significant difference in the mean distances that was consistent with the observation that negative-syndrome diagnosis is associated with a more severe form of schizophrenia. We conclude that schizophrenia is characterized by widespread changes in functional connectivity with negative syndrome patients showing a more extreme pattern of abnormality than positive syndrome patients.

Causality and synchronisation in complex systems with applications to neuroscience

This thesis presents an investigation, of synchronisation and causality, motivated by problems in computational neuroscience. The thesis addresses both theoretical and practical signal processing issues regarding the estimation of interdependence from a set of multivariate data generated by a complex underlying dynamical system. This topic is driven by a series of problems in neuroscience, which represents the principal background motive behind the material in this work. The underlying system is the human brain and the generative process of the data is based on modern electromagnetic neuroimaging methods . In this thesis, the underlying functional of the brain mechanisms are derived from the recent mathematical formalism of dynamical systems in complex networks. This is justified principally on the grounds of the complex hierarchical and multiscale nature of the brain and it offers new methods of analysis to model its emergent phenomena. A fundamental approach to study the neural activity is to investigate the connectivity pattern developed by the brain’s complex network. Three types of connectivity are important to study: 1) anatomical connectivity refering to the physical links forming the topology of the brain network; 2) effective connectivity concerning with the way the neural elements communicate with each other using the brain’s anatomical structure, through phenomena of synchronisation and information transfer; 3) functional connectivity, presenting an epistemic concept which alludes to the interdependence between data measured from the brain network. The main contribution of this thesis is to present, apply and discuss novel algorithms of functional connectivities, which are designed to extract different specific aspects of interaction between the underlying generators of the data. Firstly, a univariate statistic is developed to allow for indirect assessment of synchronisation in the local network from a single time series. This approach is useful in inferring the coupling as in a local cortical area as observed by a single measurement electrode. Secondly, different existing methods of phase synchronisation are considered from the perspective of experimental data analysis and inference of coupling from observed data. These methods are designed to address the estimation of medium to long range connectivity and their differences are particularly relevant in the context of volume conduction, that is known to produce spurious detections of connectivity. Finally, an asymmetric temporal metric is introduced in order to detect the direction of the coupling between different regions of the brain. The method developed in this thesis is based on a machine learning extensions of the well known concept of Granger causality. The thesis discussion is developed alongside examples of synthetic and experimental real data. The synthetic data are simulations of complex dynamical systems with the intention to mimic the behaviour of simple cortical neural assemblies. They are helpful to test the techniques developed in this thesis. The real datasets are provided to illustrate the problem of brain connectivity in the case of important neurological disorders such as Epilepsy and Parkinson’s disease. The methods of functional connectivity in this thesis are applied to intracranial EEG recordings in order to extract features, which characterize underlying spatiotemporal dynamics before during and after an epileptic seizure and predict seizure location and onset prior to conventional electrographic signs. The methodology is also applied to a MEG dataset containing healthy, Parkinson’s and dementia subjects with the scope of distinguishing patterns of pathological from physiological connectivity.

Effect of dopamine on synchronous neuronal oscillations in the globus pallidus-subthalamic nucleus network

Changes in the pattern of activity of neurones within the basal ganglia are relevant in the pathophysiology and symptoms of Parkinson’s disease. The globus pallidus (GP) – subthalamic nucleus (STN) network has been proposed to form a pacemaker driving regenerative synchronous bursting activity. In order to test whether this activity can be sustained in vitro a 20o parasagittal slice of mouse midbrain was developed which preserved functional connectivity between the STN and GP. Mouse STN and GP cells were characterised electrophysiologically by the presence or absence of a voltage sag in response to hyperpolarising current steps indicative of Ih and the presence of rebound depolarisations. The presence of evoked and spontaneous post-synaptic GABA and glutamatergic currents indicated functional connectivity between the STN and GP. In control slices, STN cells fired action potentials at a regular rate, activity which was unaffected by bath application of the GABAA receptor antagonist picrotoxin (50 μM) or the glutamate receptor antagonist CNQX (10 μM). Paired extracellular recordings of STN cells showed uncorrelated firing. Oscillatory burst activity was induced pharmacologically using the glutamate receptor agonist, NMDA (20 μM), in combination with the potassium channel blocker apamin (50 -100 nM). The burst activity was unaffected by bath application of picrotoxin or CNQX while paired STN recordings showed uncorrelated activity indicating that the activity is not produced by the neuronal network. Thus, no regenerative activity is evident in this mouse brain preparation, either in control slices or when bursting is pharmacologically induced, suggesting the requirement of other afferent inputs that are not present in the slice. Using single-unit extracellular recording, dopamine (30 μM) produced an excitation of STN cells. This excitation was independent of synaptic transmission and was mimicked by both the Dl-like receptor agonist SKF38393 (10 μM) and the D2-like receptor agonist quinpirole (10 μM). However, the excitation was partially reduced by the D1-like antagonist SCH23390 (2 μM) but not by the D2-like antagonists sulpiride (10 μM) and eticlopride (10 μM). Using whole-recordings, dopamine was shown to induce membrane depolarisation. This depolarisation was caused either by a D1-like receptor mediated increase in a conductance which reversed at -34 mV, consistent with a non-specific cation conductance, or a D2-like receptor mediated decrease in conductance which reversed around -100 mV, consistent with a potassium conductance. Bath application of dopamine altered the pattern of the burst-firing produced by NMDA an apamin towards a more regular pattern. This effect was associated with a decrease in amplitude and ll1crease in frequency of TTX-resistant plateau potentials which underlie the burst activity.

Changes in neural complexity during the perception of 3D images using random dot stereograms

One developing theme in consciousness research is that consciousness is not the product of any specific component of the brain, rather it is an emergent property of the changing patterns of connectivity between different specialised functional components. For example, the dynamic core hypothesis proposes that conscious experience requires high levels of neural complexity, where complexity is defined in terms of functional connectivity. To test this hypothesis, electroencephalography was recorded while participants were shown random dot-stereograms. Consistent with the dynamic core hypothesis, neural complexity increased as the participants changed from simply viewing the stereogram to consciously perceiving the hidden 3D image.