Next Article in Journal
Impact of Exposure to Chronic Light–Dark Phase Shifting Circadian Rhythm Disruption on Muscle Proteome in Periparturient Dairy Cows
Next Article in Special Issue
Proteomes Are of Proteoforms: Embracing the Complexity
Previous Article in Journal
Downregulation of Salivary Proteins, Protective against Dental Caries, in Type 1 Diabetes
Previous Article in Special Issue
Proteome Discoverer—A Community Enhanced Data Processing Suite for Protein Informatics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Evaluation of the Sensitivity of Proteomics Methods Using the Absolute Copy Number of Proteins in a Single Cell as a Metric

Department of Pharmacology and Molecular Sciences, Johns Hopkins University, Baltimore, MD 21205, USA
Submission received: 2 July 2021 / Revised: 13 July 2021 / Accepted: 15 July 2021 / Published: 20 July 2021
(This article belongs to the Special Issue Feature Review Papers in Proteomes)

Abstract

:
Proteomic technology has improved at a staggering pace in recent years, with even practitioners challenged to keep up with new methods and hardware. The most common metric used for method performance is the number of peptides and proteins identified. While this metric may be helpful for proteomics researchers shopping for new hardware, this is often not the most biologically relevant metric. Biologists often utilize proteomics in the search for protein regulators that are of a lower relative copy number in the cell. In this review, I re-evaluate untargeted proteomics data using a simple graphical representation of the absolute copy number of proteins present in a single cancer cell as a metric. By comparing single-shot proteomics data to the coverage of the most in-depth proteomic analysis of that cell line acquired to date, we can obtain a rapid metric of method performance. Using a simple copy number metric allows visualization of how proteomics has developed in both sensitivity and overall dynamic range when using both relatively long and short acquisition times. To enable reanalysis beyond what is presented here, two available web applications have been developed for single- and multi-experiment comparisons with reference protein copy number data for multiple cell lines and organisms.

1. Introduction

Researchers who are interested in analyzing the global expression of protein have more options than ever before, due to a flurry of developments in proteomics technologies over the last 20 years [1]. Today, most proteomics work is performed using liquid chromatography-coupled tandem mass spectrometry (LCMS). While a few groups use LCMS successfully to analyze intact proteins via top-down proteomics, most work today is LCMS of proteolytically digested proteins, which is often referred to as shotgun proteomics [2]. With a dizzying number of hardware platforms, reagents, and methodologies to choose from, it is natural that researchers promote their favorite technology. Metrics for the performance of different methods do exist, with relative numbers of peptide and protein identifications per unit time being a metric of choice. A challenge in evaluating peptide and protein counts, as an objective metric for overall method performance, is in the number of variables that can be altered in the data processing pipelines that can affect these results. For example, utilizing a larger potential database to compare shotgun proteomics data to invariably increases the number of peptide identifications [3]. Increasing the search space further, to evaluate an increasing number of biologically likely post-translational modifications, will have a similar effect [4,5,6].
One metric of note is the “proteomics ruler”, developed by Wizniewski et al., and enabled as an add-in feature in the Andromeda post-search analysis bioinformatics package [7,8]. The proteomics ruler uses the relative quantification data of proteins in a sample and normalizes these abundances to that of the major histones that are present in mammalian cells. Histone proteins exist in a tightly conserved ratio in mammalian systems and directly proportional to the DNA present in a cell. Furthermore, the amount of DNA in a cell is an extremely consistent value. The proteomics ruler leverages these values as constants and can generate a remarkably accurate estimate of the absolute concentration of each protein in the cells analyzed. By comparing this concentration to the protein size, it is possible to estimate the number of copies of each individual protein within a single cell of the samples being analyzed. The proteomic ruler has been successfully applied to assess a wide array of samples, from human cancer cell lines through to mammalian organs, and is a valuable tool in shotgun proteomics today.
While the proteomics ruler has been employed in multiple studies, by the teams involved in the development of the method, it has not been more widely employed to date. This is likely due, at least in part, to the fact that it is currently only available through the Perseus program. In this review, I will more widely apply the proteomics ruler as a comparison between different proteomic technologies. Despite the ethical concerns regarding the use of HeLa cells in -omics studies, it continues to be actively used by proteomics labs worldwide, and is the most obvious cell line from which to perform this re-analysis exercise [9]. It is worth noting, however, that a recently preprinted interlaboratory study has described a large degree of variation in HeLa cell lines. Gene and protein expression profiles, as well as protein copy numbers, were shown to differ between the 14 HeLa stock samples that were obtained from 13 labs around the world [10]. The results herein should be treated with extreme caution.
To date, the most comprehensive shotgun proteomic analysis of HeLa cells was performed by Bekker-Jennsen et al. [11]. In that study, the proteomic ruler was used, and the protein copy numbers were obtained for over 12,000 distinct protein groups in the cell line. Through use of a simple R scripts and complementary web-based Shiny Apps, developed for this in this work, I will use this heavily fractionated proteome as the base metric. By applying the protein copy numbers derived in this study to the proteins identified by other proteomics techniques, we can obtain a simple visualization of the relative biological sensitivity of that method compared to others.

2. Materials and Methods

2.1. Obtaining Data

The copy number data for the HeLa cell lines were obtained from the processed output of the original studies obtained from ProteomeXchange partners. A table of the files presented here with references, if publicly available, is shown in Table 1. The identifiers from the processed data from the original studies were used when possible. When processed data were not directly available for comparison, the original vendor files were processed in-house. In-house searching was performed against the UniProt SwissProt database using the appropriate alkylation modification and the oxidation of methionine as the only variable modification. All data dependent files were searched with Proteome Discoverer 2.4 using the MSAmanda 2.0 search engine [12] and Percolator for false discovery rate (FDR) estimation. All default parameters in Proteome Discoverer using the vendor provided workflow templates “PWF_QE_Basic_Percolator” and “CWS_Basic” were utilized unless otherwise noted here. Orbitrap data were searched with a 10 ppm MS1 tolerance and a 0.02 Da MS/MS tolerance if high resolution and 0.6 Da MS/MS if ion trap. Data from TOF instruments were converted to MGF with ProteoWizard and searched using a 50 ppm MS1 and 0.1Da MS/MS tolerance. All data-independent acquisition (DIA) data presented herein are based on the results from the original studies. The protein lists used for the analysis of SOMASCAN data were obtained from published studies [13]. When necessary for previously published data the UniProt identifiers were extracted by pulling the list into R using the TidyVerse and Tabulizer packages [14,15].

