MSstatsBig is designed to overcome challenges when analyzing very large mass spectrometry (MS)-based proteomics experiments. These experiments are generally (but not always) acquired with DIA and include a large number of MS runs. MSstatsBig leverages software that can work on datasets without loading them into memory. This avoids a major problem where a dataset cannot be loaded into a standard computers RAM.
MSstatsBig includes functions which are designed to replace the converters included in the MSstats package. The goal of these converters is to perform filtering on the PSM files of identified and quantified data to remove data that is not required for differential analysis. Once this data is filtered down it should be able to be loaded into your computer’s memory. After the converters are run, the standard MSstats workflow can be followed.
MSstatsBig currently includes converters for Spectronaut and FragPipe. Beyond these converters, users can manually use the underlying functions by putting their data into MSstats format and running the underlying MSstatsPreprocessBig function. This way, either by using native export format of signal processing tools or by converting raw data chunk by chunk (for example with the readr::read_delim_chunked function), MSstatsBig can be used with other popular tools such as DIA-NN.
The dataset included in this package is a small subset of a work by Clark et al. [1]. It is a large DIA dataset that includes over 100 runs. The experimental data was identified and quantified by FragPipe and the included dataset is the msstats.csv output for FragPipe.
## ProteinName PeptideSequence PrecursorCharge FragmentIon
## 1 Q86U42 (UniMod:1)AAAAAAAAAAGAAGGR 1 b4
## 2 Q86U42 (UniMod:1)AAAAAAAAAAGAAGGR 1 y6
## 3 Q86U42 (UniMod:1)AAAAAAAAAAGAAGGR 1 y7
## 4 Q86U42 (UniMod:1)AAAAAAAAAAGAAGGR 1 b5
## 5 Q86U42 (UniMod:1)AAAAAAAAAAGAAGGR 1 y8
## 6 Q86U42 (UniMod:1)AAAAAAAAAAGAAGGR 1 y9
## ProductCharge IsotopeLabelType Condition BioReplicate
## 1 1 L T 1522
## 2 1 L T 1522
## 3 1 L T 1522
## 4 1 L T 1522
## 5 1 L T 1522
## 6 1 L T 1522
## Run Intensity
## 1 CPTAC_CCRCC_W_JHU_20190112_LUMOS_C3N-01522_T 3747562.00
## 2 CPTAC_CCRCC_W_JHU_20190112_LUMOS_C3N-01522_T 770585.44
## 3 CPTAC_CCRCC_W_JHU_20190112_LUMOS_C3N-01522_T 1379359.12
## 4 CPTAC_CCRCC_W_JHU_20190112_LUMOS_C3N-01522_T 3706759.50
## 5 CPTAC_CCRCC_W_JHU_20190112_LUMOS_C3N-01522_T 77361.97
## 6 CPTAC_CCRCC_W_JHU_20190112_LUMOS_C3N-01522_T 338863.41First we run the MSstatsBig converter. The converter will save the dataset to a place on your computer, and will return an arrow object. Once then, you can read the data from text file or load the arrow data.frame into memory by using the dplyr::collect function. The “collected” data can then be treated as a standard R data.frame.
setwd(tempdir())
converted_data = bigFragPipetoMSstatsFormat(
system.file("extdata", "fgexample.csv", package = "MSstatsBig"),
"output_file.csv",
backend="arrow",
max_feature_count = 20)
# The returned arrow object needs to be collected for the remaining workflow
converted_data = as.data.frame(dplyr::collect(converted_data))Once the converter is run the standard MSstats workflow can be followed:
dataProcess function summarizes peptide-level data into protein-level estimates.groupComparison function fits a linear model to protein-level data and performs statistical inference for comparisons of selected groups.Details of the MSstats workflow can be found in [2].
##
## Attaching package: 'MSstats'## The following object is masked from 'package:grDevices':
##
## savePlot# converted_data = read.csv("output_file.csv")
summarized_data = dataProcess(converted_data,
use_log_file = FALSE)## INFO [2024-04-05 18:21:22] ** Features with one or two measurements across runs are removed.
## INFO [2024-04-05 18:21:22] ** Fractionation handled.
## INFO [2024-04-05 18:21:22] ** Updated quantification data to make balanced design. Missing values are marked by NA
## INFO [2024-04-05 18:21:22] ** Log2 intensities under cutoff = 12.59 were considered as censored missing values.
## INFO [2024-04-05 18:21:22] ** Log2 intensities = NA were considered as censored missing values.
## INFO [2024-04-05 18:21:22] ** Use all features that the dataset originally has.
## INFO [2024-04-05 18:21:22]
## # proteins: 1
## # peptides per protein: 1-1
## # features per peptide: 11-11
## INFO [2024-04-05 18:21:22]
## NAT T
## # runs 56 50
## # bioreplicates 56 50
## # tech. replicates 1 1
## INFO [2024-04-05 18:21:22] == Start the summarization per subplot...
## | | | 0% | |======================================================================| 100%
## INFO [2024-04-05 18:21:22] == Summarization is done.# Build contrast matrix
comparison = matrix(c(-1, 1),
nrow=1, byrow=TRUE)
row.names(comparison) <- c("T-NAT")
colnames(comparison) <- c("NAT", "T")
model_results = groupComparison(contrast.matrix = comparison,
data = summarized_data,
use_log_file = FALSE)## INFO [2024-04-05 18:21:22] == Start to test and get inference in whole plot ...
## | | | 0% | |======================================================================| 100%
## INFO [2024-04-05 18:21:22] == Comparisons for all proteins are done.## R version 4.3.3 (2024-02-29)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 22.04.4 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] MSstats_4.10.1 MSstatsBig_1.0.1
##
## loaded via a namespace (and not attached):
## [1] gtable_0.3.4 xfun_0.43 bslib_0.7.0
## [4] ggplot2_3.5.0 htmlwidgets_1.6.4 caTools_1.18.2
## [7] ggrepel_0.9.5 lattice_0.22-6 vctrs_0.6.5
## [10] tools_4.3.3 bitops_1.0-7 generics_0.1.3
## [13] tibble_3.2.1 fansi_1.0.6 pkgconfig_2.0.3
## [16] Matrix_1.6-5 KernSmooth_2.23-22 data.table_1.15.4
## [19] checkmate_2.3.1 assertthat_0.2.1 lifecycle_1.0.4
## [22] compiler_4.3.3 gplots_3.1.3.1 statmod_1.5.0
## [25] munsell_0.5.1 htmltools_0.5.8.1 sass_0.4.9
## [28] yaml_2.3.8 lazyeval_0.2.2 preprocessCore_1.64.0
## [31] marray_1.80.0 plotly_4.10.4 pillar_1.9.0
## [34] nloptr_2.0.3 jquerylib_0.1.4 tidyr_1.3.1
## [37] MASS_7.3-60.0.1 cachem_1.0.8 limma_3.58.1
## [40] boot_1.3-30 nlme_3.1-164 gtools_3.9.5
## [43] tidyselect_1.2.1 digest_0.6.35 dplyr_1.1.4
## [46] purrr_1.0.2 arrow_15.0.1 splines_4.3.3
## [49] fastmap_1.1.1 grid_4.3.3 colorspace_2.1-0
## [52] cli_3.6.2 magrittr_2.0.3 survival_3.5-8
## [55] utf8_1.2.4 withr_3.0.0 scales_1.3.0
## [58] backports_1.4.1 bit64_4.0.5 rmarkdown_2.26
## [61] httr_1.4.7 bit_4.0.5 lme4_1.1-35.2
## [64] evaluate_0.23 knitr_1.45 log4r_0.4.3
## [67] MSstatsConvert_1.12.1 viridisLite_0.4.2 rlang_1.1.3
## [70] Rcpp_1.0.12 glue_1.7.0 minqa_1.2.6
## [73] jsonlite_1.8.8 R6_2.5.1D. J. Clark, S. M. Dhanasekaran, F. Petralia, P. Wang and H. Zhang, “Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma,” Cell, vol. 179, pp. 964-983, 2019.
D. Kohler et al., “MSstats Version 4.0: Statistical Analyses of Quantitative Mass Spectrometry-Based Proteomic Experiments with Chromatography-Based Quantification at Scale”, J. Proteome Res. 22, 5, pp. 1466–1482, J. Proteome Res. 2023, 22, 5, 1466–1482