LORETA With Cortical Constraint: Choosing an Adequate Surface Laplacian Operator

Low resolution electromagnetic tomography (LORETA) is a well-known method for the solution of the l2-based minimization problem for EEG/MEG source reconstruction. LORETA with a volume-based source space is widely used and much effort has been invested in the theory and the application of the method in an experimental context. However, it is especially interesting to use anatomical prior knowledge and constrain the LORETA's solution to the cortical surface. This strongly reduces the number of unknowns in the inverse approach. Unlike the Laplace operator in the volume case with a rectangular and regular grid, the mesh is triangulated and highly irregular in the surface case. Thus, it is not trivial to choose or construct a Laplace operator (termed Laplace-Beltrami operator when applied to surfaces) that has the desired properties and takes into account the geometry of the mesh. In this paper, the basic methodology behind cortical LORETA is discussed and the method is applied for source reconstruction of simulated data using different Laplace-Beltrami operators in the smoothing term. The results achieved with the different operators are compared with respect to their accuracy using various measures. Conclusions about the choice of an appropriate operator are deduced from the results

Quantitative Photo Activated Localization Microscopy: Unraveling the Effects of Photoblinking

In this work we discuss how to use photophysical information for improved quantitative measurements using Photo Activated Localization Microscopy (PALM) imaging. We introduce a method that reliably estimates the number of photoblinking molecules present in a biological sample and gives a robust way to quantify proteins at the single-cell level from PALM images. We apply this method to determine the amount of β2 adrenergic receptor, a prototypical G Protein Coupled Receptor, expressed on the plasma membrane of HeLa cells

Long-Range Autocorrelations of CpG Islands in the Human Genome

In this paper, we use a statistical estimator developed in astrophysics to study the distribution and organization of features of the human genome. Using the human reference sequence we quantify the global distribution of CpG islands (CGI) in each chromosome and demonstrate that the organization of the CGI across a chromosome is non-random, exhibits surprisingly long range correlations (10 Mb) and varies significantly among chromosomes. These correlations of CGI summarize functional properties of the genome that are not captured when considering variation in any particular separate (and local) feature. The demonstration of the proposed methods to quantify the organization of CGI in the human genome forms the basis of future studies. The most illuminating of these will assess the potential impact on phenotypic variation of inter-individual variation in the organization of the functional features of the genome within and among chromosomes, and among individuals for particular chromosomes

Quantitative motion analysis of subchromosomal foci in living cells using four-dimensional microscopy

The motion of subchromosomal foci and of whole chromosome territories in live human cell nuclei was investigated in four-dimensional space-time images. Visualization of subchromosomal foci was achieved by incorporating Cy3-dUTP into the nuclear DNA of two different cell types after microinjection. A subsequent segregation of the labeled cell nuclei led to the presence of only a few labeled chromosome territories on a background of nonlabeled chromatin (. Hum. Genet. 102:241-251). This procedure yielded many distinct signals in a given cell nucleus. Motion analysis in four-dimensional space-time images was performed using single-particle tracking and a statistical approach to the detection of a possible directional motion of foci relative to the center of mass of a chromosome territory. The accuracy of the analysis was tested using simulated data sets that closely mirrored the experimental setup and using microparticles of known size. Application of the analysis tools to experimental data showed that mutual diffusion-like movements between foci located on different chromosomes were more pronounced than inside the territories. In the time range observed, movements of individual foci could best be described by a random diffusion process. The statistical test for joint directed motion of several foci inside chromosome territories revealed that foci occasionally switched from random to directional motion inside the territories

The Discontinuous Galerkin Finite Element Method for Solving the MEG and the Combined MEG/EEG Forward Problem

In Electro- (EEG) and Magnetoencephalography (MEG), one important requirement of source reconstruction is the forward model. The continuous Galerkin finite element method (CG-FEM) has become one of the dominant approaches for solving the forward problem over the last decades. Recently, a discontinuous Galerkin FEM (DG-FEM) EEG forward approach has been proposed as an alternative to CG-FEM (Engwer et al., 2017). It was shown that DG-FEM preserves the property of conservation of charge and that it can, in certain situations such as the so-called skull leakages, be superior to the standard CG-FEM approach. In this paper, we developed, implemented, and evaluated two DG-FEM approaches for the MEG forward problem, namely a conservative and a non-conservative one. The subtraction approach was used as source model. The validation and evaluation work was done in statistical investigations in multi-layer homogeneous sphere models, where an analytic solution exists, and in a six-compartment realistically shaped head volume conductor model. In agreement with the theory, the conservative DG-FEM approach was found to be superior to the non-conservative DG-FEM implementation. This approach also showed convergence with increasing resolution of the hexahedral meshes. While in the EEG case, in presence of skull leakages, DG-FEM outperformed CG-FEM, in MEG, DG-FEM achieved similar numerical errors as the CG-FEM approach, i.e., skull leakages do not play a role for the MEG modality. In particular, for the finest mesh resolution of 1 mm sources with a distance of 1.59 mm from the brain-CSF surface, DG-FEM yielded mean topographical errors (relative difference measure, RDM%) of 1.5% and mean magnitude errors (MAG%) of 0.1% for the magnetic field. However, if the goal is a combined source analysis of EEG and MEG data, then it is highly desirable to employ the same forward model for both EEG and MEG data. Based on these results, we conclude that the newly presented conservative DG-FEM can at least complement and in some scenarios even outperform the established CG-FEM approaches in EEG or combined MEG/EEG source analysis scenarios, which motivates a further evaluation of DG-FEM for applications in bioelectromagnetism

The Discontinuous Galerkin Finite Element Method for Solving the MEG and the Combined MEG/EEG Forward Problem

In Electro- (EEG) and Magnetoencephalography (MEG), one important requirement of source reconstruction is the forward model. The continuous Galerkin finite element method (CG-FEM) has become one of the dominant approaches for solving the forward problem over the last decades. Recently, a discontinuous Galerkin FEM (DG-FEM) EEG forward approach has been proposed as an alternative to CG-FEM (Engwer et al., 2017). It was shown that DG-FEM preserves the property of conservation of charge and that it can, in certain situations such as the so-called skull leakages, be superior to the standard CG-FEM approach. In this paper, we developed, implemented, and evaluated two DG-FEM approaches for the MEG forward problem, namely a conservative and a non-conservative one. The subtraction approach was used as source model. The validation and evaluation work was done in statistical investigations in multi-layer homogeneous sphere models, where an analytic solution exists, and in a six-compartment realistically shaped head volume conductor model. In agreement with the theory, the conservative DG-FEM approach was found to be superior to the non-conservative DG-FEM implementation. This approach also showed convergence with increasing resolution of the hexahedral meshes. While in the EEG case, in presence of skull leakages, DG-FEM outperformed CG-FEM, in MEG, DG-FEM achieved similar numerical errors as the CG-FEM approach, i.e., skull leakages do not play a role for the MEG modality. In particular, for the finest mesh resolution of 1 mm sources with a distance of 1.59 mm from the brain-CSF surface, DG-FEM yielded mean topographical errors (relative difference measure, RDM%) of 1.5% and mean magnitude errors (MAG%) of 0.1% for the magnetic field. However, if the goal is a combined source analysis of EEG and MEG data, then it is highly desirable to employ the same forward model for both EEG and MEG data. Based on these results, we conclude that the newly presented conservative DG-FEM can at least complement and in some scenarios even outperform the established CG-FEM approaches in EEG or combined MEG/EEG source analysis scenarios, which motivates a further evaluation of DG-FEM for applications in bioelectromagnetism