Source Reconstruction Research Articles

MEG Microstates: An Investigation of Underlying Brain Sources and Potential Neurophysiological Processes

Microstates are transient scalp configurations of brain activity measured by electroencephalography (EEG). The application of microstate analysis in magnetoencephalography (MEG) data remains challenging. In one MEG dataset (N = 113), we aimed to identify MEG microstates at rest, explore their brain sources, and relate them to changes in brain activity during open-eyes (ROE) or closed-eyes resting state (RCE) and an auditory Mismatch Negativity (MMN) task. In another dataset of simultaneously recorded EEG-MEG data (N = 21), we investigated the association between MEG and EEG microstates. Six MEG microstates (mMS) provided the best clustering of resting-state activity, each linked to different brain sources: mMS 1–2: left/right occipito-parietal; mMS 3: fronto-temporal; mMS 4: centro-medial; mMS 5–6: left/right fronto-parietal. Increases in occipital alpha power in RCE relative to ROE correlated with greater mMS 1–2 time coverage (τbs < 0.20, ps > .002), while the lateralization of deviance detection in MMN was associated with mMS 5–6 time coverage (τbs < 0.16, ps > .012). No temporal correlation was found between EEG and MEG microstates (ps > .05), despite some overlap in brain sources and global explained variance between mMS 2–3 and EEG microstates B-C (rs > 0.60, ps < .002). Hence, the MEG signal can be decomposed into microstates, but mMS brain activity clustering captures phenomena different from EEG microstates. Source reconstruction and task-related modulations link mMS to large-scale networks and localized activities. Thus, mMSs offer insights into brain dynamics and task-specific processes, complementing EEG microstates in studying physiological and dysfunctional brain activity.

An Equivalent Radiation Source Reconstruction Method Based on Enhanced Artificial Neural Network

An Equivalent Radiation Source Reconstruction Method Based on Enhanced Artificial Neural Network

CutFEM-based MEG forward modeling improves source separability and sensitivity to quasi-radial sources: A somatosensory group study.

Source analysis of magnetoencephalography (MEG) data requires the computation of the magnetic fields induced by current sources in the brain. This so-called MEG forward problem includes an accurate estimation of the volume conduction effects in the human head. Here, we introduce the Cut finite element method (CutFEM) for the MEG forward problem. CutFEM's meshing process imposes fewer restrictions on tissue anatomy than tetrahedral meshes while being able to mesh curved geometries contrary to hexahedral meshing. To evaluate the new approach, we compare CutFEM with a boundary element method (BEM) that distinguishes three tissue compartments and a 6-compartment hexahedral FEM in an n = 19 group study of somatosensory evoked fields (SEF). The neural generators of the 20 ms post-stimulus SEF components (M20) are reconstructed using both an unregularized and a regularized inversion approach. Changing the forward model resulted in reconstruction differences of about 1 centimeter in location and considerable differences in orientation. The tested 6-compartment FEM approaches significantly increase the goodness of fit to the measured data compared with the 3-compartment BEM. They also demonstrate higher quasi-radial contributions for sources below the gyral crowns. Furthermore, CutFEM improves source separability compared with both other approaches. We conclude that head models with 6 compartments rather than 3 and the new CutFEM approach are valuable additions to MEG source reconstruction, in particular for sources that are predominantly radial.

EEG signature of near-death-like experiences during syncope-induced periods of unresponsiveness

