John T. Prince

LC-MS alignment in theory and practice: a comprehensive algorithmic review

Liquid chromatography-mass spectrometry is widely used for comparative replicate sample analysis in proteomics, lipidomics and metabolomics. Before statistical comparison, registration must be established to match corresponding analytes from run to run. Alignment, the most popular correspondence approach, consists of constructing a function that warps the content of runs to most closely match a given reference sample. To date, dozens of correspondence algorithms have been proposed, creating a daunting challenge for practitioners in algorithm selection. Yet, existing reviews have highlighted only a few approaches. In this review, we describe 50 correspondence algorithms to facilitate practical algorithm selection. We elucidate the motivation for correspondence and analyze the limitations of current approaches, which include prohibitive runtimes, numerous user parameters, model limitations and the need for reference samples. We suggest and describe a paradigm shift for overcoming current correspondence limitations by building on known liquid chromatography-mass spectrometry behavior.

Mitochondrial fission mediates ceramide-induced metabolic disruption in skeletal muscle.

Ceramide is a sphingolipid that serves as an important second messenger in an increasing number of stress-induced pathways. Ceramide has long been known to affect the mitochondria, altering both morphology and physiology. We sought to assess the impact of ceramide on skeletal muscle mitochondrial structure and function. A primary observation was the rapid and dramatic division of mitochondria in ceramide-treated cells. This effect is likely to be a result of increased Drp1 (dynamin-related protein 1) action, as ceramide increased Drp1 expression and Drp1 inhibition prevented ceramide-induced mitochondrial fission. Further, we found that ceramide treatment reduced mitochondrial O2 consumption (i.e. respiration) in cultured myotubes and permeabilized red gastrocnemius muscle fibre bundles. Ceramide treatment also increased H2O2 levels and reduced Akt/PKB (protein kinase B) phosphorylation in myotubes. However, inhibition of mitochondrial fission via Drp1 knockdown completely protected the myotubes and fibre bundles from ceramide-induced metabolic disruption, including maintained mitochondrial respiration, reduced H2O2 levels and unaffected insulin signalling. These data suggest that the forced and sustained mitochondrial fission that results from ceramide accrual may alter metabolic function in skeletal muscle, which is a prominent site not only of energy demand (via the mitochondria), but also of ceramide accrual with weight gain.

Proteomics, lipidomics, metabolomics: a mass spectrometry tutorial from a computer scientist's point of view

BackgroundFor decades, mass spectrometry data has been analyzed to investigate a wide array of research interests, including disease diagnostics, biological and chemical theory, genomics, and drug development. Progress towards solving any of these disparate problems depends upon overcoming the common challenge of interpreting the large data sets generated. Despite interim successes, many data interpretation problems in mass spectrometry are still challenging. Further, though these challenges are inherently interdisciplinary in nature, the significant domain-specific knowledge gap between disciplines makes interdisciplinary contributions difficult.ResultsThis paper provides an introduction to the burgeoning field of computational mass spectrometry. We illustrate key concepts, vocabulary, and open problems in MS-omics, as well as provide invaluable resources such as open data sets and key search terms and references.ConclusionsThis paper will facilitate contributions from mathematicians, computer scientists, and statisticians to MS-omics that will fundamentally improve results over existing approaches and inform novel algorithmic solutions to open problems.

Novel Algorithms and the Benefits of Comparative Validation

Contact: 2robsmith@gmail.com, jtprince@chem.byu.edu Bioinformatic research has produced a large volume of proposed algorithmic solutions to a host of problems. Whether presented as a processing step in a clinical experiment or treated in a stand-alone publication, novel bioinformatic algorithms are often not subjected to the thorough comparative evaluation endured by their counterparts in other closely related fields—such as computer science—where an algorithm unevaluated against extant methods is considered unpublishable. Two audiences are interested in algorithmic publications: the practitioner, who may use the algorithm, and the researcher, who will work to develop solutions superior to those extant. We argue that failure during the review/publication process to require comparative evaluation for novel algorithms is detrimental to both parties. To demonstrate the dilemma, we conducted a case study of novel LC-MS alignment algorithms. Of the 48 publications from 2001 to 2012 that present alignment algorithms of which we are aware, 60% include no comparison to other methods. Another 20%

