Ashoka D. Polpitiya

DAnTE: a statistical tool for quantitative analysis of -omics data

UNLABELLED Data Analysis Tool Extension (DAnTE) is a statistical tool designed to address challenges associated with quantitative bottom-up, shotgun proteomics data. This tool has also been demonstrated for microarray data and can easily be extended to other high-throughput data types. DAnTE features selected normalization methods, missing value imputation algorithms, peptide-to-protein rollup methods, an extensive array of plotting functions and a comprehensive hypothesis-testing scheme that can handle unbalanced data and random effects. The graphical user interface (GUI) is designed to be very intuitive and user friendly. AVAILABILITY DAnTE may be downloaded free of charge at http://omics.pnl.gov/software/. SUPPLEMENTARY INFORMATION An example dataset with instructions on how to perform a series of analysis steps is available at http://omics.pnl.gov/software/

Omic data from evolved E. coli are consistent with computed optimal growth from genome‐scale models

After hundreds of generations of adaptive evolution at exponential growth, Escherichia coli grows as predicted using flux balance analysis (FBA) on genome‐scale metabolic models (GEMs). However, it is not known whether the predicted pathway usage in FBA solutions is consistent with gene and protein expression in the wild‐type and evolved strains. Here, we report that >98% of active reactions from FBA optimal growth solutions are supported by transcriptomic and proteomic data. Moreover, when E. coli adapts to growth rate selective pressure, the evolved strains upregulate genes within the optimal growth predictions, and downregulate genes outside of the optimal growth solutions. In addition, bottlenecks from dosage limitations of computationally predicted essential genes are overcome in the evolved strains. We also identify regulatory processes that may contribute to the development of the optimal growth phenotype in the evolved strains, such as the downregulation of known regulons and stringent response suppression. Thus, differential gene and protein expression from wild‐type and adaptively evolved strains supports observed growth phenotype changes, and is consistent with GEM‐computed optimal growth states.

The Drosophila melanogaster sperm proteome-II (DmSP-II)

Advances in mass spectrometry technology, high-throughput proteomics and genome annotations have resulted in significant increases in our molecular understanding of sperm composition. Using improved separation and detection methods and an updated genome annotation, a re-analysis of the Drosophila melanogaster sperm proteome (DmSP) has resulted in the identification of 956 sperm proteins. Comparative analysis with our previous proteomic dataset revealed 766 new proteins and an updated sperm proteome containing a total of 1108 proteins, termed the DmSP-II. This expanded dataset includes additional proteins with predicted sperm functions and confirms previous findings concerning the genomic organization of sperm loci. Bioinformatic and protein network analyses demonstrated high quality and reproducibility of proteome coverage across three replicate mass spectrometry runs. The use of whole-cell proteomics in conjunction with characterized phenotypes, functional annotations and pathway information has advanced our systems level understanding of sperm proteome functional networks.

DanteR: an extensible R-based tool for quantitative analysis of -omics data.

MOTIVATION The size and complex nature of mass spectrometry-based proteomics datasets motivate development of specialized software for statistical data analysis and exploration. We present DanteR, a graphical R package that features extensive statistical and diagnostic functions for quantitative proteomics data analysis, including normalization, imputation, hypothesis testing, interactive visualization and peptide-to-protein rollup. More importantly, users can easily extend the existing functionality by including their own algorithms under the Add-On tab. AVAILABILITY DanteR and its associated user guide are available for download free of charge at http://omics.pnl.gov/software/. We have an updated binary source for the DanteR package up on our website together with a vignettes document. For Windows, a single click automatically installs DanteR along with the R programming environment. For Linux and Mac OS X, users must install R and then follow instructions on the DanteR website for package installation. CONTACT rds@pnnl.gov.

Plasticity of the Systemic Inflammatory Response to Acute Infection during Critical Illness: Development of the Riboleukogram

Background Diagnosis of acute infection in the critically ill remains a challenge. We hypothesized that circulating leukocyte transcriptional profiles can be used to monitor the host response to and recovery from infection complicating critical illness. Methodology/Principal Findings A translational research approach was employed. Fifteen mice underwent intratracheal injections of live P. aeruginosa, P. aeruginosa endotoxin, live S. pneumoniae, or normal saline. At 24 hours after injury, GeneChip microarray analysis of circulating buffy coat RNA identified 219 genes that distinguished between the pulmonary insults and differences in 7-day mortality. Similarly, buffy coat microarray expression profiles were generated from 27 mechanically ventilated patients every two days for up to three weeks. Significant heterogeneity of VAP microarray profiles was observed secondary to patient ethnicity, age, and gender, yet 85 genes were identified with consistent changes in abundance during the seven days bracketing the diagnosis of VAP. Principal components analysis of these 85 genes appeared to differentiate between the responses of subjects who did versus those who did not develop VAP, as defined by a general trajectory (riboleukogram) for the onset and resolution of VAP. As patients recovered from critical illness complicated by acute infection, the riboleukograms converged, consistent with an immune attractor. Conclusions/Significance Here we present the culmination of a mouse pneumonia study, demonstrating for the first time that disease trajectories derived from microarray expression profiles can be used to quantitatively track the clinical course of acute disease and identify a state of immune recovery. These data suggest that the onset of an infection-specific transcriptional program may precede the clinical diagnosis of pneumonia in patients. Moreover, riboleukograms may help explain variance in the host response due to differences in ethnic background, gender, and pathogen. Prospective clinical trials are indicated to validate our results and test the clinical utility of riboleukograms.

