Large-scale Proteomics Experiments Research Articles

Proteomics-based method to comprehensively model the removal of host cell protein impurities.

Mechanistic models mostly focus on the target protein and some selected process- or product-related impurities. For a better process understanding, however, it is advantageous to describe also reoccurring host cell protein impurities. Within the purification of biopharmaceuticals, the binding of host cell proteins to a chromatographic resin is far from being described comprehensively. For a broader coverage of the binding characteristics, large-scale proteomic data and systems level knowledge on protein interactions are key. However, a method for determining binding parameters of the entire host cell proteome to selected chromatography resins is still lacking. In this work, we have developed a method to determine binding parameters of all detected individual host cell proteins in an Escherichia coli harvest sample from large-scale proteomics experiments. The developed method was demonstrated to model abundant and problematic proteins, which are crucial impurities to be removed. For these 15 proteins covering varying concentration ranges, the model predicts the independently measured retention time during the validation gradient well. Finally, we optimized the anion exchange chromatography capture step in silico using the determined isotherm parameters of the persistent host cell protein contaminants. From these results, strategies can be developed to separate abundant and problematic impurities from the target antigen.

Evidence for widespread translation of 5' untranslated regions.

Ribosome profiling experiments support the translation of a range of novel human open reading frames. By contrast, most peptides from large-scale proteomics experiments derive from just one source, 5' untranslated regions. Across the human genome we find evidence for 192 translated upstream regions, most of which would produce protein isoforms with extended N-terminal ends. Almost all of these N-terminal extensions are from highly abundant genes, which suggests that the novel regions we detect are just the tip of the iceberg. These upstream regions have characteristics that are not typical of coding exons. Their GC-content is remarkably high, even higher than 5' regions in other genes, and a large majority have non-canonical start codons. Although some novel upstream regions have cross-species conservation - five have orthologues in invertebrates for example - the reading frames of two thirds are not conserved beyond simians. These non-conserved regions also have no evidence of purifying selection, which suggests that much of this translation is not functional. In addition, non-conserved upstream regions have significantly more peptides in cancer cell lines than would be expected, a strong indication that an aberrant or noisy translation initiation process may play an important role in translation from upstream regions.

Visualisation of proteome-wide ordered protein abundances in Trypanosoma brucei.

Background: Trypanosoma brucei is a protozoan parasite and etiological agent of human and animal African trypanosomiasis. It has a complex life cycle, but the most studied cellular types are the in vitro cultivated bloodstream- and procyclic-forms. These correspond to the replicating, mammalian host bloodstream-dwelling, slender trypomastigotes and tsetse vector midgut-dwelling procyclic lifecycle stages, respectively. Several proteomics studies have reported the differential abundance of proteins between these in vitro cultivated cell types. However, there are no datasets providing protein abundance, from most to least abundant, within and between both cell types. Methods: We used MaxQuant software 1.6.10.4 to reprocess a recent large-scale proteomics experiment from our laboratory and extracted intensity-based quantifications of the bloodstream and procyclic form proteomes. Results: We created a web interface to visually explore protein abundances within and between the in vitro cultivated T. brucei bloodstream and procyclic form proteomes. Conclusions: The protein abundance visualization tool, searchable by protein name(s) and attribute(s), is likely to be useful to the trypanosome research community. It will allow users to contextualise their proteins of interest in terms of their abundances in the T. brucei bloodstream and procyclic form proteomes.

A community resource to mass explore the wheat grain proteome and its application to the late-maturity alpha-amylase (LMA) problem.

Late-maturity alpha-amylase (LMA) is a wheat genetic defect causing the synthesis of high isoelectric point alpha-amylase following a temperature shock during mid-grain development or prolonged cold throughout grain development, both leading to starch degradation. While the physiology is well understood, the biochemical mechanisms involved in grain LMA response remain unclear. We have applied high-throughput proteomics to 4,061 wheat flours displaying a range of LMA activities. Using an array of statistical analyses to select LMA-responsive biomarkers, we have mined them using a suite of tools applicable to wheat proteins. We observed that LMA-affected grains activated their primary metabolisms such as glycolysis and gluconeogenesis; TCA cycle, along with DNA- and RNA- binding mechanisms; and protein translation. This logically transitioned to protein folding activities driven by chaperones and protein disulfide isomerase, as well as protein assembly via dimerisation and complexing. The secondary metabolism was also mobilized with the upregulation of phytohormones and chemical and defence responses. LMA further invoked cellular structures, including ribosomes, microtubules, and chromatin. Finally, and unsurprisingly, LMA expression greatly impacted grain storage proteins, as well as starch and other carbohydrates, with the upregulation of alpha-gliadins and starch metabolism, whereas LMW glutenin, stachyose, sucrose, UDP-galactose, and UDP-glucose were downregulated. To our knowledge, this is not only the first proteomics study tackling the wheat LMA issue but also the largest plant-based proteomics study published to date. Logistics, technicalities, requirements, and bottlenecks of such an ambitious large-scale high-throughput proteomics experiment along with the challenges associated with big data analyses are discussed.

High Sensitivity Limited Material Proteomics Empowered by Data-Independent Acquisition on Linear Ion Traps.

In recent years, the concept of cell heterogeneity in biology has gained increasing attention, concomitant with a push toward technologies capable of resolving such biological complexity at the molecular level. For single-cell proteomics using Mass Spectrometry (scMS) and low-input proteomics experiments, the sensitivity of an orbitrap mass analyzer can sometimes be limiting. Therefore, low-input proteomics and scMS could benefit from linear ion traps, which provide faster scanning speeds and higher sensitivity than an orbitrap mass analyzer, however at the cost of resolution. We optimized an acquisition method that combines the orbitrap and linear ion trap, as implemented on a tribrid instrument, while taking advantage of the high-field asymmetric waveform ion mobility spectrometry (FAIMS) pro interface, with a prime focus on low-input applications. First, we compared the performance of orbitrap- versus linear ion trap mass analyzers. Subsequently, we optimized critical method parameters for low-input measurement by data-independent acquisition on the linear ion trap mass analyzer. We conclude that linear ion traps mass analyzers combined with FAIMS and Whisper flow chromatography are well-tailored for low-input proteomics experiments, and can simultaneously increase the throughput and sensitivity of large-scale proteomics experiments where limited material is available, such as clinical samples and cellular subpopulations.

Visualisation of proteome-wide ordered protein abundances in Trypanosoma brucei

Background: Trypanosoma brucei is a protozoan parasite and etiological agent of human and animal African trypanosomiasis. It has a complex life cycle, but the most studied cellular types are the in vitro cultivated bloodstream- and procyclic-forms. These correspond to the replicating, mammalian host bloodstream-dwelling, slender trypomastigotes and tsetse vector midgut-dwelling procyclic lifecycle stages, respectively. Several proteomics studies have reported the differential abundance of proteins between these in vitro cultivated cell types. However, there are no datasets providing protein abundance, from most to least abundant, within and between both cell types. Methods: We used MaxQuant software 1.6.10.4 to reprocess a recent large-scale proteomics experiment from our laboratory and extracted intensity-based quantifications of the bloodstream and procyclic form proteomes. Results: We created a web interface to visually explore protein abundances within and between the in vitro cultivated T. brucei bloodstream and procyclic form proteomes. Conclusions: The protein abundance visualization tool, searchable by protein name(s) and attribute(s), is likely to be useful to the trypanosome research community. It will allow users to contextualise their proteins of interest in terms of their abundances in the T. brucei bloodstream and procyclic form proteomes.

Visualisation of proteome-wide ordered protein abundances in Trypanosoma brucei.

Background: Trypanosoma brucei is a protozoan parasite and etiological agent of human and animal African trypanosomiasis. It has a complex life cycle, but the most studied cellular types are the in vitro cultivated bloodstream- and procyclic-forms. These correspond to the replicating, mammalian host bloodstream-dwelling, slender trypomastigotes and tsetse vector midgut-dwelling procyclic lifecycle stages, respectively. Several proteomics studies have reported the differential abundance of proteins between these in vitro cultivated cell types. However, there are no datasets providing protein abundance, from most to least abundant, within and between both cell types. Methods: We used MaxQuant software 1.6.10.4 to reprocess a recent large-scale proteomics experiment from our laboratory and extracted intensity-based quantifications of the bloodstream and procyclic form proteomes. Results: We created a web interface to visually explore protein abundances within and between the in vitro cultivated T. brucei bloodstream and procyclic form proteomes. Conclusions: The protein abundance visualization tool, searchable by protein name(s) and attribute(s), is likely to be useful to the trypanosome research community. It will allow users to contextualise their proteins of interest in terms of their abundances in the T. brucei bloodstream and procyclic form proteomes.

PolySTest: Robust Statistical Testing of Proteomics Data with Missing Values Improves Detection of Biologically Relevant Features

Statistical testing remains one of the main challenges for high-confidence detection of differentially regulated proteins or peptides in large-scale quantitative proteomics experiments by mass spectrometry. Statistical tests need to be sufficiently robust to deal with experiment intrinsic data structures and variations and often also reduced feature coverage across different biological samples due to ubiquitous missing values. A robust statistical test provides accurate confidence scores of large-scale proteomics results, regardless of instrument platform, experimental protocol and software tools. However, the multitude of different combinations of experimental strategies, mass spectrometry techniques and informatics methods complicate the decision of choosing appropriate statistical approaches. We address this challenge by introducing PolySTest, a user-friendly web service for statistical testing, data browsing and data visualization. We introduce a new method, Miss test, that simultaneously tests for missingness and feature abundance, thereby complementing common statistical tests by rescuing otherwise discarded data features. We demonstrate that PolySTest with integrated Miss test achieves higher confidence and higher sensitivity for artificial and experimental proteomics data sets with known ground truth. Application of PolySTest to mass spectrometry based large-scale proteomics data obtained from differentiating muscle cells resulted in the rescue of 10-20% additional proteins in the identified molecular networks relevant to muscle differentiation. We conclude that PolySTest is a valuable addition to existing tools and instrument enhancements that improve coverage and depth of large-scale proteomics experiments. A fully functional demo version of PolySTest and Miss test is available via http://computproteomics.bmb.sdu.dk/Apps/PolySTest.

Deep profiling and custom databases improve detection of proteoforms generated by alternative splicing.

Alternative pre-mRNA splicing has long been proposed to contribute greatly to proteome complexity. However, the extent to which mature mRNA isoforms are successfully translated into protein remains controversial. Here, we used high-throughput RNA sequencing and mass spectrometry (MS)–based proteomics to better evaluate the translation of alternatively spliced mRNAs. To increase proteome coverage and improve protein quantitation, we optimized cell fractionation and sample processing steps at both the protein and peptide level. Furthermore, we generated a custom peptide database trained on analysis of RNA-seq data with MAJIQ, an algorithm optimized to detect and quantify differential and unannotated splice junction usage. We matched tandem mass spectra acquired by data-dependent acquisition (DDA) against our custom RNA-seq based database, as well as SWISS-PROT and RefSeq databases to improve identification of splicing-derived proteoforms by 28% compared with use of the SWISS-PROT database alone. Altogether, we identified peptide evidence for 554 alternate proteoforms corresponding to 274 genes. Our increased depth and detection of proteins also allowed us to track changes in the transcriptome and proteome induced by T-cell stimulation, as well as fluctuations in protein subcellular localization. In sum, our data here confirm that use of generic databases in proteomic studies underestimates the number of spliced mRNA isoforms that are translated into protein and provides a workflow that improves isoform detection in large-scale proteomic experiments.

Promyelocytic Leukemia Protein (PML) Requirement for Interferon-induced Global Cellular SUMOylation

We report that interferon (IFN) α treatment at short and long periods increases the global cellular SUMOylation and requires the presence of the SUMO E3 ligase promyelocytic leukemia protein (PML), the organizer of PML nuclear bodies (NBs). Several PML isoforms (PMLI-PMLVII) derived from a single PML gene by alternative splicing, share the same N-terminal region but differ in their C-terminal sequences. Introducing each of the human PML isoform in PML-negative cells revealed that enhanced SUMOylation in response to IFN is orchestrated by PMLIII and PMLIV. Large-scale proteomics experiments enabled the identification of 558 SUMO sites on 389 proteins, of which 172 sites showed differential regulation upon IFNα stimulation, including K49 from UBC9, the sole SUMO E2 protein. Furthermore, IFNα induces PML-dependent UBC9 transfer to the nuclear matrix where it colocalizes with PML within the NBs and enhances cellular SUMOylation levels. Our results demonstrate that SUMOylated UBC9 and PML are key players for IFN-increased cellular SUMOylation.

Updates on resources, software tools, and databases for plant proteomics in 2016-2017.

Proteomics data processing, annotation, and analysis can often lead to major hurdles in large-scale high-throughput bottom-up proteomics experiments. Given the recent rise in protein-based big datasets being generated, efforts in in silico tool development occurrences have had an unprecedented increase; so much so, that it has become increasingly difficult to keep track of all the advances in a particular academic year. However, these tools benefit the plant proteomics community in circumventing critical issues in data analysis and visualization, as these continually developing open-source and community-developed tools hold potential in future research efforts. This review will aim to introduce and summarize more than 50 software tools, databases, and resources developed and published during 2016-2017 under the following categories: tools for data pre-processing and analysis, statistical analysis tools, peptide identification tools, databases and spectral libraries, and data visualization and interpretation tools. Intended for a well-informed proteomics community, finally, efforts in data archiving and validation datasets for the community will be discussed as well. Additionally, the author delineates the current and most commonly used proteomics tools in order to introduce novice readers to this -omics discovery platform.

The program Proteocat as a tool for planning of proteomic experiments

ProteoCat is a computer program that has been designed to help researchers in the planning of large-scale proteomic experiments. The central part of this program is the unit of hydrolysis simulation that supports 4 proteases (trypsin, lysine C, endoproteinases Asp-N and GluC). For peptides obtained after virtual hydrolysis or loaded from data files a number of properties important in mass-spectrometric experiments can be calculated and predicted; the resultant data can be analyzed or filtered (to reduce a set of peptides). The program is using new and improved modifications of own earlier developed methods for pI prediction, which can be also predicted by means of popular pKa scales proposed by other reseachers. The algorithm for prediction of peptide retention time has been realized similarly to the algorithm used in the SSRCalc program. Using ProteoCat it is possible to estimate the coverage of amino acid sequences of analyzed proteins under defined limitation on peptides detection, as well as the possibility of assembly of peptide fragments with user-defined minimal sizes of “sticky” ends. The program has a graphical user interface, written on JAVA and available at http://www.ibmc.msk.ru/LPCIT/ProteoCat.

Comparing Simplification Strategies for the Skeletal Muscle Proteome.

Skeletal muscle is a complex tissue that is dominated by the presence of a few abundant proteins. This wide dynamic range can mask the presence of lower abundance proteins, which can be a confounding factor in large-scale proteomic experiments. In this study, we have investigated a number of pre-fractionation methods, at both the protein and peptide level, for the characterization of the skeletal muscle proteome. The analyses revealed that the use of OFFGEL isoelectric focusing yielded the largest number of protein identifications (>750) compared to alternative gel-based and protein equalization strategies. Further, OFFGEL led to a substantial enrichment of a different sub-population of the proteome. Filter-aided sample preparation (FASP), coupled to peptide-level OFFGEL provided more confidence in the results due to a substantial increase in the number of peptides assigned to each protein. The findings presented here support the use of a multiplexed approach to proteome characterization of skeletal muscle, which has a recognized imbalance in the dynamic range of its protein complement.

Evaluating the use of HILIC in large-scale, multi dimensional proteomics: Horses for courses?

Despite many recent advances in instrumentation, the sheer complexity of biological samples remains a major challenge in large-scale proteomics experiments, reflecting both the large number of protein isoforms and the wide dynamic range of their expression levels. However, while the dynamic range of expression levels for different components of the proteome is estimated to be ∼107–8, the equivalent dynamic range of LC–MS is currently limited to ∼106. Sample pre-fractionation has therefore become routinely used in large-scale proteomics to reduce sample complexity during MS analysis and thus alleviate the problem of ion suppression and undersampling. There is currently a wide range of chromatographic techniques that can be applied as a first dimension separation. Here, we systematically evaluated the use of hydrophilic interaction liquid chromatography (HILIC), in comparison with hSAX, as a first dimension for peptide fractionation in a bottom-up proteomics workflow. The data indicate that in addition to its role as a useful pre-enrichment method for PTM analysis, HILIC can provide a robust, orthogonal and high-resolution method for increasing the depth of proteome coverage in large-scale proteomics experiments. The data also indicate that the choice of using either HILIC, hSAX, or other methods, is best made taking into account the specific types of biological analyses being performed.

Alternatively Spliced Homologous Exons Have Ancient Origins and Are Highly Expressed at the Protein Level.

Alternative splicing of messenger RNA can generate a wide variety of mature RNA transcripts, and these transcripts may produce protein isoforms with diverse cellular functions. While there is much supporting evidence for the expression of alternative transcripts, the same is not true for the alternatively spliced protein products. Large-scale mass spectroscopy experiments have identified evidence of alternative splicing at the protein level, but with conflicting results. Here we carried out a rigorous analysis of the peptide evidence from eight large-scale proteomics experiments to assess the scale of alternative splicing that is detectable by high-resolution mass spectroscopy. We find fewer splice events than would be expected: we identified peptides for almost 64% of human protein coding genes, but detected just 282 splice events. This data suggests that most genes have a single dominant isoform at the protein level. Many of the alternative isoforms that we could identify were only subtly different from the main splice isoform. Very few of the splice events identified at the protein level disrupted functional domains, in stark contrast to the two thirds of splice events annotated in the human genome that would lead to the loss or damage of functional domains. The most striking result was that more than 20% of the splice isoforms we identified were generated by substituting one homologous exon for another. This is significantly more than would be expected from the frequency of these events in the genome. These homologous exon substitution events were remarkably conserved—all the homologous exons we identified evolved over 460 million years ago—and eight of the fourteen tissue-specific splice isoforms we identified were generated from homologous exons. The combination of proteomics evidence, ancient origin and tissue-specific splicing indicates that isoforms generated from homologous exons may have important cellular roles.

Most highly expressed protein-coding genes have a single dominant isoform.

Although eukaryotic cells express a wide range of alternatively spliced transcripts, it is not clear whether genes tend to express a range of transcripts simultaneously across cells, or produce dominant isoforms in a manner that is either tissue-specific or regardless of tissue. To date, large-scale investigations into the pattern of transcript expression across distinct tissues have produced contradictory results. Here, we attempt to determine whether genes express a dominant splice variant at the protein level. We interrogate peptides from eight large-scale human proteomics experiments and databases and find that there is a single dominant protein isoform, irrespective of tissue or cell type, for the vast majority of the protein-coding genes in these experiments, in partial agreement with the conclusions from the most recent large-scale RNAseq study. Remarkably, the dominant isoforms from the experimental proteomics analyses coincided overwhelmingly with the reference isoforms selected by two completely orthogonal sources, the consensus coding sequence variants, which are agreed upon by separate manual genome curation teams, and the principal isoforms from the APPRIS database, predicted automatically from the conservation of protein sequence, structure, and function.

ProteoСat: a tool for planning of proteomic experiments

ProteoCat is a computer program has been designed to help researchers in the planning of large-scale proteomic experiments. The central part of this program is the subprogram of hydrolysis simulation that supports 4 proteases (trypsin, lysine C, endoproteinases AspN and GluC). For the peptides obtained after virtual hydrolysis or loaded from data file a number of properties important in mass-spectrometric experiments can be calculated or predicted. The data can be analyzed or filtered to reduce a set of peptides. The program is using new and improved modification of our methods developed to predict pI and probability of peptide detection; pI can also be predicted for a number of popular pKa's scales, proposed by other investigators. The algorithm for prediction of peptide retention time was realized similar to the algorithm used in the program SSRCalc. ProteoCat can estimate the coverage of amino acid sequences of proteins under defined limitation on peptides detection, as well as the possibility of assembly of peptide fragments with user-defined size of "sticky" ends. The program has a graphical user interface, written on JAVA and available at http://www.ibmc.msk.ru/LPCIT/ProteoCat.

MixGF: Spectral Probabilities for Mixture Spectra from more than One Peptide

In large-scale proteomic experiments, multiple peptide precursors are often cofragmented simultaneously in the same mixture tandem mass (MS/MS) spectrum. These spectra tend to elude current computational tools because of the ubiquitous assumption that each spectrum is generated from only one peptide. Therefore, tools that consider multiple peptide matches to each MS/MS spectrum can potentially improve the relatively low spectrum identification rate often observed in proteomics experiments. More importantly, data independent acquisition protocols promoting the cofragmentation of multiple precursors are emerging as alternative methods that can greatly improve the throughput of peptide identifications but their success also depends on the availability of algorithms to identify multiple peptides from each MS/MS spectrum. Here we address a fundamental question in the identification of mixture MS/MS spectra: determining the statistical significance of multiple peptides matched to a given MS/MS spectrum. We propose the MixGF generating function model to rigorously compute the statistical significance of peptide identifications for mixture spectra and show that this approach improves the sensitivity of current mixture spectra database search tools by a ≈30-390%. Analysis of multiple data sets with MixGF reveals that in complex biological samples the number of identified mixture spectra can be as high as 20% of all the identified spectra and the number of unique peptides identified only in mixture spectra can be up to 35.4% of those identified in single-peptide spectra.

POEAS: Automated Plant Phenomic Analysis Using Plant Ontology

Biological enrichment analysis using gene ontology (GO) provides a global overview of the functional role of genes or proteins identified from large-scale genomic or proteomic experiments. Phenomic enrichment analysis of gene lists can provide an important layer of information as well as cellular components, molecular functions, and biological processes associated with gene lists. Plant phenomic enrichment analysis will be useful for performing new experiments to better understand plant systems and for the interpretation of gene or proteins identified from high-throughput experiments. Plant ontology (PO) is a compendium of terms to define the diverse phenotypic characteristics of plant species, including plant anatomy, morphology, and development stages. Adoption of this highly useful ontology is limited, when compared to GO, because of the lack of user-friendly tools that enable the use of PO for statistical enrichment analysis. To address this challenge, we introduce Plant Ontology Enrichment Analysis Server (POEAS) in the public domain. POEAS uses a simple list of genes as input data and performs enrichment analysis using Ontologizer 2.0 to provide results in two levels, enrichment results and visualization utilities, to generate ontological graphs that are of publication quality. POEAS also offers interactive options to identify user-defined background population sets, various multiple-testing correction methods, different enrichment calculation methods, and resampling tests to improve statistical significance. The availability of such a tool to perform phenomic enrichment analyses using plant genes as a complementary resource will permit the adoption of PO-based phenomic analysis as part of analytical workflows. POEAS can be accessed using the URL http://caps.ncbs.res.in/poeas.

Five Gene Probes Carry Most Of The Discriminatory Power Of The 70-Gene Risk Model In Multiple Myeloma

Five Gene Probes Carry Most Of The Discriminatory Power Of The 70-Gene Risk Model In Multiple Myeloma