Sanghyun Lim

Three-dimensional reconstruction of a cardiac outline by magnetocardiography.

A 3-D cardiac visualization is significantly helpful toward clinical applications of magnetocardiography (MCG), but the cardiac reconstruction requires a segmentation process using additional image modalities. This paper proposes a 3-D cardiac outline reconstruction method using only MCG measurement data without further imaging techniques. The cardiac outline was reconstructed by a combination of both spatial filtering and coherence mapping method. The strength of cardiac activities was first estimated by the array-gain constraint minimum-norm spatial filter with recursively updated gram matrix (AGMN-RUG). Then, waveforms were reconstructed at whole source grids, and the maximum source points of an atrium and ventricle were selected as a reference, respectively. Next, the coherence between each maximum source point and whole source points was compared by the coherence mapping method. A reconstructed cardiac outline was validated by comparing with an overlapped volume ratio when the reconstructed volume was identically matched with the original volume. The results obtained by the AGMN-RUG were compared to the results by other spatial filters. The accuracy of numerical simulation and phantom experiment by the AGMN-RUG was superior 10% and 8%, respectively, than the accuracy by the standardized low-resolution electromagnetic tomography. This accuracy demonstrated the efficacy of the proposed 3-D cardiac reconstruction method.

T 1 relaxation measurement of ex-vivo breast cancer tissues at ultralow magnetic fields.

We investigated T 1 relaxations of ex-vivo cancer tissues at low magnetic fields in order to check the possibility of achieving a T 1 contrast higher than those obtained at high fields. The T 1 relaxations of fifteen pairs (normal and cancerous) of breast tissue samples were measured at three magnetic fields, 37, 62, and 122 μT, using our superconducting quantum interference device-based ultralow field nuclear magnetic resonance setup, optimally developed for ex-vivo tissue studies. A signal reconstruction based on Bayesian statistics for noise reduction was exploited to overcome the low signal-to-noise ratio. The ductal and lobular-type tissues did not exhibit meaningful T 1 contrast values between normal and cancerous tissues at the three different fields. On the other hand, an enhanced T 1 contrast was obtained for the mucinous cancer tissue.

Evaluation of Phase-Amplitude Coupling in Resting State Magnetoencephalographic Signals: Effect of Surrogates and Evaluation Approach

Phase-amplitude coupling (PAC) plays an important role in neural communication and computation. Interestingly, recent studies have indicated the presence of ubiquitous PAC phenomenon even during the resting state. Despite the importance of PAC phenomenon, estimation of significant physiological PAC is challenging because of the lack of appropriate surrogate measures to control false positives caused by non-physiological PAC. Therefore, in the present study, we evaluated PAC phenomenon during resting-state magnetoencephalography (MEG) signal and considered various surrogate measures and computational approaches widely used in the literature in addition to proposing new ones. We evaluated PAC phenomenon over the entire length of the MEG signal and for multiple shorter time segments. The results indicate that the extent of PAC phenomenon mainly depends on the surrogate measures and PAC computational methods used, as well as the evaluation approach. After a careful and critical evaluation, we found that resting-state MEG signals failed to exhibit ubiquitous PAC phenomenon, contrary to what has been suggested previously.

Individual differences in oscillatory brain activity in response to varying attentional demands during a word recall and oculomotor dual task

Every day, we face situations that involve multi-tasking. How our brain utilizes cortical resources during multi-tasking is one of many interesting research topics. In this study, we tested whether a dual-task can be differentiated in the neural and behavioral responses of healthy subjects with varying degree of working memory capacity (WMC). We combined word recall and oculomotor tasks because they incorporate common neural networks including the fronto-parietal (FP) network. Three different types of oculomotor tasks (eye fixation, Fix-EM; predictive and random smooth pursuit eye movement, P-SPEM and R-SPEM) were combined with two memory load levels (low-load: five words, high-load: 10 words) for a word recall task. Each of those dual-task combinations was supposed to create varying cognitive loads on the FP network. We hypothesize that each dual-task requires different cognitive strategies for allocating the brain’s limited cortical resources and affects brain oscillation of the FP network. In addition, we hypothesized that groups with different WMC will show differential neural and behavioral responses. We measured oscillatory brain activity with simultaneous MEG and EEG recordings and behavioral performance by word recall. Prominent frontal midline (FM) theta (4–6 Hz) synchronization emerged in the EEG of the high-WMC group experiencing R-SPEM with high-load conditions during the early phase of the word maintenance period. Conversely, significant parietal upper alpha (10–12 Hz) desynchronization was observed in the EEG and MEG of the low-WMC group experiencing P-SPEM under high-load conditions during the same period. Different brain oscillatory patterns seem to depend on each individual’s WMC and varying attentional demands from different dual-task combinations. These findings suggest that specific brain oscillations may reflect different strategies for allocating cortical resources during combined word recall and oculomotor dual-tasks.

