Evolution of genuine cross-correlation strength of focal onset seizures

To quantify the evolution of genuine zero-lag cross-correlations of focal onset seizures, we apply a recently introduced multivariate measure to broad band and to narrow-band EEG data. For frequency components below 12.5 Hz, the strength of genuine cross-correlations decreases significantly during the seizure and the immediate postseizure period, while higher frequency bands show a tendency of elevated cross-correlations during the same period. We conclude that in terms of genuine zero-lag cross-correlations, the electrical brain activity as assessed by scalp electrodes shows a significant spatial fragmentation, which might promote seizure offset.

Stochastic rank correlation: a robust merit function for 2D/3D registration of image data obtained at different energies

In this article, the authors evaluate a merit function for 2D/3D registration called stochastic rank correlation (SRC). SRC is characterized by the fact that differences in image intensity do not influence the registration result; it therefore combines the numerical advantages of cross correlation (CC)-type merit functions with the flexibility of mutual-information-type merit functions. The basic idea is that registration is achieved on a random subset of the image, which allows for an efficient computation of Spearman's rank correlation coefficient. This measure is, by nature, invariant to monotonic intensity transforms in the images under comparison, which renders it an ideal solution for intramodal images acquired at different energy levels as encountered in intrafractional kV imaging in image-guided radiotherapy. Initial evaluation was undertaken using a 2D/3D registration reference image dataset of a cadaver spine. Even with no radiometric calibration, SRC shows a significant improvement in robustness and stability compared to CC. Pattern intensity, another merit function that was evaluated for comparison, gave rather poor results due to its limited convergence range. The time required for SRC with 5% image content compares well to the other merit functions; increasing the image content does not significantly influence the algorithm accuracy. The authors conclude that SRC is a promising measure for 2D/3D registration in IGRT and image-guided therapy in general.

Analyzing spatio-temporal patterns of genuine cross-correlations

In multivariate time series analysis, the equal-time cross-correlation is a classic and computationally efficient measure for quantifying linear interrelations between data channels. When the cross-correlation coefficient is estimated using a finite amount of data points, its non-random part may be strongly contaminated by a sizable random contribution, such that no reliable conclusion can be drawn about genuine mutual interdependencies. The random correlations are determined by the signals' frequency content and the amount of data points used. Here, we introduce adjusted correlation matrices that can be employed to disentangle random from non-random contributions to each matrix element independently of the signal frequencies. Extending our previous work these matrices allow analyzing spatial patterns of genuine cross-correlation in multivariate data regardless of confounding influences. The performance is illustrated by example of model systems with known interdependence patterns. Finally, we apply the methods to electroencephalographic (EEG) data with epileptic seizure activity.

MicroRNA-21 suppression impedes medulloblastoma cell migration

Medulloblastoma (MB), the most common malignant brain tumour in children, is characterised by a high risk of leptomeningeal dissemination. But little is known about the molecular mechanisms that promote cancer cell migration in MB. Aberrant expression of miR-21 is recognised to be causatively linked to metastasis in a variety of human neoplasms including brain tumours; however its function in MB is still unknown. In this study we investigated the expression level and the role of miR-21 in MB cell migration. miR-21 was found to be up-regulated, compared to normal cerebellum, in 29/29 MB primary samples and 6/6 MB-derived cell lines. Inverse correlation was observed between miR-21 expression and the metastasis suppressor PDCD4, while miR-21 repression increased the release of PDCD4 protein, suggesting negative regulation of PDCD4 by miR-21 in MB cells. Anti-miR-21 decreased protein expression of the tumour cell invasion mediators MAP4K1 and JNK, which are also known to be negatively regulated by PDCD4, and down-regulated integrin protein that is essential for MB leptomeningeal dissemination. Moreover miR-21 knockdown in MB cells increased the expression of two eminent negative modulators of cancer cell migration, E-Cadherin and TIMP2 proteins that are known to be positively regulated by PDCD4. Finally and importantly, suppression of miR-21 decreased the motility of MB cells and reduced their migration across basement membranes in vitro. Together, these compelling data propose miR-21 pathway as a novel mechanism impacting MB cell dissemination and raises the possibility that curability of selected MB may be improved by pharmaceutical strategies directed towards microRNA-21.

Genuine cross-correlations: Which surrogate based measure reproduces analytical results best?

