|
78
|
1 <tool id="MaREA RAS Generator" name="Expression2RAS" version="1.0.7">
|
|
47
|
2 <description>- Reaction Activity Scores computation</description>
|
|
45
|
3 <macros>
|
|
|
4 <import>marea_macros.xml</import>
|
|
|
5 </macros>
|
|
|
6 <requirements>
|
|
|
7 <requirement type="package" version="0.23.0">pandas</requirement>
|
|
|
8 <requirement type="package" version="1.1.0">scipy</requirement>
|
|
49
|
9 <requirement type="package" version="0.17.1">cobra</requirement>
|
|
45
|
10 <requirement type="package" version="4.2.1">lxml</requirement>
|
|
|
11 <requirement type="package" version="0.8.1">svglib</requirement>
|
|
80
|
12 <requirement type="package" version="3.5.0">reportlab</requirement>
|
|
45
|
13 </requirements>
|
|
|
14 <command detect_errors="exit_code">
|
|
|
15 <![CDATA[
|
|
47
|
16 python $__tool_directory__/ras_generator.py
|
|
|
17 --rules_selector $cond_rule.rules_selector
|
|
|
18 --input $input
|
|
|
19 --none $none
|
|
|
20 --tool_dir $__tool_directory__
|
|
|
21 --out_log $log
|
|
|
22 --ras_output $ras_output
|
|
52
|
23 #if $cond_rule.rules_selector == 'Custom'
|
|
|
24 --custom $cond_rule.Custom_rules
|
|
|
25 #end if
|
|
45
|
26 ]]>
|
|
|
27 </command>
|
|
|
28 <inputs>
|
|
|
29 <conditional name="cond_rule">
|
|
|
30 <expand macro="options"/>
|
|
|
31 <when value="HMRcore">
|
|
|
32 </when>
|
|
|
33 <when value="Recon">
|
|
|
34 </when>
|
|
|
35 <when value="Custom">
|
|
|
36 <param name="Custom_rules" type="data" format="tabular, csv, tsv, xml" label="Custom rules" />
|
|
|
37 </when>
|
|
|
38 </conditional>
|
|
47
|
39 <param name="input" argument="--input" type="data" format="tabular, csv, tsv" label="Gene Expression dataset:" />
|
|
|
40 <param name="name" argument="--name" type="text" label="Dataset's name:" value="Dataset" help="Default: Dataset" />
|
|
|
41 <param name="none" argument="--none" type="boolean" truevalue="true" falsevalue="false" checked="true" label="(A and NaN) solved as (A)?" />
|
|
45
|
42 </inputs>
|
|
|
43
|
|
|
44 <outputs>
|
|
47
|
45 <data format="txt" name="log" label="Expression2RAS - $name - Log" />
|
|
|
46 <data format="tabular" name="ras_output" label="$name RAS"/>
|
|
45
|
47 </outputs>
|
|
47
|
48
|
|
45
|
49 <help>
|
|
|
50 <![CDATA[
|
|
|
51
|
|
|
52 What it does
|
|
|
53 -------------
|
|
48
|
54
|
|
|
55 This tool computes Reaction Activity Scores from gene expression (RNA-seq) dataset(s), as described in Graudenzi et al. Integration of transcriptomic data and metabolic networks in cancer samples reveals highly significant prognostic power. Journal of Biomedical Informatics, 2018, 87: 37-49.
|
|
|
56
|
|
|
57 Accepted files:
|
|
|
58 - A gene expression dataset
|
|
|
59
|
|
|
60 Format:
|
|
|
61 Tab-separated text file reporting the normalized expression level (e.g., TPM, RPKM, ...) of each gene (row) for a given sample (column).
|
|
|
62 Column header: sample ID.
|
|
|
63 Row header: gene ID.
|
|
|
64
|
|
|
65
|
|
|
66 Optional files:
|
|
|
67 - custom GPR (Gene-Protein-Reaction) rules. Two accepted formats:
|
|
|
68
|
|
|
69 * (Cobra Toolbox and CobraPy compliant) xml of metabolic model;
|
|
|
70 * .csv file specifyig for each reaction ID (column 1) the corresponding GPR rule (column 2).
|
|
|
71
|
|
|
72 Computation option ‘(A and NaN) solved as (A)’:
|
|
|
73 In case of missing expression value, referred to as NaN (Not a Number), for a gene joined with an AND operator in a given GPR rule, the rule ‘A and NaN’
|
|
|
74
|
|
|
75 If YES is selected: the GPR will be solved as A.
|
|
|
76
|
|
|
77 If NO is selected: the GPR will be disregarded tout-court (i.e., treated as NaN).
|
|
|
78
|
|
|
79 Example input
|
|
|
80 -------------
|
|
|
81
|
|
|
82 Custom GPR rules:
|
|
|
83
|
|
|
84 +------------+--------------------------------------+
|
|
|
85 | id | rule (with entrez-id |
|
|
|
86 +============+======================================+
|
|
|
87 | r1642 | 155060 or 10357 |
|
|
|
88 +------------+--------------------------------------+
|
|
|
89 | r1643 | 155060 or 100134869 |
|
|
|
90 +------------+--------------------------------------+
|
|
|
91 | r1640 | 155060 and 100134869 or 10357 |
|
|
|
92 +------------+--------------------------------------+
|
|
|
93
|
|
|
94 RNA-seq dataset:
|
|
|
95
|
|
|
96 +------------+----------------+----------------+----------------+
|
|
|
97 | Hugo_ID | TCGAA62670 | TCGAA62671 | TCGAA62672 |
|
|
|
98 +============+================+================+================+
|
|
|
99 | HGNC:24086 | 0.523167 | 0.371355 | 0.925661 |
|
|
|
100 +------------+----------------+----------------+----------------+
|
|
|
101 | HGNC:24086 | 0.568765 | 0.765567 | 0.456789 |
|
|
|
102 +------------+----------------+----------------+----------------+
|
|
|
103 | HGNC:9876 | 0.876545 | 0.768933 | 0.987654 |
|
|
|
104 +------------+----------------+----------------+----------------+
|
|
|
105 | HGNC:9 | 0.456788 | 0.876543 | 0.876542 |
|
|
|
106 +------------+----------------+----------------+----------------+
|
|
|
107 | HGNC:23 | 0.876543 | 0.786543 | 0.897654 |
|
|
|
108 +------------+----------------+----------------+----------------+
|
|
|
109
|
|
45
|
110 ]]>
|
|
|
111 </help>
|
|
73
|
112 <expand macro="citations" />
|
|
45
|
113 </tool>
|
|
|
114
|
