TopDownApp: An open and modular platform for analysis and visualisation of top-down proteomics data

Although Top-down (TD) proteomics techniques, aimed at the analysis of intact proteins and proteoforms, are becoming increasingly popular, efforts are needed at different levels to generalise their adoption. In this context, there are numerous improvements that are possible in the area of open science practices, including a greater application of the FAIR (Findable, Accessible, Interoperable, and Reusable) data principles. These include, for example, increased data sharing practices and readily available open data standards. Additionally, the field would benefit from the development of open data analysis workflows that can enable data reuse of public datasets, something that is increasingly common in other proteomics fields.

MetaNovo: An open-source pipeline for probabilistic peptide discovery in complex metaproteomic datasets

BACKGROUND

Microbiome research is providing important new insights into the metabolic interactions of complex microbial ecosystems involved in fields as diverse as the pathogenesis of human diseases, agriculture and climate change. Poor correlations typically observed between RNA and protein expression datasets make it hard to accurately infer microbial protein synthesis from metagenomic data. Additionally, mass spectrometry-based metaproteomic analyses typically rely on focused search sequence databases based on prior knowledge for protein identification that may not represent all the proteins present in a set of samples. Metagenomic 16S rRNA sequencing only targets the bacterial component, while whole genome sequencing is at best an indirect measure of expressed proteomes. Here we describe a novel approach, MetaNovo, that combines existing open-source software tools to perform scalable de novo sequence tag matching with a novel algorithm for probabilistic optimization of the entire UniProt knowledgebase to create tailored sequence databases for target-decoy searches directly at the proteome level, enabling metaproteomic analyses without prior expectation of sample composition or metagenomic data generation and compatible with standard downstream analysis pipelines.

RESULTS

We compared MetaNovo to published results from the MetaPro-IQ pipeline on 8 human mucosal-luminal interface samples, with comparable numbers of peptide and protein identifications, many shared peptide sequences and a similar bacterial taxonomic distribution compared to that found using a matched metagenome sequence database-but simultaneously identified many more non-bacterial peptides than the previous approaches. MetaNovo was also benchmarked on samples of known microbial composition against matched metagenomic and whole genomic sequence database workflows, yielding many more MS/MS identifications for the expected taxa, with improved taxonomic representation, while also highlighting previously described genome sequencing quality concerns for one of the organisms, and identifying an experimental sample contaminant without prior expectation.

CONCLUSIONS

By estimating taxonomic and peptide level information directly on microbiome samples from tandem mass spectrometry data, MetaNovo enables the simultaneous identification of peptides from all domains of life in metaproteome samples, bypassing the need for curated sequence databases to search. We show that the MetaNovo approach to mass spectrometry metaproteomics is more accurate than current gold standard approaches of tailored or matched genomic sequence database searches, can identify sample contaminants without prior expectation and yields insights into previously unidentified metaproteomic signals, building on the potential for complex mass spectrometry metaproteomic data to speak for itself.

Comparing Top-Down Proteoform Identification: Deconvolution, PrSM Overlap, and PTM Detection

Generating top-down tandem mass spectra (MS/MS) from complex mixtures of proteoforms benefits from improvements in fractionation, separation, fragmentation, and mass analysis. The algorithms to match MS/MS to sequences have undergone a parallel evolution, with both spectral alignment and match-counting approaches producing high-quality proteoform-spectrum matches (PrSMs). This study assesses state-of-the-art algorithms for top-down identification (ProSight PD, TopPIC, MSPathFinderT, and pTop) in their yield of PrSMs while controlling false discovery rate. We evaluated deconvolution engines (ThermoFisher Xtract, Bruker AutoMSn, Matrix Science Mascot Distiller, TopFD, and FLASHDeconv) in both ThermoFisher Orbitrap-class and Bruker maXis Q-TOF data (PXD033208) to produce consistent precursor charges and mass determinations. Finally, we sought post-translational modifications (PTMs) in proteoforms from bovine milk (PXD031744) and human ovarian tissue. Contemporary identification workflows produce excellent PrSM yields, although approximately half of all identified proteoforms from these four pipelines were specific to only one workflow. Deconvolution algorithms disagree on precursor masses and charges, contributing to identification variability. Detection of PTMs is inconsistent among algorithms. In bovine milk, 18% of PrSMs produced by pTop and TopMG were singly phosphorylated, but this percentage fell to 1% for one algorithm. Applying multiple search engines produces more comprehensive assessments of experiments. Top-down algorithms would benefit from greater interoperability.

Proteomics Standards Initiative at Twenty Years: Current Activities and Future Work

The Human Proteome Organization (HUPO) Proteomics Standards Initiative (PSI) has been successfully developing guidelines, data formats, and controlled vocabularies (CVs) for the proteomics community and other fields supported by mass spectrometry since its inception 20 years ago. Here we describe the general operation of the PSI, including its leadership, working groups, yearly workshops, and the document process by which proposals are thoroughly and publicly reviewed in order to be ratified as PSI standards. We briefly describe the current state of the many existing PSI standards, some of which remain the same as when originally developed, some of which have undergone subsequent revisions, and some of which have become obsolete. Then the set of proposals currently being developed are described, with an open call to the community for participation in the forging of the next generation of standards. Finally, we describe some synergies and collaborations with other organizations and look to the future in how the PSI will continue to promote the open sharing of data and thus accelerate the progress of the field of proteomics.

Comparative Proteomics Analysis between Maize and Sorghum Uncovers Important Proteins and Metabolic Pathways Mediating Drought Tolerance

Drought severely affects crop yield and yield stability. Maize and sorghum are major crops in Africa and globally, and both are negatively impacted by drought. However, sorghum has a better ability to withstand drought than maize. Consequently, this study identifies differences between maize and sorghum grown in water deficit conditions, and identifies proteins associated with drought tolerance in these plant species. Leaf relative water content and proline content were measured, and label-free proteomics analysis was carried out to identify differences in protein expression in the two species in response to water deficit. Water deficit enhanced the proline accumulation in sorghum roots to a higher degree than in maize, and this higher accumulation was associated with enhanced water retention in sorghum. Proteomic analyses identified proteins with differing expression patterns between the two species, revealing key metabolic pathways that explain the better drought tolerance of sorghum than maize. These proteins include phenylalanine/tyrosine ammonia-lyases, indole-3-acetaldehyde oxidase, sucrose synthase and phenol/catechol oxidase. This study highlights the importance of phenylpropanoids, sucrose, melanin-related metabolites and indole acetic acid (auxin) as determinants of the differences in drought stress tolerance between maize and sorghum. The selection of maize and sorghum genotypes with enhanced expression of the genes encoding these differentially expressed proteins, or genetically engineering maize and sorghum to increase the expression of such genes, can be used as strategies for the production of maize and sorghum varieties with improved drought tolerance.

Expansion of nickel binding- and histidine-rich proteins during gastric adaptation of Helicobacter species

Acquisition and homeostasis of essential metals during host colonization by bacterial pathogens rely on metal uptake, trafficking and storage proteins. How these factors have evolved within bacterial pathogens is poorly defined. Urease, a nickel enzyme, is essential for Helicobacter pylori to colonize the acidic stomach. Our previous data suggest that acquisition of nickel transporters and a Histidine-rich protein (HRP) involved in nickel storage in H. pylori and gastric Helicobacter spp. have been essential evolutionary events for gastric colonization. Using bioinformatics, proteomics and phylogenetics, we extended this analysis to determine how evolution has framed the repertoire of HRPs among 39 Epsilonproteobacteria; 18 gastric and 11 non-gastric enterohepatic (EH) Helicobacter spp., as well as 10 other Epsilonproteobacteria. We identified a total of 213 HRPs distributed in 22 protein families named orthologous groups (OG) with His-rich domains, including 15 newly described OGs. Gastric Helicobacter spp. are enriched in HRPs (7.7 ± 1.9 HRPs/strain) as compared to EH Helicobacter spp. (1.9 ± 1.0 HRPs/strain) with a particular prevalence of HRPs with C-terminal Histidine-rich domains in gastric species. The expression and nickel-binding capacity of several HRPs was validated in five gastric Helicobacter spp. We established the evolutionary history of new HRP families, such as the periplasmic HP0721-like proteins and the HugZ-type heme-oxygenases. The expansion of Histidine-rich extensions in gastric Helicobacter spp. proteins is intriguing but can tentatively be associated with the presence of the urease nickel-enzyme. We conclude that this HRP expansion is associated with unique properties of organisms that rely on large intracellular nickel amounts for their survival.

Proteomic Identification and Meta-Analysis in Salvia hispanica RNA-Seq de novo Assemblies

While proteomics has demonstrated its value for model organisms and for organisms with mature genome sequence annotations, proteomics has been of less value in nonmodel organisms that are unaccompanied by genome sequence annotations. This project sought to determine the value of RNA-Seq experiments as a basis for establishing a set of protein sequences to represent a nonmodel organism, in this case, the pseudocereal chia. Assembling four publicly available chia RNA-Seq datasets produced transcript sequence sets with a high BUSCO completeness, though the number of transcript sequences and Trinity "genes" varied considerably among them. After six-frame translation, ProteinOrtho detected substantial numbers of orthologs among other species within the taxonomic order Lamiales. These protein sequence databases demonstrated a good identification efficiency for three different LC-MS/MS proteomics experiments, though a seed proteome showed considerable variability in the identification of peptides based on seed protein sequence inclusion. If a proteomics experiment emphasizes a particular tissue, an RNA-Seq experiment incorporating that same tissue is more likely to support a database search identification of that proteome.

Dysregulated healing responses in diabetic wounds occur in the early stages postinjury

Chronic wounds are a serious and debilitating complication of diabetes. A better understanding of the dysregulated healing responses following injury will provide insight into the optimal time frame for therapeutic intervention. In this study, a direct comparison was done between the healing dynamics and the proteome of acute and obese diabetic wounds on days 2 and 7 following injury. Full thickness excisional wounds were induced on obese diabetic (B6.Cg-lepob/J, ob/ob, n = 14) (blood glucose 423.25 ± 127.92 mg/dL) and WT control (C57BL/6J, n = 14) (blood glucose 186.67 ± 24.5 mg/dL) mice. Histological analysis showed no signs of healing in obese DM wounds whereas complete wound closure/re-epithelisation, the formation of granulation tissue and signs of re-vascularisation, was evident in acute wounds on day 7. In obese DM wounds, substance P deficiency and increased MMP-9 activity on day 2 coincided with increased cytokine/chemokine levels within wound fluid. LC-MS/MS identified 906 proteins, of which 23 (Actn3, Itga6, Epb41, Sncg, Nefm, Rsp18, Rsp19, Rpl22, Macroh2a1, Rpn1, Ppib, Snrnp70, Ddx5, Eif3g, Tpt1, FABP5, Cavin1, Stfa1, Stfa3, Cycs, Tkt, Mb, Chmp2a) were differentially expressed in wounded tissue on day 2 (P < 0.05; more than two-fold) and 6 (Cfd, Ptms, Hp, Hmga1, Cbx3, Syap1) (P < 0.05; more than two-fold) on day 7. A large number of dysregulated proteins on day 2 was associated with an inability to progress into the proliferative stage of healing and suggest that early intervention might be pivotal for effective healing outcomes. The proteomic approach highlighted the complexity of obese DM wounds in which the dysregulation involves multiple regulatory pathways and biological processes.

Microbial function and genital inflammation in young South African women at high risk of HIV infection

BACKGROUND

Female genital tract (FGT) inflammation is an important risk factor for HIV acquisition. The FGT microbiome is closely associated with inflammatory profile; however, the relative importance of microbial activities has not been established. Since proteins are key elements representing actual microbial functions, this study utilized metaproteomics to evaluate the relationship between FGT microbial function and inflammation in 113 young and adolescent South African women at high risk of HIV infection. Women were grouped as having low, medium, or high FGT inflammation by K-means clustering according to pro-inflammatory cytokine concentrations.

RESULTS

A total of 3186 microbial and human proteins were identified in lateral vaginal wall swabs using liquid chromatography-tandem mass spectrometry, while 94 microbial taxa were included in the taxonomic analysis. Both metaproteomics and 16S rRNA gene sequencing analyses showed increased non-optimal bacteria and decreased lactobacilli in women with FGT inflammatory profiles. However, differences in the predicted relative abundance of most bacteria were observed between 16S rRNA gene sequencing and metaproteomics analyses. Bacterial protein functional annotations (gene ontology) predicted inflammatory cytokine profiles more accurately than bacterial relative abundance determined by 16S rRNA gene sequence analysis, as well as functional predictions based on 16S rRNA gene sequence data (p < 0.0001). The majority of microbial biological processes were underrepresented in women with high inflammation compared to those with low inflammation, including a Lactobacillus-associated signature of reduced cell wall organization and peptidoglycan biosynthesis. This signature remained associated with high FGT inflammation in a subset of 74 women 9 weeks later, was upheld after adjusting for Lactobacillus relative abundance, and was associated with in vitro inflammatory cytokine responses to Lactobacillus isolates from the same women. Reduced cell wall organization and peptidoglycan biosynthesis were also associated with high FGT inflammation in an independent sample of ten women.

CONCLUSIONS