Dynamical structure function identifiability conditions enabling signal structure reconstruction

Networks of controlled dynamical systems exhibit a variety of interconnection patterns that could be interpreted as the structure of the system. One such interpretation of system structure is a systems signal structure, characterized as the open-loop causal dependencies among manifest variables and represented by its dynamical structure function. Although this notion of structure is among the weakest available, previous work has shown that if no a priori structural information is known about the system, not even the Boolean structure of the dynamical structure function is identifiable. Consequently, one method previously suggested for obtaining the necessary a priori structural information is to leverage knowledge about target specificity of the controlled inputs. This work extends these results to demonstrate precisely the a priori structural information that is both necessary and sufficient to reconstruct the network from input-output data. This extension is important because it significantly broadens the applicability of the identifiability conditions, enabling the design of network reconstruction experiments that were previously impossible due to practical constraints on the types of actuation mechanisms available to the engineer or scientist. The work is motivated by the proteomics problem of reconstructing the Per-Arnt-Sim Kinase pathway used in the metabolism of sugars.

Sequence and structural characterization of great salt lake bacteriophage CW02, a member of the T7-like supergroup.

ABSTRACT Halophage CW02 infects a Salinivibrio costicola-like bacterium, SA50, isolated from the Great Salt Lake. Following isolation, cultivation, and purification, CW02 was characterized by DNA sequencing, mass spectrometry, and electron microscopy. A conserved module of structural genes places CW02 in the T7 supergroup, members of which are found in diverse aquatic environments, including marine and freshwater ecosystems. CW02 has morphological similarities to viruses of the Podoviridae family. The structure of CW02, solved by cryogenic electron microscopy and three-dimensional reconstruction, enabled the fitting of a portion of the bacteriophage HK97 capsid protein into CW02 capsid density, thereby providing additional evidence that capsid proteins of tailed double-stranded DNA phages have a conserved fold. The CW02 capsid consists of bacteriophage lambda gpD-like densities that likely contribute to particle stability. Turret-like densities were found on icosahedral vertices and may represent a unique adaptation similar to what has been seen in other extremophilic viruses that infect archaea, such as Sulfolobus turreted icosahedral virus and halophage SH1.

Metriculator: quality assessment for mass spectrometry-based proteomics.

SUMMARY Quality control in mass spectrometry-based proteomics remains subjective, labor-intensive and inconsistent between laboratories. We introduce Metriculator, a software designed to facilitate long-term storage of extensive performance metrics as introduced by NIST in 2010. Metriculator features a web interface that generates interactive comparison plots for contextual understanding of metric values and an automated metric generation toolkit. The comparison plots are designed for at-a-glance determination of outliers and trends in the datasets, together with relevant statistical comparisons. Easy-to-use quantitative comparisons and a framework for integration plugins will encourage a culture of quality assurance within the proteomics community. AVAILABILITY AND IMPLEMENTATION Available under the MIT license at http://github.com/princelab/metriculator.

The Genomes, Proteomes, and Structures of Three Novel Phages That Infect the Bacillus cereus Group and Carry Putative Virulence Factors

ABSTRACT This article reports the results of studying three novel bacteriophages, JL, Shanette, and Basilisk, which infect the pathogen Bacillus cereus and carry genes that may contribute to its pathogenesis. We analyzed host range and superinfection ability, mapped their genomes, and characterized phage structure by mass spectrometry and transmission electron microscopy (TEM). The JL and Shanette genomes were 96% similar and contained 217 open reading frames (ORFs) and 220 ORFs, respectively, while Basilisk has an unrelated genome containing 138 ORFs. Mass spectrometry revealed 23 phage particle proteins for JL and 15 for Basilisk, while only 11 and 4, respectively, were predicted to be present by sequence analysis. Structural protein homology to well-characterized phages suggested that JL and Shanette were members of the family Myoviridae, which was confirmed by TEM. The third phage, Basilisk, was similar only to uncharacterized phages and is an unrelated siphovirus. Cryogenic electron microscopy of this novel phage revealed a T=9 icosahedral capsid structure with the major capsid protein (MCP) likely having the same fold as bacteriophage HK97 MCP despite the lack of sequence similarity. Several putative virulence factors were encoded by these phage genomes, including TerC and TerD involved in tellurium resistance. Host range analysis of all three phages supports genetic transfer of such factors within the B. cereus group, including B. cereus, B. anthracis, and B. thuringiensis. This study provides a basis for understanding these three phages and other related phages as well as their contributions to the pathogenicity of B. cereus group bacteria. IMPORTANCE The Bacillus cereus group of bacteria contains several human and plant pathogens, including B. cereus, B. anthracis, and B. thuringiensis. Phages are intimately linked to the evolution of their bacterial hosts and often provide virulence factors, making the study of B. cereus phages important to understanding the evolution of pathogenic strains. Herein we provide the results of detailed study of three novel B. cereus phages, two highly related myoviruses (JL and Shanette) and an unrelated siphovirus (Basilisk). The detailed characterization of host range and superinfection, together with results of genomic, proteomic, and structural analyses, reveal several putative virulence factors as well as the ability of these phages to infect different pathogenic species.