The analysis of short segments of noise-contaminated, multivariate real world data constitutes a challenge. In this paper we compare several techniques of analysis, which are supposed to correctly extract the amount of genuine cross-correlations from a multivariate data set. In order to test for the quality of their performance we derive time series from a linear test model, which allows the analytical derivation of genuine correlations. We compare the numerical estimates of the four measures with the analytical results for different correlation pattern. In the bivariate case all but one measure performs similarly well. However, in the multivariate case measures based on the eigenvalues of the equal-time cross-correlation matrix do not extract exclusively information about the amount of genuine correlations, but they rather reflect the spatial organization of the correlation pattern. This may lead to failures when interpreting the numerical results as illustrated by an application to three electroencephalographic recordings of three patients suffering from pharmacoresistent epilepsy.

Theoretical models of planetary system formation: mass vs. semi-major axis

Context. Planet formation models have been developed during the past years to try to reproduce what has been observed of both the solar system and the extrasolar planets. Some of these models have partially succeeded, but they focus on massive planets and, for the sake of simplicity, exclude planets belonging to planetary systems. However, more and more planets are now found in planetary systems. This tendency, which is a result of radial velocity, transit, and direct imaging surveys, seems to be even more pronounced for low-mass planets. These new observations require improving planet formation models, including new physics, and considering the formation of systems. Aims: In a recent series of papers, we have presented some improvements in the physics of our models, focussing in particular on the internal structure of forming planets, and on the computation of the excitation state of planetesimals and their resulting accretion rate. In this paper, we focus on the concurrent effect of the formation of more than one planet in the same protoplanetary disc and show the effect, in terms of architecture and composition of this multiplicity. Methods: We used an N-body calculation including collision detection to compute the orbital evolution of a planetary system. Moreover, we describe the effect of competition for accretion of gas and solids, as well as the effect of gravitational interactions between planets. Results: We show that the masses and semi-major axes of planets are modified by both the effect of competition and gravitational interactions. We also present the effect of the assumed number of forming planets in the same system (a free parameter of the model), as well as the effect of the inclination and eccentricity damping. We find that the fraction of ejected planets increases from nearly 0 to 8% as we change the number of embryos we seed the system with from 2 to 20 planetary embryos. Moreover, our calculations show that, when considering planets more massive than ~5 M⊕, simulations with 10 or 20 planetary embryos statistically give the same results in terms of mass function and period distribution.

Details of left ventricular radial wall motion supporting the ventricular theory of the third heart sound obtained by cardiac MR.

OBJECTIVE Obtaining new details of radial motion of left ventricular (LV) segments using velocity-encoding cardiac MRI. METHODS Cardiac MR examinations were performed on 14 healthy volunteers aged between 19 and 26 years. Cine images for navigator-gated phase contrast velocity mapping were acquired using a black blood segmented κ-space spoiled gradient echo sequence with a temporal resolution of 13.8 ms. Peak systolic and diastolic radial velocities as well as radial velocity curves were obtained for 16 ventricular segments. RESULTS Significant differences among peak radial velocities of basal and mid-ventricular segments have been recorded. Particular patterns of segmental radial velocity curves were also noted. An additional wave of outward radial movement during the phase of rapid ventricular filling, corresponding to the expected timing of the third heart sound, appeared of particular interest. CONCLUSION The technique has allowed visualization of new details of LV radial wall motion. In particular, higher peak systolic radial velocities of anterior and inferior segments are suggestive of a relatively higher dynamics of anteroposterior vs lateral radial motion in systole. Specific patterns of radial motion of other LV segments may provide additional insights into LV mechanics. ADVANCES IN KNOWLEDGE The outward radial movement of LV segments impacted by the blood flow during rapid ventricular filling provides a potential substrate for the third heart sound. A biphasic radial expansion of the basal anteroseptal segment in early diastole is likely to be related to the simultaneous longitudinal LV displacement by the stretched great vessels following repolarization and their close apposition to this segment.

Real-Time Multimodal Retinal Image Registration for a Computer Assisted Laser Photocoagulation System

An algorithm for the real-time registration of a retinal video sequence captured with a scanning digital ophthalmoscope (SDO) to a retinal composite image is presented. This method is designed for a computer-assisted retinal laser photocoagulation system to compensate for retinal motion and hence enhance the accuracy, speed, and patient safety of retinal laser treatments. The procedure combines intensity and feature-based registration techniques. For the registration of an individual frame, the translational frame-to-frame motion between preceding and current frame is detected by normalized cross correlation. Next, vessel points on the current video frame are identified and an initial transformation estimate is constructed from the calculated translation vector and the quadratic registration matrix of the previous frame. The vessel points are then iteratively matched to the segmented vessel centerline of the composite image to refine the initial transformation and register the video frame to the composite image. Criteria for image quality and algorithm convergence are introduced, which assess the exclusion of single frames from the registration process and enable a loss of tracking signal if necessary. The algorithm was successfully applied to ten different video sequences recorded from patients. It revealed an average accuracy of 2.47 ± 2.0 pixels (∼23.2 ± 18.8 μm) for 2764 evaluated video frames and demonstrated that it meets the clinical requirements.

