Reverse Engineering Gene Networks with ANN: Variability in Network Inference Algorithms

Motivation :Reconstructing the topology of a gene regulatory network is one of the key tasks in systems biology. Despite of the wide variety of proposed methods, very little work has been dedicated to the assessment of their stability properties. Here we present a methodical comparison of the performance of a novel method (RegnANN) for gene network inference based on multilayer perceptrons with three reference algorithms (ARACNE, CLR, KELLER), focussing our analysis on the prediction variability induced by both the network intrinsic structure and the available data. Results: The extensive evaluation on both synthetic data and a selection of gene modules of "Escherichia coli" indicates that all the algorithms suffer of instability and variability issues with regards to the reconstruction of the topology of the network. This instability makes objectively very hard the task of establishing which method performs best. Nevertheless, RegnANN shows MCC scores that compare very favorably with all the other inference methods tested. Availability: The software for the RegnANN inference algorithm is distributed under GPL3 and it is available at the corresponding author home page (http://mpba.fbk.eu/grimaldi/regnann-supmat

Inference algorithms for gene networks: a statistical mechanics analysis

The inference of gene regulatory networks from high throughput gene expression data is one of the major challenges in systems biology. This paper aims at analysing and comparing two different algorithmic approaches. The first approach uses pairwise correlations between regulated and regulating genes; the second one uses message-passing techniques for inferring activating and inhibiting regulatory interactions. The performance of these two algorithms can be analysed theoretically on well-defined test sets, using tools from the statistical physics of disordered systems like the replica method. We find that the second algorithm outperforms the first one since it takes into account collective effects of multiple regulators

Direct constraint on the distance of y2 Velorum from AMBER/VLTI observations

In this work, we present the first AMBER observations, of the Wolf-Rayet and O (WR+O) star binary system y2 Velorum. The AMBER instrument was used with the telescopes UT2, UT3, and UT4 on baselines ranging from 46m to 85m. It delivered spectrally dispersed visibilities, as well as differential and closure phases, with a resolution R = 1500 in the spectral band 1.95-2.17 micron. We interpret these data in the context of a binary system with unresolved components, neglecting in a first approximation the wind-wind collision zone flux contribution. We show that the AMBER observables result primarily from the contribution of the individual components of the WR+O binary system. We discuss several interpretations of the residuals, and speculate on the detection of an additional continuum component, originating from the free-free emission associated with the wind-wind collision zone (WWCZ), and contributing at most to the observed K-band flux at the 5% level. The expected absolute separation and position angle at the time of observations were 5.1±0.9mas and 66±15° respectively. However, we infer a separation of 3.62+0.11-0.30 mas and a position angle of 73+9-11°. Our analysis thus implies that the binary system lies at a distance of 368+38-13 pc, in agreement with recent spectrophotometric estimates, but significantly larger than the Hipparcos value of 258+41-31 pc

Detecting a stochastic gravitational wave background with the Laser Interferometer Space Antenna

The random superposition of many weak sources will produce a stochastic background of gravitational waves that may dominate the response of the LISA (Laser Interferometer Space Antenna) gravitational wave observatory. Unless something can be done to distinguish between a stochastic background and detector noise, the two will combine to form an effective noise floor for the detector. Two methods have been proposed to solve this problem. The first is to cross-correlate the output of two independent interferometers. The second is an ingenious scheme for monitoring the instrument noise by operating LISA as a Sagnac interferometer. Here we derive the optimal orbital alignment for cross-correlating a pair of LISA detectors, and provide the first analytic derivation of the Sagnac sensitivity curve.Comment: 9 pages, 11 figures. Significant changes to the noise estimate

Increasing the Depth of Current Understanding: Sensitivity Testing of Deep-Sea Larval Dispersal Models for Ecologists

Larval dispersal is an important ecological process of great interest to conservation and the establishment of marine protected areas. Increasing numbers of studies are turning to biophysical models to simulate dispersal patterns, including in the deep-sea, but for many ecologists unassisted by a physical oceanographer, a model can present as a black box. Sensitivity testing offers a means to test the models' abilities and limitations and is a starting point for all modelling efforts. The aim of this study is to illustrate a sensitivity testing process for the unassisted ecologist, through a deep-sea case study example, and demonstrate how sensitivity testing can be used to determine optimal model settings, assess model adequacy, and inform ecological interpretation of model outputs. Five input parameters are tested (timestep of particle simulator (TS), horizontal (HS) and vertical separation (VS) of release points, release frequency (RF), and temporal range (TR) of simulations) using a commonly employed pairing of models. The procedures used are relevant to all marine larval dispersal models. It is shown how the results of these tests can inform the future set up and interpretation of ecological studies in this area. For example, an optimal arrangement of release locations spanning a release area could be deduced; the increased depth range spanned in deep-sea studies may necessitate the stratification of dispersal simulations with different numbers of release locations at different depths; no fewer than 52 releases per year should be used unless biologically informed; three years of simulations chosen based on climatic extremes may provide results with 90% similarity to five years of simulation; and this model setup is not appropriate for simulating rare dispersal events. A step-by-step process, summarising advice on the sensitivity testing procedure, is provided to inform all future unassisted ecologists looking to run a larval dispersal simulation

Lateral distribution of muons in IceCube cosmic ray events

In cosmic ray air showers, the muon lateral separation from the center of the shower is a measure of the transverse momentum that the muon parent acquired in the cosmic ray interaction. IceCube has observed cosmic ray interactions that produce muons laterally separated by up to 400 m from the shower core, a factor of 6 larger distance than previous measurements. These muons originate in high p(T) (>2 GeV/c) interactions from the incident cosmic ray, or high-energy secondary interactions. The separation distribution shows a transition to a power law at large values, indicating the presence of a hard p(T) component that can be described by perturbative quantum chromodynamics. However, the rates and the zenith angle distributions of these events are not well reproduced with the cosmic ray models tested here, even those that include charm interactions. This discrepancy may be explained by a larger fraction of kaons and charmed particles than is currently incorporated in the simulations.R. Abbasi ... G. C. Hill ... et al. (IceCube Collaboration

IceCube sensitivity for low-energy neutrinos from nearby supernovae

This paper describes the response of the IceCube neutrino telescope located at the geographic south pole to outbursts of MeV neutrinos from the core collapse of nearby massive stars. IceCube was completed in December 2010 forming a lattice of 5160 photomultiplier tubes that monitor a volume of ∼1 km3 in the deep Antarctic ice for particle induced photons. The telescope was designed to detect neutrinos with energies greater than 100 GeV. Owing to subfreezing ice temperatures, the photomultiplier dark noise rates are particularly low. Hence IceCube can also detect large numbers of MeV neutrinos by observing a collective rise in all photomultiplier rates on top of the dark noise. With 2 ms timing resolution, IceCube can detect subtle features in the temporal development of the supernova neutrino burst. For a supernova at the galactic center, its sensitivity matches that of a background-free megaton-scale supernova search experiment. The sensitivity decreases to 20 standard deviations at the galactic edge (30 kpc) and 6 standard deviations at the Large Magellanic Cloud (50 kpc). IceCube is sending triggers from potential supernovae to the Supernova Early Warning System. The sensitivity to neutrino properties such as the neutrino hierarchy is discussed, as well as the possibility to detect the neutronization burst, a short outbreak of νe ’s released by electron capture on protons soon after collapse. Tantalizing signatures, such as the formation of a quark star or a black hole as well as the characteristics of shock waves, are investigated to illustrate IceCube’s capability for supernova detection

Use of event-level neutrino telescope data in global fits for theories of new physics

We present a fast likelihood method for including event-level neutrino telescope data in parameter explorations of theories for new physics, and announce its public release as part of DarkSUSY 5.0.6. Our construction includes both angular and spectral information about neutrino events, as well as their total number. We also present a corresponding measure for simple model exclusion, which can be used for single models without reference to the rest of a parameter space. We perform a number of supersymmetric parameter scans with IceCube data to illustrate the utility of the method: example global fits and a signal recovery in the constrained minimal supersymmetric standard model (CMSSM), and a model exclusion exercise in a 7-parameter phenomenological version of the MSSM. The final IceCube detector con figuration will probe almost the entire focus-point region of the CMSSM, as well as a number of MSSM-7 models that will not otherwise be accessible to e. g. direct detection. Our method accurately recovers the mock signal, and provides tight constraints on model parameters and derived quantities. We show that the inclusion of spectral information significantly improves the accuracy of the recovery, providing motivation for its use in future IceCube analyses

Search for time-independent neutrino emission from astrophysical sources with 3 yr of IceCube data

We present the results of a search for neutrino point sources using the IceCube data collected between 2008 April and 2011 May with three partially completed configurations of the detector: the 40-, 59-, and 79-string configurations. The live-time of this data set is 1040 days. An unbinned maximum likelihood ratio test was used to search for an excess of neutrinos above the atmospheric background at any given direction in the sky. By adding two more years of data with improved event selection and reconstruction techniques, the sensitivity was improved by a factor of 3.5 or more with respect to the previously published results obtained with the 40-string configuration of IceCube. We performed an all-sky survey and a dedicated search using a catalog of a priori selected objects observed by other telescopes. In both searches, the data are compatible with the background-only hypothesis. In the absence of evidence for a signal, we set upper limits on the flux of muon neutrinos. For an E-2 neutrino spectrum, the observed limits are (0.9-5) x 10(-12) TeV-1 cm(-2) s(-1) for energies between 1 TeV and 1 PeV in the northern sky and (0.9-23.2) x 10(-12) TeV-1 cm(-2) s(-1) for energies between 10(2) TeV and 10(2) PeV in the southern sky. We also report upper limits for neutrino emission from groups of sources that were selected according to theoretical models or observational parameters and analyzed with a stacking approach. Some of the limits presented already reach the level necessary to quantitatively test current models of neutrino emission