The impact of MEG source reconstruction method on source-space connectivity estimation: A comparison between minimum-norm solution and beamforming.

Despite numerous important contributions, the investigation of brain connectivity with magnetoencephalography (MEG) still faces multiple challenges. One critical aspect of source-level connectivity, largely overlooked in the literature, is the putative effect of the choice of the inverse method on the subsequent cortico-cortical coupling analysis. We set out to investigate the impact of three inverse methods on source coherence detection using simulated MEG data. To this end, thousands of randomly located pairs of sources were created. Several parameters were manipulated, including inter- and intra-source correlation strength, source size and spatial configuration. The simulated pairs of sources were then used to generate sensor-level MEG measurements at varying signal-to-noise ratios (SNR). Next, the source level power and coherence maps were calculated using three methods (a) L2-Minimum-Norm Estimate (MNE), (b) Linearly Constrained Minimum Variance (LCMV) beamforming, and (c) Dynamic Imaging of Coherent Sources (DICS) beamforming. The performances of the methods were evaluated using Receiver Operating Characteristic (ROC) curves. The results indicate that beamformers perform better than MNE for coherence reconstructions if the interacting cortical sources consist of point-like sources. On the other hand, MNE provides better connectivity estimation than beamformers, if the interacting sources are simulated as extended cortical patches, where each patch consists of dipoles with identical time series (high intra-patch coherence). However, the performance of the beamformers for interacting patches improves substantially if each patch of active cortex is simulated with only partly coherent time series (partial intra-patch coherence). These results demonstrate that the choice of the inverse method impacts the results of MEG source-space coherence analysis, and that the optimal choice of the inverse solution depends on the spatial and synchronization profile of the interacting cortical sources. The insights revealed here can guide method selection and help improve data interpretation regarding MEG connectivity estimation

Mutual Information in Frequency and its Application to Measure Cross-Frequency Coupling in Epilepsy

We define a metric, mutual information in frequency (MI-in-frequency), to detect and quantify the statistical dependence between different frequency components in the data, referred to as cross-frequency coupling and apply it to electrophysiological recordings from the brain to infer cross-frequency coupling. The current metrics used to quantify the cross-frequency coupling in neuroscience cannot detect if two frequency components in non-Gaussian brain recordings are statistically independent or not. Our MI-in-frequency metric, based on Shannon's mutual information between the Cramer's representation of stochastic processes, overcomes this shortcoming and can detect statistical dependence in frequency between non-Gaussian signals. We then describe two data-driven estimators of MI-in-frequency: one based on kernel density estimation and the other based on the nearest neighbor algorithm and validate their performance on simulated data. We then use MI-in-frequency to estimate mutual information between two data streams that are dependent across time, without making any parametric model assumptions. Finally, we use the MI-in- frequency metric to investigate the cross-frequency coupling in seizure onset zone from electrocorticographic recordings during seizures. The inferred cross-frequency coupling characteristics are essential to optimize the spatial and spectral parameters of electrical stimulation based treatments of epilepsy.Comment: This paper is accepted for publication in IEEE Transactions on Signal Processing and contains 15 pages, 9 figures and 1 tabl

Neural population coding: combining insights from microscopic and mass signals

Behavior relies on the distributed and coordinated activity of neural populations. Population activity can be measured using multi-neuron recordings and neuroimaging. Neural recordings reveal how the heterogeneity, sparseness, timing, and correlation of population activity shape information processing in local networks, whereas neuroimaging shows how long-range coupling and brain states impact on local activity and perception. To obtain an integrated perspective on neural information processing we need to combine knowledge from both levels of investigation. We review recent progress of how neural recordings, neuroimaging, and computational approaches begin to elucidate how interactions between local neural population activity and large-scale dynamics shape the structure and coding capacity of local information representations, make them state-dependent, and control distributed populations that collectively shape behavior

Spatial Characterization and Estimation of Intracardiac Propagation Patterns During Atrial Fibrillation

This doctoral thesis is in the field of biomedical signal processing with focus on methods for the analysis of atrial fibrillation (AF). Paper I of the present thesis addresses the challenge of extracting spatial properties of AF from body surface signals. Different parameters are extracted to estimate the preferred direction of atrial activation and the complexity of the atrial activation pattern. In addition, the relation of the spatial properties to AF organization, which is quantified by AF frequency, is evaluated. While no significant correlation between the preferred direction of atrial activation and AF frequency could be observed, the complexity of the atrial activation pattern was found to increase with AF frequency. The remaining three papers deal with the analysis of the propagation of the electrical activity in the atria during AF based on intracardiac signals. In Paper II, a time-domain method to quantify propagation patterns along a linear catheter based on the detected atrial activation times is developed. Taking aspects on intra-atrial signal organization into account, the detected activation times are combined into wavefronts, and parameters related to the consistency of the wavefronts over time and the activation order along the catheter are extracted. Furthermore, the potential relationship of the extracted parameters to established measures from body surface signals is investigated. While the degree of wavefront consistency was not reflected by the applied body surface measures, AF frequency could distinguish between recordings with different degrees of intra-atrial signal organization. This supports the role of AF frequency as an organization measure of AF. In Paper III, a novel method to analyze intracardiac propagation patterns based on causality analysis in the frequency domain is introduced. In particular, the approach is based on the partial directed coherence (PDC), which evaluates directional coupling between multiple signals in the frequency domain. The potential of the method is illustrated with simulation scenarios based on a detailed ionic model of the human atrial cell as well as with real data recordings, selected to present typical propagation mechanisms and recording situations in atrial tachyarrhythmias. For simulated data, the PDC is correctly reflecting the direction of coupling and thus the propagation between all recording sites. For real data, clear propagation patterns are identified which agree with previous clinical observations. Thus, the results illustrate the ability of the novel approach to identify propagation patterns from intracardiac signals during AF which can provide important information about the underlying AF mechanisms, potentially improving the planning and outcome of ablation. However, spurious couplings over long distances can be observed when analyzing real data comprised by a large number of simultaneously recorded signals, which gives room for further improvement of the method. The derivation of the PDC is entirely based on the fit of a multivariate autoregressive (MVAR) model, commonly estimated by the least-squares (LS) method. In Paper IV, the adaptive group least absolute selection and shrinkage operator (LASSO) is introduced in order to avoid overfitting of the MVAR model and to incorporate prior information such as sparsity of the solution. The sparsity can be motivated by the observation that direct couplings over longer distances are likely to be zero during AF; an information which has been further incorporated by proposing distance-adaptive group LASSO. In simulations, adaptive and distance-adaptive group LASSO are found to be superior to LS estimation in terms of both detection and estimation accuracy. In addition, the results of both simulations and real data analysis indicate that further improvements can be achieved when the distance between the recording sites is known or can be estimated. This further promotes the PDC as a method for analysis of AF propagation patterns, which may contribute to a better understanding of AF mechanisms as well as improved AF treatment

Connectivity Analysis in EEG Data: A Tutorial Review of the State of the Art and Emerging Trends

Understanding how different areas of the human brain communicate with each other is a crucial issue in neuroscience. The concepts of structural, functional and effective connectivity have been widely exploited to describe the human connectome, consisting of brain networks, their structural connections and functional interactions. Despite high-spatial-resolution imaging techniques such as functional magnetic resonance imaging (fMRI) being widely used to map this complex network of multiple interactions, electroencephalographic (EEG) recordings claim high temporal resolution and are thus perfectly suitable to describe either spatially distributed and temporally dynamic patterns of neural activation and connectivity. In this work, we provide a technical account and a categorization of the most-used data-driven approaches to assess brain-functional connectivity, intended as the study of the statistical dependencies between the recorded EEG signals. Different pairwise and multivariate, as well as directed and non-directed connectivity metrics are discussed with a pros-cons approach, in the time, frequency, and information-theoretic domains. The establishment of conceptual and mathematical relationships between metrics from these three frameworks, and the discussion of novel methodological approaches, will allow the reader to go deep into the problem of inferring functional connectivity in complex networks. Furthermore, emerging trends for the description of extended forms of connectivity (e.g., high-order interactions) are also discussed, along with graph-theory tools exploring the topological properties of the network of connections provided by the proposed metrics. Applications to EEG data are reviewed. In addition, the importance of source localization, and the impacts of signal acquisition and pre-processing techniques (e.g., filtering, source localization, and artifact rejection) on the connectivity estimates are recognized and discussed. By going through this review, the reader could delve deeply into the entire process of EEG pre-processing and analysis for the study of brain functional connectivity and learning, thereby exploiting novel methodologies and approaches to the problem of inferring connectivity within complex networks

Electrophysiological characterization of the hyperdirect pathway and its functional relevance for subthalamic deep brain stimulation

The subthalamic nucleus (STN) receives input from various cortical areas via hyperdirect pathway (HDP) which bypasses the basal-ganglia loop. Recently, the HDP has gained increasing interest, because of its relevance for STN deep brain stimulation (DBS). To understand the HDP's role cortical responses evoked by STN-DBS have been investigated. These responses have short (<2 ms), medium (2–15 ms), and long (20–70 ms) latencies. Medium-latency responses are supposed to represent antidromic cortical activations via HDP. Together with long-latency responses the medium responses can potentially be used as biomarker of DBS efficacy as well as side effects. We here propose that the activation sequence of the cortical evoked responses can be conceptualized as high frequency oscillations (HFO) for signal analysis. HFO might therefore serve as marker for antidromic activation. Using existing knowledge on HFO recordings, this approach allows data analyses and physiological modeling to advance the pathophysiological understanding of cortical DBS-evoked high-frequency activity