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.

Doppler Orbit Determination of Deep Space Probes by the Bernese GNSS Software: First Results of the Combined Orbit Determination from DSN and Inter-Satellite Ka-Band Data from the Grail Mission

Navigation of deep space probes is most commonly operated using the spacecraft Doppler tracking technique. Orbital parameters are determined from a series of repeated measurements of the frequency shift of a microwave carrier over a given integration time. Currently, both ESA and NASA operate antennas at several sites around the world to ensure the tracking of deep space probes. Just a small number of software packages are nowadays used to process Doppler observations. The Astronomical Institute of the University of Bern (AIUB) has recently started the development of Doppler data processing capabilities within the Bernese GNSS Software. This software has been extensively used for Precise Orbit Determination of Earth orbiting satellites using GPS data collected by on-board receivers and for subsequent determination of the Earth gravity field. In this paper, we present the currently achieved status of the Doppler data modeling and orbit determination capabilities in the Bernese GNSS Software using GRAIL data. In particular we will focus on the implemented orbit determination procedure used for the combined analysis of Doppler and intersatellite Ka-band data. We show that even at this earlier stage of the development we can achieve an accuracy of few mHz on two-way S-band Doppler observation and of 2 µm/s on KBRR data from the GRAIL primary mission phase.

Resolving Atomic Connectivity in Graphene Nanostructure Junctions

We report on the structural characterization of junctions between atomically well-defined graphene nanoribbons (GNRs) by means of low-temperature, noncontact scanning probe microscopy. We show that the combination of simultaneously acquired frequency shift and tunneling current maps with tight binding (TB) simulations allows a comprehensive characterization of the atomic connectivity in the GNR junctions. The proposed approach can be generally applied to the investigation of graphene nanomaterials and their interconnections and is thus expected to become an important tool in the development of graphene-based circuitry.

Peri-ictal correlation dynamics of high-frequency (80-200 Hz) intracranial EEG

To assess (1) how large-scale correlation of intracranial EEG signals in the high-frequency range (80-200Hz) evolves from the pre-ictal, through the ictal into the postictal state and (2) whether the contribution of local neuronal activity to large-scale EEG correlation differentiates epileptogenic from non-epileptogenic brain tissue.

The impact of input fluctuations on the frequency-current relationships of layer 5 pyramidal neurons in the rat medial prefrontal cortex

The role of irregular cortical firing in neuronal computation is still debated, and it is unclear how signals carried by fluctuating synaptic potentials are decoded by downstream neurons. We examined in vitro frequency versus current (f-I) relationships of layer 5 (L5) pyramidal cells of the rat medial prefrontal cortex (mPFC) using fluctuating stimuli. Studies in the somatosensory cortex show that L5 neurons become insensitive to input fluctuations as input mean increases and that their f-I response becomes linear. In contrast, our results show that mPFC L5 pyramidal neurons retain an increased sensitivity to input fluctuations, whereas their sensitivity to the input mean diminishes to near zero. This implies that the discharge properties of L5 mPFC neurons are well suited to encode input fluctuations rather than input mean in their firing rates, with important consequences for information processing and stability of persistent activity at the network level.

Automatic Reduction of Artifacts in EEG-Signals

Electroencephalograms (EEG) are often contaminated with high amplitude artifacts limiting the usability of data. Methods that reduce these artifacts are often restricted to certain types of artifacts, require manual interaction or large training data sets. Within this paper we introduce a novel method, which is able to eliminate many different types of artifacts without manual intervention. The algorithm first decomposes the signal into different sub-band signals in order to isolate different types of artifacts into specific frequency bands. After signal decomposition with principal component analysis (PCA) an adaptive threshold is applied to eliminate components with high variance corresponding to the dominant artifact activity. Our results show that the algorithm is able to significantly reduce artifacts while preserving the EEG activity. Parameters for the algorithm do not have to be identified for every patient individually making the method a good candidate for preprocessing in automatic seizure detection and prediction algorithms.

