Sashimi plots: Quantitative visualization of RNA sequencing read alignments

We introduce Sashimi plots, a quantitative multi-sample visualization of mRNA sequencing reads aligned to gene annotations. Sashimi plots are made using alignments (stored in the SAM/BAM format) and gene model annotations (in GFF format), which can be custom-made by the user or obtained from databases such as Ensembl or UCSC. We describe two implementations of Sashimi plots: (1) a stand-alone command line implementation aimed at making customizable publication quality figures, and (2) an implementation built into the Integrated Genome Viewer (IGV) browser, which enables rapid and dynamic creation of Sashimi plots for any genomic region of interest, suitable for exploratory analysis of alternatively spliced regions of the transcriptome. Isoform expression estimates outputted by the MISO program can be optionally plotted along with Sashimi plots. Sashimi plots can be used to quickly screen differentially spliced exons along genomic regions of interest and can be used in publication quality figures. The Sashimi plot software and documentation is available from: http://genes.mit.edu/burgelab/miso/docs/sashimi.htmlComment: 2 figure

A genomic and evolutionary approach reveals non-genetic drug resistance in malaria

Background: Drug resistance remains a major public health challenge for malaria treatment and eradication. Individual loci associated with drug resistance to many antimalarials have been identified, but their epistasis with other resistance mechanisms has not yet been elucidated. Results: We previously described two mutations in the cytoplasmic prolyl-tRNA synthetase (cPRS) gene that confer resistance to halofuginone. We describe here the evolutionary trajectory of halofuginone resistance of two independent drug resistance selections in Plasmodium falciparum. Using this novel methodology, we discover an unexpected non-genetic drug resistance mechanism that P. falciparum utilizes before genetic modification of the cPRS. P. falciparum first upregulates its proline amino acid homeostasis in response to halofuginone pressure. We show that this non-genetic adaptation to halofuginone is not likely mediated by differential RNA expression and precedes mutation or amplification of the cPRS gene. By tracking the evolution of the two drug resistance selections with whole genome sequencing, we further demonstrate that the cPRS locus accounts for the majority of genetic adaptation to halofuginone in P. falciparum. We further validate that copy-number variations at the cPRS locus also contribute to halofuginone resistance. Conclusions: We provide a three-step model for multi-locus evolution of halofuginone drug resistance in P. falciparum. Informed by genomic approaches, our results provide the first comprehensive view of the evolutionary trajectory malaria parasites take to achieve drug resistance. Our understanding of the multiple genetic and non-genetic mechanisms of drug resistance informs how we will design and pair future anti-malarials for clinical use. Electronic supplementary material The online version of this article (doi:10.1186/s13059-014-0511-2) contains supplementary material, which is available to authorized users

Responses to Bacteria, Virus, and Malaria Distinguish the Etiology of Pediatric Clinical Pneumonia

RationalePlasma-detectable biomarkers that rapidly and accurately diagnose bacterial infections in children with suspected pneumonia could reduce the morbidity of respiratory disease and decrease the unnecessary use of antibiotic therapy.ObjectivesUsing 56 markers measured in a multiplexed immunoassay, we sought to identify proteins and protein combinations that could discriminate bacterial from viral or malarial diagnoses.MethodsWe selected 80 patients with clinically diagnosed pneumonia (as defined by the World Health Organization) who also met criteria for bacterial, viral, or malarial infection based on clinical, radiographic, and laboratory results. Ten healthy community control subjects were enrolled to assess marker reliability. Patients were subdivided into two sets: one for identifying potential markers and another for validating them.Measurements and main resultsThree proteins (haptoglobin, tumor necrosis factor receptor 2 or IL-10, and tissue inhibitor of metalloproteinases 1) were identified that, when combined through a classification tree signature, accurately classified patients into bacterial, malarial, and viral etiologies and misclassified only one patient with bacterial pneumonia from the validation set. The overall sensitivity and specificity of this signature for the bacterial diagnosis were 96 and 86%, respectively. Alternative combinations of markers with comparable accuracy were selected by support vector machine and regression models and included haptoglobin, IL-10, and creatine kinase-MB.ConclusionsCombinations of plasma proteins accurately identified children with a respiratory syndrome who were likely to have bacterial infections and who would benefit from antibiotic therapy. When used in conjunction with malaria diagnostic tests, they may improve diagnostic specificity and simplify treatment decisions for clinicians

Transcriptional categorization of the etiology of pneumonia syndrome in pediatric patients in malaria endemic areas

<p><strong><em>Background.</em></strong> Pediatric acute respiratory distress in tropical settings is very common. Bacterial pneumonia is a major contributor to morbidity and mortality and requires adequate diagnosis for correct treatment. A rapid test that could identify bacterial (vs. other) infections would have great clinical utility.</p> <p><strong><em>Methods and Findings.</em></strong> We performed RNA-Seq and analyzed the transcriptomes of 68 pediatric patients with well-characterized clinical phenotype to identify transcriptional features associated with each disease class. We refined the features to predictive models (support vector machine, elastic net) and validated those models on an independent test set of 37 patients (80 – 85% accuracy).</p> <p><strong><em>Conclusions.</em></strong> We have identified sets of genes that are differentially expressed in pediatric patients with pneumonia syndrome attributable to different infections and requiring different therapeutic interventions. This study demonstrates that human transcription signatures in infected patients recapitulate the underlying biology and provide models for predicting a bacterial diagnosis to inform treatment</p> <ul> </ul

Proteomics, Post-translational Modifications, and Integrative Analyses Reveal Molecular Heterogeneity within Medulloblastoma Subgroups

There is a pressing need to identify therapeutic targets in tumors with low mutation rates such as the malignant pediatric brain tumor medulloblastoma. To address this challenge, we quantitatively profiled global proteomes and phospho-proteomes of 45 medulloblastoma samples. Integrated analyses revealed that tumors with similar RNA expression vary extensively at the post-transcriptional and post-translational levels. We identified distinct pathways associated with two subsets of SHH tumors, and found post-translational modifications of MYC that are associated with poor outcomes in group 3 tumors. We found kinases associated with subtypes and showed that inhibiting PRKDC sensitizes MYC-driven cells to radiation. Our study shows that proteomics enables a more comprehensive, functional readout, providing a foundation for future therapeutic strategies