2.2. Compiling the Absolute Copy Numbers

The UniProt identifier for the best protein identifier was removed from the final protein report from the reference data deposited. The majority protein and protein group identifiers were used from MaxQuant and Proteome Discoverer, respectively. MaxQuant assigned the majority protein to the accession that possesses greater than 50% of the peptides from all proteins with equivalent evidence [16]. Proteome Discoverer utilizes different logic depending on the version in use. In versions 1.0–2.0, the protein group is assigned to the protein with the highest total percentage coverage. In versions 2.1 onward, the largest protein in a group bearing equivalent evidence to all other proteins is assigned [17]. A recent analysis of 22 search engines demonstrated that resulting protein lists change little from the same MS/MS evidence [18], and therefore the two will be considered as equivalent in this light reanalysis.
The processed data from the Bekker-Jennsen et al. study identified 14,238 distinct proteins, of which 14.178 were assigned a copy number estimate based on IBAQ and the Proteomics Ruler Perseus accessory program [11]. I will refer to this as the HeLaHF dataset for the remainder of this work. Table 1 is a summary of the files used in this review.

2.3. Visualization of Copy Number Distribution in R/Shiny

The copy numbers of the Proteome Ruler and those applied to the identifiers are plotted using the base R histogram functions utilizing 30 bins and plotted versus number of total protein counts as well as normalized to density [27]. All work was performed in R studio. The R script is publicly available at https://github.com/orsburn/copynumbeR (accessed on 14 July 2021). The Shiny Apps files presented in this work and additional Proteomic Ruler base datasets can be found at https://www.lcmsmethods.org/methodtesting (accessed on 15 April 2021). All files utilized in this review were obtained from ProteomeXchange partner repositories as referenced in the Supplemental Information [28].

3. Results

3.1. Generational Improvements in Proteomics Hardware for Data Dependent Analysis

One of the most powerful forces driving the growth of proteomics as a field has been the increase in LCMS hardware performance over time. For a more thorough review of this topic, please see “The One Hour Yeast Proteome”, which thoroughly covers this topic up to the date of its publication [29]. A clear example of this increase in performance is the step from the original Q-Exactive system to the Q-Exactive high-field (HF) system. Several improvements in the architecture of the HF system exist over the “Classic” system, and these have been thoroughly described by others [20,30]. A more refined series of lenses, segmented quadrupoles with more symmetrical isolation efficiency and a lower instrument overhead, undoubtedly have effects on instrument sensitivity. However, these changes were implemented in the Q-Exactive “Plus” system as well. In this author’s hands, the “Classic” and “Plus” have a similar performance for global proteomics (data not shown). However, replacing the larger D30 Orbitrap system for the smaller diameter D20 Orbitrap has a marked change in the overall instrument performance. By increasing the curvature of the electric field in the D20 system, spectra of the same resolution can be obtained in approximately half the time of the earlier design, effectively doubling the spectral acquisition rate [31]. Typical results from these two systems, utilizing a 15 cm PepMap column and identical HPLC systems for 200 ng of HeLa digest standard, are 2000 proteins for 120 min for the classic and 3400 proteins for the HF (Table 1). Figure 1A,C represent the distribution of the absolute copy numbers of proteins found in HeLa cells in grey, with the distribution of the proteins from the classic and HF devices, respectively, in blue. A vertical line was added for reference, to flag the apex of the copy number counts for the “Classic” system. The 1400 additional proteins that were identified by the HF, using this identical sample and chromatography system, predominantly belong to lower copy number protein groups. Figure 1B,D represent the normalized density distribution of these same numbers, with the vertical bar again indicating the apex reference point of the “Classic” system.

3.2. Use of Absolute Copy Number for Optimization of Chromatographic Conditions

A comprehensive optimization study of the Orbitrap Fusion II “Lumos” instruments, by Espadas et al., included work with gradient optimization on 50 cm columns using 1000 ng injections of HeLa digest. Table 1 includes a summary of the protein identifications that were obtained in the work, using 60, 90, 120, and 240 min gradients. At 60 min, the instrument identified an impressive 4475 protein groups and 240 min increased this number to 5604, with other times falling in-between these two points. As shown in Figure 2, increasing the length of the gradient did, in every case, increase the number of proteins that were identified, as well as leading to the discovery of lower copy number proteins. The increase in protein identifications is clearer when taking the protein numbers into account, rather than by visualizing in this manner, because the distribution of 10% of the proteome across 30 bins appears slight to the eye, demonstrating a limitation of copy number visualization. A better metric can be obtained by considering the increase in identified proteins per unit time, as shown in Figure 2E. By increasing the gradient time, more proteins are identified, but with decreasing returns per unit time. While 300 protein groups may be added by extending the run time from 60 to 90 min, to identify an additional 600 protein groups to the number acquired with a 120 min gradient, the total run time must double. While there are limited studies in the literature to observe in this manner, this does appear to be a trend that is instrument- and chromatography-specific. Table 1 includes the summary of three files from an Orbitrap Elite system generated by the authors. In our hands, the Orbitrap Elite appears to have a maximum coverage of approximately 4000 protein groups, almost regardless of the gradient length and scan type utilized. The only single-shot runs that have exceeded 4200 protein groups utilized a 75 cm column and a 540 min gradient, more than 2× the amount of time necessary to identify 3900 unique protein groups from 200 ng of HeLa digest sample.

3.3. Rapid Proteomics Methods