Programmed Cell Death Protein 5 Interacts with the Cytosolic Chaperonin Containing Tailless Complex Polypeptide 1 (CCT) to Regulate β-Tubulin Folding

Background: Programmed cell death protein 5 (PDCD5) has been proposed to act as a pro-apoptotic factor with tumor suppressor capabilities. Results: PDCD5 forms a complex with the cytosolic chaperonin CCT and inhibits β-tubulin folding. Conclusion: PDCD5 functions as a modulator of CCT to regulate β-tubulin folding. Significance: PDCD5 may exert its pro-apoptotic function by blocking β-tubulin folding. Programmed cell death protein 5 (PDCD5) has been proposed to act as a pro-apoptotic factor and tumor suppressor. However, the mechanisms underlying its apoptotic function are largely unknown. A proteomics search for binding partners of phosducin-like protein, a co-chaperone for the cytosolic chaperonin containing tailless complex polypeptide 1 (CCT), revealed a robust interaction between PDCD5 and CCT. PDCD5 formed a complex with CCT and β-tubulin, a key CCT-folding substrate, and specifically inhibited β-tubulin folding. Cryo-electron microscopy studies of the PDCD5·CCT complex suggested a possible mechanism of inhibition of β-tubulin folding. PDCD5 bound the apical domain of the CCTβ subunit, projecting above the folding cavity without entering it. Like PDCD5, β-tubulin also interacts with the CCTβ apical domain, but a second site is found at the sensor loop deep within the folding cavity. These orientations of PDCD5 and β-tubulin suggest that PDCD5 sterically interferes with β-tubulin binding to the CCTβ apical domain and inhibits β-tubulin folding. Given the importance of tubulins in cell division and proliferation, PDCD5 might exert its apoptotic function at least in part through inhibition of β-tubulin folding.

Massifquant: open-source Kalman filter-based XC-MS isotope trace feature detection

MOTIVATION Isotope trace (IT) detection is a fundamental step for liquid or gas chromatography mass spectrometry (XC-MS) data analysis that faces a multitude of technical challenges on complex samples. The Kalman filter (KF) application to IT detection addresses some of these challenges; it discriminates closely eluting ITs in the m/z dimension, flexibly handles heteroscedastic m/z variances and does not bin the m/z axis. Yet, the behavior of this KF application has not been fully characterized, as no cost-free open-source implementation exists and incomplete evaluation standards for IT detection persist. RESULTS Massifquant is an open-source solution for KF IT detection that has been subjected to novel and rigorous methods of performance evaluation. The presented evaluation with accompanying annotations and optimization guide sets a new standard for comparative IT detection. Compared with centWave, matchedFilter and MZMine2-alternative IT detection engines-Massifquant detected more true ITs in a real LC-MS complex sample, especially low-intensity ITs. It also offers competitive specificity and equally effective quantitation accuracy. AVAILABILITY AND IMPLEMENTATION Massifquant is integrated into XCMS with GPL license ≥ 2.0 and hosted by Bioconductor: http://bioconductor.org. Annotation data are archived at http://hdl.lib.byu.edu/1877/3232. Parameter optimization code and documentation is hosted at https://github.com/topherconley/optimize-it.