Uniform approach to linear and nonlinear interrelation patterns in multivariate time series

Currently, a variety of linear and nonlinear measures is in use to investigate spatiotemporal interrelation patterns of multivariate time series. Whereas the former are by definition insensitive to nonlinear effects, the latter detect both nonlinear and linear interrelation. In the present contribution we employ a uniform surrogate-based approach, which is capable of disentangling interrelations that significantly exceed random effects and interrelations that significantly exceed linear correlation. The bivariate version of the proposed framework is explored using a simple model allowing for separate tuning of coupling and nonlinearity of interrelation. To demonstrate applicability of the approach to multivariate real-world time series we investigate resting state functional magnetic resonance imaging (rsfMRI) data of two healthy subjects as well as intracranial electroencephalograms (iEEG) of two epilepsy patients with focal onset seizures. The main findings are that for our rsfMRI data interrelations can be described by linear cross-correlation. Rejection of the null hypothesis of linear iEEG interrelation occurs predominantly for epileptogenic tissue as well as during epileptic seizures.

Observation of Associated Near-Side and Away-Side Long-Range Correlations in √sNN=5.02 TeV Proton-Lead Collisions with the ATLAS Detector

Two-particle correlations in relative azimuthal angle (Delta phi) and pseudorapidity (Delta eta) are measured in root S-NN = 5.02 TeV p + Pb collisions using the ATLAS detector at the LHC. The measurements are performed using approximately 1 mu b(-1) of data as a function of transverse momentum (p(T)) and the transverse energy (Sigma E-T(Pb)) summed over 3.1 < eta < 4.9 in the direction of the Pb beam. The correlation function, constructed from charged particles, exhibits a long-range (2 < vertical bar Delta eta vertical bar < 5) "near-side" (Delta phi similar to 0) correlation that grows rapidly with increasing Sigma E-T(Pb). A long-range "away-side" (Delta phi similar to pi) correlation, obtained by subtracting the expected contributions from recoiling dijets and other sources estimated using events with small Sigma E-T(Pb), is found to match the near-side correlation in magnitude, shape (in Delta eta and Delta phi) and Sigma E-T(Pb) dependence. The resultant Delta phi correlation is approximately symmetric about pi/2, and is consistent with a dominant cos2 Delta phi modulation for all Sigma E-T(Pb) ranges and particle p(T).

Transient rhythmic network activity in the somatosensory cortex evoked by distributed input in vitro

The initiation and maintenance of physiological and pathophysiological oscillatory activity depends on the synaptic interactions within neuronal networks. We studied the mechanisms underlying evoked transient network oscillation in acute slices of the adolescent rat somatosensory cortex and modeled its underpinning mechanisms. Oscillations were evoked by brief spatially distributed noisy extracellular stimulation, delivered via bipolar electrodes. Evoked transient network oscillation was detected with multi-neuron patch-clamp recordings under different pharmacological conditions. The observed oscillations are in the frequency range of 2-5 Hz and consist of 4-12 mV large, 40-150 ms wide compound synaptic events with rare overlying action potentials. This evoked transient network oscillation is only weakly expressed in the somatosensory cortex and requires increased [K+]o of 6.25 mM and decreased [Ca2+]o of 1.5 mM and [Mg2+]o of 0.5 mM. A peak in the cross-correlation among membrane potential in layers II/III, IV and V neurons reflects the underlying network-driven basis of the evoked transient network oscillation. The initiation of the evoked transient network oscillation is accompanied by an increased [K+]o and can be prevented by the K+ channel blocker quinidine. In addition, a shift of the chloride reversal potential takes place during stimulation, resulting in a depolarizing type A GABA (GABAA) receptor response. Blockade of alpha-amino-3-hydroxy-5-methyl-4-isoxazole-proprionate (AMPA), N-methyl-D-aspartate (NMDA), or GABA(A) receptors as well as gap junctions prevents evoked transient network oscillation while a reduction of AMPA or GABA(A) receptor desensitization increases its duration and amplitude. The apparent reversal potential of -27 mV of the evoked transient network oscillation, its pharmacological profile, as well as the modeling results suggest a mixed contribution of glutamatergic, excitatory GABAergic, and gap junctional conductances in initiation and maintenance of this oscillatory activity. With these properties, evoked transient network oscillation resembles epileptic afterdischarges more than any other form of physiological or pathophysiological neocortical oscillatory activity.

Theta burst transcranial magnetic stimulation is associated with increased EEG synchronization in the stimulated relative to unstimulated cerebral hemisphere

