Empirical Bayes models for multiple probe type microarrays at the probe level

<p>Abstract</p> <p>Background</p> <p>When analyzing microarray data a primary objective is often to find differentially expressed genes. With empirical Bayes and penalized t-tests the sample variances are adjusted towards a global estimate, producing more stable results compared to ordinary t-tests. However, for Affymetrix type data a clear dependency between variability and intensity-level generally exists, even for logged intensities, most clearly for data at the probe level but also for probe-set summarizes such as the MAS5 expression index. As a consequence, adjustment towards a global estimate results in an intensity-level dependent false positive rate.</p> <p>Results</p> <p>We propose two new methods for finding differentially expressed genes, Probe level Locally moderated Weighted median-t (PLW) and Locally Moderated Weighted-t (LMW). Both methods use an empirical Bayes model taking the dependency between variability and intensity-level into account. A global covariance matrix is also used allowing for differing variances between arrays as well as array-to-array correlations. PLW is specially designed for Affymetrix type arrays (or other multiple-probe arrays). Instead of making inference on probe-set summaries, comparisons are made separately for each perfect-match probe and are then summarized into one score for the probe-set.</p> <p>Conclusion</p> <p>The proposed methods are compared to 14 existing methods using five spike-in data sets. For RMA and GCRMA processed data, PLW has the most accurate ranking of regulated genes in four out of the five data sets, and LMW consistently performs better than all examined moderated t-tests when used on RMA, GCRMA, and MAS5 expression indexes.</p

Metarhizium brunneum Blastospore Pathogenesis in Aedes aegypti Larvae: Attack on Several Fronts Accelerates Mortality

Aedes aegypti is the vector of a wide range of diseases (e.g. yellow fever, dengue, Chikungunya and Zika) which impact on over half the world's population. Entomopathogenic fungi such as Metarhizium anisopliae and Beauveria bassiana have been found to be highly efficacious in killing mosquito larvae but only now are the underlying mechanisms for pathogenesis being elucidated. Recently it was shown that conidia of M. anisopliae caused stress induced mortality in Ae. aegypti larvae, a different mode of pathogenicity to that normally seen in terrestrial hosts. Blastospores constitute a different form of inoculum produced by this fungus when cultured in liquid media and although blastospores are generally considered to be more virulent than conidia no evidence has been presented to explain why. In our study, using a range of biochemical, molecular and microscopy methods, the infection process of Metarhizium brunneum (formerly M. anisopliae) ARSEF 4556 blastospores was investigated. It appears that the blastospores, unlike conidia, readily adhere to and penetrate mosquito larval cuticle. The blastospores are readily ingested by the larvae but unlike the conidia are able infect the insect through the gut and rapidly invade the haemocoel. The fact that pathogenicity related genes were upregulated in blastospores exposed to larvae prior to invasion, suggests the fungus was detecting host derived cues. Similarly, immune and defence genes were upregulated in the host prior to infection suggesting mosquitoes were also able to detect pathogen-derived cues. The hydrophilic blastospores produce copious mucilage, which probably facilitates adhesion to the host but do not appear to depend on production of Pr1, a cuticle degrading subtilisin protease, for penetration since protease inhibitors did not significantly alter blastospore virulence. The fact the blastospores have multiple routes of entry (cuticle and gut) may explain why this form of the inoculum killed Ae. aegypti larvae in a relatively short time (12-24hrs), significantly quicker than when larvae were exposed to conidia. This study shows that selecting the appropriate form of inoculum is important for efficacious control of disease vectors such as Ae. aegypti

Disease Expression in Autosomal Recessive Retinal Dystrophy Associated With Mutations in the DRAM2 Gene

