Although originally developed for liquid chromatography-tandem mass spectrometry (LC-MS/MS) data, CluMSID can also be used with direct infusion-tandem mass spectrometry (DI-MS/MS) data.
Generally, the missing retention time dimension makes feature annotation in metabolomics harder but if only direct infusion data is at hand, CluMSID can help to get an overview of the chemodiversity of a sample measured by DI-MS/MS.
In this example, we will use a similar sample (1uL Pseudomonas aeruginosa PA14 cell extract) as in the General Tutorial, measured on the same machine, a Bruker maxisHD qTOF operated in ESI-(+) mode with auto-MS/MS but without chromatographic separation.
We load the file from the
The extraction of spectra works the same way as with LC-MS/MS data:
Merging of redundant spectra is less straightforward when retention time is not available. Depending on the MS/MS method it can be next to impossible to decide whether two spectra with the same precursor m/z and similar fragmentation patterns derive from the same analyte or from two different but structurally similar ones.
In this example, we would like to merge spectra with identical precursor ions only if they were recorded one right after another. We can do so by setting
rt_tolerance to 1 second:
We see that we have hardly reduced the number of spectra in the list. If we would decide to merge all spectra with identical precursor m/z from the entire run, we could do so by setting
rt_tolerance to the duration of the run, in this case approx. 250 seconds:
The resulting number of spectra is drastically lower but the danger of merging spectra that do not actually derive from the same analyte is also very big.
In this very explorative example, we skip the integration of previous knowledge on feature identities and generate a distance matrix right away:
Starting from this distance matrix, we can use all the data exploration functions that
CluMSID offers. In this example workflow, we look at a cluster dendrogram:
It is directly obvious that we have some spectra that are nearly identical and thus most likely derive from the same analyte, e.g. the many spectra with a precursor m/z of 270.19. But we still see nice clustering of similar spectra with different precursor m/z, e.g. the huge gray cluster that contains a lot of different alkylquinolone type metabolites (see General Tutorial).
In conclusion, CluMSID is very useful to provide an overview of spectral similarities within DI-MS/MS runs but wherever annotation is in the focus, one should not do without the additional layer of information created by chromatographic separation.
sessionInfo() #> R version 4.0.0 (2020-04-24) #> Platform: x86_64-pc-linux-gnu (64-bit) #> Running under: Ubuntu 18.04.4 LTS #> #> Matrix products: default #> BLAS: /home/biocbuild/bbs-3.11-bioc/R/lib/libRblas.so #> LAPACK: /home/biocbuild/bbs-3.11-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: #>  stats graphics grDevices utils datasets methods base #> #> other attached packages: #>  CluMSIDdata_1.3.0 CluMSID_1.4.0 #> #> loaded via a namespace (and not attached): #>  httr_1.4.1 Biobase_2.48.0 tidyr_1.0.2 #>  vsn_3.56.0 viridisLite_0.3.0 jsonlite_1.6.1 #>  foreach_1.5.0 network_1.16.0 gtools_3.8.2 #>  assertthat_0.2.1 highr_0.8 BiocManager_1.30.10 #>  affy_1.66.0 stats4_4.0.0 yaml_2.2.1 #>  impute_1.62.0 pillar_1.4.3 lattice_0.20-41 #>  glue_1.4.0 limma_3.44.0 digest_0.6.25 #>  RColorBrewer_1.1-2 colorspace_1.4-1 htmltools_0.4.0 #>  preprocessCore_1.50.0 plyr_1.8.6 MALDIquant_1.19.3 #>  XML_3.99-0.3 pkgconfig_2.0.3 zlibbioc_1.34.0 #>  purrr_0.3.4 scales_1.1.0 gdata_2.18.0 #>  sna_2.5 affyio_1.58.0 BiocParallel_1.22.0 #>  tibble_3.0.1 dbscan_1.1-5 IRanges_2.22.0 #>  ggplot2_3.3.0 ellipsis_0.3.0 lazyeval_0.2.2 #>  BiocGenerics_0.34.0 magrittr_1.5 crayon_1.3.4 #>  statnet.common_4.3.0 evaluate_0.14 GGally_1.5.0 #>  ncdf4_1.17 doParallel_1.0.15 nlme_3.1-147 #>  MASS_7.3-51.6 gplots_3.0.3 mzR_2.22.0 #>  tools_4.0.0 data.table_1.12.8 lifecycle_0.2.0 #>  stringr_1.4.0 plotly_18.104.22.168 MSnbase_2.14.0 #>  S4Vectors_0.26.0 munsell_0.5.0 pcaMethods_1.80.0 #>  compiler_4.0.0 mzID_1.26.0 caTools_1.18.0 #>  rlang_0.4.5 grid_4.0.0 iterators_1.0.12 #>  htmlwidgets_1.5.1 bitops_1.0-6 rmarkdown_2.1 #>  gtable_0.3.0 codetools_0.2-16 reshape_0.8.8 #>  R6_2.4.1 knitr_1.28 dplyr_0.8.5 #>  ProtGenerics_1.20.0 KernSmooth_2.23-17 ape_5.3 #>  stringi_1.4.6 parallel_4.0.0 Rcpp_22.214.171.124 #>  vctrs_0.2.4 tidyselect_1.0.0 xfun_0.13 #>  coda_0.19-3