In taxonomicallyinformedannotation/tima-r: Taxonomically Informed Metabolite Annotation

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)
knitr::opts_knit$set(root.dir = "..")

This vignette describes the main steps of the annotation process.

Structural annotations of your features

For the moment, we support 3 different types of annotations:

• Internal MS^1^ exact mass-based library search
• Internal MS^2^ library search (experimental and in silico)
• SIRIUS

MS^1^-based

These annotations are of the lowest possible quality. However, they allow to annotate unusual adducts, in-source fragments thanks to different small tricks implemented. Try to really restrict the adduct list and structure-organism pairs you want to consider as possibilities explode rapidly.

library("timaR")

source(file = "inst/scripts/annotate_masses.R")

Spectral

source(file = "inst/scripts/annotate_spectra.R")

We use the spectral entropy from https://doi.org/10.1038/s41592-021-01331-z for matching.

In case, a python implementation of the spectral matching steps is also available at: https://github.com/mandelbrot-project/spectral_lib_matcher. The python version also includes other similarity measures.

source(file = "inst/scripts/prepare_annotations_spectra.R")

source(file = "inst/scripts/prepare_annotations_gnps.R")

Fingerprint-based

Sirius

As SIRIUS jobs are long to perform, we provide example SIRIUS workspaces (both SIRIUS 5 and 6). Note that spectral matches from SIRIUS are not supported for now. They have been generated on the 20 first lines of the example MGF with the following command:

```{bash eval=FALSE, include=TRUE}

this is run on SIRIUS 6

sirius \ --noCite \ --input=data/source/example_spectra_mini.mgf \ --output=data/interim/annotations/example_sirius.sirius/ \ --maxmz=800 \ config \ --AlgorithmProfile=orbitrap \ --StructureSearchDB=BIO \ --Timeout.secondsPerTree=10 \ --Timeout.secondsPerInstance=10 \ formulas \ zodiac \ fingerprints \ classes \ structures \ denovo-structures \ summaries \ --full-summary

this is run on SIRIUS 5

sirius \ --noCite \ --input data/source/example_spectra_mini.mgf \ --output data/interim/annotations/example_sirius/ \ --maxmz 800 \ config \ --AlgorithmProfile orbitrap \ --StructureSearchDB BIO \ --Timeout.secondsPerTree 10 \ --Timeout.secondsPerInstance 10 \ formula \ zodiac \ fingerprint \ structure \ compound-classes \ write-summaries \ --full-summary

These parameters were not optimized and were only used to give an example output.
If you are using the cli, do not forget to generate the summaries with the `--full-summary` option, or if you use the gui, generate them by clicking the corresponding icon.
You can get an example running:

```r
source(file = "inst/scripts/get_example_sirius.R")

The sirius workspace should have yourPattern_sirius as name and be placed in data/interim/annotations (else it will not be found).

source(file = "inst/scripts/prepare_annotations_sirius.R")

targets::tar_visnetwork(
  names = starts_with("ann_s"),
  exclude = c(
    contains("benchmark"),
    contains("par_"),
    contains("paths")
  ),
  targets_only = TRUE,
  degree_from = 8
)

Annotations are now prepared and can be used for further processing. Your features are not only informed with structural information but also, chemical class information. The latter might be corresponding or not to the chemical class of your annotated structure, depending on the consistency of your annotations.

Chemical class annotation of your features

Within our workflow, we offer a new way to attribute chemical classes to your features. It is analog to Network Annotation Propagation, but uses the edges of your network instead of the clusters. This makes more sense in our view, as also recently illustrated by CANOPUS.

We are currently also working on CANOPUS integration for chemical class annotation but this implies way heavier computations and we want to offer our users a fast solution.

Generating a network

A network is generated during the process. The edges are created based on the spectral entropy similarity calculated between your spectra (see https://doi.org/10.1038/s41592-021-01331-z).

source(file = "inst/scripts/create_edges_spectra.R")

If needed, you can get an example of what your minimal feature table should look like by running (no parameters needed):

source(file = "inst/scripts/get_example_features.R")

source(file = "inst/scripts/prepare_features_tables.R")

source(file = "inst/scripts/prepare_features_edges.R")

source(file = "inst/scripts/create_components.R")

source(file = "inst/scripts/prepare_features_components.R")

Biological source annotation

This step allows you to attribute biological source information to your features. If all your features come from a single extract, it will attribute the biological source of your extract to all your features. If you have multiple extracts aligned, it will take the n (according to your parameters) highest intensities of your aligned feature table and attribute the biological source of corresponding extracts.

source(file = "inst/scripts/prepare_taxa.R")

Filter annotations (based on retention time)

This step allows you to filter out the annotation of all the tools used, based on your own internal (experimental or predicted) retention times library. It is optional. If you do not have one, it will simply group the annotations of all tools.

source(file = "inst/scripts/filter_annotations.R")

You are almost there! See already all the steps accomplished!

targets::tar_visnetwork(
  names = matches("^ann"),
  exclude = c(
    contains("benchmark"),
    contains("par_"),
    contains("paths")
  ),
  targets_only = TRUE,
  degree_from = 8
)

We now recommend you to read the next vignette.

taxonomicallyinformedannotation/tima-r documentation built on July 17, 2024, 7:11 p.m.