Both the presence of specific microbial taxa in the FGT and their properties and activities are critical determinants of FGT inflammation. Our findings support those of previous studies suggesting that peptidoglycan is directly immunosuppressive, and identify a possible avenue for biotherapeutic development to reduce inflammation in the FGT. To facilitate further investigations of microbial activities, we have developed the FGT-DB application that is available at http://fgtdb.org/ . Video Abstract.

Interrogating Fractionation and Other Sources of Variability in Shotgun Proteomes Using Quality Metrics

The increasing amount of publicly available proteomics data creates opportunities for data scientists to investigate quality metrics in novel ways. QuaMeter IDFree is used to generate quality metrics from 665 RAW files and 97 WIFF files representing publicly available "shotgun" mass spectrometry datasets. These experiments are selected to represent Mycobacterium tuberculosis lysates, mouse MDSCs, and exosomes derived from human cell lines. Machine learning techniques are demonstrated to detect outliers within experiments and it is shown that quality metrics may be used to distinguish sources of variability among these experiments. In particular, the findings demonstrate that according to nested ANOVA performed on an SDS-PAGE shotgun principal component analysis, runs of fractions from the same gel regions cluster together rather than technical replicates, close temporal proximity, or even biological samples. This indicates that the individual fraction may have had a higher impact on the quality metrics than other factors. In addition, sample type, instrument type, mass analyzer, fragmentation technique, and digestion enzyme are identified as sources of variability. From a quality control perspective, the importance of study design and in particular, the run order, is illustrated in seeking ways to limit the impact of technical variability.

Targeted next-generation sequencing identifies novel variants in candidate genes for Parkinson's disease in Black South African and Nigerian patients

BACKGROUND

The prevalence of Parkinson's disease (PD) is increasing in sub-Saharan Africa, but little is known about the genetics of PD in these populations. Due to their unique ancestry and diversity, sub-Saharan African populations have the potential to reveal novel insights into the pathobiology of PD. In this study, we aimed to characterise the genetic variation in known and novel PD genes in a group of Black South African and Nigerian patients.

METHODS

We recruited 33 Black South African and 14 Nigerian PD patients, and screened them for sequence variants in 751 genes using an Ion AmpliSeq™ Neurological Research panel. We used bcftools to filter variants and annovar software for the annotation. Rare variants were prioritised using MetaLR and MetaSVM prediction scores. The effect of a variant on ATP13A2's protein structure was investigated by molecular modelling.

RESULTS

We identified 14,655 rare variants with a minor allele frequency ≤ 0.01, which included 2448 missense variants. Notably, no common pathogenic mutations were identified in these patients. Also, none of the known PD-associated mutations were found highlighting the need for more studies in African populations. Altogether, 54 rare variants in 42 genes were considered deleterious and were prioritized, based on MetaLR and MetaSVM scores, for follow-up studies. Protein modelling showed that the S1004R variant in ATP13A2 possibly alters the conformation of the protein.

CONCLUSIONS

We identified several rare variants predicted to be deleterious in sub-Saharan Africa PD patients; however, further studies are required to determine the biological effects of these variants and their possible role in PD. Studies such as these are important to elucidate the genetic aetiology of this disorder in patients of African ancestry.

SWATH-MS based proteomic profiling of pancreatic ductal adenocarcinoma tumours reveals the interplay between the extracellular matrix and related intracellular pathways

Pancreatic cancer accounts for 2.8% of new cancer cases worldwide and is projected to become the second leading cause of cancer-related deaths by 2030. Patients of African ancestry appear to be at an increased risk for pancreatic ductal adenocarcinoma (PDAC), with more severe disease and outcomes. The purpose of this study was to map the proteomic and genomic landscape of a cohort of PDAC patients of African ancestry. Thirty tissues (15 tumours and 15 normal adjacent tissues) were obtained from consenting South African PDAC patients. Optimisation of the sample preparation method allowed for the simultaneous extraction of high-purity protein and DNA for SWATH-MS and OncoArray SNV analyses. We quantified 3402 proteins with 49 upregulated and 35 downregulated proteins at a minimum 2.1 fold change and FDR adjusted p-value (q-value) ≤ 0.01 when comparing tumour to normal adjacent tissue. Many of the upregulated proteins in the tumour samples are involved in extracellular matrix formation (ECM) and related intracellular pathways. In addition, proteins such as EMIL1, KBTB2, and ZCCHV involved in the regulation of ECM proteins were observed to be dysregulated in pancreatic tumours. Downregulation of pathways involved in oxygen and carbon dioxide transport were observed. Genotype data showed missense mutations in some upregulated proteins, such as MYPN, ESTY2 and SERPINB8. Approximately 11% of the dysregulated proteins, including ISLR, BP1, PTK7 and OLFL3, were predicted to be secretory proteins. These findings help in further elucidating the biology of PDAC and may aid in identifying future plausible markers for the disease.

Immunometabolic Signatures Predict Risk of Progression to Active Tuberculosis and Disease Outcome

There remains a pressing need for biomarkers that can predict who will progress to active tuberculosis (TB) after exposure to Mycobacterium tuberculosis (MTB) bacterium. By analyzing cohorts of household contacts of TB index cases (HHCs) and a stringent non-human primate (NHP) challenge model, we evaluated whether integration of blood transcriptional profiling with serum metabolomic profiling can provide new understanding of disease processes and enable improved prediction of TB progression. Compared to either alone, the combined application of pre-existing transcriptome- and metabolome-based signatures more accurately predicted TB progression in the HHC cohorts and more accurately predicted disease severity in the NHPs. Pathway and data-driven correlation analyses of the integrated transcriptional and metabolomic datasets further identified novel immunometabolomic signatures significantly associated with TB progression in HHCs and NHPs, implicating cortisol, tryptophan, glutathione, and tRNA acylation networks. These results demonstrate the power of multi-omics analysis to provide new insights into complex disease processes.

Sipros Ensemble improves database searching and filtering for complex metaproteomics

Motivation

Complex microbial communities can be characterized by metagenomics and metaproteomics. However, metagenome assemblies often generate enormous, and yet incomplete, protein databases, which undermines the identification of peptides and proteins in metaproteomics. This challenge calls for increased discrimination of true identifications from false identifications by database searching and filtering algorithms in metaproteomics.

Results

Sipros Ensemble was developed here for metaproteomics using an ensemble approach. Three diverse scoring functions from MyriMatch, Comet and the original Sipros were incorporated within a single database searching engine. Supervised classification with logistic regression was used to filter database searching results. Benchmarking with soil and marine microbial communities demonstrated a higher number of peptide and protein identifications by Sipros Ensemble than MyriMatch/Percolator, Comet/Percolator, MS-GF+/Percolator, Comet & MyriMatch/iProphet and Comet & MyriMatch & MS-GF+/iProphet. Sipros Ensemble was computationally efficient and scalable on supercomputers.

Availability and implementation

Freely available under the GNU GPL license at http://sipros.omicsbio.org.

Supplementary information

Supplementary data are available at Bioinformatics online.

Quality control in mass spectrometry-based proteomics

Mass spectrometry is a highly complex analytical technique and mass spectrometry-based proteomics experiments can be subject to a large variability, which forms an obstacle to obtaining accurate and reproducible results. Therefore, a comprehensive and systematic approach to quality control is an essential requirement to inspire confidence in the generated results. A typical mass spectrometry experiment consists of multiple different phases including the sample preparation, liquid chromatography, mass spectrometry, and bioinformatics stages. We review potential sources of variability that can impact the results of a mass spectrometry experiment occurring in all of these steps, and we discuss how to monitor and remedy the negative influences on the experimental results. Furthermore, we describe how specialized quality control samples of varying sample complexity can be incorporated into the experimental workflow and how they can be used to rigorously assess detailed aspects of the instrument performance.

A Serum Circulating miRNA Signature for Short-Term Risk of Progression to Active Tuberculosis Among Household Contacts

Biomarkers that predict who among recently Mycobacterium tuberculosis (MTB)-exposed individuals will progress to active tuberculosis are urgently needed. Intracellular microRNAs (miRNAs) regulate the host response to MTB and circulating miRNAs (c-miRNAs) have been developed as biomarkers for other diseases. We performed machine-learning analysis of c-miRNA measurements in the serum of adult household contacts (HHCs) of TB index cases from South Africa and Uganda and developed a c-miRNA-based signature of risk for progression to active TB. This c-miRNA-based signature significantly discriminated HHCs within 6 months of progression to active disease from HHCs that remained healthy in an independent test set [ROC area under the ROC curve (AUC) 0.74, progressors < 6 Mo to active TB and ROC AUC 0.66, up to 24 Mo to active TB], and complements the predictions of a previous cellular mRNA-based signature of TB risk.

Challenges and Opportunities for Biological Mass Spectrometry Core Facilities in the Developing World

The developing world is seeing rapid growth in the availability of biological mass spectrometry (MS), particularly through core facilities. As proteomics and metabolomics becomes locally feasible for investigators in these nations, application areas associated with high burden in these nations, such as infectious disease, will see greatly increased research output. This article evaluates the rapid growth of MS in South Africa (currently approaching 20 laboratories) as a model for establishing MS core facilities in other nations of the developing world. Facilities should emphasize new services rather than new instruments. The reduction of the delays associated with reagent and other supply acquisition would benefit both facilities and the users who make use of their services. Instrument maintenance and repair, often mediated by an in-country business for an international vendor, is also likely to operate on a slower schedule than in the wealthiest nations. A key challenge to facilities in the developing world is educating potential facility users in how best to design experiments for proteomics and metabolomics, what reagents are most likely to introduce problematic artifacts, and how to interpret results from the facility. Here, we summarize the experience of 6 different institutions to raise the level of biological MS available to researchers in South Africa.

Characterization of Plant Glycoproteins: Analysis of Plant Glycopeptide Mass Spectrometry Data with plantGlycoMS, a Package in the R Statistical Computing Environment

plantGlycoMS is a set of tools, implemented in R, which is used to assess and validate glycopeptide spectrum matches (gPSMs). Validity of gPSMs is based on characteristic fragmentation patterns of glycopeptides (gPSMvalidator), adherence of the glycan moiety to the known N-glycan biosynthesis pathway in plants (pGlycoFilter), and elution of the glycopeptide within the observed retention time window of other glycopeptides sharing the same peptide backbone (rt.Restrict). plantGlycoMS also contains a tool for relative quantitation of glycoforms based on selected ion chromatograms of glycopeptide ion precursors in the mass spectrometry level 1 data (glycoRQ). This protocol walks the user through this workflow with example mass spectrometry data obtained for a plant glycoprotein.

Quality Assessments of Long-Term Quantitative Proteomic Analysis of Breast Cancer Xenograft Tissues

Clinical proteomics requires large-scale analysis of human specimens to achieve statistical significance. We evaluated the long-term reproducibility of an iTRAQ (isobaric tags for relative and absolute quantification)-based quantitative proteomics strategy using one channel for reference across all samples in different iTRAQ sets. A total of 148 liquid chromatography tandem mass spectrometric (LC-MS/MS) analyses were completed, generating six 2D LC-MS/MS data sets for human-in-mouse breast cancer xenograft tissues representative of basal and luminal subtypes. Such large-scale studies require the implementation of robust metrics to assess the contributions of technical and biological variability in the qualitative and quantitative data. Accordingly, we derived a quantification confidence score based on the quality of each peptide-spectrum match to remove quantification outliers from each analysis. After combining confidence score filtering and statistical analysis, reproducible protein identification and quantitative results were achieved from LC-MS/MS data sets collected over a 7-month period. This study provides the first quality assessment on long-term stability and technical considerations for study design of a large-scale clinical proteomics project.

Proteomics Standards Initiative: Fifteen Years of Progress and Future Work

The Proteomics Standards Initiative (PSI) of the Human Proteome Organization (HUPO) has now been developing and promoting open community standards and software tools in the field of proteomics for 15 years. Under the guidance of the chair, cochairs, and other leadership positions, the PSI working groups are tasked with the development and maintenance of community standards via special workshops and ongoing work. Among the existing ratified standards, the PSI working groups continue to update PSI-MI XML, MITAB, mzML, mzIdentML, mzQuantML, mzTab, and the MIAPE (Minimum Information About a Proteomics Experiment) guidelines with the advance of new technologies and techniques. Furthermore, new standards are currently either in the final stages of completion (proBed and proBAM for proteogenomics results as well as PEFF) or in early stages of design (a spectral library standard format, a universal spectrum identifier, the qcML quality control format, and the Protein Expression Interface (PROXI) web services Application Programming Interface). In this work we review the current status of all of these aspects of the PSI, describe synergies with other efforts such as the ProteomeXchange Consortium, the Human Proteome Project, and the metabolomics community, and provide a look at future directions of the PSI.

Proteogenomic Investigation of Strain Variation in Clinical Mycobacterium tuberculosis Isolates

Mycobacterium tuberculosis consists of a large number of different strains that display unique virulence characteristics. Whole-genome sequencing has revealed substantial genetic diversity among clinical M. tuberculosis isolates, and elucidating the phenotypic variation encoded by this genetic diversity will be of the utmost importance to fully understand M. tuberculosis biology and pathogenicity. In this study, we integrated whole-genome sequencing and mass spectrometry (GeLC-MS/MS) to reveal strain-specific characteristics in the proteomes of two clinical M. tuberculosis Latin American-Mediterranean isolates. Using this approach, we identified 59 peptides containing single amino acid variants, which covered ∼9% of all coding nonsynonymous single nucleotide variants detected by whole-genome sequencing. Furthermore, we identified 29 distinct peptides that mapped to a hypothetical protein not present in the M. tuberculosis H37Rv reference proteome. Here, we provide evidence for the expression of this protein in the clinical M. tuberculosis SAWC3651 isolate. The strain-specific databases enabled confirmation of genomic differences (i.e., large genomic regions of difference and nonsynonymous single nucleotide variants) in these two clinical M. tuberculosis isolates and allowed strain differentiation at the proteome level. Our results contribute to the growing field of clinical microbial proteogenomics and can improve our understanding of phenotypic variation in clinical M. tuberculosis isolates.

Host blood RNA signatures predict the outcome of tuberculosis treatment

Biomarkers for tuberculosis treatment outcome will assist in guiding individualized treatment and evaluation of new therapies. To identify candidate biomarkers, RNA sequencing of whole blood from a well-characterized TB treatment cohort was performed. Application of a validated transcriptional correlate of risk for TB revealed symmetry in host gene expression during progression from latent TB infection to active TB disease and resolution of disease during treatment, including return to control levels after drug therapy. The symmetry was also seen in a TB disease signature, constructed from the TB treatment cohort, that also functioned as a strong correlate of risk. Both signatures identified patients at risk of treatment failure 1–4 weeks after start of therapy. Further mining of the transcriptomes revealed an association between treatment failure and suppressed expression of mitochondrial genes before treatment initiation, leading to development of a novel baseline (pre-treatment) signature of treatment failure. These novel host responses to TB treatment were integrated into a five-gene real-time PCR-based signature that captures the clinically relevant responses to TB treatment and provides a convenient platform for stratifying patients according to their risk of treatment failure. Furthermore, this 5-gene signature is shown to correlate with the pulmonary inflammatory state (as measured by PET-CT) and can complement sputum-based Gene Xpert for patient stratification, providing a rapid and accurate alternative to current methods.

The Human Proteome Organization-Proteomics Standards Initiative Quality Control Working Group: Making Quality Control More Accessible for Biological Mass Spectrometry

To have confidence in results acquired during biological mass spectrometry experiments, a systematic approach to quality control is of vital importance. Nonetheless, until now, only scattered initiatives have been undertaken to this end, and these individual efforts have often not been complementary. To address this issue, the Human Proteome Organization-Proteomics Standards Initiative has established a new working group on quality control at its meeting in the spring of 2016. The goal of this working group is to provide a unifying framework for quality control data. The initial focus will be on providing a community-driven standardized file format for quality control. For this purpose, the previously proposed qcML format will be adapted to support a variety of use cases for both proteomics and metabolomics applications, and it will be established as an official PSI format. An important consideration is to avoid enforcing restrictive requirements on quality control but instead provide the basic technical necessities required to support extensive quality control for any type of mass spectrometry-based workflow. We want to emphasize that this is an open community effort, and we seek participation from all scientists with an interest in this field.

Chronic intermittent alcohol disrupts the GluN2B-associated proteome and specifically regulates group I mGlu receptor-dependent long-term depression

N-Methyl-d-aspartate receptors (NMDARs) are major targets of both acute and chronic alcohol, as well as regulators of plasticity in a number of brain regions. Aberrant plasticity may contribute to the treatment resistance and high relapse rates observed in alcoholics. Recent work suggests that chronic alcohol treatment preferentially modulates both the expression and subcellular localization of NMDARs containing the GluN2B subunit. Signaling through synaptic and extrasynaptic GluN2B-NMDARs has already been implicated in the pathophysiology of various other neurological disorders. NMDARs interact with a large number of proteins at the glutamate synapse, and a better understanding of how alcohol modulates this proteome is needed. We employed a discovery-based proteomic approach in subcellular fractions of hippocampal tissue from chronic intermittent alcohol (CIE)-exposed C57Bl/6J mice to gain insight into alcohol-induced changes in GluN2B signaling complexes. Protein enrichment analyses revealed changes in the association of post-synaptic proteins, including scaffolding, glutamate receptor and PDZ-domain binding proteins with GluN2B. In particular, GluN2B interaction with metabotropic glutamate (mGlu)1/5 receptor-dependent long-term depression (LTD)-associated proteins such as Arc and Homer 1 was increased, while GluA2 was decreased. Accordingly, we found a lack of mGlu1/5 -induced LTD while α1 -adrenergic receptor-induced LTD remained intact in hippocampal CA1 following CIE. These data suggest that CIE specifically disrupts mGlu1/5 -LTD, representing a possible connection between NMDAR and mGlu receptor signaling. These studies not only demonstrate a new way in which alcohol can modulate plasticity in the hippocampus but also emphasize the utility of this discovery-based proteomic approach to generate new hypotheses regarding alcohol-related mechanisms.

Recognizing millions of consistently unidentified spectra across hundreds of shotgun proteomics datasets

Mass spectrometry (MS) is the main technology used in proteomics approaches. However, on average 75% of spectra analysed in an MS experiment remain unidentified. We propose to use spectrum clustering at a large-scale to shed a light on these unidentified spectra. PRoteomics IDEntifications database (PRIDE) Archive is one of the largest MS proteomics public data repositories worldwide. By clustering all tandem MS spectra publicly available in PRIDE Archive, coming from hundreds of datasets, we were able to consistently characterize three distinct groups of spectra: 1) incorrectly identified spectra, 2) spectra correctly identified but below the set scoring threshold, and 3) truly unidentified spectra. Using a multitude of complementary analysis approaches, we were able to identify less than 20% of the consistently unidentified spectra. The complete spectrum clustering results are available through the new version of the PRIDE Cluster resource (http://www.ebi.ac.uk/pride/cluster). This resource is intended, among other aims, to encourage and simplify further investigation into these unidentified spectra.

Greazy: Open-Source Software for Automated Phospholipid Tandem Mass Spectrometry Identification

Lipid identification from data produced with high-throughput technologies is essential to the elucidation of the roles played by lipids in cellular function and disease. Software tools for identifying lipids from tandem mass (MS/MS) spectra have been developed, but they are often costly or lack the sophistication of their proteomics counterparts. We have developed Greazy, an open source tool for the automated identification of phospholipids from MS/MS spectra, that utilizes methods similar to those developed for proteomics. From user-supplied parameters, Greazy builds a phospholipid search space and associated theoretical MS/MS spectra. Experimental spectra are scored against search space lipids with similar precursor masses using a peak score based on the hypergeometric distribution and an intensity score utilizing the percentage of total ion intensity residing in matching peaks. The LipidLama component filters the results via mixture modeling and density estimation. We assess Greazy's performance against the NIST 2014 metabolomics library, observing high accuracy in a search of multiple lipid classes. We compare Greazy/LipidLama against the commercial lipid identification software LipidSearch and show that the two platforms differ considerably in the sets of identified spectra while showing good agreement on those spectra identified by both. Lastly, we demonstrate the utility of Greazy/LipidLama with different instruments. We searched data from replicates of alveolar type 2 epithelial cells obtained with an Orbitrap and from human serum replicates generated on a quadrupole-time-of-flight (Q-TOF). These findings substantiate the application of proteomics derived methods to the identification of lipids. The software is available from the ProteoWizard repository: http://tiny.cc/bumbershoot-vc12-bin64 .

Reproducibility of Differential Proteomic Technologies in CPTAC Fractionated Xenografts

The NCI Clinical Proteomic Tumor Analysis Consortium (CPTAC) employed a pair of reference xenograft proteomes for initial platform validation and ongoing quality control of its data collection for The Cancer Genome Atlas (TCGA) tumors. These two xenografts, representing basal and luminal-B human breast cancer, were fractionated and analyzed on six mass spectrometers in a total of 46 replicates divided between iTRAQ and label-free technologies, spanning a total of 1095 LC-MS/MS experiments. These data represent a unique opportunity to evaluate the stability of proteomic differentiation by mass spectrometry over many months of time for individual instruments or across instruments running dissimilar workflows. We evaluated iTRAQ reporter ions, label-free spectral counts, and label-free extracted ion chromatograms as strategies for data interpretation (source code is available from http://homepages.uc.edu/~wang2x7/Research.htm ). From these assessments, we found that differential genes from a single replicate were confirmed by other replicates on the same instrument from 61 to 93% of the time. When comparing across different instruments and quantitative technologies, using multiple replicates, differential genes were reproduced by other data sets from 67 to 99% of the time. Projecting gene differences to biological pathways and networks increased the degree of similarity. These overlaps send an encouraging message about the maturity of technologies for proteomic differentiation.

proBAMsuite, a Bioinformatics Framework for Genome-Based Representation and Analysis of Proteomics Data

To facilitate genome-based representation and analysis of proteomics data, we developed a new bioinformatics framework, proBAMsuite, in which a central component is the protein BAM (proBAM) file format for organizing peptide spectrum matches (PSMs)¹ within the context of the genome. proBAMsuite also includes two R packages, proBAMr and proBAMtools, for generating and analyzing proBAM files, respectively. Applying proBAMsuite to three recently published proteomics datasets, we demonstrated its utility in facilitating efficient genome-based sharing, interpretation, and integration of proteomics data. First, the interpretation of proteomics data is significantly enhanced with the rich genomic annotation information. Second, PSMs can be easily reannotated using user-specified gene annotation schemes and assembled into both protein and gene identifications. Third, using the genome as a common reference, proBAMsuite facilitates seamless proteomics and proteogenomics data integration. Finally, proBAM files can be readily visualized in genome browsers and thus bring proteomics data analysis to a general audience beyond the proteomics community. Results from this study establish proBAMsuite as a useful bioinformatics framework for proteomics and proteogenomics research.

Comprehensive Characterization of Glycosylation and Hydroxylation of Basement Membrane Collagen IV by High-Resolution Mass Spectrometry

Collagen IV is the main structural protein that provides a scaffold for assembly of basement membrane proteins. Posttranslational modifications such as hydroxylation of proline and lysine and glycosylation of lysine are essential for the functioning of collagen IV triple-helical molecules. These modifications are highly abundant posing a difficult challenge for in-depth characterization of collagen IV using conventional proteomics approaches. Herein, we implemented an integrated pipeline combining high-resolution mass spectrometry with different fragmentation techniques and an optimized bioinformatics workflow to study posttranslational modifications in mouse collagen IV. We achieved 82% sequence coverage for the α1 chain, mapping 39 glycosylated hydroxylysine, 148 4-hydroxyproline, and seven 3-hydroxyproline residues. Further, we employed our pipeline to map the modifications on human collagen IV and achieved 85% sequence coverage for the α1 chain, mapping 35 glycosylated hydroxylysine, 163 4-hydroxyproline, and 14 3-hydroxyproline residues. Although lysine glycosylation heterogeneity was observed in both mouse and human, 21 conserved sites were identified. Likewise, five 3-hydroxyproline residues were conserved between mouse and human, suggesting that these modification sites are important for collagen IV function. Collectively, these are the first comprehensive maps of hydroxylation and glycosylation sites in collagen IV, which lay the foundation for dissecting the key role of these modifications in health and disease.

An Analysis of the Sensitivity of Proteogenomic Mapping of Somatic Mutations and Novel Splicing Events in Cancer

Improvements in mass spectrometry (MS)-based peptide sequencing provide a new opportunity to determine whether polymorphisms, mutations, and splice variants identified in cancer cells are translated. Herein, we apply a proteogenomic data integration tool (QUILTS) to illustrate protein variant discovery using whole genome, whole transcriptome, and global proteome datasets generated from a pair of luminal and basal-like breast-cancer-patient-derived xenografts (PDX). The sensitivity of proteogenomic analysis for singe nucleotide variant (SNV) expression and novel splice junction (NSJ) detection was probed using multiple MS/MS sample process replicates defined here as an independent tandem MS experiment using identical sample material. Despite analysis of over 30 sample process replicates, only about 10% of SNVs (somatic and germline) detected by both DNA and RNA sequencing were observed as peptides. An even smaller proportion of peptides corresponding to NSJ observed by RNA sequencing were detected (<0.1%). Peptides mapping to DNA-detected SNVs without a detectable mRNA transcript were also observed, suggesting that transcriptome coverage was incomplete (∼80%). In contrast to germline variants, somatic variants were less likely to be detected at the peptide level in the basal-like tumor than in the luminal tumor, raising the possibility of differential translation or protein degradation effects. In conclusion, this large-scale proteogenomic integration allowed us to determine the degree to which mutations are translated and identify gaps in sequence coverage, thereby benchmarking current technology and progress toward whole cancer proteome and transcriptome analysis.

The SEQUEST family tree

Since its introduction in 1994, SEQUEST has gained many important new capabilities, and a host of successor algorithms have built upon its successes. This Account and Perspective maps the evolution of this important tool and charts the relationships among contributions to the SEQUEST legacy. Many of the changes represented improvements in computing speed by clusters and graphics cards. Mass spectrometry innovations in mass accuracy and activation methods led to shifts in fragment modeling and scoring strategies. These changes, as well as the movement of laboratories and lab members, have led to great diversity among the members of the SEQUEST family. Graphical Abstract ᅟ.

The 2012/2013 ABRF Proteomic Research Group Study: Assessing Longitudinal Intralaboratory Variability in Routine Peptide Liquid Chromatography Tandem Mass Spectrometry Analyses

Questions concerning longitudinal data quality and reproducibility of proteomic laboratories spurred the Protein Research Group of the Association of Biomolecular Resource Facilities (ABRF-PRG) to design a study to systematically assess the reproducibility of proteomic laboratories over an extended period of time. Developed as an open study, initially 64 participants were recruited from the broader mass spectrometry community to analyze provided aliquots of a six bovine protein tryptic digest mixture every month for a period of nine months. Data were uploaded to a central repository, and the operators answered an accompanying survey. Ultimately, 45 laboratories submitted a minimum of eight LC-MSMS raw data files collected in data-dependent acquisition (DDA) mode. No standard operating procedures were enforced; rather the participants were encouraged to analyze the samples according to usual practices in the laboratory. Unlike previous studies, this investigation was not designed to compare laboratories or instrument configuration, but rather to assess the temporal intralaboratory reproducibility. The outcome of the study was reassuring with 80% of the participating laboratories performing analyses at a medium to high level of reproducibility and quality over the 9-month period. For the groups that had one or more outlying experiments, the major contributing factor that correlated to the survey data was the performance of preventative maintenance prior to the LC-MSMS analyses. Thus, the Protein Research Group of the Association of Biomolecular Resource Facilities recommends that laboratories closely scrutinize the quality control data following such events. Additionally, improved quality control recording is imperative. This longitudinal study provides evidence that mass spectrometry-based proteomics is reproducible. When quality control measures are strictly adhered to, such reproducibility is comparable among many disparate groups. Data from the study are available via ProteomeXchange under the accession code PXD002114.

Correcting systematic bias and instrument measurement drift with mzRefinery

MOTIVATION

Systematic bias in mass measurement adversely affects data quality and negates the advantages of high precision instruments.

RESULTS

We introduce the mzRefinery tool for calibration of mass spectrometry data files. Using confident peptide spectrum matches, three different calibration methods are explored and the optimal transform function is chosen. After calibration, systematic bias is removed and the mass measurement errors are centered at 0 ppm. Because it is part of the ProteoWizard package, mzRefinery can read and write a wide variety of file formats.

AVAILABILITY AND IMPLEMENTATION

The mzRefinery tool is part of msConvert, available with the ProteoWizard open source package at http://proteowizard.sourceforge.net/

CONTACT

samuel.payne@pnnl.gov.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Advanced Proteogenomic Analysis Reveals Multiple Peptide Mutations and Complex Immunoglobulin Peptides in Colon Cancer

Aiming toward an improved understanding of the regulation of proteins in cancer, recent studies from the Clinical Proteomic Tumor Analysis Consortium (CPTAC) have focused on analyzing cancer tissue using proteomic technologies and workflows. Although many proteogenomics approaches for the study of cancer samples have been proposed, serious methodological challenges remain, especially in the identification of multiple mutational variants or structural variations such as fusion gene events. In addition, although immune system genes play an important role in cancer, identification of IgG peptides remains challenging in proteomic data sets. Here, we describe an integrative proteogenomic method that extends the limit of proteogenomic searches to identify multiple variant peptides as well as immunoglobulin gene variations/rearrangements using customized mining of RNA-seq data. Our results also provide the first extensive characterization of tumor immune response and demonstrate the potential of this method to improve the molecular characterization of tumor subtypes.

Proteomic analysis of colon and rectal carcinoma using standard and customized databases

Understanding proteomic differences underlying the different phenotypic classes of colon and rectal carcinoma is important and may eventually lead to a better assessment of clinical behavior of these cancers. We here present a comprehensive description of the proteomic data obtained from 90 colon and rectal carcinomas previously subjected to genomic analysis by The Cancer Genome Atlas (TCGA). Here, the primary instrument files and derived secondary data files are compiled and presented in forms that will allow further analyses of the biology of colon and rectal carcinoma. We also discuss new challenges in processing these large proteomic datasets for relevant proteins and protein variants.

Wavelet-based peak detection and a new charge inference procedure for MS/MS implemented in ProteoWizard's msConvert

We report the implementation of high-quality signal processing algorithms into ProteoWizard, an efficient, open-source software package designed for analyzing proteomics tandem mass spectrometry data. Specifically, a new wavelet-based peak-picker (CantWaiT) and a precursor charge determination algorithm (Turbocharger) have been implemented. These additions into ProteoWizard provide universal tools that are independent of vendor platform for tandem mass spectrometry analyses and have particular utility for intralaboratory studies requiring the advantages of different platforms convergent on a particular workflow or for interlaboratory investigations spanning multiple platforms. We compared results from these tools to those obtained using vendor and commercial software, finding that in all cases our algorithms resulted in a comparable number of identified peptides for simple and complex samples measured on Waters, Agilent, and AB SCIEX quadrupole time-of-flight and Thermo Q-Exactive mass spectrometers. The mass accuracy of matched precursor ions also compared favorably with vendor and commercial tools. Additionally, typical analysis runtimes (∼1–100 ms per MS/MS spectrum) were short enough to enable the practical use of these high-quality signal processing tools for large clinical and research data sets.

A standardized framing for reporting protein identifications in mzIdentML 1.2

Inferring which protein species have been detected in bottom-up proteomics experiments has been a challenging problem for which solutions have been maturing over the past decade. While many inference approaches now function well in isolation, comparing and reconciling the results generated across different tools remains difficult. It presently stands as one of the greatest barriers in collaborative efforts such as the Human Proteome Project and public repositories such as the PRoteomics IDEntifications (PRIDE) database. Here we present a framework for reporting protein identifications that seeks to improve capabilities for comparing results generated by different inference tools. This framework standardizes the terminology for describing protein identification results, associated with the HUPO-Proteomics Standards Initiative (PSI) mzIdentML standard, while still allowing for differing methodologies to reach that final state. It is proposed that developers of software for reporting identification results will adopt this terminology in their outputs. While the new terminology does not require any changes to the core mzIdentML model, it represents a significant change in practice, and, as such, the rules will be released via a new version of the mzIdentML specification (version 1.2) so that consumers of files are able to determine whether the new guidelines have been adopted by export software.

Employing ProteoWizard to Convert Raw Mass Spectrometry Data

After raw data have been captured by mass spectrometers in biological LC-MS/MS experiments, they must be converted from vendor-specific binary files to open-format files for manipulation by most software. This protocol details the use of ProteoWizard software for this conversion, taking format features, coding options, and vendor particularities into account. This protocol will aid researchers in preparing their data for analysis by database search engines and other bioinformatics tools.

QC metrics from CPTAC raw LC-MS/MS data interpreted through multivariate statistics

Shotgun proteomics experiments integrate a complex sequence of processes, any of which can introduce variability. Quality metrics computed from LC-MS/MS data have relied upon identifying MS/MS scans, but a new mode for the QuaMeter software produces metrics that are independent of identifications. Rather than evaluating each metric independently, we have created a robust multivariate statistical toolkit that accommodates the correlation structure of these metrics and allows for hierarchical relationships among data sets. The framework enables visualization and structural assessment of variability. Study 1 for the Clinical Proteomics Technology Assessment for Cancer (CPTAC), which analyzed three replicates of two common samples at each of two time points among 23 mass spectrometers in nine laboratories, provided the data to demonstrate this framework, and CPTAC Study 5 provided data from complex lysates under Standard Operating Procedures (SOPs) to complement these findings. Identification-independent quality metrics enabled the differentiation of sites and run-times through robust principal components analysis and subsequent factor analysis. Dissimilarity metrics revealed outliers in performance, and a nested ANOVA model revealed the extent to which all metrics or individual metrics were impacted by mass spectrometer and run time. Study 5 data revealed that even when SOPs have been applied, instrument-dependent variability remains prominent, although it may be reduced, while within-site variability is reduced significantly. Finally, identification-independent quality metrics were shown to be predictive of identification sensitivity in these data sets. QuaMeter and the associated multivariate framework are available from http://fenchurch.mc.vanderbilt.edu and http://homepages.uc.edu/~wang2x7/, respectively.

Neurovascular unit on a chip: implications for translational applications

The blood-brain barrier (BBB) dynamically controls exchange between the brain and the body, but this interaction cannot be studied directly in the intact human brain or sufficiently represented by animal models. Most existing in vitro BBB models do not include neurons and glia with other BBB elements and do not adequately predict drug efficacy and toxicity. Under the National Institutes of Health Microtissue Initiative, we are developing a three-dimensional, multicompartment, organotypic microphysiological system representative of a neurovascular unit of the brain. The neurovascular unit system will serve as a model to study interactions between the central nervous system neurons and the cerebral spinal fluid (CSF) compartment, all coupled to a realistic blood-surrogate supply and venous return system that also incorporates circulating immune cells and the choroid plexus. Hence all three critical brain barriers will be recapitulated: blood-brain, brain-CSF, and blood-CSF. Primary and stem cell-derived human cells will interact with a variety of agents to produce critical chemical communications across the BBB and between brain regions. Cytomegalovirus, a common herpesvirus, will be used as an initial model of infections regulated by the BBB. This novel technological platform, which combines innovative microfluidics, cell culture, analytical instruments, bioinformatics, control theory, neuroscience, and drug discovery, will replicate chemical communication, molecular trafficking, and inflammation in the brain. The platform will enable targeted and clinically relevant nutritional and pharmacologic interventions for or prevention of such chronic diseases as obesity and acute injury such as stroke, and will uncover potential adverse effects of drugs. If successful, this project will produce clinically useful technologies and reveal new insights into how the brain receives, modifies, and is affected by drugs, other neurotropic agents, and diseases.

Identification of proteins at active, stalled, and collapsed replication forks using isolation of proteins on nascent DNA (iPOND) coupled with mass spectrometry

Both DNA and chromatin need to be duplicated during each cell division cycle. Replication happens in the context of defects in the DNA template and other forms of replication stress that present challenges to both genetic and epigenetic inheritance. The replication machinery is highly regulated by replication stress responses to accomplish this goal. To identify important replication and stress response proteins, we combined isolation of proteins on nascent DNA (iPOND) with quantitative mass spectrometry. We identified 290 proteins enriched on newly replicated DNA at active, stalled, and collapsed replication forks. Approximately 16% of these proteins are known replication or DNA damage response proteins. Genetic analysis indicates that several of the newly identified proteins are needed to facilitate DNA replication, especially under stressed conditions. Our data provide a useful resource for investigators studying DNA replication and the replication stress response and validate the use of iPOND combined with mass spectrometry as a discovery tool.

IDPQuantify: combining precursor intensity with spectral counts for protein and peptide quantification

Differentiating and quantifying protein differences in complex samples produces significant challenges in sensitivity and specificity. Label-free quantification can draw from two different information sources: precursor intensities and spectral counts. Intensities are accurate for calculating protein relative abundance, but values are often missing due to peptides that are identified sporadically. Spectral counting can reliably reproduce difference lists, but differentiating peptides or quantifying all but the most concentrated protein changes is usually beyond its abilities. Here we developed new software, IDPQuantify, to align multiple replicates using principal component analysis, extract accurate precursor intensities from MS data, and combine intensities with spectral counts for significant gains in differentiation and quantification. We have applied IDPQuantify to three comparative proteomic data sets featuring gold standard protein differences spiked in complicated backgrounds. The software is able to associate peptides with peaks that are otherwise left unidentified to increase the efficiency of protein quantification, especially for low-abundance proteins. By combing intensities with spectral counts from IDPicker, it gains an average of 30% more true positive differences among top differential proteins. IDPQuantify quantifies protein relative abundance accurately in these test data sets to produce good correlations between known and measured concentrations.

HI-bone: a scoring system for identifying phenylisothiocyanate-derivatized peptides based on precursor mass and high intensity fragment ions

Peptide sequence matching algorithms used for peptide identification by tandem mass spectrometry (MS/MS) enumerate theoretical peptides from the database, predict their fragment ions, and match them to the experimental MS/MS spectra. Here, we present an approach for scoring MS/MS identifications based on the high mass accuracy matching of precursor ions, the identification of a high intensity b1 fragment ion, and partial sequence tags from phenylthiocarbamoyl-derivatized peptides. This derivatization process boosts the b1 fragment ion signal, which turns it into a powerful feature for peptide identification. We demonstrate the effectiveness of our scoring system by implementing it on a computational tool called "HI-bone" and by identifying mass spectra of an Escherichia coli sample acquired on an Orbitrap Velos instrument using Higher-energy C-trap dissociation. Following this strategy, we identified 1614 peptide spectrum matches with a peptide false discovery rate (FDR) below 1%. These results were significantly higher than those from Mascot and SEQUEST using a similar FDR.

Basophile: accurate fragment charge state prediction improves peptide identification rates

In shotgun proteomics, database search algorithms rely on fragmentation models to predict fragment ions that should be observed for a given peptide sequence. The most widely used strategy (Naive model) is oversimplified, cleaving all peptide bonds with equal probability to produce fragments of all charges below that of the precursor ion. More accurate models, based on fragmentation simulation, are too computationally intensive for on-the-fly use in database search algorithms. We have created an ordinal-regression-based model called Basophile that takes fragment size and basic residue distribution into account when determining the charge retention during CID/higher-energy collision induced dissociation (HCD) of charged peptides. This model improves the accuracy of predictions by reducing the number of unnecessary fragments that are routinely predicted for highly-charged precursors. Basophile increased the identification rates by 26% (on average) over the Naive model, when analyzing triply-charged precursors from ion trap data. Basophile achieves simplicity and speed by solving the prediction problem with an ordinal regression equation, which can be incorporated into any database search software for shotgun proteomic identification.