Purpose: To determine the disease course of retinal dystrophy caused by recessive variants in the DRAM2 (damage-regulated autophagy modulator 2) gene. Methods: Sixteen individuals with DRAM2-retinopathy were examined (six families; age range, 19–56 years, includes one pre-symptomatic case). The change in visual acuity over time was studied, and electrophysiology (n = 6), retina-tracking perimetry (n = 1), fundus autofluorescence (FAF) imaging (n = 6), and optical coherence tomography (OCT; n = 12) were performed. Results: All symptomatic patients presented with central visual loss (15/15) unaccompanied either by nyctalopia or light-hypersensitivity; most (11/15) developed symptoms in the third decade of life. A granular macular appearance, often with associated white/yellow dots, was an early fundoscopic feature. There was an ill-defined ring of hyperautofluorescence on FAF. Optical coherence tomography revealed loss of the ellipsoid zone perifoveally in a 19-year-old pre-symptomatic individual. The central atrophic area enlarged over time and fundoscopy showed peripheral degeneration in seven of the nine individuals that were examined ≥10 years after becoming symptomatic; some of these subjects developed nyctalopia and light hypersensitivity. Electrophysiology revealed generalized retinal dysfunction in three of the five individuals that were tested ≥10 years after becoming symptomatic. Conclusions: Patients with DRAM2-retinopathy are typically asymptomatic in the first two decades of life and present with central visual loss and a maculopathy. A faint hyperautofluorescent ring on FAF can be a suggestive feature. The retinal periphery is frequently affected later in the disease process. Photoreceptor degeneration is likely to be the primary event and future studies on DRAM2-retinopathy are expected to provide important insights into retinal autophagy

Formulation, characterisation and flexographic printing of novel Boger fluids to assess the effects of ink elasticity on print uniformity

Model elastic inks were formulated, rheologically characterised in shear and extension, and printed via flexography to assess the impact of ink elasticity on print uniformity. Flexography is a roll-to-roll printing process with great potential in the mass production of printed electronics for which understanding layer uniformity and the influence of rheology is of critical importance. A new set of flexo-printable Boger fluids was formulated by blending polyvinyl alcohol and high molecular weight polyacrylamide to provide inks of varying elasticity. During print trials, the phenomenon of viscous fingering was observed in all prints, with those of the Newtonian ink exhibiting a continuous striping in the printing direction. Increasing elasticity significantly influenced this continuity, disrupting it and leading to a quantifiable decrease in the overall relative size of the printed finger features. As such, ink elasticity was seen to have a profound effect on flexographic printing uniformity, showing the rheological tuning of inks may be a route to obtaining specific printed features

Jet energy measurement with the ATLAS detector in proton-proton collisions at root s=7 TeV

The jet energy scale and its systematic uncertainty are determined for jets measured with the ATLAS detector at the LHC in proton-proton collision data at a centre-of-mass energy of √s = 7TeV corresponding to an integrated luminosity of 38 pb-1. Jets are reconstructed with the anti-kt algorithm with distance parameters R=0. 4 or R=0. 6. Jet energy and angle corrections are determined from Monte Carlo simulations to calibrate jets with transverse momenta pT≥20 GeV and pseudorapidities {pipe}η{pipe}<4. 5. The jet energy systematic uncertainty is estimated using the single isolated hadron response measured in situ and in test-beams, exploiting the transverse momentum balance between central and forward jets in events with dijet topologies and studying systematic variations in Monte Carlo simulations. The jet energy uncertainty is less than 2. 5 % in the central calorimeter region ({pipe}η{pipe}<0. 8) for jets with 60≤pT<800 GeV, and is maximally 14 % for pT<30 GeV in the most forward region 3. 2≤{pipe}η{pipe}<4. 5. The jet energy is validated for jet transverse momenta up to 1 TeV to the level of a few percent using several in situ techniques by comparing a well-known reference such as the recoiling photon pT, the sum of the transverse momenta of tracks associated to the jet, or a system of low-pT jets recoiling against a high-pT jet. More sophisticated jet calibration schemes are presented based on calorimeter cell energy density weighting or hadronic properties of jets, aiming for an improved jet energy resolution and a reduced flavour dependence of the jet response. The systematic uncertainty of the jet energy determined from a combination of in situ techniques is consistent with the one derived from single hadron response measurements over a wide kinematic range. The nominal corrections and uncertainties are derived for isolated jets in an inclusive sample of high-pT jets. Special cases such as event topologies with close-by jets, or selections of samples with an enhanced content of jets originating from light quarks, heavy quarks or gluons are also discussed and the corresponding uncertainties are determined. © 2013 CERN for the benefit of the ATLAS collaboration

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

Search for R-parity-violating supersymmetry in events with four or more leptons in sqrt(s) =7 TeV pp collisions with the ATLAS detector