Diurnal rhythms result in significant changes in the cellular protein complement in the cyanobacterium Cyanothece 51142.

Cyanothece sp. ATCC 51142 is a diazotrophic cyanobacterium notable for its ability to perform oxygenic photosynthesis and dinitrogen fixation in the same single cell. Previous transcriptional analysis revealed that the existence of these incompatible cellular processes largely depends on tightly synchronized expression programs involving ∼30% of genes in the genome. To expand upon current knowledge, we have utilized sensitive proteomic approaches to examine the impact of diurnal rhythms on the protein complement in Cyanothece 51142. We found that 250 proteins accounting for ∼5% of the predicted ORFs from the Cyanothece 51142 genome and 20% of proteins detected under alternating light/dark conditions exhibited periodic oscillations in their abundances. Our results suggest that altered enzyme activities at different phases during the diurnal cycle can be attributed to changes in the abundance of related proteins and key compounds. The integration of global proteomics and transcriptomic data further revealed that post-transcriptional events are important for temporal regulation of processes such as photosynthesis in Cyanothece 51142. This analysis is the first comprehensive report on global quantitative proteomics in a unicellular diazotrophic cyanobacterium and uncovers novel findings about diurnal rhythms.

Geometry and Control of Human Eye Movements

In this paper, we study the human oculomotor system as a simple mechanical control system. It is a well known physiological fact that all eye movements obey Listings law, which states that eye orientations form a subset consisting of rotation matrices for which the axes are orthogonal to the normal gaze direction. First, we discuss the geometry of this restricted configuration space (referred to as the Listing space). Then we formulate the system as a simple mechanical control system with a holonomic constraint. We propose a realistic model with musculotendon complexes and address the question of controlling the gaze. As an example, an optimal energy control problem is formulated and numerically solved

DtaRefinery, a Software Tool for Elimination of Systematic Errors from Parent Ion Mass Measurements in Tandem Mass Spectra Data Sets

Hybrid two-stage mass spectrometers capable of both highly accurate mass measurement and high throughput MS/MS fragmentation have become widely available in recent years, allowing for significantly better discrimination between true and false MS/MS peptide identifications by the application of a relatively narrow window for maximum allowable deviations of measured parent ion masses. To fully gain the advantage of highly accurate parent ion mass measurements, it is important to limit systematic mass measurement errors. Based on our previous studies of systematic biases in mass measurement errors, here, we have designed an algorithm and software tool that eliminates the systematic errors from the peptide ion masses in MS/MS data. We demonstrate that the elimination of the systematic mass measurement errors allows for the use of tighter criteria on the deviation of measured mass from theoretical monoisotopic peptide mass, resulting in a reduction of both false discovery and false negative rates of peptide identification. A software implementation of this algorithm called DtaRefinery reads a set of fragmentation spectra, searches for MS/MS peptide identifications using a FASTA file containing expected protein sequences, fits a regression model that can estimate systematic errors, and then corrects the parent ion mass entries by removing the estimated systematic error components. The output is a new file with fragmentation spectra with updated parent ion masses. The software is freely available.

Combined pulsed-Q dissociation and electron transfer dissociation for identification and quantification of iTRAQ-labeled phosphopeptides.

Here, we report a new approach that integrates pulsed Q dissociation (PQD) and electron transfer dissociation (ETD) techniques for confident and quantitative identification of iTRAQ-labeled phosphopeptides. The use of isobaric tags for relative and absolute quantification enables a high-throughput quantification of peptides via reporter ion signals in the low m/z range of tandem mass spectra. PQD, a form of ion trap collision activated dissociation, allows for detection of low mass-to-charge fragment ions, and electron transfer dissociation is especially useful for sequencing peptides that contain post-translational modifications. Analysis of the phosphoproteome of human fibroblast cells using a sensitive linear ion trap mass spectrometer demonstrated that this hybrid approach improves both identification and quantification of phosphopeptides. ETD improved phosphopeptide identification, while PQD provides improved quantification of iTRAQ-labeled phosphopeptides.

Using systems biology to simplify complex disease: Immune cartography

What if there was a rapid, inexpensive, and accurate blood diagnostic that could determine which patients were infected, identify the organism(s) responsible, and identify patients who were not responding to therapy? We hypothesized that systems analysis of the transcriptional activity of circulating immune effector cells could be used to identify conserved elements in the host response to systemic inflammation, and furthermore, to discriminate between sterile and infectious etiologies. We review herein a validated, systems biology approach demonstrating that 1) abdominal and pulmonary sepsis diagnoses can be made in mouse models using microarray (RNA) data from circulating blood, 2) blood microarray data can be used to differentiate between the host response to Gram-negative and Gram-positive pneumonia, 3) the endotoxin response of normal human volunteers can be mapped at the level of gene expression, and 4) a similar strategy can be used in the critically ill to follow septic patients and quantitatively determine immune recovery. These findings provide the foundation of immune cartography and demonstrate the potential of this approach for rapidly diagnosing sepsis and identifying pathogens. Further, our data suggest a new approach to determine how specific pathogens perturb the physiology of circulating leukocytes in a cell-specific manner. Large, prospective clinical trails are needed to validate the clinical utility of leukocyte RNA diagnostics (e.g., the riboleukogram).