As described in the GC-EI-MS tutorial, CluMSID can also be used to analyse low resolution data – although using low resolution data comes at a cost.
In this example, we will use a similar sample (1uL Pseudomonas aeruginosa PA14 cell extract) as in the General Tutorial, measured with similar chromatography but on a different mass spectrometer, a Bruker amaZon ion trap instrument operated in ESI-(+) mode with auto-MS/MS. In addition to introducing a workflow for low resolution LC-MS/MS data, this example also demonstrates that CluMSID can work with data from different types of mass spectrometers.
We load the file from the
The extraction of spectra works the same way as with high resolution LC-MS/MS data:
Like in the GC-EI-MS example, we have to adjust
mz_tolerance to a much higher value compared to high resolution data, while the retention time tolerance can remain unchanged.
We see that we have similar numbers of spectra as in the General Tutorial, because we tried to keep all parameters except for the mass spectrometer type comparable.
As we do not have low resolution spectral libraries at hand, we skip the integration of previous knowledge on feature identities in this example and generate a distance matrix right away:
Starting from this distance matrix, we can use all the data exploration functions that
When we now make an MDS plot, we learn that the similarity data is very different from the comparable high resolution data:
To get a better overview of the data and the general similarity behaviour, we create a heat map of the distance matrix:
We clearly see that the heat map is generally a lot “warmer” than in the General Tutorial (an intuition that is supported by the histogram in the top-left corner), i.e. we have a higher general degree of similarity between spectra. That is not surprising as the m/z information has much fewer levels than in high resolution data and fragments of different sum formula are more likely to have indistinguishable mass-to-charge ratios.
We also see that some more or less compact clusters can be identified. This is easier to inspect in the dendrogram visualisation:
In conclusion, CluMSID is capable of processing low resolution LC-MS/MS data and if high resolution data is not available, it can be very useful to provide an overview of spectral similarities in low resolution data, thereby helping metabolite annotation if some individual metabolites can be identified by comparison to authentic standards. However, concerning feature annotation, high resolution methods should always be favoured for the many benefits they provide.
sessionInfo() #> R version 3.6.3 (2020-02-29) #> Platform: x86_64-pc-linux-gnu (64-bit) #> Running under: Ubuntu 18.04.4 LTS #> #> Matrix products: default #> BLAS: /home/biocbuild/bbs-3.10-bioc/R/lib/libRblas.so #> LAPACK: /home/biocbuild/bbs-3.10-bioc/R/lib/libRlapack.so #> #> locale: #>  LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C #>  LC_TIME=en_US.UTF-8 LC_COLLATE=C #>  LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 #>  LC_PAPER=en_US.UTF-8 LC_NAME=C #>  LC_ADDRESS=C LC_TELEPHONE=C #>  LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C #> #> attached base packages: #>  stats4 parallel stats graphics grDevices utils datasets #>  methods base #> #> other attached packages: #>  magrittr_1.5 metaMSdata_1.22.0 metaMS_1.22.0 #>  CAMERA_1.42.0 xcms_3.8.2 MSnbase_2.12.0 #>  ProtGenerics_1.18.0 S4Vectors_0.24.4 mzR_2.20.0 #>  Rcpp_188.8.131.52 BiocParallel_1.20.1 Biobase_2.46.0 #>  BiocGenerics_0.32.0 CluMSIDdata_1.2.0 CluMSID_1.2.1 #> #> loaded via a namespace (and not attached): #>  colorspace_1.4-1 ellipsis_0.3.0 htmlTable_1.13.3 #>  base64enc_0.1-3 rstudioapi_0.11 farver_2.0.3 #>  affyio_1.56.0 fansi_0.4.1 codetools_0.2-16 #>  splines_3.6.3 ncdf4_1.17 doParallel_1.0.15 #>  impute_1.60.0 robustbase_0.93-6 knitr_1.28 #>  Formula_1.2-3 jsonlite_1.6.1 cluster_2.1.0 #>  vsn_3.54.0 png_0.1-7 graph_1.64.0 #>  readr_1.3.1 BiocManager_1.30.10 compiler_3.6.3 #>  httr_1.4.1 backports_1.1.6 assertthat_0.2.1 #>  Matrix_1.2-18 lazyeval_0.2.2 limma_3.42.2 #>  cli_2.0.2 acepack_1.4.1 htmltools_0.4.0 #>  tools_3.6.3 igraph_1.2.5 coda_0.19-3 #>  gtable_0.3.0 glue_1.4.0 affy_1.64.0 #>  RANN_2.6.1 dplyr_0.8.5 MALDIquant_1.19.3 #>  statnet.common_4.3.0 vctrs_0.2.4 multtest_2.42.0 #>  gdata_2.18.0 ape_5.3 preprocessCore_1.48.0 #>  nlme_3.1-147 iterators_1.0.12 xfun_0.13 #>  stringr_1.4.0 network_1.16.0 lifecycle_0.2.0 #>  gtools_3.8.2 XML_3.99-0.3 DEoptimR_1.0-8 #>  zlibbioc_1.32.0 MASS_7.3-51.5 scales_1.1.0 #>  hms_0.5.3 pcaMethods_1.78.0 RBGL_1.62.1 #>  MassSpecWavelet_1.52.0 RColorBrewer_1.1-2 yaml_2.2.1 #>  gridExtra_2.3 ggplot2_3.3.0 rpart_4.1-15 #>  reshape_0.8.8 latticeExtra_0.6-29 stringi_1.4.6 #>  highr_0.8 foreach_1.5.0 checkmate_2.0.0 #>  caTools_1.18.0 rlang_0.4.5 pkgconfig_2.0.3 #>  bitops_1.0-6 mzID_1.24.0 evaluate_0.14 #>  lattice_0.20-41 purrr_0.3.3 labeling_0.3 #>  htmlwidgets_1.5.1 tidyselect_1.0.0 GGally_1.5.0 #>  plyr_1.8.6 R6_2.4.1 IRanges_2.20.2 #>  gplots_3.0.3 Hmisc_4.4-0 sna_2.5 #>  pillar_1.4.3 foreign_0.8-76 survival_3.1-12 #>  nnet_7.3-13 tibble_3.0.0 crayon_1.3.4 #>  KernSmooth_2.23-16 plotly_184.108.40.206 rmarkdown_2.1 #>  jpeg_0.1-8.1 grid_3.6.3 data.table_1.12.8 #>  digest_0.6.25 tidyr_1.0.2 dbscan_1.1-5 #>  munsell_0.5.0 viridisLite_0.3.0