Integration of EEG source imaging and fMRI during continuous viewing of natural movies

Electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) are noninvasive neuroimaging tools which can be used to measure brain activity with excellent temporal and spatial resolution, respectively. By combining the neural and hemodynamic recordings from these modalities, we can gain better insight into how and where the brain processes complex stimuli, which may be especially useful in patients with different neural diseases. However, due to their vastly different spatial and temporal resolutions, the integration of EEG and fMRI recordings is not always straightforward. One fundamental obstacle has been that paradigms used for EEG experiments usually rely on event-related paradigms, while fMRI is not limited in this regard. Therefore, here we ask whether one can reliably localize stimulus-driven EEG activity using the continuously varying feature intensities occurring in natural movie stimuli presented over relatively long periods of time. Specifically, we asked whether stimulus-driven aspects in the EEG signal would be co-localized with the corresponding stimulus-driven BOLD signal during free viewing of a movie. Secondly, we wanted to integrate the EEG signal directly with the BOLD signal, by estimating the underlying impulse response function (IRF) that relates the BOLD signal to the underlying current density in the primary visual area (V1). We made sequential fMRI and 64-channel EEG recordings in seven subjects who passively watched 2-min-long segments of a James Bond movie. To analyze EEG data in this natural setting, we developed a method based on independent component analysis (ICA) to reject EEG artifacts due to blinks, subject movement, etc., in a way unbiased by human judgment. We then calculated the EEG source strength of this artifact-free data at each time point of the movie within the entire brain volume using low-resolution electromagnetic tomography (LORETA). This provided for every voxel in the brain (i.e., in 3D space) an estimate of the current density at every time point. We then carried out a correlation between the time series of visual contrast changes in the movie with that of EEG voxels. We found the most significant correlations in visual area V1, just as seen in previous fMRI studies (Bartels A, Zeki, S, Logothetis NK. Natural vision reveals regional specialization to local motion and to contrast-invariant, global flow in the human brain. Cereb Cortex 2008;18(3):705–717), but on the time scale of milliseconds rather than of seconds. To obtain an estimate of how the EEG signal relates to the BOLD signal, we calculated the IRF between the BOLD signal and the estimated current density in area V1. We found that this IRF was very similar to that observed using combined intracortical recordings and fMRI experiments in nonhuman primates. Taken together, these findings open a new approach to noninvasive mapping of the brain. It allows, firstly, the localization of feature-selective brain areas during natural viewing conditions with the temporal resolution of EEG. Secondly, it provides a tool to assess EEG/BOLD transfer functions during processing of more natural stimuli. This is especially useful in combined EEG/fMRI experiments, where one can now potentially study neural-hemodynamic relationships across the whole brain volume in a noninvasive manner.     

Evaluation of Traumatic Spinal Cord Injury in a Rat Model Using 99mTc-GA-5 as a Potential In Vivo Tracer

Spinal cord injury (SCI) refers to the damage suffered in the spinal cord by any trauma or pathology. The purpose of this work was to determine whether ^99mTc-GA-5, a radiotracer targeting Glial Fibrillary Acidic Protein (GFAP), can reveal in vivo the reactivation of astrocytes in a murine model with SCI. A method for the ^99mTc radiolabeling of the mouse anti-GFAP monoclonal antibody GA-5 was implemented. Radiochemical characterization was performed, and radioimmunohistochemistry assays were used to evaluate the integrity of ^99mTc-GA-5. MicroSPECT/CT was used for in vivo imaging to trace SCI in the rats. No alterations in the GA-5’s recognition/specificity ability were observed after the radiolabeling. The GA-5’s radiolabeling procedure implemented in this work offers a practical method to allow the in vivo following of this monoclonal antibody to evaluate its biodistribution and specificity for GFAP receptors using SPECT/CT molecular imaging.     

Not So Innocent After All: Interfacial Chemistry Determines Charge-Transport Efficiency in Single-Molecule Junctions

Most studies in molecular electronics focus on altering the molecular wire backbone to tune the electrical properties of the whole junction. However, it is often overlooked that the chemical structure of the groups anchoring the molecule to the metallic electrodes influences the electronic structure of the whole system and, therefore, its conductance. We synthesised electron-accepting dithienophosphole oxide derivatives and fabricated their single-molecule junctions. We found that the anchor group has a dramatic effect on charge-transport efficiency: in our case, electron-deficient 4-pyridyl contacts suppress conductance, while electron-rich 4-thioanisole termini promote efficient transport. Our calculations show that this is due to minute changes in charge distribution, probed at the electrode interface. Our findings provide a framework for efficient molecular junction design, especially valuable for compounds with strong electron withdrawing/donating backbones.     

Bismuth-Catalyzed Amide Reduction

In this article we report that a cationic version of Akiba's Bi^III complex catalyzes the reduction of amides to amines using silane as hydride donor. The catalytic system features low catalyst loadings and mild conditions, en route to secondary and tertiary aryl- and alkylamines. The system tolerates functional groups such as alkene, ester, nitrile, furan and thiophene. Kinetic studies on the reaction mechanism result in the identification of a reaction network with an important product inhibition that is in agreement with the experimental reaction profiles.     

Numerical analysis for a system coupling curve evolution attached orthogonally to a fixed boundary,to a reaction–diffusion equation on the curve

We consider a semidiscrete finite element approximation for a system consisting of the evolution of a planar curve evolving by forced curve shortening flow inside a given bounded domain , such that the curve meets the boundary orthogonally, and the forcing is a function of the solution of a reaction–diffusion equation that holds on the evolving curve. We prove optimal order error bounds for the resulting approximation and present numerical experiments.     

A new ϕ-FEM approach for problems with natural boundary conditions

We present a new finite element method, called $\varphi$-FEM, to solve numerically elliptic partial differential equations with natural (Neumann or Robin) boundary conditions using simple computational grids, not fitted to the boundary of the physical domain. The boundary data are taken into account using a level-set function, which is a popular tool to deal with complicated or evolving domains. Our approach belongs to the family of fictitious domain methods (or immersed boundary methods) and is close to recent methods of CutFEM/XFEM type. Contrary to the latter, $\varphi$-FEM does not need any nonstandard numerical integration on cut mesh elements or on the actual boundary, while assuring the optimal convergence orders with finite elements of any degree and providing reasonably well conditioned discrete problems. In the first version of $\varphi$-FEM, only essential (Dirichlet) boundary conditions was considered. Here, to deal with natural boundary conditions, we introduce the gradient of the primary solution as an auxiliary variable. This is done only on the mesh cells cut by the boundary, so that the size of the numerical system is only slightly increased. We prove theoretically the optimal convergence of our scheme and a bound on the discrete problem conditioning, independent of the mesh cuts. The numerical experiments confirm these results.