This document reproduces the Figures presented in (2016). For a description of the theory behind applications shown here we refer to the original manuscript. The results differ slightly due to technical changes or bugfixes in msqc1 that have been implemented after the manuscript was printed.
msqc1 1.30.0
Here we provide an experiment data package containing all the spread sheets used for figures and supplemental figures published in Kockmann et al. (2016). The raw instrument files are accessible for registered users through https://fgcz-bfabric.uzh.ch (2022) as project 1959. Raw LC-MS data from all five platforms were imported into Skyline 3.1, see MacLean et al. (2010), processed, and are available though https://panoramaweb.org/labkey/MSQC1.url when published. This package contains data.frame objects exported by Skyline making the data available for the R world.
The scatter-plot displays the reference Light:Heavy ratio versus the on-column amount of heavy peptide of the MSQC1 peptide mix. Note, x and y axis are drawn in log scale.
show the dilution series
table(msqc1_dil$relative.amount)
##
## 0.025 0.05 0.2 1 2 5
## 2550 2550 2550 2550 2550 2550
show the peptide frequency
table(msqc1_dil$Peptide.Sequence)
##
## ADVTPADFSEWSK ALIVLAHSER AVQQPDGLAVLGIFLK DGLDAASYYAPVR
## 378 864 828 378
## EGHLSPDIVAEQK ESDTSYVSLK FEDENFILK FSTVAGESGSADTVR
## 1170 864 684 864
## GAGAFGYFEVTHDITK GAGSSEPVTGLDAK GGPFSDSYR GTFIIDPAAVIR
## 792 378 864 378
## GTFIIDPGGVIR GYSIFSYATK LFLQFGAQGSPFLK LGGNEQVTR
## 378 648 378 378
## NLSVEDAAR SADFTNFDPR TAENFR TPVISGGPYEYR
## 828 864 144 378
## TPVITGAPYEYR VEATFGVDESNAK VLDALQAIK VSFELFADK
## 378 378 864 864
## YILAGVENSK
## 378
show ion types
table(msqc1_dil$Fragment.Ion)
##
## b8 precursor precursor [M+1] precursor [M+2] y10
## 72 2106 2106 2106 720
## y11 y12 y4 y5 y6
## 144 144 396 990 1620
## y7 y8 y9
## 2106 1638 1152
show instruments
table(msqc1_dil$instrument)
##
## QExactive QExactiveHF QTRAP TRIPLETOF TSQVantage
## 3996 4032 1242 4032 1998
Sigma mix peptide level signals - The graph displays the log2 L:H area ratios of eight technical replicates over 14 peptides from five MS platforms. The 14 panels were ordered by the on column amount of heavy peptide per vial (0.8, 4, 20, 40, 80, 200, 500, 1000 fmol). The black line indicates the theoretical L:H ratio as reported in the QC certificate by Sigma-Aldrich. In a perfect setting all data points would be located close to the black line indicating a perfect match between theoretical and measured log2 L:H ratios. The dark gray boxes correspond to a 0.5 and the light grey boxes to a deviation of 1 from the expected value (black line).
S <- .shape_for_volcano(S=msqc1_8rep, peptides)
msqc1:::.figure_setup()
xyplot(-log(p.value, 10) ~ log2FC | instrument, data=S, group=Peptide.Sequence,
panel = function(...){
panel.abline(v=c(-1,1), col='lightgray')
panel.abline(h=-log(0.05,10), col='lightgray')
panel.xyplot(...)
},
ylab='-log10 of p-value',
xlab='log2 fold change',
layout=c(1,5),
auto.key = list(space = "right", points = TRUE, lines = FALSE, cex=1))
Ratio stability upon analyte dilution - Each scatter-plot panel displays the experimental derived log2 L:H ratios versus the relative amount. The panels are ordered by the SIL on column amount (lower left to upper right). Color grouping was done by instrument. The loess fit curves were added for visualizing the trend. The SIL value given in each panel legend is valid for the relative amount of 1. The horizontal black line indicates the theoretical log2 L:H ratio.
Accuracy - The graph displays in each panel a sensitivity curves for a given relative amount.
sessionInfo()
## R version 4.3.1 (2023-06-16)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 22.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.18-bioc/R/lib/libRblas.so
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_GB LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## time zone: America/New_York
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] msqc1_1.30.0 lattice_0.22-5 BiocStyle_2.30.0
##
## loaded via a namespace (and not attached):
## [1] digest_0.6.33 R6_2.5.1 bookdown_0.36
## [4] fastmap_1.1.1 xfun_0.40 magrittr_2.0.3
## [7] cachem_1.0.8 knitr_1.44 htmltools_0.5.6.1
## [10] rmarkdown_2.25 cli_3.6.1 grid_4.3.1
## [13] sass_0.4.7 jquerylib_0.1.4 compiler_4.3.1
## [16] tools_4.3.1 evaluate_0.22 bslib_0.5.1
## [19] Rcpp_1.0.11 magick_2.8.1 yaml_2.3.7
## [22] BiocManager_1.30.22 jsonlite_1.8.7 rlang_1.1.1
Kockmann, Tobias, Christian Trachsel, Christian Panse, Åsa Wåhlander, Nathalie Selevsek, Jonas Grossmann, Witold E. Wolski, and Ralph Schlapbach. 2016. “Targeted proteomics coming of age - SRM, PRM and DIA performance evaluated from a core facility perspective.” PROTEOMICS. http://onlinelibrary.wiley.com/doi/10.1002/pmic.201500502/full.
MacLean, B., D. M. Tomazela, N. Shulman, M. Chambers, G. L. Finney, B. Frewen, R. Kern, D. L. Tabb, D. C. Liebler, and M. J. MacCoss. 2010. “Skyline: an open source document editor for creating and analyzing targeted proteomics experiments.” Bioinformatics 26 (7): 966–68.
Panse, Christian, Christian Trachsel, and Can Türker. 2022. “Bridging Data Management Platforms and Visualization Tools to Enable Ad-Hoc and Smart Analytics in Life Sciences.” Journal of Integrative Bioinformatics. https://doi.org/10.1515/jib-2022-0031.