1
|
1 <tool id="isoDE" name="IsoDE2">
|
|
2 <description>Compute gene Differential Expression based on IsoEM2 output </description>
|
|
3 <command interpreter="bash">isoDE2.sh
|
|
4 -c1
|
|
5 $condition1
|
|
6 #for $r in $condition1replicates
|
|
7 ${r.c1Rep}
|
|
8 #end for
|
|
9 -c2
|
|
10 $condition2
|
|
11 #for $r in $condition2replicates
|
|
12 ${r.c2Rep}
|
|
13 #end for
|
|
14 -pval $pval
|
|
15 -geneFPKMout $geneFPKM
|
|
16 -geneTPMout $geneTPM
|
|
17 -isoFPKMout $isoformFPKM
|
|
18 -isoTPMout $isoformTPM
|
|
19 </command>
|
|
20
|
|
21 <inputs>
|
|
22 <param name="condition1" type="data" label="Select data for Condition 1" format="gz" help="Condition 1 isoEM2 compressed output file"/>
|
|
23 <!--
|
|
24 <param format="toolshed.gz" name="condition1" type="data" label="Select data for Condition 1" help="Condition 1 isoEM2 compressed output file"/>
|
|
25 -->
|
|
26 <repeat name="condition1replicates" title="Replicates for Condition 1">
|
|
27 <param name="c1Rep" label="Add replicate" type="data" format="gz" data_ref="condtion1" />
|
|
28 </repeat>
|
|
29
|
|
30 <param format="gz" name="condition2" type="data" label="Select data for Condition 2" help="Condition 2 isoEM2 compressed output file"/>
|
|
31 <!--
|
|
32 <param format="toolshed.gz" name="condition2" type="data" label="Select data for Condition 2" help="Condition 2 IsoEM2 compressed output file"/>
|
|
33 -->
|
|
34 <repeat name="condition2replicates" title="Replicates for Condition 2">
|
|
35 <param format="gz" name="c2Rep" label="Add replicate" type="data" data_ref="condtion2" />
|
|
36 </repeat>
|
|
37
|
|
38
|
|
39 <param name="pval" label="p-value" type="float" value="0.05" help="Desired p-value to for which a reliable fold change level will be reported" />
|
|
40
|
|
41 </inputs>
|
|
42 <outputs>
|
|
43 <data format="tabular" name="geneFPKM" label="isoDE gene fpkm" />
|
|
44 <data format="tabular" name="isoformFPKM" label="isoDE isoform fpkm" />
|
|
45 <data format="tabular" name="geneTPM" label="isoDE gene tpm" />
|
|
46 <data format="tabular" name="isoformTPM" label="isoDE isoform tpm" />
|
|
47 </outputs>
|
|
48
|
|
49 <help>
|
|
50 **What it does**
|
|
51
|
|
52 Computes gene and isoform differential expression between two conditions (example tumor and normal) for both Fragment per Kilobase of transcript length per Million
|
|
53 bases (FPKM) and Transcripts per Million (TPM) values. The computation is based on the boostraping output generated by IsoEM2. The number of bootstrap iterations for
|
|
54 IsoEM2 should be >= 20 (suggested 200).
|
|
55
|
|
56 **Input**
|
|
57
|
|
58 * - One or more IsoEM output files (compressed tar files) for each of the two conditions. More than one file can be used if there are replicated for either condition
|
|
59 * - Desired p-value to for which a reliable fold change level will be reported
|
|
60 *
|
|
61
|
|
62
|
|
63
|
|
64
|
|
65 **Output**
|
|
66
|
|
67 * four output files containinag results for Gene FPKM DE, Gene TPM DE, Isoform FPKM DE, and Isoform TPM DE. The four files have identical format with the following fields
|
|
68 * 1- Gene/isoform ID
|
|
69 * 2- Reliable log_2(FC) : conservative estimate of fold change in log base 2.
|
|
70 * For the confidence level specified by the input p-value, fold change of gene/isoform abundance (FPKM/TPM) in condition 2 compared condition 1 is
|
|
71 * at least 2 ^ absoulte value of this field.The sign indicates the direction, +ve means over expressed in condition 2, -ve means underexpressed in
|
|
72 * condition 1. NDE indicates that no change was detected.
|
|
73 * 3- log_2(condition 1 FPKM (or TPM)/condition 1 FPKM(or TPM)) based on IsoEM2 run without bootstrapping
|
|
74 * 4- condition 1 FPKM (or TPM) based on IsoEM2 run without bootstrapping
|
|
75 * 5- condition 2 FPKM (or TPM) based on IsoEM2 run without bootstrapping
|
|
76
|
|
77
|
|
78 </help>
|
|
79 </tool>
|