A common and accurate criticism of proteomics technology is the speed at which data could be obtained [32,33]. With improvements in mass spectrometry, chromatography material, and sample preparation techniques, progress has been made toward shortening the total time between the samples being obtained and analyzed. In both the commercial sector and in clinical sciences, both high-flow proteomics and capillary separation are utilized to increase both the up-time and speed of analysis [34,35,36,37]. Nanoflow is still the most utilized separation method in proteomics, but efforts are underway to decrease nanoLC run times [38]. As a summary of recent developments, Figure 3 is a visualization of three generations of Orbitrap instruments, and the results obtained when analyzing samples with gradients less than 30 min in length. The Orbitrap Velos file used nanoflow liquid chromatography and parallelized use of the Orbitrap for MS1 and ion trap for MS/MS, and achieved identification of 1171 protein groups in this reanalysis. Recent work using shorter gradients on high-field Orbitrap systems, demonstrates the progress in hardware performance over the last 10 years. The use of a chromatography system with rapid pre-formed gradients, the EvoSep1, demonstrated remarkable coverage when used for acquisition on a Q-Exactive HF system (Figure 3C) [38]. Building on the EvoSep technology, a high-field asymmetric waveform ion mobility spectrometer (FAIMS), equipped Exploris 480 system, identified 3182 proteins with high-resolution MS/MS spectra and a single FAIMS compensating voltage of −70 EV (Figure 3C) [23]. While alterations at this level are reasonably easy to visualize, it is important to consider how scaling and binning can affect any visualization. As an example, Figure 3C demonstrates an overlay of these same results when using twice the number of bins.

3.4. Absolute Sensitivity in Single-Shot Proteomics Today

The field of proteomics has been almost impossible for insiders to keep up with, as new methods appear even more frequently than new hardware [39]. Recently, new hardware designs that are leveraging sophisticated ion mobility devices have appeared, which have challenged the status quo of Orbitrap dominance in proteomics. Figure 4 is a comparison of the highest coverage single-shot analyses that the authors have seen to date. The files shown are from published studies using 2 h gradients, where A is a 1000 ng injection of HeLa digest ion an Orbitrap Fusion 2 “Lumos” system. Figure 4B is a file from a TIMSTOF Pro system, and Figure 4C is from the recently published results on the TIMSTOF system operated in pasefDIA. The Lumos achieves an impressive 5098 protein groups when reprocessed with MSAmanda 2.0 for this comparison. The TIMSTOF Pro system, when analyzed in the same manner, achieves 5970 protein groups using the same software [22]. The recently published data using pasefDIA achieves a remarkable 7699 protein groups when processed by the authors in the original study, which is a number that is over 54% of the total proteome reference numbers [40].

3.5. Match between Runs

Match between runs (MBR) allows proteomics results to be additive in nature, with identifications made by MS/MS in one run to be applied to a second run if the chromatographic features and isotopic profile match within set parameters [41]. A recent method, called BoxCar, leans heavily on MBR. In BoxCar, multiple MS1 scans are gas-phase fractionated to obtain a more democratic distribution of MS1 signal. By collecting multiple fractions and capping each gas-phase fraction at a set limit, it is more difficult for high-abundance ions to fully fill the Orbitrap, and therefore suppress the signal of other coeluting ions. BoxCar results in an increase in the signal-to-noise ratio of lower abundance ions, at a cost of cycle time, due to the amount of time spent acquiring MS1 scans [25]. Figure 5 is an example of the results obtained from a 60 min BoxCar injection of 500 ng of HeLa lysate. When MBR is employed, BoxCar on Q-Exactive HF system can identify 6479 protein groups when matched against a highly fractionated DDA library that is generated on the same system. Figure 5B shows the number of those proteins that are identified in the absence of an MS1 library or MBR. Methods that are derived from BoxCar that utilizes parallelization in the ion trap on a Tribrid device and BoxCarDIA have demonstrated promise in alleviating the relative cycle time hits from the original method [42,43]. In addition, recent work has demonstrated a mechanism for estimating false discovery rates in MBR, which could go a long way toward realizing the potential of this strategy [44,45].

3.6. Single-Cell Proteomics

During the review of this text, a preprint from Hartlmayr and Ctortecka et al. demonstrated single-cell proteomics, utilizing a novel semi-automated platform coupled to a FAIMS-equipped Exploris 480 system [26]. The authors used this platform to analyze single cells, including those prepared from a HeLa cell line. This timely study allows a more direct analysis of the relative accuracy of the proteomic ruler data itself, as well as captures an understanding of the current state of single-cell proteomics technology today. Both label-free proteomics and TMT-labeled proteomics, utilizing carrier channels in a similar manner to SCoPE, were performed. As shown in Figure 6, both the methods obtained similar coverage, with a single unlabeled cell resulting in 608 protein groups, while single cells that were loaded with a 20 cell carrier channels resulted in 769 protein groups in our reanalysis. While altogether remarkable that such a depth of coverage can be obtained from single human cells, the fact that the distribution does not appear fully biased to the single highest calculated copy number proteins hints at some level of uncertainty in the proteomic ruler data. As single cells are prepped in a much different manner than bulk cell homogenates, and with the recently described variability in HeLa cell cultures globally, this may be altogether unsurprising.

3.7. Additional Methods

SomaScan is a commercially available alternative to LCMS proteomics workflows which utilized nucleic acid aptamers.
SomaScan aptamers are arranged in multiple configurations, with the SomaScan 1300 containing the highest relative number of targets. Figure 7 is a comparison of the 1300 kit and targets against the HeLa HF copy number library [19]. Despite the relatively small number of targets compared to any modern LCMS-based proteomics method described in this work, this alternative technology targets proteins across the entire dynamic range. The median log copy number of proteins in the SOMASCAN 1300 kit is 4.47, giving it one of the deepest dynamic range distributions of any technology reviewed here. However, as shown in Figure 7B, the small number of relative targets is dwarfed by even a 21 min LCMS method, utilizing some of today’s best hardware, such as a FAIMS-equipped Exploris 480.
The Supplemental Material contains copy number distribution data from all the files described in Table S1.

4. Conclusions

Today, we still do not have a complete picture of what the proteome is, and it seems likely that shotgun proteomics may never be able to fully answer critical questions, such as “how many human proteoforms are there?” [46]. With proteoforms being the next currency in proteomics, top-down technology will need to continue to develop at a rocket pace to eventually pick up the slack [47]. Shotgun proteomics will not be going away anytime soon, and the tools have continued to mature to a point where we are competitive in coverage, time, and sensitivity with today’s genomics and transcriptomics techniques [48]. Proteomics, as a field, still has plenty to overcome, most notably standardizing sample preparation and methodology, and the maturation of informatics. When we control for sample preparation, LCMS-based proteomics has demonstrated remarkable intra- and inter-lab reproducibility [49]. As a growing and maturing field, with diverse biological problems to confront, method development and optimization will need to continue until all organisms and organelles have been successfully tackled. The goal of this review was to take a snapshot of where we are today, and to use visualization of protein copy number depth as an additional tool when making inevitable decisions.
New hardware is released by instrument vendors nearly every year. If the goal of a lab is to obtain deeper proteomic depth with no further alteration in workflows, moving up to the newer generation of hardware may be the best solution for that task. On the same hardware, increasing the LCMS gradient time should almost always lead to increases in the peptide and protein numbers, but, as shown in Figure 2, this may be a slope of steeply diminishing returns. For some projects, increasing the total acquisition time from 60 min to 240 min may be an acceptable solution. For commercial labs or those with more users than hardware, that may be too steep of a price for an additional 1000 protein identifications. A common criticism of proteomics has always been the relatively low number of samples in each study, relative to genomics or transcriptomic studies [50]. As such, some groups are under pressure to use shorter acquisition times to address clinical cohorts, and several of today’s hardware advances appear well-suited to these tasks (Figure 3). Today’s current-generation instruments can identify over 3000 proteins from a HeLa standard in as little as 21 min of total run time, nearly three times the number of slower instruments of previous generations (Figure 4).
We must also consider the quality of the evidence that we consider acceptable for making a peptide or protein identification. While some groups have argued in the past, [33] and some more recently [51], that a high-resolution mass and retention time is sufficient to assign a peptide identification, it is fair to say that this is not universally accepted today. However, peptide identification that is supplemented with algorithms such as match between runs (MBR) is becoming increasingly adopted, and may be essential in some cases to meet the increasing expectations of collaborators. With deep analyses of the accuracy of these tools [44], providing reasons for optimism, and the recent description of a method for false discovery rate estimation for MBR, MS1-based peptide identification may soon experience a renaissance [45]. The recently described BoxCar method for quadrupole Orbitrap systems leans more heavily on MBR than any previously described technique. By sacrificing the number of MS/MS scans per run, to obtain more MS1 scans with improved signal-to-noise ratios, peptides identified by the former decrease markedly (Figure 5). Without the use of tools that can perform MBR, or without MS1 libraries, BoxCar can appear to be a waste of time and effort. The increase in the signal-to-noise ratio is too tempting of a target, however, and work continues to develop methods based on these innovative methods [42,43]. BoxCar is not the only method today that sacrifices the number of MS/MS scans per LCMS run to obtain alternative data. The FAIMS front end for current-generation Orbitrap instruments can be operated at multiple compensation voltages (CV) in each run. Although each CV requires an MS1 scan and a corresponding decrease in the available MS/MS acquisition time, this gas-phase fractionation leads to overall increases in protein identifications [23].
One of the fastest growing areas in proteomics today is the application to single cells [52,53]. The LCMS community appears to be divided into two distinct camps, those leveraging the newest hardware advances to increase sequencing depth in unlabeled samples [54] and those utilizing chemical tags to amplify peptide signals [55,56,57]. Today, both techniques appear to be able to identify a few hundred proteins per cell, with each innovation adding just a bit more to the overall depth [58]. With a better understanding of today’s limits, such as the maximum loading channels that should be used in amplification-based experiments, further advances will continue to chip away at these limits [59]. One key limitation is the sample preparation, and two preprints that were posted in early 2021 have demonstrated nanoliter semi-automated workflows to address these limitations [26,60]. These are reaping obvious dividends, resulting in some of the most comprehensive single-cell proteomes described to date (Figure 6). With a technology of such focus that can realistically only obtain data on the highest copy number proteins, every advance should move these histograms toward distributions further to the left.
Finally, although LCMS has had a monopoly on proteomics for decades, this is clearly no longer the case. New technologies are appearing today to challenge the status quo in more directly measuring organism phenotypes. SOMASCAN is one early example that has continued to gain ground. The use of aptamer technology appears to be less biased by the absolute protein copy numbers in a cell than LCMS technology (Figure 7). As these arrays continue to increase in the number of probes that may be utilized per sample, SOMASCAN may provide a complementary technology for the identification of proteins that are the most difficult for LCMS.
Today, we can choose from a variety of tools for proteomics, each with its own strengths and weaknesses. Our hope is that this review provides some insight into where proteomics is today. Furthermore, I hope that the tools created in the construction of this manuscript can be helpful to other researchers as they make inevitable decisions and compromises.
I hope that this review helps provide some biological perspective of the sensitivity of proteomics in use today. Further methods are compared in the Supplemental Information. To help facilitate further investigation in this regard, the simple tools used in this work have been made publicly available at www.lcmsmethods.org/methodtesting (accessed on 14 July 2021).

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/proteomes9030034/s1, Supplemental Table S1: the peptide and protein counts of other files analyzed in this review. Supplemental Figures: plotted copy numbers of other files in this review.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

I would like to thank Alexis Norris for extensive assistance with the FlexDashboard and Conor Jenkins for verifying data accuracy with Python tools. In addition, I would like to thank the reviewers of this manuscript for the most constructive review process this author has experienced to date.

Conflicts of Interest

No conflicts of interest to report.

