Package: peakPantheR
Authors: Arnaud Wolfer, Goncalo Correia
The peakPantheR
package is designed for the detection, integration and reporting of pre-defined features in MS files (e.g. compounds, fragments, adducts, …).
The graphical user interface implements all of peakPantheR
’s functionalities and can be preferred to understand the methodology, select the best parameters on a subset of the samples before running the command line, or to visually explore results.
Using the faahKO raw MS dataset as an example, this vignette will:
This vignette employ the .csv
or .RData
files generated from faahKO in the vignette Getting Started with peakPantheR.
The graphical user interface is started as follow:
The graphical interface is divided in 5 main tabs, Import Data, Run annotation, Diagnostic: plot & update, View results and Export results.
\(~\)
The first input format is using a .RData
file containing a peakPantheRAnnotation named annotationObject
. This object can be annotated or not, for example loading a previously run annotation (see the Export section for more details).
\(~\)
The second input format consists of multiple .csv
files describing the targeted features, spectra to process and corresponding metadata (optional). Spectra can also be directly selected on disk.
\(~\)
With the targeted features and spectra defined, Run annotation handles the fitting parameter selection as well as downstream computation. First the use of updated regions of interest (uROI
) and fallback integration regions (FIR
) can be selected if available. If uROI
haven’t been previously defined, the option will be crossed out. Secondly the curve fitting model to use can be selected from the interface. Finally Parallelisation
enables the selection of the number of CPU cores to employ for parallel file access and processing.
\(~\)
Note:
The targeted regions of interest (
ROI
) should represent a good starting point for feature integration, however it might be necessary to refine these boundary box to the specific analytical run considered. This ensures a successful integration over all the spectra irrespective of potential chromatographic equilibration differences or retention time drift.Updated regions of interest (
uROI
) can be defined and will supplantROI
.uROI
can for example be manually defined to “tighten” or correct theROI
and avoid erroneous integration. Another use ofuROI
is to encompass the integration region in each sample throughout the run without targeting any excess spectral region that would interfere with the correct analysis.Fallback integration regions (
FIR
) are defined as spectral regions that will be integrated (i.e. integrating the baseline signal) when no successful chromatographic peak could be detected in a sample.FIR
shouldn’t reasonably stretch further than the minimum and maximum bound (RT / m/z) of all found peaks across all samples for a given feature: this way no excess signal will be considered.
\(~\)
With all features integrated in all samples, Diagnostic provide tools to assess the quality of the peak integration and refine integration boundaries by setting uROI
and FIR
adapted to the specific chromatographic run being processed.
\(~\)
Annotation statistics
summarises the success in integrating each targeted feature. The ratio of peaks found (%)
, ratio of peaks filled (%)
and the average ppm error
and RT deviation (s)
will highlight a feature that wasn’t reliably integrated over a large number of samples. Visual evaluation (see below) and the adjustment of uROI
or FIR
might assist in tuning the integration of said feature.
\(~\)
Update uROI/FIR
automatically sets uROI
and FIR
for each feature based on the RT / m/z boundaries of the peaks successfully integrated.
\(~\)
Diagnostic plot
offer a visualisation of a selected feature across all samples in order of analysis. This visualisation highlights the fitting of the feature in each sample, as well as the change in RT / m/z (of the peak apex) and peak area through time. Samples can be automatically coloured based on a sample metadata column.
\(~\)
Once uROI
and FIR
are successfully set, it is possible to go back to the Run annotation tab and refit all features in all samples (Note: this will overwrite the current results).
\(~\)
If the features integration are satisfactory, View results regroups all the integration results
\(~\)
Overall results
displays a fitting property for all targeted features (as columns) and all spectra (as rows).
\(~\)
Results per targeted feature
displays all fitting properties (as columns) for all samples (as rows) for a selected targeted feature.
\(~\)
Results per sample
displays all fitting properties (as columns) for all targeted features (as rows) for a selected sample.
\(~\)
The Export tab manages the saving of input parameters, annotation results and automated reporting.
The peakPantheRAnnotation
in it’s current state can be saved as a .RData
file which can subsequently be reloaded. The .csv
files defining the current analysis can also be exported to reproduce the current processing.
All diagnostic plots from the Diagnostic
tab can be automatically saved to disk for rapid evaluation. This can be executed in parallel if a large number of plots have to be generated.
Finally each fitting property can be saved as a .csv
file with all samples as rows and all targeted features as columns. Additionally a summary table will present the integration success rate for each targeted feature.
\(~\)
If a very high number of targeted features and samples are to be processed, peakPantheR
’s command line functions are more efficient, as they can be integrated in scripts and the reporting automated.
\(~\)
#> ─ Session info ───────────────────────────────────────────────────────────────
#> setting value
#> version R version 4.2.2 (2022-10-31)
#> os Ubuntu 20.04.5 LTS
#> system x86_64, linux-gnu
#> ui X11
#> language (EN)
#> collate C
#> ctype en_US.UTF-8
#> tz America/New_York
#> date 2022-12-12
#> pandoc 2.5 @ /usr/bin/ (via rmarkdown)
#>
#> ─ Packages ───────────────────────────────────────────────────────────────────
#> package * version date (UTC) lib source
#> affy 1.76.0 2022-12-12 [2] Bioconductor
#> affyio 1.68.0 2022-12-12 [2] Bioconductor
#> assertthat 0.2.1 2019-03-21 [2] CRAN (R 4.2.2)
#> Biobase * 2.58.0 2022-12-12 [2] Bioconductor
#> BiocGenerics * 0.44.0 2022-12-12 [2] Bioconductor
#> BiocManager 1.30.19 2022-10-25 [2] CRAN (R 4.2.2)
#> BiocParallel * 1.32.4 2022-12-12 [2] Bioconductor
#> BiocStyle * 2.26.0 2022-12-12 [2] Bioconductor
#> bitops 1.0-7 2021-04-24 [2] CRAN (R 4.2.2)
#> bookdown 0.30 2022-11-09 [2] CRAN (R 4.2.2)
#> bslib 0.4.1 2022-11-02 [2] CRAN (R 4.2.2)
#> cachem 1.0.6 2021-08-19 [2] CRAN (R 4.2.2)
#> callr 3.7.3 2022-11-02 [2] CRAN (R 4.2.2)
#> cli 3.4.1 2022-09-23 [2] CRAN (R 4.2.2)
#> clue 0.3-63 2022-11-19 [2] CRAN (R 4.2.2)
#> cluster 2.1.4 2022-08-22 [2] CRAN (R 4.2.2)
#> codetools 0.2-18 2020-11-04 [2] CRAN (R 4.2.2)
#> colorspace 2.0-3 2022-02-21 [2] CRAN (R 4.2.2)
#> crayon 1.5.2 2022-09-29 [2] CRAN (R 4.2.2)
#> DBI 1.1.3 2022-06-18 [2] CRAN (R 4.2.2)
#> DelayedArray 0.24.0 2022-12-12 [2] Bioconductor
#> DEoptimR 1.0-11 2022-04-03 [2] CRAN (R 4.2.2)
#> devtools 2.4.5 2022-10-11 [2] CRAN (R 4.2.2)
#> digest 0.6.31 2022-12-11 [2] CRAN (R 4.2.2)
#> doParallel * 1.0.17 2022-02-07 [2] CRAN (R 4.2.2)
#> dplyr 1.0.10 2022-09-01 [2] CRAN (R 4.2.2)
#> DT 0.26 2022-10-19 [2] CRAN (R 4.2.2)
#> ellipsis 0.3.2 2021-04-29 [2] CRAN (R 4.2.2)
#> evaluate 0.18 2022-11-07 [2] CRAN (R 4.2.2)
#> faahKO * 1.38.0 2022-12-08 [2] Bioconductor
#> fansi 1.0.3 2022-03-24 [2] CRAN (R 4.2.2)
#> farver 2.1.1 2022-07-06 [2] CRAN (R 4.2.2)
#> fastmap 1.1.0 2021-01-25 [2] CRAN (R 4.2.2)
#> foreach * 1.5.2 2022-02-02 [2] CRAN (R 4.2.2)
#> fs 1.5.2 2021-12-08 [2] CRAN (R 4.2.2)
#> generics 0.1.3 2022-07-05 [2] CRAN (R 4.2.2)
#> GenomeInfoDb 1.34.4 2022-12-12 [2] Bioconductor
#> GenomeInfoDbData 1.2.9 2022-11-08 [2] Bioconductor
#> GenomicRanges 1.50.1 2022-12-12 [2] Bioconductor
#> ggplot2 3.4.0 2022-11-04 [2] CRAN (R 4.2.2)
#> glue 1.6.2 2022-02-24 [2] CRAN (R 4.2.2)
#> gtable 0.3.1 2022-09-01 [2] CRAN (R 4.2.2)
#> highr 0.9 2021-04-16 [2] CRAN (R 4.2.2)
#> htmltools 0.5.4 2022-12-07 [2] CRAN (R 4.2.2)
#> htmlwidgets 1.5.4 2021-09-08 [2] CRAN (R 4.2.2)
#> httpuv 1.6.6 2022-09-08 [2] CRAN (R 4.2.2)
#> impute 1.72.1 2022-12-12 [2] Bioconductor
#> IRanges 2.32.0 2022-12-12 [2] Bioconductor
#> iterators * 1.0.14 2022-02-05 [2] CRAN (R 4.2.2)
#> jquerylib 0.1.4 2021-04-26 [2] CRAN (R 4.2.2)
#> jsonlite 1.8.4 2022-12-06 [2] CRAN (R 4.2.2)
#> knitr 1.41 2022-11-18 [2] CRAN (R 4.2.2)
#> labeling 0.4.2 2020-10-20 [2] CRAN (R 4.2.2)
#> later 1.3.0 2021-08-18 [2] CRAN (R 4.2.2)
#> lattice 0.20-45 2021-09-22 [2] CRAN (R 4.2.2)
#> lifecycle 1.0.3 2022-10-07 [2] CRAN (R 4.2.2)
#> limma 3.54.0 2022-12-12 [2] Bioconductor
#> magick 2.7.3 2021-08-18 [2] CRAN (R 4.2.2)
#> magrittr 2.0.3 2022-03-30 [2] CRAN (R 4.2.2)
#> MALDIquant 1.22 2022-11-10 [2] CRAN (R 4.2.2)
#> MASS 7.3-58.1 2022-08-03 [2] CRAN (R 4.2.2)
#> MassSpecWavelet 1.64.0 2022-12-12 [2] Bioconductor
#> Matrix 1.5-3 2022-11-11 [2] CRAN (R 4.2.2)
#> MatrixGenerics 1.10.0 2022-12-12 [2] Bioconductor
#> matrixStats 0.63.0 2022-11-18 [2] CRAN (R 4.2.2)
#> memoise 2.0.1 2021-11-26 [2] CRAN (R 4.2.2)
#> mime 0.12 2021-09-28 [2] CRAN (R 4.2.2)
#> miniUI 0.1.1.1 2018-05-18 [2] CRAN (R 4.2.2)
#> minpack.lm 1.2-2 2022-04-13 [2] CRAN (R 4.2.2)
#> MsCoreUtils 1.10.0 2022-12-12 [2] Bioconductor
#> MsFeatures 1.6.0 2022-12-12 [2] Bioconductor
#> MSnbase * 2.24.0 2022-12-12 [2] Bioconductor
#> munsell 0.5.0 2018-06-12 [2] CRAN (R 4.2.2)
#> mzID 1.36.0 2022-12-12 [2] Bioconductor
#> mzR * 2.32.0 2022-12-12 [2] Bioconductor
#> ncdf4 1.20 2022-12-03 [2] CRAN (R 4.2.2)
#> pander * 0.6.5 2022-03-18 [2] CRAN (R 4.2.2)
#> pcaMethods 1.90.0 2022-12-12 [2] Bioconductor
#> peakPantheR * 1.12.2 2022-12-12 [1] Bioconductor
#> pillar 1.8.1 2022-08-19 [2] CRAN (R 4.2.2)
#> pkgbuild 1.4.0 2022-11-27 [2] CRAN (R 4.2.2)
#> pkgconfig 2.0.3 2019-09-22 [2] CRAN (R 4.2.2)
#> pkgload 1.3.2 2022-11-16 [2] CRAN (R 4.2.2)
#> plyr 1.8.8 2022-11-11 [2] CRAN (R 4.2.2)
#> pracma 2.4.2 2022-09-22 [2] CRAN (R 4.2.2)
#> preprocessCore 1.60.0 2022-12-12 [2] Bioconductor
#> prettyunits 1.1.1 2020-01-24 [2] CRAN (R 4.2.2)
#> processx 3.8.0 2022-10-26 [2] CRAN (R 4.2.2)
#> profvis 0.3.7 2020-11-02 [2] CRAN (R 4.2.2)
#> promises 1.2.0.1 2021-02-11 [2] CRAN (R 4.2.2)
#> ProtGenerics * 1.30.0 2022-12-12 [2] Bioconductor
#> ps 1.7.2 2022-10-26 [2] CRAN (R 4.2.2)
#> purrr 0.3.5 2022-10-06 [2] CRAN (R 4.2.2)
#> R6 2.5.1 2021-08-19 [2] CRAN (R 4.2.2)
#> RANN 2.6.1 2019-01-08 [2] CRAN (R 4.2.2)
#> RColorBrewer 1.1-3 2022-04-03 [2] CRAN (R 4.2.2)
#> Rcpp * 1.0.9 2022-07-08 [2] CRAN (R 4.2.2)
#> RCurl 1.98-1.9 2022-10-03 [2] CRAN (R 4.2.2)
#> remotes 2.4.2 2021-11-30 [2] CRAN (R 4.2.2)
#> rlang 1.0.6 2022-09-24 [2] CRAN (R 4.2.2)
#> rmarkdown 2.18 2022-11-09 [2] CRAN (R 4.2.2)
#> robustbase 0.95-0 2022-04-02 [2] CRAN (R 4.2.2)
#> S4Vectors * 0.36.1 2022-12-12 [2] Bioconductor
#> sass 0.4.4 2022-11-24 [2] CRAN (R 4.2.2)
#> scales 1.2.1 2022-08-20 [2] CRAN (R 4.2.2)
#> sessioninfo 1.2.2 2021-12-06 [2] CRAN (R 4.2.2)
#> shiny 1.7.3 2022-10-25 [2] CRAN (R 4.2.2)
#> shinycssloaders 1.0.0 2020-07-28 [2] CRAN (R 4.2.2)
#> stringi 1.7.8 2022-07-11 [2] CRAN (R 4.2.2)
#> stringr 1.5.0 2022-12-02 [2] CRAN (R 4.2.2)
#> SummarizedExperiment 1.28.0 2022-12-12 [2] Bioconductor
#> tibble 3.1.8 2022-07-22 [2] CRAN (R 4.2.2)
#> tidyselect 1.2.0 2022-10-10 [2] CRAN (R 4.2.2)
#> urlchecker 1.0.1 2021-11-30 [2] CRAN (R 4.2.2)
#> usethis 2.1.6 2022-05-25 [2] CRAN (R 4.2.2)
#> utf8 1.2.2 2021-07-24 [2] CRAN (R 4.2.2)
#> vctrs 0.5.1 2022-11-16 [2] CRAN (R 4.2.2)
#> vsn 3.66.0 2022-12-12 [2] Bioconductor
#> withr 2.5.0 2022-03-03 [2] CRAN (R 4.2.2)
#> xcms * 3.20.0 2022-12-12 [2] Bioconductor
#> xfun 0.35 2022-11-16 [2] CRAN (R 4.2.2)
#> XML 3.99-0.13 2022-12-04 [2] CRAN (R 4.2.2)
#> xtable 1.8-4 2019-04-21 [2] CRAN (R 4.2.2)
#> XVector 0.38.0 2022-12-12 [2] Bioconductor
#> yaml 2.3.6 2022-10-18 [2] CRAN (R 4.2.2)
#> zlibbioc 1.44.0 2022-12-12 [2] Bioconductor
#>
#> [1] /tmp/Rtmpo1JCkD/Rinst5a7014f3ee241
#> [2] /home/biocbuild/bbs-3.16-bioc/R/library
#>
#> ──────────────────────────────────────────────────────────────────────────────