During fainting, disconnected consciousness may emerge in the form of dream-like experiences. Characterized by extra-ordinary and mystical features, these subjective experiences have been associated to near-death-like experiences (NDEs-like). We here aim to assess brain activity during syncope-induced disconnected consciousness by means of high-density EEG monitoring. Transient loss of consciousness and unresponsiveness were induced in 27 healthy volunteers through hyperventilation, orthostasis, and Valsalva maneuvers. Upon awakening, subjects were asked to report memories, if any. The Greyson NDE scale was used to evaluate the potential phenomenological content experienced during the syncope-induced periods of unresponsiveness. EEG source reconstruction assessed cortical activations during fainting, which were regressed out with subjective reports collected upon recovery of normal consciousness. We also conducted functional connectivity, graph-theoretic and complexity analyses. High quality high-density EEG data were obtained in 22 volunteers during syncope and unresponsiveness (lasting 22±8 s). NDE-like features (Greyson NDE scale total score ≥7/32) were apparent for eight volunteers and characterized by higher activity in delta, theta and beta2 bands in temporal and frontal regions. The richness of the NDE-like content was associated with delta, theta and beta2 bands cortical current densities, in temporal, parietal and frontal lobes, including insula, right temporoparietal junction, and cingulate cortex. Our analyses also revealed a higher complexity and that networks related to delta, theta, and beta2 bands were characterized by a higher overall connectivity paralleled by a higher segregation (i.e., local efficiency) and a higher integration (i.e., global efficiency) for the NDE-like group compared to the non-NDE-like group. Fainting-induced NDE-like episodes seem to be sustained by surges of neural activity representing promising markers of disconnected consciousness.

Neural mechanisms for secondary suppression of emotional distractors: Evidence from concurrent electroencephalography-magnetoencephalography data.

Distractor suppression allows us to remain on-task in the presence of distractions by filtering task-irrelevant information from ongoing cognitive processing and responding. Electrophysiological studies have revealed that this key feature of selective attention is a dynamic process that involves at least two distinct stages of processing. Two important aspects of these processing stages remain unclear: Whether the processing of emotional distractors at an earlier stage is automatic, as reflected in the N2/early posterior negativity (EPN) component; and what functional-anatomical brain systems are recruited in each stage. The present study addresses these issues by measuring brain activity with concurrent electroencephalography-magnetoencephalography (MEG) recordings while participants performed a combined rapid serial visual presentation and motion tracking task. Event-related potentials (ERP) showed significant effects of attentional capture and attentional modulation during two time windows marked by the N2/EPN and P3b ERP components. Source reconstruction of concurrent MEG measurements revealed activation of the left visual association cortex and anterior cingulate cortex during the N2/EPN time window, activation of the insula during the early phase of the P3b and anterior cingulate cortex activity during the later phase of the P3b. The findings provide novel evidence establishing a connection between the increased N2 response to negative pictures and the activation of the cingulate gyrus, which facilitates the suppression of distractions during demanding cognitive tasks. In addition, distinct activation patterns were observed in the insula and anterior cingulate cortex during the P3b time window, indicating that attentional control mediated by the anterior cingulate cortex operates to suppress the processing of distracting emotional stimuli. (PsycInfo Database Record (c) 2024 APA, all rights reserved).

Noninvasive reconstruction of internal heat source in biological tissue using adaptive simulated annealing algorithm

The heat distribution information of human lesions is of great value for disease analysis, diagnosis, and treatment. It is a typical inverse problem of heat conduction that deriving the distribution of internal heat sources from the temperature distribution on the body surface. This paper transforms such an inverse problem of bio-heat transfer into a direct one, thereby avoiding complex boundary conditions and regularization processes. To noninvasively reconstruct the internal heat source and its corresponding 3D temperature field in biological tissue, the adaptive simulated annealing (ASA) algorithm is used in the simulation module, where the position P(x, y, z) of point heat source in biological tissue and its corresponding temperature T are set as the optimization variables. Under a certain optimized sample, one can obtain the simulated temperature distributing on the surface of the module, then subtract the simulated temperature from the measured temperature of the same surface which was measured using a thermal infrared imager. If the sum of absolute values of the difference is smaller, it indicates that the current sample is closer to the true location and temperature of the heat source. When the values of optimization variables are determined, the corresponding 3D temperature field is also confirmed. The simulation results show the simulated position and temperature of the heat source are very approximate with those of the real experimental module. The method presented in this paper has enormous potential and promising prospects in clinical research and application, such as tumor hyperthermia, disease thermal diagnosis technology, etc.

Cortical Generators and Connections Underlying Phoneme Perception: A Mismatch Negativity and P300 Investigation.

