Differential expression analysis for sequence count data

*Motivation:* High-throughput nucleotide sequencing provides quantitative readouts in assays for RNA expression (RNA-Seq), protein-DNA binding (ChIP-Seq) or cell counting (barcode sequencing). Statistical inference of differential signal in such data requires estimation of their variability throughout the dynamic range. When the number of replicates is small, error modelling is needed to achieve statistical power.

*Results:* We propose an error model that uses the negative binomial distribution, with variance and mean linked by local regression, to model the null distribution of the count data. The method controls type-I error and provides good detection power. 

*Availability:* A free open-source R software package, _DESeq_, is available from the Bioconductor project and from "http://www-huber.embl.de/users/anders/DESeq":http://www-huber.embl.de/users/anders/DESeq

Genes Suggest Ancestral Colour Polymorphisms Are Shared across Morphologically Cryptic Species in Arctic Bumblebees

email Suzanne orcd idCopyright: © 2015 Williams et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Understanding valuation of travel time changes: are preferences different under different stated choice design settings?

Stated choice (SC) experiments are the most popular method to estimate the value of travel time changes (VTTC) of a population. In the simplest VTTC experiment, the SC design variables are time changes and cost changes. The levels of these variables create a particular setting from which preferences are inferred. This paper tries to answer the question “do preferences vary with SC settings?”. For this, we investigate the role of the variables used in the SC experiment on the estimation of the set of VTTC (i.e. mean and covariates). Ideally, one would like to observe the same individuals completing different SC experiments. Since that option is not available, an alternative approach is to use a large dataset of responses, and split it according to different levels of the variable of interest. We refer to this as partial data analysis. The estimation of the same model on each sub-sample provides insights into potential effects of the variable of interest. This approach is applied in relation to three design variables on the data for the last national VTTC study in the UK, using state-of-the-art model specifications. The results show several ways in which the estimated set of VTTC can be affected by the levels of SC design variables. We conclude that model estimates (including the VTTC and covariates) are different in different settings. Hence by focussing the survey on specific settings, sample level results will be affected accordingly. Our findings have implications for appraisal and can inform the construction of future SC experiments

Fermi Gamma-ray Imaging of a Radio Galaxy

The Fermi Gamma-ray Space Telescope has detected the gamma-ray glow emanating from the giant radio lobes of the radio galaxy Centaurus A. The resolved gamma-ray image shows the lobes clearly separated from the central active source. In contrast to all other active galaxies detected so far in high-energy gamma-rays, the lobe flux constitutes a considerable portion (>1/2) of the total source emission. The gamma-ray emission from the lobes is interpreted as inverse Compton scattered relic radiation from the cosmic microwave background (CMB), with additional contribution at higher energies from the infrared-to-optical extragalactic background light (EBL). These measurements provide gamma-ray constraints on the magnetic field and particle energy content in radio galaxy lobes, and a promising method to probe the cosmic relic photon fields.Comment: 27 pages, includes Supplementary Online Material; corresponding authors: C.C. Cheung, Y. Fukazawa, J. Knodlseder, L. Stawar

Terrestrial implications for the maritime geoarchaeological resource: A view from the Lower Palaeolithic

Stone tools and faunal remains have been recovered from the English Channel and the North Sea through trawling, dredging for aggregates, channel clearance, and coring. These finds highlight the potential for a maritime Lower Palaeolithic archaeological resource. It is proposed here that any Lower Palaeolithic artefacts, faunal remains, and sediments deposited in the maritime zone during dry, low-stand phases were once (and may still be) contextually similar to their counterparts in the terrestrial Lower Palaeolithic records of north-western Europe. Given these similarities, can interpretive models and analytical frameworks developed for terrestrial archaeology be profitably applied to an assessment of the potential value of any maritime resource? The terrestrial geoarchaeological resource for the Lower Palaeolithic is dominated by artefacts and ecofacts that have been fluvially reworked. The spatio-temporal resolution of these data varies from entire river valleys and marine isotope stages to river channel gravel bar surfaces and decadal timescales, thus supporting a variety of questions and approaches. However, the structure of the terrestrial resource also highlights two fundamental limitations in current maritime knowledge that can restrict the application of terrestrial approaches to any potential maritime resource: (i) how have the repetitive transgressions and regressions of the Middle and Late Pleistocene modified the terrace landforms and sediments associated with the river systems of the English Channel and southern North Sea basins?; and (ii) do the surviving submerged terrace landforms and fluvial sedimentary deposits support robust geochronological models, as is the case with the classical terrestrial terrace sequences? This paper highlights potential approaches to these questions, and concludes that the fluvial palaeogeography, Pleistocene fossils, and potential Lower Palaeolithic artefacts of the maritime geoarchaeological resource can be profitably investigated in future as derived, low-resolution data sets, facilitating questions of colonisation, occupation, demography, and material culture

Detection of Gamma-Ray Emission from the Starburst Galaxies M82 and NGC 253 with the Large Area Telescope on Fermi

We report the detection of high-energy gamma-ray emission from two starburst galaxies using data obtained with the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. Steady point-like emission above 200 MeV has been detected at significance levels of 6.8 sigma and 4.8 sigma respectively, from sources positionally coincident with locations of the starburst galaxies M82 and NGC 253. The total fluxes of the sources are consistent with gamma-ray emission originating from the interaction of cosmic rays with local interstellar gas and radiation fields and constitute evidence for a link between massive star formation and gamma-ray emission in star-forming galaxies.Comment: Submitted to ApJ Letter

Extreme Food-Plant Specialisation in Megabombus Bumblebees as a Product of Long Tongues Combined with Short Nesting Seasons

© 2015 Huang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. http://creativecommons.org/licenses/by/4.0/ The attached file is the published version of the article

Evo-devo of human adolescence: beyond disease models of early puberty

Despite substantial heritability in pubertal development, much variation remains to be explained, leaving room for the influence of environmental factors to adjust its phenotypic trajectory in the service of fitness goals. Utilizing evolutionary development biology (evo-devo), we examine adolescence as an evolutionary life-history stage in its developmental context. We show that the transition from the preceding stage of juvenility entails adaptive plasticity in response to energy resources, other environmental cues, social needs of adolescence and maturation toward youth and adulthood. Using the evolutionary theory of socialization, we show that familial psychosocial stress fosters a fast life history and reproductive strategy rather than early maturation being just a risk factor for aggression and delinquency. Here we explore implications of an evolutionary-developmental-endocrinological-anthropological framework for theory building, while illuminating new directions for research

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal