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

The Caenorhabditis chemoreceptor gene families

Background: Chemoreceptor proteins mediate the first step in the transduction of environmental chemical stimuli, defining the breadth of detection and conferring stimulus specificity. Animal genomes contain families of genes encoding chemoreceptors that mediate taste, olfaction, and pheromone responses. The size and diversity of these families reflect the biology of chemoperception in specific species. Results: Based on manual curation and sequence comparisons among putative G-protein-coupled chemoreceptor genes in the nematode Caenorhabditis elegans, we identified approximately 1300 genes and 400 pseudogenes in the 19 largest gene families, most of which fall into larger superfamilies. In the related species C. briggsae and C. remanei, we identified most or all genes in each of the 19 families. For most families, C. elegans has the largest number of genes and C. briggsae the smallest number, suggesting changes in the importance of chemoperception among the species. Protein trees reveal family-specific and species-specific patterns of gene duplication and gene loss. The frequency of strict orthologs varies among the families, from just over 50% in two families to less than 5% in three families. Several families include large species-specific expansions, mostly in C. elegans and C. remanei. Conclusion: Chemoreceptor gene families in Caenorhabditis species are large and evolutionarily dynamic as a result of gene duplication and gene loss. These dynamics shape the chemoreceptor gene complements in Caenorhabditis species and define the receptor space available for chemosensory responses. To explain these patterns, we propose the gray pawn hypothesis: individual genes are of little significance, but the aggregate of a large number of diverse genes is required to cover a large phenotype space.JHT was supported by NIH grant RO1GM48700 and HMR by R01AI56081

Glial Hsp70 Protects K+ Homeostasis in the Drosophila Brain during Repetitive Anoxic Depolarization

Neural tissue is particularly vulnerable to metabolic stress and loss of ion homeostasis. Repetitive stress generally leads to more permanent dysfunction but the mechanisms underlying this progression are poorly understood. We investigated the effects of energetic compromise in Drosophila by targeting the Na+/K+-ATPase. Acute ouabain treatment of intact flies resulted in subsequent repetitive comas that led to death and were associated with transient loss of K+ homeostasis in the brain. Heat shock pre-conditioned flies were resistant to ouabain treatment. To control the timing of repeated loss of ion homeostasis we subjected flies to repetitive anoxia while recording extracellular [K+] in the brain. We show that targeted expression of the chaperone protein Hsp70 in glial cells delays a permanent loss of ion homeostasis associated with repetitive anoxic stress and suggest that this is a useful model for investigating molecular mechanisms of neuroprotection

Recurrent Modification of a Conserved Cis-Regulatory Element Underlies Fruit Fly Pigmentation Diversity

The development of morphological traits occurs through the collective action of networks of genes connected at the level of gene expression. As any node in a network may be a target of evolutionary change, the recurrent targeting of the same node would indicate that the path of evolution is biased for the relevant trait and network. Although examples of parallel evolution have implicated recurrent modification of the same gene and cis-regulatory element (CRE), little is known about the mutational and molecular paths of parallel CRE evolution. In Drosophila melanogaster fruit flies, the Bric-à-brac (Bab) transcription factors control the development of a suite of sexually dimorphic traits on the posterior abdomen. Female-specific Bab expression is regulated by the dimorphic element, a CRE that possesses direct inputs from body plan (ABD-B) and sex-determination (DSX) transcription factors. Here, we find that the recurrent evolutionary modification of this CRE underlies both intraspecific and interspecific variation in female pigmentation in the melanogaster species group. By reconstructing the sequence and regulatory activity of the ancestral Drosophila melanogaster dimorphic element, we demonstrate that a handful of mutations were sufficient to create independent CRE alleles with differing activities. Moreover, intraspecific and interspecific dimorphic element evolution proceeded with little to no alterations to the known body plan and sex-determination regulatory linkages. Collectively, our findings represent an example where the paths of evolution appear biased to a specific CRE, and drastic changes in function were accompanied by deep conservation of key regulatory linkages. © 2013 Rogers et al

Transforming growth factor beta signaling: The master sculptor of fingers

Transforming growth factor beta (TGF?) constitutes a large and evolutionarily conserved superfamily of secreted factors that play essential roles in embryonic development, cancer, tissue regeneration, and human degenerative pathology. Studies of this signaling cascade in the regulation of cellular and tissue changes in the three-dimensional context of a developing embryo have notably advanced in the understanding of the action mechanism of these growth factors. In this review, we address the role of TGF? signaling in the developing limb, focusing on its essential function in the morphogenesis of the autopod. As we discuss in this work, modern mouse genetic experiments together with more classical embryological approaches in chick embryos, provided very valuable information concerning the role of TGF? and Activin family members in the morphogenesis of the digits of tetrapods, including the formation of phalanxes, digital tendons, and interphalangeal joints. We emphasize the importance of the Activin and TGF? proteins as digit inducing factors and their critical interaction with the BMP signaling to sculpt the hand and foot morphology

“Even if You Know Everything You Can Forget”: Health Worker Perceptions of Mobile Phone Text-Messaging to Improve Malaria Case-Management in Kenya

This paper presents the results of a qualitative study to investigate the perceptions and experiences of health workers involved in a a cluster-randomized controlled trial of a novel intervention to improve health worker malaria case-management in 107 government health facilities in Kenya. The intervention involved sending text-messages about paediatric outpatient malaria case-management accompanied by “motivating” quotes to health workers’ mobile phones. Ten malaria messages were developed reflecting recommendations from the Kenyan national guidelines. Two messages were delivered per day for 5 working days and the process was repeated for 26 weeks (May to October 2009). The accompanying quotes were unique to each message. The intervention was delivered to 119 health workers and there were significant improvements in correct artemether-lumefantrine (AL) management both immediately after the intervention (November 2009) and 6 months later (May 2010). In-depth interviews with 24 health workers were undertaken to investigate the possible drivers of this change. The results suggest high acceptance of all components of the intervention, with the active delivery of information in an on the job setting, the ready availability of new and stored text messages and the perception of being kept ‘up to date’ as important factors influencing practice. Applying the construct of stages of change we infer that in this intervention the SMS messages were operating primarily at the action and maintenance stages of behaviour change achieving their effect by creating an enabling environment and providing a prompt to action for the implementation of case management practices that had already been accepted as the clinical norm by the health workers. Future trials testing the effectiveness of SMS reminders in creating an enabling environment for the establishment of new norms in clinical practice as well as in providing a prompt to action for the implementation of the new case-management guidelines are justified

A Family of Chemoreceptors in Tribolium castaneum (Tenebrionidae: Coleoptera)

Chemoperception in invertebrates is mediated by a family of G-protein-coupled receptors (GPCR). To date nothing is known about the molecular mechanisms of chemoperception in coleopteran species. Recently the genome of Tribolium castaneum was sequenced for use as a model species for the Coleoptera. Using blast searches analyses of the T. castaneum genome with previously predicted amino acid sequences of insect chemoreceptor genes, a putative chemoreceptor family consisting of 62 gustatory receptors (Grs) and 26 olfactory receptors (Ors) was identified. The receptors have seven transmembrane domains (7TMs) and all belong to the GPCR receptor family. The expression of the T. castaneum chemoreceptor genes was investigated using quantification real- time RT-PCR and in situ whole mount RT-PCR analysis in the antennae, mouth parts, and prolegs of the adults and larvae. All of the predicted TcasGrs were expressed in the labium, maxillae, and prolegs of the adults but TcasGr13, 19, 28, 47, 62, 98, and 61 were not expressed in the prolegs. The TcasOrs were localized only in the antennae and not in any of the beetles gustatory organs with one exception; the TcasOr16 (like DmelOr83b), which was localized in the antennae, labium, and prolegs of the beetles. A group of six TcasGrs that presents a lineage with the sugar receptors subfamily in Drosophila melanogaster were localized in the lacinia of the Tribolium larvae. TcasGr1, 3, and 39, presented an ortholog to CO2 receptors in D. melanogaster and Anopheles gambiae was recorded. Low expression of almost all of the predicted chemoreceptor genes was observed in the head tissues that contain the brains and suboesophageal ganglion (SOG). These findings demonstrate the identification of a chemoreceptor family in Tribolium, which is evolutionarily related to other insect species

Prescription Medicines and the Risk of Road Traffic Crashes: A French Registry-Based Study

Using three nationwide databases in France, Ludivine Orriols, Emmanuel Lagarde, and colleagues provide evidence that prescribed medicines contribute to the risk of experiencing a road traffic crash

Diffractive Dijet Production at sqrt(s)=630 and 1800 GeV at the Fermilab Tevatron

We report a measurement of the diffractive structure function $F_{jj}^D$ of the antiproton obtained from a study of dijet events produced in association with a leading antiproton in $\bar pp$ collisions at $\sqrt s=630$ GeV at the Fermilab Tevatron. The ratio of $F_{jj}^D$ at $\sqrt s=630$ GeV to $F_{jj}^D$ obtained from a similar measurement at $\sqrt s=1800$ GeV is compared with expectations from QCD factorization and with theoretical predictions. We also report a measurement of the $\xi$ ($x$-Pomeron) and $\beta$ ($x$ of parton in Pomeron) dependence of $F_{jj}^D$ at $\sqrt s=1800$ GeV. In the region $0.035<\xi<0.095$, $|t|<1$ GeV$^2$ and $\beta<0.5$, $F_{jj}^D(\beta,\xi)$ is found to be of the form $\beta^{-1.0\pm 0.1} \xi^{-0.9\pm 0.1}$, which obeys $\beta$-$\xi$ factorization.Comment: LaTeX, 9 pages, Submitted to Phys. Rev. Letter

Measurement of the B0 anti-B0 oscillation frequency using l- D*+ pairs and lepton flavor tags

The oscillation frequency Delta-md of B0 anti-B0 mixing is measured using the partially reconstructed semileptonic decay anti-B0 -> l- nubar D*+ X. The data sample was collected with the CDF detector at the Fermilab Tevatron collider during 1992 - 1995 by triggering on the existence of two lepton candidates in an event, and corresponds to about 110 pb-1 of pbar p collisions at sqrt(s) = 1.8 TeV. We estimate the proper decay time of the anti-B0 meson from the measured decay length and reconstructed momentum of the l- D*+ system. The charge of the lepton in the final state identifies the flavor of the anti-B0 meson at its decay. The second lepton in the event is used to infer the flavor of the anti-B0 meson at production. We measure the oscillation frequency to be Delta-md = 0.516 +/- 0.099 +0.029 -0.035 ps-1, where the first uncertainty is statistical and the second is systematic.Comment: 30 pages, 7 figures. Submitted to Physical Review