The cortical generators of the pure tone MMN and P300 have been thoroughly studied. Their nature and interaction with respect to phoneme perception, however, is poorly understood. Accordingly, the cortical sources and functional connections that underlie the MMN and P300 in relation to passive and active speech sound perception were identified.An inattentive and attentive phonemic oddball paradigm, eliciting a MMN and P300 respectively, were administered in 60 healthy adults during simultaneous high-density EEG recording. For both the MMN and P300, eLORETA source reconstruction was performed. The maximal cross-correlation was calculated between ROI-pairs to investigate inter-regional functional connectivity specific to passive and active deviant processing.MMN activation clusters were identified in the temporal (insula, superior temporal gyrus and temporal pole), frontal (rostral middle frontal and pars opercularis) and parietal (postcentral and supramarginal gyrus) cortex. Passive discrimination of deviant phonemes was aided by a network connecting right temporoparietal cortices to left frontal areas. For the P300, clusters with significantly higher activity were found in the frontal (caudal middle frontal and precentral), parietal (precuneus) and cingulate (posterior and isthmus) cortex. Significant intra- and interhemispheric connections between parietal, cingulate and occipital regions constituted the network governing active phonemic target detection. A predominantly bilateral network was found to underly both the MMN and P300.While passive phoneme discrimination is aided by a fronto-temporo-parietal network, active categorization calls on a network entailing fronto-parieto-cingulate cortices. Neural processing of phonemic contrasts, as reflected by the MMN and P300, does not appear to show pronounced lateralization to the language-dominant hemisphere.

Numerical Reconstruction of Heat Source on a Semi-Infinite Rod

Inverse heat conduction on an unbounded domain has always been the focus and difficulty in engineering and mechanical fields. In this paper, we consider the inverse problem of identifying the heat source by using the terminal measured temperature data. Being different from other ordinary heat conduction problems, our mathematical model is defined in a semi-infinite region, which may lead to many numerical methods that cannot be used directly. Based on the artificial boundary method, we first introduce the truncated boundary and derive the exact boundary conditions on it. Then we use the finite difference method to establish an approximate difference scheme for the forward problem. For the inverse problem, we design the Landweber and conjugate gradient iterative algorithm and do numerical experiments. Numerical results verify the correctness and effectiveness of our algorithm.

Far-field pattern accuracy improvement in planar near-field measurement using infinitesimal dipole modeling

This paper discusses the challenges in accurately measuring the far-field performance of array antennas, which have become more complex and larger due to recent trends in increasing the number of elements and the antenna size. The far-field measurement method, which is a general method for measuring a radiation pattern in a far-field, is difficult to be utilized at a high frequency due to limitations in physical chamber size and difficulty in diagnosing an antenna. In order to improve these problems, a near-field measurement method is used as an alternative. Even though using a near-field measurement, the accuracy of the derived far-field pattern is reduced at low elevation angles due to zero-padding of data outside the measurement area. The integration of the near-field measurement with the source reconstruction is a promising approach to enhance the measurement limits of the far-field pattern analysis. This paper aims to extend the previous work to accurately predict the far-field radiation pattern, especially in the low elevation angles. The verification of the proposed method is performed for a standard gain horn and 4 by 4 patch array antennas, with and without beam steering.

Unsupervised separation of nonlinearly mixed event-related potentials using manifold clustering and non-negative matrix factorization

Event-related potentials (ERPs) can quantify brain responses to reveal the neural mechanisms of sensory perception. However, ERPs often reflect nonlinear mixture responses to multiple sources of sensory stimuli, and an accurate separation of the response to each stimulus remains a challenge. This study aimed to separate the ERP into nonlinearly mixed source components specific to individual stimuli. We developed an unsupervised learning method based on clustering of manifold structures of mixture signals combined with channel optimization for signal source reconstruction using non-negative matrix factorization (NMF). Specifically, we first implemented manifold learning based on Local Tangent Space Alignment (LTSA) to extract the spatial manifold structure of multi-resolution sub-signals separated via wavelet packet transform. We then used fuzzy entropy to extract the dynamical process of the manifold structures and performed a k-means clustering to separate different sources. Lastly, we used NMF to obtain the optimal contributions of multiple channels to ensure accurate source reconstructions. We evaluated our developed approach using a simulated ERP dataset with known ground truth of two components of ERP mixture signals. Our results show that the correlation coefficient between the reconstructed source signal and the true source signal was 92.8 % and that the separation accuracy in ERP amplitude was 91.6 %. The results show that our unsupervised separation approach can accurately separate ERP signals from nonlinear mixture source components. The outcomes provide a promising way to isolate brain responses to multiple stimulus sources during multisensory perception.