References

  1. Yates, J.R.; Ruse, C.I.; Nakorchevsky, A. Proteomics by mass spectrometry: Approaches, advances, and applications. Annu. Rev. Biomed. Eng. 2009, 11, 49–79. [Google Scholar] [CrossRef] [Green Version]
  2. Lössl, P.; van de Waterbeemd, M.; Heck, A.J. The diverse and expanding role of mass spectrometry in structural and molecular biology. EMBO J. 2016, 35, 2634–2657. [Google Scholar] [CrossRef] [PubMed]
  3. Prakash, A.; Majumder, S.; Ahmad, S.; Varkey, M.; Anish, T.A.; Jenkins, C.; Rigby, M.; Orsburn, B. Detection and verification of 2.3 million cancer mutations in NCI60 cancer cell lines with a cloud search engine. J. Proteom. 2019, 209, 103488. [Google Scholar] [CrossRef]
  4. Prakash, A.; Ahmad, S.; Majumder, S.; Jenkins, C.; Orsburn, B. Bolt: A New Age Peptide Search Engine for Comprehensive MS/MS Sequencing through Vast Protein Databases in Minutes. J. Am. Soc. Mass Spectrom. 2019, 30, 2408–2418. [Google Scholar] [CrossRef] [PubMed]
  5. Kong, A.T.; Leprevost, F.V.; Avtonomov, D.M.; Mellacheruvu, D.; Nesvizhskii, A.I. MSFragger: Ultrafast and comprehensive peptide identification in mass spectrometry-based proteomics. Nat. Methods 2017, 14, 513–520. [Google Scholar] [CrossRef] [Green Version]
  6. Solntsev, S.K.; Shortreed, M.R.; Frey, B.L.; Smith, L.M. Enhanced Global Post-translational Modification Discovery with MetaMorpheus. J. Proteome Res. 2018, 17, 1844–1851. [Google Scholar] [CrossRef]
  7. Wiśniewski, J.R.; Hein, M.Y.; Cox, J.; Mann, M. A “proteomic ruler” for protein copy number and concentration estimation without spike-in standards. Mol. Cell. Proteom. 2014, 13, 3497–3506. [Google Scholar] [CrossRef] [Green Version]
  8. Wiśniewski, J.R. Label-Free and Standard-Free Absolute Quantitative Proteomics Using the “Total Protein” and “Proteomic Ruler” Approaches. Methods Enzymol. 2017, 585, 49–60. [Google Scholar]
  9. Köcher, T.; Pichler, P.; De Pra, M.; Rieux, L.; Swart, R.; Mechtler, K. Development and performance evaluation of an ultralow flow nanoliquid chromatography-tandem mass spectrometry set-up. Proteomics 2014, 14, 1999–2007. [Google Scholar] [CrossRef]
  10. Liu, Y.; Mi, Y.; Mueller, T.; Kreibich, S.; Williams, E.G.; Van Drogen, A.; Borel, C.; Germain, P.-L.; Frank, M.; Bludau, I.; et al. Genomic, Proteomic and Phenotypic Heterogeneity in HeLa Cells across Laboratories: Implications for Reproducibility of Research Results. bioRxiv 2018. [Google Scholar] [CrossRef]
  11. Bekker-Jensen, D.B.; Kelstrup, C.D.; Batth, T.S.; Larsen, S.C.; Haldrup, C.; Bramsen, J.B.; Sorensen, K.D.; Hoyer, S.; Orntoft, T.F.; Andersen, C.L.; et al. An Optimized Shotgun Strategy for the Rapid Generation of Comprehensive Human Proteomes. Cell Syst. 2017, 4, 587–599.e4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Dorfer, V.; Pichler, P.; Stranzl, T.; Stadlmann, J.; Taus, T.; Winkler, S.; Mechtler, K. MS Amanda, a universal identification algorithm optimized for high accuracy tandem mass spectra. J. Proteome Res. 2014, 13, 3679–3684. [Google Scholar] [CrossRef]
  13. Coombs, K.M.; Simon, P.F.; McLeish, N.J.; Zahedi-Amiri, A.; Kobasa, D. Aptamer profiling of A549 cells infected with low-pathogenicity and high-pathogenicity influenza viruses. Viruses 2019, 11, 1028. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the Tidyverse. J. Open Source Softw. 2019, 4, 1686. [Google Scholar] [CrossRef]
  15. Leeper, T.J.; Tabulizer: Bindings for Tabula PDF Table Extractor Library. R Packag Version 0.2.2. 2018. Available online: https://rdrr.io/cran/tabulizer/ (accessed on 15 July 2021).
  16. Cox, J.; Hein, M.Y.; Luber, C.A.; Paron, I.; Nagaraj, N.; Mann, M. Accurate Proteome-wide Label-free Quantification by Delayed Normalization and Maximal Peptide Ratio Extraction, Termed MaxLFQ. Mol. Cell. Proteom. 2014, 13, 2513–2526. [Google Scholar] [CrossRef] [Green Version]
  17. Orsburn, B.C. Proteome Discoverer—A Community Enhanced Data Processing Suite for Protein Informatics. Proteomes 2021, 9, 15. [Google Scholar] [CrossRef]
  18. Matthiesen, R.; Prieto, G.; Beck, H.C. Comparing peptide spectra matches across search engines. Methods Mol. Biol. 2020, 2051, 133–143. [Google Scholar] [PubMed]
  19. Xiong, H.; Yan, J.; Cai, S.; He, Q.; Peng, D.; Liu, Z.; Liu, Y. Cancer protein biomarker discovery based on nucleic acid aptamers. Int. J. Biol. Macromol. 2019, 132, 190–202. [Google Scholar] [CrossRef]
  20. Scheltema, R.A.; Hauschild, J.-P.; Lange, O.; Hornburg, D.; Denisov, E.; Damoc, E.; Kuehn, A.; Makarov, A.; Mann, M. The Q Exactive HF, a Benchtop Mass Spectrometer with a Pre-filter, High-performance Quadrupole and an Ultra-high-field Orbitrap Analyzer. Mol. Cell. Proteom. 2014, 13, 3698–3708. [Google Scholar] [CrossRef] [Green Version]
  21. Espadas, G.; Borràs, E.; Chiva, C.; Sabidó, E. Evaluation of different peptide fragmentation types and mass analyzers in data-dependent methods using an Orbitrap Fusion Lumos Tribrid mass spectrometer. Proteomics 2017, 17, 1600416. [Google Scholar] [CrossRef]
  22. Meier, F.; Brunner, A.D.; Koch, S.; Koch, H.; Lubeck, M.; Krause, M.; Goedecke, N.; Decker, J.; Kosinski, T.; Park, M.A.; et al. Online parallel accumulation–serial fragmentation (PASEF) with a novel trapped ion mobility mass spectrometer. Mol. Cell. Proteom. 2018, 17, 2534–2545. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Bekker-Jensen, D.B.; Martínez-Val, A.; Steigerwald, S.; Rüther, P.; Fort, K.L.; Arrey, T.N.; Harder, A.; Makarov, A.; Olsen, J.V. A compact quadrupole-orbitrap mass spectrometer with FAIMS interface improves proteome coverage in short LC gradients. Mol. Cell. Proteom. 2020, 19, 716–729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Meier, F.; Brunner, A.-D.; Frank, M.; Ha, A.; Voytik, E.; Kaspar-Schoenefeld, S.; Lubeck, M.; Raether, O.; Aebersold, R.; Collins, B.C.; et al. Parallel accumulation—Serial fragmentation combined with data-independent acquisition (diaPASEF): Bottom-up proteomics with near optimal ion usage. bioRxiv 2019. [Google Scholar] [CrossRef]
  25. Meier, F.; Geyer, P.E.; Virreira Winter, S.; Cox, J.; Mann, M. BoxCar acquisition method enables single-shot proteomics at a depth of 10,000 proteins in 100 minutes. Nat. Methods 2018, 15, 440–448. [Google Scholar] [CrossRef]
  26. Hartlmayr, D.; Ctortecka, C.; Seth, A.; Mendjan, S.; Tourniaire, G.; Mechtler, K. An automated workflow for label-free and multiplexed single cell proteomics sample preparation at unprecedented sensitivity. bioRxiv 2021. [Google Scholar] [CrossRef]
  27. Moraga, P.; Moraga, P. Interactive dashboards with flexdashboard and Shiny. In Geospatial Health Data; Chapman and Hall/CRC: Boca Raton, FL, USA, 2019. [Google Scholar]
  28. Vizcaíno, J.A.; Deutsch, E.W.; Wang, R.; Csordas, A.; Reisinger, F.; Ríos, D.; Dianes, J.A.; Sun, Z.; Farrah, T.; Bandeira, N.; et al. ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nat. Biotechnol. 2014, 32, 223–226. [Google Scholar] [CrossRef]
  29. Hebert, A.S.; Richards, A.L.; Bailey, D.J.; Ulbrich, A.; Coughlin, E.E.; Westphall, M.S.; Coon, J.J. The One Hour Yeast Proteome. Mol. Cell. Proteom. 2014, 13, 339–347. [Google Scholar] [CrossRef] [Green Version]
  30. Kelstrup, C.D.; Jersie-Christensen, R.R.; Batth, T.S.; Arrey, T.N.; Kuehn, A.; Kellmann, M.; Olsen, J.V. Rapid and deep proteomes by faster sequencing on a benchtop quadrupole ultra-high-field orbitrap mass spectrometer. J. Proteome Res. 2014, 13, 6187–6195. [Google Scholar] [CrossRef]
  31. Kanawati, B.; Schmitt-Kopplin, P. Fundamentals and Applications of Fourier Transform Mass Spectrometry; Elsevier: Amsterdam, The Netherlands, 2019; ISBN 9780128140147. [Google Scholar]
  32. Zubarev, R.A. The challenge of the proteome dynamic range and its implications for in-depth proteomics. Proteomics 2013, 13, 723–726. [Google Scholar] [CrossRef]
  33. Baker, E.S.; Liu, T.; Petyuk, V.A.; Burnum-Johnson, K.E.; Ibrahim, Y.M.; Anderson, G.A.; Smith, R.D. Mass spectrometry for translational proteomics: Progress and clinical implications. Genome Med. 2012, 4, 63. [Google Scholar] [CrossRef] [Green Version]
  34. Rinas, A.; Jenkins, C.; Orsburn, B. Assessing a commercial capillary electrophoresis interface (ZipChip) for shotgun proteomi applications. bioRxiv 2019. [Google Scholar] [CrossRef]
  35. Jenkins, C.; Orsburn, B. Standard Flow Multiplexed Proteomics (SFloMPro)—An Accessible and Cost-Effective Alternative to NanoLC Workflows. bioRxiv 2020. [Google Scholar] [CrossRef]
  36. Grebe, S.K.G.; Singh, R.J. Clinical peptide and protein quantification by mass spectrometry (MS). TrAC Trends Anal. Chem. 2016, 84, 131–143. [Google Scholar] [CrossRef]
  37. Fernández-Niño, S.M.G.; Smith-Moritz, A.M.; Chan, L.J.G.; Adams, P.D.; Heazlewood, J.L.; Petzold, C.J. Standard flow liquid chromatography for shotgun proteomics in bioenergy research. Front. Bioeng. Biotechnol. 2015, 3, 44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Bache, N.; Geyer, P.E.; Bekker-Jensen, D.B.; Hoerning, O.; Falkenby, L.; Treit, P.V.; Doll, S.; Paron, I.; Müller, J.B.; Meier, F.; et al. A novel LC system embeds analytes in pre-formed gradients for rapid, ultra-robust proteomics. Mol. Cell. Proteom. 2018, 17, 2284–2296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Eliuk, S.; Makarov, A. Evolution of Orbitrap Mass Spectrometry Instrumentation. Annu. Rev. Anal. Chem. 2015, 8, 61–80. [Google Scholar] [CrossRef] [PubMed]
  40. Meier, F.; Brunner, A.D.; Frank, M.; Ha, A.; Bludau, I.; Voytik, E.; Kaspar-Schoenefeld, S.; Lubeck, M.; Raether, O.; Bache, N.; et al. diaPASEF: Parallel accumulation–serial fragmentation combined with data-independent acquisition. Nat. Methods 2020, 17, 1229–1236. [Google Scholar] [CrossRef]
  41. Tyanova, S.; Temu, T.; Cox, J. The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. Nat. Protoc. 2016, 11, 2301–2319. [Google Scholar] [CrossRef]
  42. Jenkins, C.; Orsburn, B. BoxCar assisted MS fragmentation (BAMF). bioRxiv 2019. [Google Scholar] [CrossRef] [Green Version]
  43. Mehta, D.; Scandola, S.; Uhrig, R.G. BoxCar and library-free data-independent acquisition substantially improve the depth, range, and completeness of label-free quantitative proteomics in Arabidopsis. bioRxiv 2021. [Google Scholar] [CrossRef]
  44. Lim, M.Y.; Paulo, J.A.; Gygi, S.P. Evaluating False Transfer Rates from the Match-between-Runs Algorithm with a Two-Proteome Model. J. Proteome Res. 2019, 18, 4020–4026. [Google Scholar] [CrossRef]
  45. Yu, F.; Haynes, S.E.; Nesvizhskii, A.I. Label-free quantification with FDR-controlled match-between-runs. bioRxiv 2020. [Google Scholar] [CrossRef]
  46. Aebersold, R.; Agar, J.N.; Amster, I.J.; Baker, M.S.; Bertozzi, C.R.; Boja, E.S.; Costello, C.E.; Cravatt, B.F.; Fenselau, C.; Garcia, B.A.; et al. How many human proteoforms are there? Nat. Chem. Biol. 2018, 14, 206–214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Smith, L.M.; Kelleher, N.L. Proteoforms as the next proteomics currency. Science 2018, 359, 1106–1107. [Google Scholar] [CrossRef] [PubMed]
  48. Specht, H.; Emmott, E.; Koller, T.; Slavov, N. High-throughput single-cell proteomics quantifies the emergence of macrophage heterogeneity. bioRxiv 2019. [Google Scholar] [CrossRef] [Green Version]
  49. Abbatiello, S.E.; Schilling, B.; Mani, D.R.; Zimmerman, L.J.; Hall, S.C.; MacLean, B.; Albertolle, M.; Allen, S.; Burgess, M.; Cusack, M.P.; et al. Large-scale interlaboratory study to develop, analytically validate and apply highly multiplexed, quantitative peptide assays to measure cancer-relevant proteins in plasma. Mol. Cell. Proteom. 2015, 14, 2357–2374. [Google Scholar] [CrossRef] [Green Version]
  50. Levin, Y. The role of statistical power analysis in quantitative proteomics. Proteomics 2011, 11, 2565–2567. [Google Scholar] [CrossRef]
  51. Ivanov, M.V.; Bubis, J.A.; Gorshkov, V.; Tarasova, I.A.; Levitsky, L.I.; Lobas, A.A.; Solovyeva, E.M.; Pridatchenko, M.L.; Kjeldsen, F.; Gorshkov, M.V. DirectMS1: MS/MS-Free Identification of 1000 Proteins of Cellular Proteomes in 5 Minutes. Anal. Chem. 2020, 92, 4326–4333. [Google Scholar] [CrossRef] [Green Version]
  52. Slavov, N. Unpicking the proteome in single cells. Science 2020, 367, 512–513. [Google Scholar] [CrossRef]
  53. Ctortecka, C.; Mechtler, K. The rise of single-cell proteomics. Anal. Sci. Adv. 2021, 2, 84–89. [Google Scholar] [CrossRef]
  54. Kelly, R.; Zhu, Y.; Liang, Y.; Cong, Y.; Piehowski, P.; Dou, M.; Zhao, R.; Qian, W.-J.; Burnum-Johnson, K.; Ansong, C. Single Cell Proteome Mapping of Tissue Heterogeneity Using Microfluidic Nanodroplet Sample Processing and Ultrasensitive LC-MS. J. Biomol. Tech. 2019, 30, S61. [Google Scholar]
  55. Specht, H.; Harmange, G.; Perlman, D.H.; Emmott, E.; Niziolek, Z.; Budnik, B.; Slavov, N. Automated sample preparation for high-throughput single-cell proteomics. bioRxiv 2018. [Google Scholar] [CrossRef]
  56. Huffman, R.G.; Chen, A.; Specht, H.; Slavov, N. DO-MS: Data-Driven Optimization of Mass Spectrometry Methods. J. Proteome Res. 2019, 18, 2493–2500. [Google Scholar] [CrossRef]
  57. Budnik, B.; Levy, E.; Harmange, G.; Slavov, N. Mass-spectrometry of single mammalian cells quantifies proteome heterogeneity during cell differentiation. bioRxiv 2017. [Google Scholar] [CrossRef] [Green Version]
  58. Stejskal, K.; Op de Beeck, J.; Dürnberger, G.; Jacobs, P.; Mechtler, K. Ultra-sensitive nanoLC-MS using second generation micro pillar array LC technology with Orbitrap Exploris 480 and FAIMS PRO to enable single cell proteomics. bioRxiv 2021. [Google Scholar] [CrossRef]
  59. Cheung, T.K.; Lee, C.Y.; Bayer, F.P.; McCoy, A.; Kuster, B.; Rose, C.M. Defining the carrier proteome limit for single-cell proteomics. Nat. Methods 2021, 18, 76–83. [Google Scholar] [CrossRef]
  60. Leduc, A.; Huffman, R.G.; Slavov, N. Droplet sample preparation for single-cell proteomics applied to the cell cycle. bioRxiv 2021. [Google Scholar] [CrossRef]
Figure 1. The copy number distribution of two identical sample injections and chromatography systems on a Q-Exactive classic by protein counts (A) and density (B), compared to a Q-Exactive HF system (C,D). The vertical line represents the apex copy number on the classic system.
Figure 1. The copy number distribution of two identical sample injections and chromatography systems on a Q-Exactive classic by protein counts (A) and density (B), compared to a Q-Exactive HF system (C,D). The vertical line represents the apex copy number on the classic system.
Proteomes 09 00034 g001
Figure 2. Visualizing diminishing returns in gradient extension. The distribution of protein copy numbers on an Orbitrap Fusion Lumos with the same chromatography conditions and utilizing the instrument’s highest scan acquisition rates with a total gradient length of 60, 90, 120, and 240 min (AD, respectively). (E) A plot of the number of proteins identified versus total gradient time demonstrating the exponential increase in run time required to improve coverage under these conditions.
Figure 2. Visualizing diminishing returns in gradient extension. The distribution of protein copy numbers on an Orbitrap Fusion Lumos with the same chromatography conditions and utilizing the instrument’s highest scan acquisition rates with a total gradient length of 60, 90, 120, and 240 min (AD, respectively). (E) A plot of the number of proteins identified versus total gradient time demonstrating the exponential increase in run time required to improve coverage under these conditions.
Proteomes 09 00034 g002
Figure 3. The distribution of protein absolute copy numbers in three methods of less than 30 min in total length. (A) An Orbitrap Velos 30 min gradient. (B) Results from a 21 min analysis on a Q-Exactive high-field system. (C) A 21 min gradient utilizing a FAIMS-equipped Exploris 480 system with a single compensating voltage of −70 EV. (D) An overlay of these data using a 60 bins with (A) in red, (B) in green and (C) in orange.
Figure 3. The distribution of protein absolute copy numbers in three methods of less than 30 min in total length. (A) An Orbitrap Velos 30 min gradient. (B) Results from a 21 min analysis on a Q-Exactive high-field system. (C) A 21 min gradient utilizing a FAIMS-equipped Exploris 480 system with a single compensating voltage of −70 EV. (D) An overlay of these data using a 60 bins with (A) in red, (B) in green and (C) in orange.
Proteomes 09 00034 g003
Figure 4. A comparison of three of high-coverage proteomics methods utilizing a 2 h total run time. The vertical line is used here as a reference point to help visualize the shift in copy number distribution toward the apex of the deepest coverage method. (A) Lumos system with 1000 ng injection and 50 cm column. (B) TIMSTOF Pro system operating in pasefDDA mode. (C) A modified TIMSTOF operating in pasefDIA mode.
Figure 4. A comparison of three of high-coverage proteomics methods utilizing a 2 h total run time. The vertical line is used here as a reference point to help visualize the shift in copy number distribution toward the apex of the deepest coverage method. (A) Lumos system with 1000 ng injection and 50 cm column. (B) TIMSTOF Pro system operating in pasefDDA mode. (C) A modified TIMSTOF operating in pasefDIA mode.
Proteomes 09 00034 g004
Figure 5. A comparison between (A) protein groups identified in a 60 min BoxCar run when match between runs is used on a comprehensive MS1 library. (B) The same file when match between runs is not employed.
Figure 5. A comparison between (A) protein groups identified in a 60 min BoxCar run when match between runs is used on a comprehensive MS1 library. (B) The same file when match between runs is not employed.
Proteomes 09 00034 g005
Figure 6. The copy number distribution of protein copy numbers identified in two single-cell methods. A TMT experiment consisting of single cells with a carrier channel of 20 cells is shown in blue. A label-free analysis of a single HeLa cell is overlaid in red. Apexes are defined by color and the mean log copy for each single-cell run is provided.
Figure 6. The copy number distribution of protein copy numbers identified in two single-cell methods. A TMT experiment consisting of single cells with a carrier channel of 20 cells is shown in blue. A label-free analysis of a single HeLa cell is overlaid in red. Apexes are defined by color and the mean log copy for each single-cell run is provided.
Proteomes 09 00034 g006
Figure 7. (A) The distribution of SomaScan 1300 targets by copy number by count. (B) Comparison of this distribution against a recently published 21 min method utilizing an Exploris 480 equipped with FAIMS with 60 bins utilized to increase the granularity of visualization.
Figure 7. (A) The distribution of SomaScan 1300 targets by copy number by count. (B) Comparison of this distribution against a recently published 21 min method utilizing an Exploris 480 equipped with FAIMS with 60 bins utilized to increase the granularity of visualization.
Proteomes 09 00034 g007
Table 1. A summary of select files described in this text with references to the original study. A more complete table of studies and files analyzed is available as supplemental information.
Table 1. A summary of select files described in this text with references to the original study. A more complete table of studies and files analyzed is available as supplemental information.
File DescriptionNumber of ProteinsMean Log Copy
Number
HeLa HF 2018 (23 h) [11]14,1794.2
SomaScan 1300 [19]13084.47
QE Classic 200 ng 120 min [20]20165.72
QE HF 200 ng 120 min [20]34875.4
Lumos HCD-IT 60 min [21]44355.25
Lumos HCD-IT 90 min47705.21
Lumos HCD-IT 120 min50985.17
Lumos HCD-IT 240 min56045.09
Velos OT-IT 30 min (PRIDE PXD011070)11715.98
TIMSTOF Pro pasefDDA 120 min [22]59705.04
Exploris 480 FAIMS 21 min [23]31825.32
pasefDIA 120 min [24]76994.77
QE HF BoxCar 1 ug 60 min (MBR) [25]64795.16
QE HF BoxCar 1 ug 60 min (MS/MS) [25]25055.77
Exploris 480 Single-Cell TMT 20× Carrier [26]7696.15
Exploris 480 Single-Cell LFQ [26]6086.17
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Orsburn, B.C. Evaluation of the Sensitivity of Proteomics Methods Using the Absolute Copy Number of Proteins in a Single Cell as a Metric. Proteomes 2021, 9, 34. https://0-doi-org.brum.beds.ac.uk/10.3390/proteomes9030034

AMA Style

Orsburn BC. Evaluation of the Sensitivity of Proteomics Methods Using the Absolute Copy Number of Proteins in a Single Cell as a Metric. Proteomes. 2021; 9(3):34. https://0-doi-org.brum.beds.ac.uk/10.3390/proteomes9030034

Chicago/Turabian Style

Orsburn, Benjamin C. 2021. "Evaluation of the Sensitivity of Proteomics Methods Using the Absolute Copy Number of Proteins in a Single Cell as a Metric" Proteomes 9, no. 3: 34. https://0-doi-org.brum.beds.ac.uk/10.3390/proteomes9030034

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop