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

Paroxysmal autonomic instability with dystonia in a patient with tuberculous meningitis: a case report

<p>Abstract</p> <p>Introduction</p> <p>This case report describes an extremely rare combination of paroxysmal autonomic instability with dystonia and tuberculous meningitis. Paroxysmal autonomic instability with dystonia is normally associated with severe traumatic brain injury.</p> <p>Case presentation</p> <p>A 69-year-old man of Indonesian descent was initially suspected of having a community-acquired pneumonia, which was seen on chest X-ray and computed tomography of the chest. However, a bronchoscopy showed no abnormalities. He was treated with amoxicillin-clavulanic acid in combination with ciprofloxacin. However, nine days after admission he was disorientated and complained of headache. Neurological examination revealed no further abnormalities. A lumbar puncture revealed no evidence of meningitis. He was then transferred to our hospital. At that time, initial cultures of bronchial fluid for <it>Mycobacterium tuberculosis </it>turned positive, as well as polymerase chain reaction for <it>Mycobacterium tuberculosis</it>. Later, during his stay in our intensive care unit, he developed periods with hypertension, sinus tachycardia, excessive transpiration, decreased oxygen saturation with tachypnea, pink foamy sputum, and high fever. This constellation of symptoms was accompanied by dystonia in the first days. These episodes lasted approximately 30 minutes and improved after administration of morphine, benzodiazepines or clonidine. Magnetic resonance imaging showed an abnormal signal in the region of the hippocampus, thalamus and the anterior parts of the lentiform nucleus and caudate nucleus.</p> <p>Conclusions</p> <p>In patients with (tuberculous) meningitis and episodes of extreme hypertension and fever, paroxysmal autonomic instability with dystonia should be considered.</p

Constraints on Nucleon Decay via "Invisible" Modes from the Sudbury Neutrino Observatory

Data from the Sudbury Neutrino Observatory have been used to constrain the lifetime for nucleon decay to ``invisible'' modes, such as n -> 3 nu. The analysis was based on a search for gamma-rays from the de-excitation of the residual nucleus that would result from the disappearance of either a proton or neutron from O16. A limit of tau_inv > 2 x 10^{29} years is obtained at 90% confidence for either neutron or proton decay modes. This is about an order of magnitude more stringent than previous constraints on invisible proton decay modes and 400 times more stringent than similar neutron modes.Comment: Update includes missing efficiency factor (limits change by factor of 2) Submitted to Physical Review Letter

First Neutrino Observations from the Sudbury Neutrino Observatory

The first neutrino observations from the Sudbury Neutrino Observatory are presented from preliminary analyses. Based on energy, direction and location, the data in the region of interest appear to be dominated by 8B solar neutrinos, detected by the charged current reaction on deuterium and elastic scattering from electrons, with very little background. Measurements of radioactive backgrounds indicate that the measurement of all active neutrino types via the neutral current reaction on deuterium will be possible with small systematic uncertainties. Quantitative results for the fluxes observed with these reactions will be provided when further calibrations have been completed.Comment: Latex, 7 pages, 10 figures, Invited paper at Neutrino 2000 Conference, Sudbury, Canada, June 16-21, 2000 to be published in the Proceeding

Isoform Diversity and Regulation in Peripheral and Central Neurons Revealed through RNA-Seq

To fully understand cell type identity and function in the nervous system there is a need to understand neuronal gene expression at the level of isoform diversity. Here we applied Next Generation Sequencing of the transcriptome (RNA-Seq) to purified sensory neurons and cerebellar granular neurons (CGNs) grown on an axonal growth permissive substrate. The goal of the analysis was to uncover neuronal type specific isoforms as a prelude to understanding patterns of gene expression underlying their intrinsic growth abilities. Global gene expression patterns were comparable to those found for other cell types, in that a vast majority of genes were expressed at low abundance. Nearly 18% of gene loci produced more than one transcript. More than 8000 isoforms were differentially expressed, either to different degrees in different neuronal types or uniquely expressed in one or the other. Sensory neurons expressed a larger number of genes and gene isoforms than did CGNs. To begin to understand the mechanisms responsible for the differential gene/isoform expression we identified transcription factor binding sites present specifically in the upstream genomic sequences of differentially expressed isoforms, and analyzed the 3′ untranslated regions (3′ UTRs) for microRNA (miRNA) target sites. Our analysis defines isoform diversity for two neuronal types with diverse axon growth capabilities and begins to elucidate the complex transcriptional landscape in two neuronal populations

Nephele: genotyping via complete composition vectors and MapReduce

<p>Abstract</p> <p>Background</p> <p>Current sequencing technology makes it practical to sequence many samples of a given organism, raising new challenges for the processing and interpretation of large genomics data sets with associated metadata. Traditional computational phylogenetic methods are ideal for studying the evolution of gene/protein families and using those to infer the evolution of an organism, but are less than ideal for the study of the whole organism mainly due to the presence of insertions/deletions/rearrangements. These methods provide the researcher with the ability to group a set of samples into distinct genotypic groups based on sequence similarity, which can then be associated with metadata, such as host information, pathogenicity, and time or location of occurrence. Genotyping is critical to understanding, at a genomic level, the origin and spread of infectious diseases. Increasingly, genotyping is coming into use for disease surveillance activities, as well as for microbial forensics. The classic genotyping approach has been based on phylogenetic analysis, starting with a multiple sequence alignment. Genotypes are then established by expert examination of phylogenetic trees. However, these traditional single-processor methods are suboptimal for rapidly growing sequence datasets being generated by next-generation DNA sequencing machines, because they increase in computational complexity quickly with the number of sequences.</p> <p>Results</p> <p>Nephele is a suite of tools that uses the complete composition vector algorithm to represent each sequence in the dataset as a vector derived from its constituent k-mers by passing the need for multiple sequence alignment, and affinity propagation clustering to group the sequences into genotypes based on a distance measure over the vectors. Our methods produce results that correlate well with expert-defined clades or genotypes, at a fraction of the computational cost of traditional phylogenetic methods run on traditional hardware. Nephele can use the open-source Hadoop implementation of MapReduce to parallelize execution using multiple compute nodes. We were able to generate a neighbour-joined tree of over 10,000 16S samples in less than 2 hours.</p> <p>Conclusions</p> <p>We conclude that using Nephele can substantially decrease the processing time required for generating genotype trees of tens to hundreds of organisms at genome scale sequence coverage.</p

Experimental Incubations Elicit Profound Changes in Community Transcription in OMZ Bacterioplankton

Sequencing of microbial community RNA (metatranscriptome) is a useful approach for assessing gene expression in microorganisms from the natural environment. This method has revealed transcriptional patterns in situ, but can also be used to detect transcriptional cascades in microcosms following experimental perturbation. Unambiguously identifying differential transcription between control and experimental treatments requires constraining effects that are simply due to sampling and bottle enclosure. These effects remain largely uncharacterized for “challenging” microbial samples, such as those from anoxic regions that require special handling to maintain in situ conditions. Here, we demonstrate substantial changes in microbial transcription induced by sample collection and incubation in experimental bioreactors. Microbial communities were sampled from the water column of a marine oxygen minimum zone by a pump system that introduced minimal oxygen contamination and subsequently incubated in bioreactors under near in situ oxygen and temperature conditions. Relative to the source water, experimental samples became dominated by transcripts suggestive of cell stress, including chaperone, protease, and RNA degradation genes from diverse taxa, with strong representation from SAR11-like alphaproteobacteria. In tandem, transcripts matching facultative anaerobic gammaproteobacteria of the Alteromonadales (e.g., Colwellia) increased 4–13 fold up to 43% of coding transcripts, and encoded a diverse gene set suggestive of protein synthesis and cell growth. We interpret these patterns as taxon-specific responses to combined environmental changes in the bioreactors, including shifts in substrate or oxygen availability, and minor temperature and pressure changes during sampling with the pump system. Whether such changes confound analysis of transcriptional patterns may vary based on the design of the experiment, the taxonomic composition of the source community, and on the metabolic linkages between community members. These data highlight the impressive capacity for transcriptional changes within complex microbial communities, underscoring the need for caution when inferring in situ metabolism based on transcript abundances in experimental incubations

How Linear Tension Converts to Curvature: Geometric Control of Bone Tissue Growth

This study investigated how substrate geometry influences in-vitro tissue formation at length scales much larger than a single cell. Two-millimetre thick hydroxyapatite plates containing circular pores and semi-circular channels of 0.5 mm radius, mimicking osteons and hemi-osteons respectively, were incubated with MC3T3-E1 cells for 4 weeks. The amount and shape of the tissue formed in the pores, as measured using phase contrast microscopy, depended on the substrate geometry. It was further demonstrated, using a simple geometric model, that the observed curvature-controlled growth can be derived from the assembly of tensile elements on a curved substrate. These tensile elements are cells anchored on distant points of the curved surface, thus creating an actin “chord” by generating tension between the adhesion sites. Such a chord model was used to link the shape of the substrate to cell organisation and tissue patterning. In a pore with a circular cross-section, tissue growth increases the average curvature of the surface, whereas a semi-circular channel tends to be flattened out. Thereby, a single mechanism could describe new tissue growth in both cortical and trabecular bone after resorption due to remodelling. These similarities between in-vitro and in-vivo patterns suggest geometry as an important signal for bone remodelling

Update of complications and functional outcome of the ileo-pouch anal anastomosis: overview of evidence and meta-analysis of 96 observational studies

Item does not contain fulltextOBJECTIVE: The objective of this study is to provide a comprehensive update of the outcome of the ileo-pouch anal anastomosis (IPAA). DATA SOURCES: An extensive search in PubMed, EMBASE, and The Cochrane Library was conducted. STUDY SELECTION AND DATA EXTRACTION: All studies published after 2000 reporting on complications or functional outcome after a primary open IPAA procedure for UC or FAP were selected. Study characteristics, functional outcome, and complications were extracted. DATA SYNTHESIS: A review with similar methodology conducted 10 years earlier was used to evaluate developments in outcome over time. Pooled estimates were compared using a random-effects logistic meta-analyzing technique. Analyses focusing on the effect of time of study conductance, centralization, and variation in surgical techniques were performed. RESULTS: Fifty-three studies including 14,966 patients were included. Pooled rates of pouch failure and pelvic sepsis were 4.3% (95% CI, 3.5-6.3) and 7.5% (95% CI 6.1-9.1), respectively. Compared to studies published before 2000, a reduction of 2.5% was observed in the pouch failure rate (p = 0.0038). Analysis on the effect of the time of study conductance confirmed a decline in pouch failure. Functional outcome remained stable over time, with a 24-h defecation frequency of 5.9 (95% CI, 5.0-6.9). Technical surgery aspects did not have an important effect on outcome. CONCLUSION: This review provides up to date outcome estimates of the IPAA procedure that can be useful as reference values for practice and research. It is also shows a reduction in pouch failure over time.1 juli 201