Theta burst transcranial magnetic stimulation (TBS) may induce behavioural changes that outlast the stimulation period. The neurophysiological basis of these behavioural changes are currently under investigation. Given the evidence that cortical information processing relies on transient synchronization and desynchronization of neuronal assemblies, we set out to test whether TBS is associated with changes of neuronal synchronization as assessed by surface EEG. In four healthy subjects one TBS train of 600 pulses (200 bursts, each burst consisting of 3 pulses at 30 Hz, repeated at intervals of 100 ms) was applied over the right frontal eye field and EEG synchronization was assessed in a time-resolved manner over 60 min by using a non-overlapping moving window. For each time step the linear cross-correlation matrix for six EEG channels of the right and for the six homotopic EEG channels of the left hemisphere were computed and their largest eigenvalues used to assess changes of synchronization. Synchronization was computed for broadband EEG and for the delta, theta, alpha, beta and gamma frequency bands. In all subjects EEG synchronization of the stimulated hemisphere was significantly and persistently increased relative to EEG synchronization of the unstimulated hemisphere. This effect occurred immediately after TBS for the theta, alpha, beta and gamma frequency bands and 10-20 min after TBS for broadband and delta frequency band EEG. Our results demonstrate that TBS is associated with increased neuronal synchronization of the cerebral hemisphere ipsilateral to the stimulation site relative to the unstimulated hemisphere. We speculate that enhanced synchronization interferes with cortical information processing and thus may be a neurophysiological correlate of the impaired behavioural performance detected previously.

Skeletal benefits after long-term retirement in former elite female gymnasts

Bone strength benefits after long-term retirement from elite gymnastics in terms of bone geometry and volumetric BMD were studied by comparing retired female gymnasts to moderately active age-matched women. In a cross-sectional study, 30 retired female gymnasts were compared with 30 age-matched moderately active controls. Bone geometric and densitometric parameters were measured by pQCT at the distal epiphyses and shafts of the tibia, femur, radius, and humerus. Muscle cross-sectional areas were assessed from the shaft scans. Independent t-tests were conducted on bone and muscle variables to detect differences between the two groups. The gymnasts had retired for a mean of 6.1 +/- 0.4 yr and were engaged in radial and humeral shafts, cortical cross-sectional area (CSA), total CSA, BMC, and strength strain index (SSIpol) were significantly greater (13-38%, p

Parameter-free extraction of the functional network topology from intracranial EEG recordings: a time-resolved study of graph properties in focal onset seizures

Rationale: Focal onset epileptic seizures are due to abnormal interactions between distributed brain areas. By estimating the cross-correlation matrix of multi-site intra-cerebral EEG recordings (iEEG), one can quantify these interactions. To assess the topology of the underlying functional network, the binary connectivity matrix has to be derived from the cross-correlation matrix by use of a threshold. Classically, a unique threshold is used that constrains the topology [1]. Our method aims to set the threshold in a data-driven way by separating genuine from random cross-correlation. We compare our approach to the fixed threshold method and study the dynamics of the functional topology. Methods: We investigate the iEEG of patients suffering from focal onset seizures who underwent evaluation for the possibility of surgery. The equal-time cross-correlation matrices are evaluated using a sliding time window. We then compare 3 approaches assessing the corresponding binary networks. For each time window: * Our parameter-free method derives from the cross-correlation strength matrix (CCS)[2]. It aims at disentangling genuine from random correlations (due to finite length and varying frequency content of the signals). In practice, a threshold is evaluated for each pair of channels independently, in a data-driven way. * The fixed mean degree (FMD) uses a unique threshold on the whole connectivity matrix so as to ensure a user defined mean degree. * The varying mean degree (VMD) uses the mean degree of the CCS network to set a unique threshold for the entire connectivity matrix. * Finally, the connectivity (c), connectedness (given by k, the number of disconnected sub-networks), mean global and local efficiencies (Eg, El, resp.) are computed from FMD, CCS, VMD, and their corresponding random and lattice networks. Results: Compared to FMD and VMD, CCS networks present: *topologies that are different in terms of c, k, Eg and El. *from the pre-ictal to the ictal and then post-ictal period, topological features time courses that are more stable within a period, and more contrasted from one period to the next. For CCS, pre-ictal connectivity is low, increases to a high level during the seizure, then decreases at offset. k shows a ‘‘U-curve’’ underlining the synchronization of all electrodes during the seizure. Eg and El time courses fluctuate between the corresponding random and lattice networks values in a reproducible manner. Conclusions: The definition of a data-driven threshold provides new insights into the topology of the epileptic functional networks.