A shift from oral to intravenous iron supplementation therapy is observed over time in a large swiss cohort of patients with inflammatory bowel disease

BACKGROUND In 2007, leading international experts in the field of inflammatory bowel disease (IBD) recommended intravenous (IV) iron supplements over oral (PO) ones because of superior effectiveness and better tolerance. We aimed to determine the percentage of patients with IBD undergoing iron therapy and to assess the dynamics of iron prescription habits (IV versus PO). METHODS We analyzed anonymized data on patients with Crohn's disease and ulcerative colitis extracted from the Helsana database. Helsana is a Swiss health insurance company providing coverage for 18% of the Swiss population (1.2 million individuals). RESULTS In total, 629 patients with Crohn's disease (61% female) and 398 patients with ulcerative colitis (57% female) were identified; mean observation time was 31.8 months for Crohn's disease and 31.0 months for ulcerative colitis patients. Of all patients with IBD, 27.1% were prescribed iron (21.1% in males; 31.1% in females). Patients treated with steroids, immunomodulators, and/or anti-tumor necrosis factor drugs were more frequently treated with iron supplements when compared with those not treated with any medications (35.0% versus 20.9%, odds ratio, 1.94; P < 0.001). The frequency of IV iron prescriptions increased significantly from 2006 to 2009 for both genders (males: from 2.6% to 10.1%, odds ratio = 3.84, P < 0.001; females: from 5.3% to 12.1%, odds ratio = 2.26, P = 0.002), whereas the percentage of PO iron prescriptions did not change. CONCLUSIONS Twenty-seven percent of patients with IBD were treated with iron supplements. Iron supplements administered IV were prescribed more frequently over time. These prescription habits are consistent with the implementation of guidelines on the management of iron deficiency in IBD.

Are Greenhouse Gas Signals of Northern Hemisphere winter extra-tropical cyclone activity dependent on the identification and tracking algorithm?

For Northern Hemisphere extra-tropical cyclone activity, the dependency of a potential anthropogenic climate change signal on the identification method applied is analysed. This study investigates the impact of the used algorithm on the changing signal, not the robustness of the climate change signal itself. Using one single transient AOGCM simulation as standard input for eleven state-of-the-art identification methods, the patterns of model simulated present day climatologies are found to be close to those computed from re-analysis, independent of the method applied. Although differences in the total number of cyclones identified exist, the climate change signals (IPCC SRES A1B) in the model run considered are largely similar between methods for all cyclones. Taking into account all tracks, decreasing numbers are found in the Mediterranean, the Arctic in the Barents and Greenland Seas, the mid-latitude Pacific and North America. Changing patterns are even more similar, if only the most severe systems are considered: the methods reveal a coherent statistically significant increase in frequency over the eastern North Atlantic and North Pacific. We found that the differences between the methods considered are largely due to the different role of weaker systems in the specific methods.

Vegetation shift and laurophyllisation: the possible role of forest fires

Southern Switzerland is a fire prone area where fire has to be considered as a natural environmental factor. In the past decades, fire frequency has tended to increase due to changes in landscape management. The most common type of fire is surface fire which normally breaks out during the vegetation resting period. Usually this type of fire shows short residence time (rapid spread), low to medium fire intensity and limited size. South-facing slopes are particularly fire-prone, so that very high fire frequency is possible: under these conditions passive resistant species and postfire resprouting species are favoured, usually leading to a reduction in the number of surviving species to a few fire adapted sprouters. Evergreen broadleaves are extremely sensitive to repeated fires. A simulation of the potential vegetation of southern Switzerland under climatic changed conditions evidenced the coincidence of the potential area of spreading forests rich in evergreen broad-leaved species with the most fire-prone area of the region. Therefore, in future, wildfires could play an important regulating role: most probably they will not stop the large-scale laurophyllisation of the thermophilous forests of southern Switzerland, but at sites with high fire frequency the vegetation shift could be slowed or even prevented by fire-disturbances.