|
1
|
1 <tool id="abims_CAMERA_annotateDiffreport" name="CAMERA.annotate" version="2.1.2">
|
|
0
|
2
|
|
|
3 <description>CAMERA annotate function. Returns annotation results (isotope peaks, adducts and fragments) and a diffreport if more than one condition.</description>
|
|
|
4
|
|
|
5 <requirements>
|
|
|
6 <requirement type="package" version="3.1.2">R</requirement>
|
|
|
7 <requirement type="binary">Rscript</requirement>
|
|
|
8 <requirement type="package" version="1.44.0">xcms</requirement>
|
|
|
9 <requirement type="package" version="1.22.0">camera</requirement>
|
|
|
10 <requirement type="package" version="2.1">camera_w4m_script</requirement>
|
|
|
11 </requirements>
|
|
|
12
|
|
|
13 <stdio>
|
|
|
14 <exit_code range="1:" level="fatal" />
|
|
|
15 </stdio>
|
|
|
16
|
|
1
|
17 <command><![CDATA[
|
|
0
|
18 CAMERA.r
|
|
1
|
19 xfunction annotatediff
|
|
|
20 image $image
|
|
|
21
|
|
|
22 variableMetadataOutput $variableMetadata
|
|
|
23 dataMatrixOutput $datamatrix
|
|
|
24 xsetRdataOutput $rdata
|
|
|
25
|
|
|
26 convert_param $convert_param num_digits $num_digits
|
|
|
27 nSlaves \${GALAXY_SLOTS:-1}
|
|
0
|
28 sigma $sigma perfwhm $perfwhm
|
|
|
29 maxcharge $maxcharge maxiso $maxiso minfrac $minfrac
|
|
|
30 ppm $ppm mzabs $mzabs intval $intval
|
|
|
31 max_peaks $max_peaks
|
|
|
32 quick $quick_block.quick
|
|
|
33 #if $quick_block.quick == "FALSE"
|
|
|
34 polarity $quick_block.polarity
|
|
|
35 cor_eic_th $quick_block.cor_eic_th
|
|
|
36 graphMethod $quick_block.graphMethod
|
|
|
37 pval $quick_block.pval
|
|
|
38 calcCiS $quick_block.calcCiS
|
|
|
39 calcIso $quick_block.calcIso
|
|
|
40 calcCaS $quick_block.calcCaS
|
|
|
41 multiplier $quick_block.multiplier
|
|
|
42 #end if
|
|
|
43
|
|
|
44 #if $options.option == "show":
|
|
|
45 eicmax $options.eicmax
|
|
|
46 eicwidth $options.eicwidth
|
|
|
47 value $options.value
|
|
|
48 sortpval $options.sortpval
|
|
|
49 h $options.h
|
|
|
50 w $options.w
|
|
|
51 mzdec $options.mzdec
|
|
1
|
52 new_file_path ${__new_file_path__}/primary_${output_diffreport.id}_
|
|
0
|
53 #else
|
|
|
54 eicmax 0
|
|
|
55 eicwidth 200
|
|
|
56 value "into"
|
|
|
57 sortpval TRUE
|
|
|
58 h 480
|
|
|
59 w 640
|
|
|
60 mzdec 2
|
|
1
|
61 new_file_path ${__new_file_path__}/primary_${output_diffreport.id}_
|
|
0
|
62 #end if
|
|
|
63
|
|
|
64
|
|
1
|
65 ]]></command>
|
|
0
|
66
|
|
|
67 <inputs>
|
|
|
68
|
|
|
69
|
|
|
70 <param name="image" type="data" label="RData file" format="rdata.xcms.fillpeaks,rdata" help="output file from another function xcms (fillPeaks)" />
|
|
|
71
|
|
|
72 <param name="convert_param" type="boolean" checked="false" truevalue="TRUE" falsevalue="FALSE" label="Convert retention time (seconds) into minutes" help="Convert the columns rtmed, rtmin and rtmax into minutes"/>
|
|
1
|
73 <param name="num_digits" type="integer" value="4" label="" help="Number of decimal places for mass values reported in ions identifiers. A minimum of 4 decimal places should be set" />
|
|
0
|
74
|
|
|
75 <param name="sigma" type="integer" value="6" label="groupFWHM: multiplier of the standard deviation" help="[sigma]" />
|
|
|
76 <param name="perfwhm" type="float" value="0.6" max="1" min="0" label="groupFWHM: percentage of FWHM width" help="[perfwhm]" />
|
|
|
77 <param name="maxcharge" type="integer" value="3" label="findIsotopes: max. ion charge" help="[maxcharge]" />
|
|
|
78 <param name="maxiso" type="integer" value="4" label="findIsotopes: max. number of expected isotopes" help="[maxiso]" />
|
|
|
79 <param name="minfrac" type="float" value="0.5" max="1" min="0" label="findIsotopes: The percentage number of samples, which must satisfy the C12/C13 rule for isotope annotation" help="[minfrac]" />
|
|
|
80 <param name="ppm" type="integer" value="5" label="General ppm error" help="[ppm]" />
|
|
|
81 <param name="mzabs" type="float" value="0.015" label="mzabs" help="General absolut error in m/z" />
|
|
|
82 <param name="intval" type="select" label="General used intensity value" help="[intval] See the help section below">
|
|
|
83 <option value="into" selected="true">into</option>
|
|
|
84 <option value="maxo" >maxo</option>
|
|
|
85 <option value="intb">intb</option>
|
|
|
86 </param>
|
|
|
87
|
|
|
88
|
|
|
89 <param name="max_peaks" type="integer" value="100" label="How much peaks will be calculated in every thread using the parallel mode" help="[max_peaks]" />
|
|
|
90 <conditional name="quick_block">
|
|
|
91 <param name="quick" type="select" label="Quick mode" help="[quick] If TRUE, use only groupFWHM and findIsotopes functions. Else if FALSE, use also groupCorr and findAdducts">
|
|
|
92 <option value="TRUE">TRUE</option>
|
|
|
93 <option value="FALSE" selected="true">FALSE</option>
|
|
|
94 </param>
|
|
|
95 <when value="FALSE">
|
|
|
96 <param name="polarity" type="select" label="findAdducts: Which polarity mode was used for measuring of the ms sample" help="polarity">
|
|
|
97 <option value="positive" >positive</option>
|
|
|
98 <option value="negative" selected="true">negative</option>
|
|
|
99 </param>
|
|
|
100 <param name="cor_eic_th" type="float" value="0.75" max="1" min="0" label="groupCorr: correlation threshold (0..1)" help="[cor_eic_th]" />
|
|
|
101 <param name="graphMethod" type="select" label="groupCorr: Method selection for grouping peaks after correlation analysis into pseudospectra" help="[graphMethod]">
|
|
|
102 <option value="hcs" selected="true">hcs</option>
|
|
|
103 <option value="lpc">lpc</option>
|
|
|
104 </param>
|
|
|
105 <param name="pval" type="float" value="0.05" max="1" min="0" label="groupCorr: significant correlation threshold" help="[pval]" />
|
|
|
106 <param name="calcCiS" type="boolean" checked="true" truevalue="TRUE" falsevalue="FALSE" label="groupCorr: Use correlation inside samples for peak grouping" help="[calcCiS]"/>
|
|
|
107 <param name="calcIso" type="boolean" checked="false" truevalue="TRUE" falsevalue="FALSE" label="groupCorr: Use isotopic relationship for peak grouping" help="[calcIso]"/>
|
|
|
108 <param name="calcCaS" type="boolean" checked="false" truevalue="TRUE" falsevalue="FALSE" label="groupCorr: Use correlation across samples for peak grouping" help="[calcCaS]"/>
|
|
|
109 <param name="multiplier" type="integer" value="3" label="findAdducts: If no ruleset is provided, calculate ruleset with max. number n of [nM+x] clusterions" help="[multiplier]" />
|
|
|
110 </when>
|
|
|
111 <when value="TRUE">
|
|
|
112 <param name="polarity" type="hidden" value="quick" label="for the output label" help="for the output label" />
|
|
|
113 </when>
|
|
|
114 </conditional>
|
|
|
115
|
|
|
116
|
|
|
117 <!-- Annotatediffreport specific parameters -->
|
|
|
118 <conditional name="options">
|
|
|
119 <param name="option" type="select" label="Number of condition">
|
|
|
120 <option value="hide" selected="true">One condition</option>
|
|
|
121 <option value="show">Two or more conditions</option>
|
|
|
122 </param>
|
|
|
123 <when value="show">
|
|
|
124 <param name="eicmax" type="integer" value="0" label="diffreport: number of the most significantly different analytes to create EICs for" help="[eicmax]" />
|
|
|
125 <param name="eicwidth" type="integer" value="200" label="diffreport: width (in seconds) of EICs produced" help="[eicwidth]" />
|
|
|
126 <param name="value" type="select" label="diffreport: Intensity values to be used for the diffreport" help="[value]">
|
|
|
127 <option value="into" selected="true">into</option>
|
|
|
128 <option value="maxo" >maxo</option>
|
|
|
129 <option value="intb">intb</option>
|
|
|
130 </param>
|
|
|
131 <param name="h" type="integer" value="480" label="diffreport: Numeric variable for the height of the eic and boxplots that are printed out" help="[height]" />
|
|
|
132 <param name="w" type="integer" value="640" label="diffreport: Numeric variable for the width of the eic and boxplots print out made" help="[width]" />
|
|
|
133 <param name="mzdec" type="integer" value="2" label="diffreport: Number of decimal places of title m/z values in the eic plot" help="[mzdec]" />
|
|
|
134 <param name="sortpval" type="boolean" checked="false" truevalue="TRUE" falsevalue="FALSE" label="diffreport: logical indicating whether the reports should be sorted by p-value" help="[sortpval]"/>
|
|
|
135 </when>
|
|
|
136 <when value="hide">
|
|
|
137 </when>
|
|
|
138 </conditional>
|
|
|
139
|
|
|
140 </inputs>
|
|
|
141
|
|
|
142 <outputs>
|
|
1
|
143 <data name="variableMetadata" format="tabular" label="${image.name[:-6]}.annotate.variableMetadata.tsv" />
|
|
0
|
144 <data name="datamatrix" format="tabular" label="${image.name[:-6]}.annotate.dataMatrix.tsv" />
|
|
|
145 <data name="rdata" format="rdata.camera.quick" label="${image.name[:-6]}.annotate.${$quick_block.polarity}.Rdata" >
|
|
|
146 <change_format>
|
|
|
147 <when input="quick_block.polarity" value="positive" format="rdata.camera.positive" />
|
|
|
148 <when input="quick_block.polarity" value="negative" format="rdata.camera.negative" />
|
|
|
149 </change_format>
|
|
1
|
150 </data>
|
|
|
151 <data name="output_diffreport" format="text" label="${image.name[:-6]}.annotateDiffreport" hidden="True">
|
|
|
152 <filter>(options['option'] == 'show')</filter>
|
|
|
153 </data>
|
|
0
|
154 </outputs>
|
|
|
155
|
|
|
156 <tests>
|
|
|
157 <test>
|
|
|
158 <param name="image" value="xset.group.retcor.group.fillPeaks.RData"/>
|
|
|
159 <param name="convert_param" value="False"/>
|
|
|
160 <param name="num_digits" value="0"/>
|
|
|
161 <param name="sigma" value="6"/>
|
|
|
162 <param name="perfwhm" value="0.6"/>
|
|
|
163 <param name="maxcharge" value="3"/>
|
|
|
164 <param name="maxiso" value="4"/>
|
|
|
165 <param name="minfrac" value="0.5"/>
|
|
|
166 <param name="ppm" value="5"/>
|
|
|
167 <param name="mzabs" value="0.015"/>
|
|
|
168 <param name="intval" value="into"/>
|
|
|
169 <param name="max_peaks" value="100"/>
|
|
|
170 <param name="quick_block.quick" value="FALSE"/>
|
|
|
171 <param name="quick_block.polarity" value="negative"/>
|
|
|
172 <param name="quick_block.cor_eic_th" value="0.75"/>
|
|
|
173 <param name="quick_block.graphMethod" value="hcs"/>
|
|
|
174 <param name="quick_block.pval" value="0.05"/>
|
|
|
175 <param name="quick_block.calcCiS" value="True"/>
|
|
|
176 <param name="quick_block.calcIso" value="False"/>
|
|
|
177 <param name="quick_block.calcCaS" value="False"/>
|
|
|
178 <param name="quick_block.multiplier" value="3"/>
|
|
|
179 <param name="options.option" value="show"/>
|
|
|
180 <param name="options.eicmax" value="200"/>
|
|
|
181 <param name="options.eicwidth" value="200"/>
|
|
|
182 <param name="options.value" value="into"/>
|
|
|
183 <param name="options.h" value="480"/>
|
|
|
184 <param name="options.w" value="640"/>
|
|
|
185 <param name="options.mzdec" value="2"/>
|
|
|
186 <param name="options.sortpval" value="False"/>
|
|
|
187 <output name="variableMetadata" file="xset.group.retcor.group.fillPeaks.annotate.variableMetadata.tsv" />
|
|
|
188 <output name="datamatrix" file="xset.group.retcor.group.fillPeaks.annotate.dataMatrix.tsv" />
|
|
|
189 <output name="rdata" file="xset.group.retcor.group.fillPeaks.annotate.negative.Rdata" />
|
|
|
190 </test>
|
|
|
191 </tests>
|
|
|
192
|
|
|
193
|
|
|
194
|
|
1
|
195 <help><![CDATA[
|
|
0
|
196
|
|
|
197 .. class:: infomark
|
|
|
198
|
|
|
199 **Authors** Carsten Kuhl ckuhl@ipb-halle.de, Ralf Tautenhahn rtautenh@scripps.edu, Steffen Neumann sneumann@@ipb-halle.de
|
|
|
200
|
|
|
201 .. class:: infomark
|
|
|
202
|
|
|
203 **Galaxy integration** ABiMS TEAM - UPMC/CNRS - Station biologique de Roscoff and Yann Guitton yann.guitton@univ-nantes.fr - part of Workflow4Metabolomics.org [W4M]
|
|
|
204
|
|
|
205 | Contact support@workflow4metabolomics.org for any questions or concerns about the Galaxy implementation of this tool.
|
|
|
206
|
|
|
207 ---------------------------------------------------
|
|
|
208
|
|
|
209
|
|
|
210 ================
|
|
|
211 CAMERA.annotate
|
|
|
212 ================
|
|
|
213
|
|
|
214 -----------
|
|
|
215 Description
|
|
|
216 -----------
|
|
|
217
|
|
|
218 The R-package CAMERA is a Collection of Algorithms for MEtabolite
|
|
|
219 pRofile Annotation. Its primary purpose is the annotation and evaluation of
|
|
|
220 LC-MS data. It includes algorithms for annotation of isotope peaks, adducts
|
|
|
221 and fragments in peak lists. Additional methods cluster mass signals that
|
|
|
222 originate from a single metabolite, based on rules for mass differences and
|
|
|
223 peak shape comparison. To use the strength of already existing programs,
|
|
|
224 CAMERA is designed to interact directly with processed peak data from the
|
|
|
225 R-package **xcms**.
|
|
|
226
|
|
|
227 **What it does?**
|
|
|
228
|
|
|
229 The CAMERA annotation procedure can be split into two parts: We want to answer the questions which peaks occur from the same molecule and secondly compute its exact mass and annotate the ion species. Therefore CAMERA annotation workflow contains following primary functions: 1. peak grouping after retention time (**groupFWHM**) 2. peak group verification with peakshape correlation (**groupCorr**) Both methods separate peaks into different groups, which we define as ”pseu- dospectra”. Those pseudospectra can consists from one up to 100 ions, de- pending on the molecules amount and ionizability. Afterwards the exposure of the ion species can be performed with: 2 1. annotation of possible isotopes (**findIsotopes**) 2. annotation of adducts and calculating hypothetical masses for the group (**findAdducts**) This workflow results in a data-frame similar to a xcms peak table, that can be easily stored in a comma separated table .csv (Excel-readable).
|
|
|
230
|
|
|
231 If you have two or more conditions, it will return a diffreport result within the annotation results.
|
|
|
232 The diffreport result shows the most significant differences between two sets of samples. Optionally
|
|
|
233 create extracted ion chromatograms for the most significant differences.
|
|
|
234
|
|
|
235
|
|
|
236 -----------------
|
|
|
237 Workflow position
|
|
|
238 -----------------
|
|
|
239
|
|
|
240
|
|
|
241 **Upstream tools**
|
|
|
242
|
|
|
243 ========================= ==================== ====================== ==========
|
|
|
244 Name output file format parameter
|
|
|
245 ========================= ==================== ====================== ==========
|
|
|
246 xcms.fillPeaks xset.fillPeaks.RData rdata.xcms.fillpeaks RData file
|
|
|
247 ========================= ==================== ====================== ==========
|
|
|
248
|
|
|
249
|
|
|
250 **Downstream tools**
|
|
|
251
|
|
|
252 +---------------------------+---------------------------------------+------------------------------------------------+
|
|
|
253 | Name | Output file | Format |
|
|
|
254 +===========================+=======================================+================================================+
|
|
|
255 |CAMERA_combinexsAnnot |xset.annotate.dataMatrix.tsv | rdata.camera.positive or rdata.camera.negative |
|
|
|
256 +---------------------------+---------------------------------------+------------------------------------------------+
|
|
|
257 |Determine Vdk or Lowess |xset.annotate.dataMatrix.tsv | Tabular |
|
|
|
258 +---------------------------+---------------------------------------+------------------------------------------------+
|
|
|
259 |Normalization Vdk/Lowess |xset.annotate.dataMatrix.tsv | Tabular |
|
|
|
260 +---------------------------+---------------------------------------+------------------------------------------------+
|
|
|
261 |Anova |xset.annotate.dataMatrix.tsv | Tabular |
|
|
|
262 +---------------------------+---------------------------------------+------------------------------------------------+
|
|
|
263 |PCA |xset.annotate.dataMatrix.tsv | Tabular |
|
|
|
264 +---------------------------+---------------------------------------+------------------------------------------------+
|
|
|
265 |Hierarchical Clustering |xset.annotate.dataMatrix.tsv | Tabular |
|
|
|
266 +---------------------------+---------------------------------------+------------------------------------------------+
|
|
|
267
|
|
|
268
|
|
|
269 The output file **xset.annotate.dataMatrix.tsv** is an tabular file. You can continue your analysis using it in the following tools:
|
|
|
270
|
|
|
271 | Determine Vdk or Lowess
|
|
|
272 | Normalization Vdk/Lowess
|
|
|
273 | Anova
|
|
|
274 | PCA
|
|
|
275 | Hierarchical Clustering
|
|
|
276
|
|
|
277
|
|
|
278 **General schema of the metabolomic workflow**
|
|
|
279
|
|
|
280 .. image:: annotate_workflow.png
|
|
|
281
|
|
|
282 -----------
|
|
|
283 Input files
|
|
|
284 -----------
|
|
|
285
|
|
|
286 +---------------------------+---------------------------+
|
|
|
287 | Parameter : num + label | Format |
|
|
|
288 +===========================+===========================+
|
|
|
289 | 1 : RData file | rdata.xcms.fillpeaks |
|
|
|
290 +---------------------------+---------------------------+
|
|
|
291
|
|
|
292
|
|
|
293 ----------
|
|
|
294 Parameters
|
|
|
295 ----------
|
|
|
296
|
|
|
297 diffreport: Intensity values to be used for the diffreport
|
|
|
298 ----------------------------------------------------------
|
|
|
299
|
|
|
300 | If **value="into"**, integrated peak intensities are used.
|
|
|
301 | If **value="maxo"**, maximum peak intensities are used.
|
|
|
302 | If **value="intb"**, baseline corrected integrated peak intensities are used (only available if peak detection was done by ‘findPeaks.centWave’).
|
|
|
303
|
|
|
304
|
|
|
305
|
|
|
306 ------------
|
|
|
307 Output files
|
|
|
308 ------------
|
|
|
309
|
|
|
310 xset.annotate.variableMetadata.tsv
|
|
|
311
|
|
|
312 | For each metabolite (row) :
|
|
|
313 | the value of the intensity in each sample, fold, anova, mzmed, mzmin, mzmax, rtmed, rtmin, rtmax, npeaks, isotopes, adduct and pcgroup
|
|
|
314
|
|
|
315 xset.annotate.dataMatrix.tsv
|
|
|
316
|
|
|
317 | A tabular file which represents for each metabolite (row), the value of the intensity in each sample (column).
|
|
|
318
|
|
|
319 xset.annotate.zip
|
|
|
320
|
|
|
321 | It contains filebase_eic, filebase_box and filebase.tsv for one conditon vs another (Anova analysis).
|
|
|
322
|
|
|
323 xset.annotate.Rdata rdata.camera.quick or rdata.camera.positive or rdata.camera.negative
|
|
|
324
|
|
|
325 | Rdata file, that be used outside Galaxy in R.
|
|
|
326
|
|
|
327 ------
|
|
|
328
|
|
|
329 .. class:: infomark
|
|
|
330
|
|
|
331 The output **"xset.annotate.dataMatrix.tsv"** is a tabular file. You can continue your analysis using it in the following tools of the workflow:
|
|
|
332
|
|
|
333 | Determine Vdk or Lowess (Data correction)
|
|
|
334 | Normalization Vdk/Lowess (Data correction)
|
|
|
335 | Anova (Statistical analysis)
|
|
|
336 | PCA (Statistical analysis)
|
|
|
337 | Hierarchical Clustering (Statistical analysis)
|
|
|
338
|
|
|
339
|
|
|
340 ---------------------------------------------------
|
|
|
341
|
|
|
342 ---------------
|
|
|
343 Working example
|
|
|
344 ---------------
|
|
|
345
|
|
|
346 Input files
|
|
|
347 -----------
|
|
|
348
|
|
|
349 | RData file -> **xset.fillPeaks.RData**
|
|
|
350
|
|
|
351 Parameters
|
|
|
352 ----------
|
|
|
353
|
|
|
354 | sortpval -> **false**
|
|
|
355 | sigma -> **6 (default)**
|
|
|
356 | perfwhm -> **0.6 (default)**
|
|
|
357 | maxcharge -> **3 (default)**
|
|
|
358 | maxiso -> **4 (default)**
|
|
|
359 | minfrac -> **0.5 (default)**
|
|
|
360 | ppm -> **500**
|
|
|
361 | mzabs -> **0.015(default)**
|
|
|
362 | Advanced options -> **hide**
|
|
|
363 | Numver of condition -> Two or more conditioons
|
|
|
364 | eicmax -> **200**
|
|
|
365 | eicwidth -> **200 (default)**
|
|
|
366
|
|
|
367
|
|
|
368
|
|
|
369 Output files
|
|
|
370 ------------
|
|
|
371
|
|
|
372 **Example of an xset.annotate.dataMatrix.tsv output**
|
|
|
373
|
|
|
374 .. image:: annotatediffreport_dataMatrix.png
|
|
|
375
|
|
|
376 **Example of a part of xset.annotate.variableMetadata.tsv output**
|
|
|
377
|
|
|
378 .. image:: annotatediffreport_variableMetadata.png
|
|
|
379
|
|
1
|
380 ]]></help>
|
|
0
|
381
|
|
|
382 <citations>
|
|
|
383 <citation type="doi">10.1021/ac202450g</citation>
|
|
|
384 <citation type="doi">10.1093/bioinformatics/btu813</citation>
|
|
|
385 </citations>
|
|
|
386
|
|
|
387
|
|
|
388 </tool>
|
|
|
389
|
|
|
390
|