A search for new phenomena in final states with four or more leptons (electrons or muons) is presented. The analysis is based on 4.7 fb−1 of $\sqrt{s}=7\;\mathrm{TeV}$ proton-proton collisions delivered by the Large Hadron Collider and recorded with the ATLAS detector. Observations are consistent with Standard Model expectations in two signal regions: one that requires moderate values of missing transverse momentum and another that requires large effective mass. The results are interpreted in a simplified model of R-parity-violating supersymmetry in which a 95% CL exclusion region is set for charged wino masses up to 540 GeV. In an R-parity-violating MSUGRA/CMSSM model, values of m 1/2 up to 820 GeV are excluded for 10 < tan β < 40

Search for high-mass resonances decaying to dilepton final states in pp collisions at s√=7 TeV with the ATLAS detector

The ATLAS detector at the Large Hadron Collider is used to search for high-mass resonances decaying to an electron-positron pair or a muon-antimuon pair. The search is sensitive to heavy neutral Z′ gauge bosons, Randall-Sundrum gravitons, Z * bosons, techni-mesons, Kaluza-Klein Z/γ bosons, and bosons predicted by Torsion models. Results are presented based on an analysis of pp collisions at a center-of-mass energy of 7 TeV corresponding to an integrated luminosity of 4.9 fb−1 in the e + e − channel and 5.0 fb−1 in the μ + μ −channel. A Z ′ boson with Standard Model-like couplings is excluded at 95 % confidence level for masses below 2.22 TeV. A Randall-Sundrum graviton with coupling k/MPl=0.1 is excluded at 95 % confidence level for masses below 2.16 TeV. Limits on the other models are also presented, including Technicolor and Minimal Z′ Models

Long-Branch Attraction Bias and Inconsistency in Bayesian Phylogenetics

Bayesian inference (BI) of phylogenetic relationships uses the same probabilistic models of evolution as its precursor maximum likelihood (ML), so BI has generally been assumed to share ML's desirable statistical properties, such as largely unbiased inference of topology given an accurate model and increasingly reliable inferences as the amount of data increases. Here we show that BI, unlike ML, is biased in favor of topologies that group long branches together, even when the true model and prior distributions of evolutionary parameters over a group of phylogenies are known. Using experimental simulation studies and numerical and mathematical analyses, we show that this bias becomes more severe as more data are analyzed, causing BI to infer an incorrect tree as the maximum a posteriori phylogeny with asymptotically high support as sequence length approaches infinity. BI's long branch attraction bias is relatively weak when the true model is simple but becomes pronounced when sequence sites evolve heterogeneously, even when this complexity is incorporated in the model. This bias—which is apparent under both controlled simulation conditions and in analyses of empirical sequence data—also makes BI less efficient and less robust to the use of an incorrect evolutionary model than ML. Surprisingly, BI's bias is caused by one of the method's stated advantages—that it incorporates uncertainty about branch lengths by integrating over a distribution of possible values instead of estimating them from the data, as ML does. Our findings suggest that trees inferred using BI should be interpreted with caution and that ML may be a more reliable framework for modern phylogenetic analysis

Multiple Chromosomal Rearrangements Structured the Ancestral Vertebrate Hox-Bearing Protochromosomes

While the proposal that large-scale genome expansions occurred early in vertebrate evolution is widely accepted, the exact mechanisms of the expansion—such as a single or multiple rounds of whole genome duplication, bloc chromosome duplications, large-scale individual gene duplications, or some combination of these—is unclear. Gene families with a single invertebrate member but four vertebrate members, such as the Hox clusters, provided early support for Ohno's hypothesis that two rounds of genome duplication (the 2R-model) occurred in the stem lineage of extant vertebrates. However, despite extensive study, the duplication history of the Hox clusters has remained unclear, calling into question its usefulness in resolving the role of large-scale gene or genome duplications in early vertebrates. Here, we present a phylogenetic analysis of the vertebrate Hox clusters and several linked genes (the Hox “paralogon”) and show that different phylogenies are obtained for Dlx and Col genes than for Hox and ErbB genes. We show that these results are robust to errors in phylogenetic inference and suggest that these competing phylogenies can be resolved if two chromosomal crossover events occurred in the ancestral vertebrate. These results resolve conflicting data on the order of Hox gene duplications and the role of genome duplication in vertebrate evolution and suggest that a period of genome reorganization occurred after genome duplications in early vertebrates