Early effect of onabotulinumtoxinA on EEG-based functional connectivity in patients with chronic migraine: A pilot study.

In this pilot prospective cohort study, we aimed to evaluate, using high-density electroencephalography (HD-EEG), the longitudinal changes in functional connectivity (FC) in patients with chronic migraine (CM) treated with onabotulinumtoxinA (OBTA). OBTA is a treatment for CM. Several studies have shown the modulatory action of OBTA on the central nervous system; however, research on migraine is limited. This study was conducted at the Neurology Unit of "Policlinico Tor Vergata," Rome, Italy, and included 12 adult patients with CM treated with OBTA and 15 healthy controls (HC). Patients underwent clinical scales at enrollment (T0) and 3 months (T1) from the start of treatment. HD-EEG was recorded using a 64-channel system in patients with CM at T0 and T1. A source reconstruction method was used to identify brain activity. FC in δ-θ-α-β-low-γ bands was analyzed using the weighted phase-lag index. FC changes between HCs and CM at T0 and T1 were assessed using cross-validation methods to estimate the results' reliability. Compared to HCs at T0, patients with CM showed hyperconnected networks in δ (p = 0.046, area under the receiver operating characteristic curve [AUC: 0.76-0.98], Cohen's κ [0.65-0.93]) and β (p = 0.031, AUC [0.68-0.95], Cohen's κ [0.51-0.84]), mainly involving orbitofrontal, occipital, temporal pole and orbitofrontal, superior temporal, occipital, cingulate areas, and hypoconnected networks in α band (p = 0.029, AUC [0.80-0.99], Cohen's κ [0.42-0.77]), predominantly involving cingulate, temporal pole, and precuneus. Patients with CM at T1, compared to T0, showed hypoconnected networks in δ band (p = 0.032, AUC [0.73-0.99], Cohen's κ [0.53-0.90]) and hyperconnected networks in α band (p = 0.048, AUC [0.58-0.93], Cohen's κ [0.37-0.78]), involving the sensorimotor, orbitofrontal, cingulate, and temporal cortex. These preliminary results showed that patients with CM presented disrupted EEG-FC compared to controls restored by a single session of OBTA treatment, suggesting a primary central modulatory action of OBTA.

The effects of spatial leakage correction on the reliability of EEG-based functional connectivity networks.

Electroencephalography (EEG) functional connectivity (FC) estimates are confounded by the volume conduction problem. This effect can be greatly reduced by applying FC measures insensitive to instantaneous, zero-lag dependencies (corrected measures). However, numerous studies showed that FC measures sensitive to volume conduction (uncorrected measures) exhibit higher reliability and higher subject-level identifiability. We tested how source reconstruction contributed to the reliability difference of EEG FC measures on a large (n = 201) resting-state data set testing eight FC measures (including corrected and uncorrected measures). We showed that the high reliability of uncorrected FC measures in resting state partly stems from source reconstruction: idiosyncratic noise patterns define a baseline resting-state functional network that explains a significant portion of the reliability of uncorrected FC measures. This effect remained valid for template head model-based, as well as individual head model-based source reconstruction. Based on our findings we made suggestions how to best use spatial leakage corrected and uncorrected FC measures depending on the main goals of the study.

Phase space representation of sound field excited by a noise source in underwater acoustic waveguide.

The analysis of the field excited in a waveguide by a point noise source is performed using the phase space representation of this field given by the distribution of its amplitude in the depth-angle-time space. The transition from the traditional description of the field amplitude as a function of depth and time to phase space representation is performed using the coherent state expansion developed in quantum mechanics. In this paper, the correlation function of noise signals arriving at different points of the phase plane depth-angle is investigated. Numerical simulation data show that measurements of signal correlations in phase space, performed with the help of a receiving vertical antenna, can be used as input data in solving the problem of source localization and reconstruction of unknown parameters of the sound speed profile. It is shown that in phase space there is an analog of the classical interference pattern observed in the distribution of sound intensity in the distance-frequency plane. The slopes of striations in this interference pattern, as in the conventional one, are given by the Chuprov waveguide invariant.

PDE-constrained Optimization for Electroencephalographic Source Reconstruction

PDE-constrained Optimization for Electroencephalographic Source Reconstruction

Accuracy of dipole source reconstruction in the 3-layer BEM model against the 5-layer BEM-FMM model.

To compare cortical dipole fitting spatial accuracy between the widely used yet highly simplified 3-layer and modern more realistic 5-layer BEM-FMM models with and without adaptive mesh refinement (AMR) methods. We generate simulated noiseless 256-channel EEG data from 5-layer (7-compartment) meshes of 15 subjects from the Connectome Young Adult dataset. For each subject, we test four dipole positions, three sets of conductivity values, and two types of head segmentation. We use the boundary element method (BEM) with fast multipole method (FMM) acceleration, with or without (AMR), for forward modeling. Dipole fitting is carried out with the FieldTrip MATLAB toolbox. The average position error (across all tested dipoles, subjects, and models) is ~4 mm, with a standard deviation of ~2 mm. The orientation error is ~20° on average, with a standard deviation of ~15°. Without AMR, the numerical inaccuracies produce a larger disagreement between the 3- and 5-layer models, with an average position error of ~8 mm (6 mm standard deviation), and an orientation error of 28° (28° standard deviation). The low-resolution 3-layer models provide excellent accuracy in dipole localization. On the other hand, dipole orientation is retrieved less accurately. Therefore, certain applications may require more realistic models for practical source reconstruction. AMR is a critical component for improving the accuracy of forward EEG computations using a high-resolution 5-layer volume conduction model. Improving EEG source reconstruction accuracy is important for several clinical applications, including epilepsy and other seizure-inducing conditions.

EEG dynamic source imaging using a regularized optimization with spatio-temporal constraints.

One of the most important needs in neuroimaging is brain dynamic source imaging with high spatial and temporal resolution. EEG source imaging estimates the underlying sources from EEG recordings, which provides enhanced spatial resolution with intrinsically high temporal resolution. To ensure identifiability in the underdetermined source reconstruction problem, constraints on EEG sources are essential. This paper introduces a novel method for estimating source activities based on spatio-temporal constraints and a dynamic source imaging algorithm. The method enhances time resolution by incorporating temporal evolution of neural activity into a regularization function. Additionally, two spatial regularization constraints based on and norms are applied in the transformed domain to address both focal and spread neural activities, achieved through spatial gradient and Laplacian transform. Performance evaluation, conducted quantitatively using synthetic datasets, discusses the influence of parameters such as source extent, number of sources, correlation level, and SNR level on temporal and spatial metrics. Results demonstrate that the proposed method provides superior spatial and temporal reconstructions compared to state-of-the-art inverse solutions including STRAPS, sLORETA, SBL, dSPM, and MxNE. This improvement is attributed to the simultaneous integration of transformed spatial and temporal constraints. When applied to a real auditory ERP dataset, our algorithm accurately reconstructs brain source time series and locations, effectively identifying the origins of auditory evoked potentials. In conclusion, our proposed method with spatio-temporal constraints outperforms the state-of-the-art algorithms in estimating source distribution and time courses.

Spelling principles matter: An ERP study investigating the processing of different types of pseudohomophones

Spelling in any writing system is governed by fundamental principles. We examined the processing of two types of pseudohomophones constructed from words whose spellings are based on different principles – on the traditional principle of writing, requiring memorization of their spelling, and on the morphological principle, allowing the determination of their spelling from another word with the same morpheme (root) to examine the dependence of the occurrence of orthography–phonology conflict on spelling principles. Event-related potentials were recorded from 22 volunteers during silent reading. Pseudohomophones based on the morphological principle increased the N400 amplitude, emphasizing semantic and morphological processing importance. The P600 component showed significant effects for differentiating words and pseudohomophones based on the traditional principle, predominantly indicating the involvement of memory and reanalysis processes. Source reconstruction demonstrates that both pseudohomophones activate the left inferior frontal gyrus. However, pseudohomophones based on the traditional principle additionally activate the right and left postcentral gyrus, indicating the involvement of additional areas in the differentiation process. The earlier differences for stimuli based on the morphological principle indicate access to smaller units (morphemes), whereas stimuli based on the traditional principle require whole word processing. Our findings underscore the significant role of spelling principles in orthographic processing.

Mu-Suppression Neurofeedback Training Targeting the Mirror Neuron System: A Pilot Study

Neurofeedback training (NFT) is a promising adjuvant intervention method. The desynchronization of mu rhythm (8–13 Hz) in the electroencephalogram (EEG) over centro-parietal areas is known as a valid indicator of mirror neuron system (MNS) activation, which has been associated with social skills. Still, the effect of neurofeedback training on the MNS requires to be well investigated. The present study examined the possible impact of NFT with a mu suppression training protocol encompassing 15 NFT sessions (45 min each) on 16 healthy neurotypical participants. In separate pre- and post-training sessions, 64-channel EEG was recorded while participants (1) observed videos with various types of movements (including complex goal-directed hand movements and social interaction scenes) and (2) performed the "Reading the Mind in the Eyes Test" (RMET). EEG source reconstruction analysis revealed statistically significant mu suppression during hand movement observation across MNS-attributed fronto-parietal areas after NFT. The frequency analysis showed no significant mu suppression after NFT, despite the fact that numerical mu suppression appeared to be visible in a majority of participants during goal-directed hand movement observation. At the behavioral level, RMET accuracy scores did not suggest an effect of NFT on the ability to interpret subtle emotional expressions, although RMET response times were reduced after NFT. In conclusion, the present study exhibited preliminary and partial evidence that mu suppression NFT can induce mu suppression in MNS-attributed areas. More powerful experimental designs and longer training may be necessary to induce substantial and consistent mu suppression, particularly while observing social scenarios.

A general method to reconstruct strong gravitational lenses based on the singular perturbative approach

ABSTRACT The number of gravitational arcs systems detected is increasing quickly and should even increase at a faster rate in the near future. This wealth of new gravitational arcs requires the development of a purely automated method to reconstruct the lens and source. A general reconstruction method based on the singular perturbative approach is proposed in this paper. This method generates a lens and source reconstruction directly from the gravitational arc image. The method is fully automated and works in two steps. The first step is to generate a guess solution based on the circular solution in the singular perturbative approach. The second step is to break the sign degeneracy and to refine the solution by using a general source model. The refinement of the solution is conducted step-by-step to avoid the source-lens degeneracy issue. One important asset of this automated method is that the lens solution is written in universal terms which allows the computation of statistics. Considering the large number of lenses which should be available in the near future this ability to compute unbiased statistics is an important asset.

The auditory P2 is influenced by pitch changes but not pitch strength and consists of two separate subcomponents

Abstract The P2 component of the auditory evoked potential has previously been shown to depend on the acoustic stimulus properties and prior exposure to the materials. Here, we show that it is also affected by acoustic changes, as P2 amplitudes were strongly enhanced in response to voice pitch changes with a stepwise pattern compared to dynamic pitch changes typical for natural speech, and also reflected the magnitude of these pitch changes. Furthermore, it is demonstrated that neither the P2 nor any other component is affected by the harmonicity of the materials. Despite no prior exposure and a weaker pitch, artificially created inharmonic versions of the materials elicited similar activity throughout the auditory cortex. This suggests that so-called harmonic template neurons observed in animal studies are either absent or do not exist in sufficient number in the human auditory cortex to detect their activity extracranially. Crucially, morphology as well as scalp maps and source reconstructions of the EEG data showed that the P2 appears to consist of two separate subcomponents. While the “P2a” was localised to the auditory cortex, the subsequent “P2b” included generators spread across the auditory cortex and association areas. The two subcomponents thus likely reflect processing at different stages of the auditory pathway.