RMassBank: Non-standard usage
27 October 2020
1 Introduction
##
## Attaching package: 'gplots'## The following object is masked from 'package:stats':
##
## lowessThis vignette assumes you are familiar with the standard usage of RMassBank, which is documented in
vignette("RMassBank")2 Skipping recalibration
For instances where recalibration is not wanted, e.g. there is
not enough data, or the user wants to use non-recalibrated
data, recalibration can be deactivated. To do this, the recalibrator
entry in the settings must be set to recalibrate.identity. This can be
done in the settings file directly (preferred):
recalibrator:
MS1: recalibrate.identity
MS2: recalibrate.identityOr, alternatively, the settings can be adapted directly via R code.
RmbDefaultSettings()
rmbo <- getOption("RMassBank")
rmbo$recalibrator <- list(
"MS1" = "recalibrate.identity",
"MS2" = "recalibrate.identity"
)
options("RMassBank" = rmbo)To show the results of using a non-recalibrated workflow, we load a workspace with pre-processed data:
w <- loadMsmsWorkspace(system.file("results/pH_narcotics_RF.RData",
package="RMassBankData"))## Warning in .updateObject.RmbWorkspace.1to2(w, ..., verbose): You are
## loading an archive from an old RMassBank version. The aggregate tables are
## not loaded from the original object, but recomputed.## Warning in .updateObject.RmbWorkspace.1to2(w, ..., verbose): If you hand-
## edited any aggregate table, the information might not be retained in the
## new object.## Warning in .updateObject.RmbWorkspace.1to2(w, ..., verbose): You are
## loading an archive from an old RMassBank version. The aggregate tables are
## not loaded from the original object, but recomputed.## Warning in .updateObject.RmbWorkspace.1to2(w, ..., verbose): If you hand-
## edited any aggregate table, the information might not be retained in the
## new object.## Warning in .updateObject.RmbWorkspace.1to2(w, ..., verbose): You are
## loading an archive from an old RMassBank version. The multiplicity
## filtering results are not loaded from the original object, but recomputed.## Warning in .updateObject.RmbWorkspace.1to2(w, ..., verbose): If you hand-
## edited any multiplicity filtering results, the information might not be
## retained in the new object.The recalibration curve:
recal <- makeRecalibration(w@parent, "pH",
recalibrateBy = rmbo$recalibrateBy,
recalibrateMS1 = rmbo$recalibrateMS1,
recalibrator = list(MS1="recalibrate.loess",MS2="recalibrate.loess"),
recalibrateMS1Window = 15)
w@rc <- recal$rc
w@rc.ms1 <- recal$rc.ms1
w@parent <- w
plotRecalibration(w)
Some example peaks to show the effect of recalibration:
w@spectra[[1]]@parent@mz[30:32]## [1] 133.0346 133.0381 133.9772w@spectra[[1]]@children[[1]]@mz[15:17]## [1] 87.00286 88.13592 97.10767Now reprocess the recalibration step with the
above set recalibration.identity:
w <- msmsWorkflow(w, steps=4)## msmsWorkflow: Step 4. Recalibrate m/z values in raw spectra## msmsWorkflow: Done.
The recalibration graph shows that the recalibration “curve” will do no recalibration. To verify, we can show the same peaks as before:
w@spectra[[1]]@parent@mz[30:32]## [1] 133.0346 133.0381 133.9772w@spectra[[1]]@children[[1]]@mz[15:17]## [1] 87.00286 88.13592 97.107673 Combining multiplicities
Standard multiplicity filtering, which is configurable in the settings, eliminates peaks which are observed only once for a compound. This eliminates spurious formula matches for random noise efficiently. It works well if either many spectra are recorded per compound, or if the same collision energy is present twice (e.g. with different resolutions). It sometimes fails for spectra on the “outer end” of the recorded collision energies when that spectrum is only present once – peaks which appear only in the highest or only in the lowest recorded energy can be erroneously deleted. To prevent this, one can re-run the workflow, read a second set of spectra for every compound (the second most intense) and combine the peak multiplicities of the two analyzed runs. (Mutiplicity filtering can also be switched off completely.)
Example:
RmbDefaultSettings()
getOption("RMassBank")$multiplicityFilter## [1] 2# to make processing faster, we only use 3 spectra per compound
rmbo <- getOption("RMassBank")
rmbo$spectraList <- list(
list(mode="CID", ces = "35%", ce = "35 % (nominal)", res = 7500),
list(mode="HCD", ces = "15%", ce = "15 % (nominal)", res = 7500),
list(mode="HCD", ces = "30%", ce = "30 % (nominal)", res = 7500)
)
options(RMassBank = rmbo)
loadList(system.file("list/NarcoticsDataset.csv",
package="RMassBankData"))## Loaded compoundlist successfullyw <- newMsmsWorkspace()
files <- list.files(system.file("spectra", package="RMassBankData"),
".mzML", full.names = TRUE)
w@files <- files[1:2]First, the spectra are read and processed until reanalysis (step 7) normally:
w1 <- msmsWorkflow(w, mode="pH", steps=c(1))## msmsWorkflow: Step 1. Acquire all MSMS spectra from files## msmsWorkflow: Done.# Here we artificially cut spectra out to make the workflow run faster for the vignette:
w1@spectra <- as(lapply(w1@spectra, function(s)
{
s@children <- s@children[1:3]
s
}),"SimpleList")
w1 <- msmsWorkflow(w1, mode="pH", steps=c(2:7))## msmsWorkflow: Step 2. First analysis pre recalibration## msmsWorkflow: Step 3. Aggregate all spectra## msmsWorkflow: Step 4. Recalibrate m/z values in raw spectra## msmsWorkflow: Step 5. Reanalyze recalibrated spectra## msmsWorkflow: Step 6. Aggregate recalibrated results## msmsWorkflow: Step 7. Reanalyze fail peaks for N2 + O## msmsWorkflow: Done.
Subsequently, we re-read and process the “confirmation spectra”, i.e. the second-best spectra from the files. Therefore, we will have two sets of spectra for each compound and every real peak should in theory occur twice.
w2 <- msmsWorkflow(w, mode="pH", steps=c(1), confirmMode = 1)## msmsWorkflow: Step 1. Acquire all MSMS spectra from files## msmsWorkflow: Done.# Here we artificially cut spectra out to make the workflow run faster for the vignette:
w2@spectra <- as(lapply(w2@spectra, function(s)
{
s@children <- s@children[1:3]
s
}),"SimpleList")
w2 <- msmsWorkflow(w2, mode="pH", steps=c(2:7))## msmsWorkflow: Step 2. First analysis pre recalibration## msmsWorkflow: Step 3. Aggregate all spectra## msmsWorkflow: Step 4. Recalibrate m/z values in raw spectra## msmsWorkflow: Step 5. Reanalyze recalibrated spectra## msmsWorkflow: Step 6. Aggregate recalibrated results## msmsWorkflow: Step 7. Reanalyze fail peaks for N2 + O## msmsWorkflow: Done.
Finally, we combine the two workspaces for multiplicity filtering, and apply the last step in the workflow (multiplicity filtering).
wTotal <- combineMultiplicities(c(w1, w2))
wTotal <- msmsWorkflow(wTotal, steps=8, mode="pH", archivename = "output")## msmsWorkflow: Step 8. Peak multiplicity filtering## msmsWorkflow: Done.Subsequently, we can proceed as usual with mbWorkflow:
mb <- newMbWorkspace(wTotal)
# [...] load lists, execute workflow etc.4 Session information
sessionInfo()## R version 4.0.3 (2020-10-10)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 18.04.5 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.12-bioc/R/lib/libRblas.so
## LAPACK: /home/biocbuild/bbs-3.12-bioc/R/lib/libRlapack.so
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_US.UTF-8 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=en_US.UTF-8
## [9] LC_ADDRESS=en_US.UTF-8 LC_TELEPHONE=en_US.UTF-8
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=en_US.UTF-8
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] gplots_3.1.0 RMassBankData_1.27.0 RMassBank_3.0.0
## [4] Rcpp_1.0.5 BiocStyle_2.18.0
##
## loaded via a namespace (and not attached):
## [1] lattice_0.20-41 gtools_3.8.2 png_0.1-7
## [4] assertthat_0.2.1 digest_0.6.27 foreach_1.5.1
## [7] R6_2.4.1 plyr_1.8.6 mzID_1.28.0
## [10] stats4_4.0.3 evaluate_0.14 httr_1.4.2
## [13] ggplot2_3.3.2 pillar_1.4.6 itertools_0.1-3
## [16] zlibbioc_1.36.0 rlang_0.4.8 curl_4.3
## [19] magick_2.5.0 S4Vectors_0.28.0 preprocessCore_1.52.0
## [22] rmarkdown_2.5 mzR_2.24.0 BiocParallel_1.24.0
## [25] stringr_1.4.0 ProtGenerics_1.22.0 munsell_0.5.0
## [28] fingerprint_3.5.7 compiler_4.0.3 xfun_0.18
## [31] pkgconfig_2.0.3 BiocGenerics_0.36.0 pcaMethods_1.82.0
## [34] htmltools_0.5.0 tidyselect_1.1.0 tibble_3.0.4
## [37] bookdown_0.21 IRanges_2.24.0 codetools_0.2-16
## [40] XML_3.99-0.5 crayon_1.3.4 dplyr_1.0.2
## [43] bitops_1.0-6 MASS_7.3-53 grid_4.0.3
## [46] gtable_0.3.0 lifecycle_0.2.0 affy_1.68.0
## [49] magrittr_1.5 rcdk_3.5.0 scales_1.1.1
## [52] KernSmooth_2.23-17 ncdf4_1.17 stringi_1.5.3
## [55] impute_1.64.0 affyio_1.60.0 doParallel_1.0.16
## [58] limma_3.46.0 ellipsis_0.3.1 generics_0.0.2
## [61] vctrs_0.3.4 rjson_0.2.20 iterators_1.0.13
## [64] tools_4.0.3 Biobase_2.50.0 MSnbase_2.16.0
## [67] glue_1.4.2 purrr_0.3.4 rcdklibs_2.3
## [70] parallel_4.0.3 yaml_2.2.1 colorspace_1.4-1
## [73] BiocManager_1.30.10 vsn_3.58.0 caTools_1.18.0
## [76] MALDIquant_1.19.3 rJava_0.9-13 knitr_1.30