Dynamic pattern decoding of source-reconstructed MEG or EEG data: Perspective of multivariate pattern analysis and signal leakage

Recently, an increasing number of studies have employed multivariate pattern analysis (MVPA) rather than univariate analysis for the dynamic pattern decoding of event-related responses recorded with a MEG/EEG sensor. The use of the MVPA approach for source-reconstructed MEG/EEG data is uncommon. For these data, we need to consider the source orientation information and the signal leakage among brain regions. In the present study, we evaluate the perspective of the MVPA approach in the context of source orientation information and signal leakage in source-reconstructed MEG data. We perform face vs. tool object category decoding (FvsT-OCD) of event-related responses from single or multiple voxels from a brain region using a univariate analysis approach and/or the MVPA approach. We also propose and perform symmetric signal leakage correction of source-reconstructed data using an independent component analysis-based approach. FvsT-OCD using single voxel information shows higher sensitivity for the MVPA approach than univariate analysis, as the MVPA approach efficiently utilizes information on all three dipole orientations and is less affected by inter-subject variability. The MVPA approach shows higher sensitivity for FvsT-OCD when considering information from multiple voxels than for a single voxel in a brain region. This finding suggests that the MVPA approach captures the latent multivoxel distributed pattern. However, the results may be partly or entirely attributable to signal leakage between brain regions, as the sensitivity is substantially reduced after signal leakage correction. A consideration of signal leakage is therefore essential during the evaluation of MVPA outcomes.

Frontoparietal EEG alpha-phase synchrony reflects differential attentional demands during word recall and oculomotor dual-tasks.

To study the relationship between the varying degrees of cognitive load and long-range synchronization among neural networks, we utilized a dual-task paradigm combining concurrent word recall working memory tasks and oculomotor tasks that differentially activate the common frontoparietal (FP) network. We hypothesized that each dual-task combination would generate differential neuronal activation patterns among long-range connection during word retention period. Given that the FP alpha-phase synchronization is involved in attentional top-down processes, one would expect that the long-range synchronization pattern is affected by the degrees of dual-task demand. We measured a single-trial phase locking value in the alpha frequency (8–12 Hz) with electroencephalography in healthy participants. Single-trial phase locking value characterized the synchronization between two brain signals. Our results revealed that different amounts of FP alpha-phase synchronization were produced by different dual-task combinations, particularly during the early phase of the word retention period. These differences were dependent on the individual’s working memory capacity and memory load. Our study shows that during dual-task, each oculomotor task, which is subserved by distinct neural network, generates different modulation patterns on long-range neuronal activation and FP alpha-phase synchronization seems to reflect these differential cognitive loads.

Heartfelt Face Perception via the Interoceptive Pathway – a MEG study

The somatic marker hypothesis proposes that the cortical representation of visceral signals is a crucial component of emotion processing. No previous study has investigated the causal relationships among brain regions that process visceral information during emotional perception. In this magnetoencephalography study of 32 healthy subjects, heartbeat-evoked potentials (HEPs), which reflect the cortical processing of heartbeats, were modulated by the perception of a sad face. The modulation effect was localized to the prefrontal cortices, the globus pallidus, and an interoceptive network including the right anterior insula (RAI) and anterior cingulate cortex (RACC). Importantly, our Granger causality analysis provides the first evidence for increased causal flow of heartbeat information from the RAI to the RACC during sad face perception. Moreover, this HEP modulation effect was neither an artefact nor an effect of visual-evoked potentials. These findings provide important progress in the understanding of brain-body interactions during emotion processing.

Magnetoencephalographic study of event‐related fields and cortical oscillatory changes during cutaneous warmth processing

Thermoreception is an important cutaneous sense, which plays a role in the maintenance of our body temperature and in the detection of potential noxious heat stimulation. In this study, we investigated event‐related fields (ERFs) and neural oscillatory activities, which were modulated by warmth stimulation. We developed a warmth stimulator that could elicit a warmth sensation, without pain or tactile sensation, by using a deep‐penetrating 980‐nm diode laser. The index finger of each participant (n = 24) was irradiated with the laser warmth stimulus, and the cortical responses were measured using magnetoencephalography (MEG). The ERFs and oscillatory responses had late latencies (∼1.3 s and 1.0–1.5 s for ERFs and oscillatory responses, respectively), which could be explained by a slow conduction velocity of warmth‐specific C‐fibers. Cortical sources of warmth‐related ERFs were seen in the bilateral primary and secondary somatosensory cortices (SI and SII), posterior part of the anterior cingulate cortex (pACC), ipsilateral primary motor, and premotor cortex. Thus, we suggested that SI, SII, and pACC play a role in processing the warmth sensation. Time–frequency analysis demonstrated the suppression of the alpha (8–13 Hz) and beta (18–23 Hz) band power in the bilateral sensorimotor cortex. We proposed that the suppressions in alpha and beta band power are involved in the automatic response to the input of warmth stimulation and sensorimotor interactions. The delta band power (1–4 Hz) increased in the frontal, temporal, and cingulate cortices. The power changes in delta band might be related with the attentional processes during the warmth stimulation.

Development of a bio-magnetic measurement system and sensor configuration analysis for rats

Magnetoencephalography (MEG) based on superconducting quantum interference devices enables the measurement of very weak magnetic fields (10-1000 fT) generated from the human or animal brain. In this article, we introduce a small MEG system that we developed specifically for use with rats. Our system has the following characteristics: (1) variable distance between the pick-up coil and outer Dewar bottom (∼5 mm), (2) small pick-up coil (4 mm) for high spatial resolution, (3) good field sensitivity (45∼ 80fT/cm/Hz), (4) the sensor interval satisfies the Nyquist spatial sampling theorem, and (5) small source localization error for the region to be investigated. To reduce source localization error, it is necessary to establish an optimal sensor layout. To this end, we simulated confidence volumes at each point on a grid on the surface of a virtual rat head. In this simulation, we used locally fitted spheres as model rat heads. This enabled us to consider more realistic volume currents. We constrained the model such that the dipoles could have only four possible orientations: the x- and y-axes from the original coordinates, and two tangentially layered dipoles (local x- and y-axes) in the locally fitted spheres. We considered the confidence volumes according to the sensor layout and dipole orientation and positions. We then conducted a preliminary test with a 4-channel MEG system prior to manufacturing the multi-channel system. Using the 4-channel MEG system, we measured rat magnetocardiograms. We obtained well defined P-, QRS-, and T-waves in rats with a maximum value of 15 pT/cm. Finally, we measured auditory evoked fields and steady state auditory evoked fields with maximum values 400 fT/cm and 250 fT/cm, respectively.

Approximate Subject Specific Pseudo MRI from an Available MRI Dataset for MEG Source Imaging

Computation of headmodel and sourcemodel from the subjects MRI scan is an essential step for source localization of magnetoencephalography (MEG) (or EEG) sensor signals. In the absence of a real MRI scan, pseudo MRI (i.e., associated headmodel and sourcemodel) is often approximated from an available standard MRI template or pool of MRI scans considering the subjects digitized head surface. In the present study, we approximated two types of pseudo MRI (i.e., associated headmodel and sourcemodel) using an available pool of MRI scans with the focus on MEG source imaging. The first was the first rank pseudo MRI; that is, the MRI scan in the dataset having the lowest objective registration error (ORE) after being registered (rigid body transformation with isotropic scaling) to the subjects digitized head surface. The second was the averaged rank pseudo MRI that is generated by averaging of headmodels and sourcemodels from multiple MRI scans respectively, after being registered to the subjects digitized head surface. Subject level analysis showed that the mean upper bound of source location error for the approximated sourcemodel in reference to the real one was 10 ± 3 mm for the averaged rank pseudo MRI, which was significantly lower than the first rank pseudo MRI approach. Functional group source response in the brain to visual stimulation in the form of event-related power (ERP) at the time latency of peak amplitude showed noticeably identical source distribution for first rank pseudo MRI, averaged rank pseudo MRI, and real MRI. The source localization error for functional peak response was significantly lower for averaged rank pseudo MRI compared to first rank pseudo MRI. We conclude that it is feasible to use approximated pseudo MRI, particularly the averaged rank pseudo MRI, as a substitute for real MRI without losing the generality of the functional group source response.