|
6
|
1 <tool id="GetDifferentialExpression" name="get differential expression">
|
|
|
2 <description>Get the differential expression between 2 conditions using Fisher's exact test, on regions defined by a third file.</description>
|
|
|
3 <command interpreter="python">
|
|
|
4 ../Java/Python/GetDifferentialExpression.py -i $formatType.inputFileName1
|
|
|
5 #if $formatType.FormatInputFileName1 == 'bed':
|
|
|
6 -f bed
|
|
|
7 #elif $formatType.FormatInputFileName1 == 'gff':
|
|
|
8 -f gff
|
|
|
9 #elif $formatType.FormatInputFileName1 == 'gff2':
|
|
|
10 -f gff2
|
|
|
11 #elif $formatType.FormatInputFileName1 == 'gff3':
|
|
|
12 -f gff3
|
|
|
13 #elif $formatType.FormatInputFileName1 == 'sam':
|
|
|
14 -f sam
|
|
|
15 #elif $formatType.FormatInputFileName1 == 'gtf':
|
|
|
16 -f gtf
|
|
|
17 #end if
|
|
|
18
|
|
|
19 -j $formatType2.inputFileName2
|
|
|
20 #if $formatType2.FormatInputFileName2 == 'bed':
|
|
|
21 -g bed
|
|
|
22 #elif $formatType2.FormatInputFileName2 == 'gff':
|
|
|
23 -g gff
|
|
|
24 #elif $formatType2.FormatInputFileName2 == 'gff2':
|
|
|
25 -g gff2
|
|
|
26 #elif $formatType2.FormatInputFileName2 == 'gff3':
|
|
|
27 -g gff3
|
|
|
28 #elif $formatType2.FormatInputFileName2 == 'sam':
|
|
|
29 -g sam
|
|
|
30 #elif $formatType2.FormatInputFileName2 == 'gtf':
|
|
|
31 -g gtf
|
|
|
32 #end if
|
|
|
33
|
|
|
34 -k $formatTypeRef.inputFileNameRef
|
|
|
35 #if $formatTypeRef.FormatInputFileNameRef == 'bed':
|
|
|
36 -l bed
|
|
|
37 #elif $formatTypeRef.FormatInputFileNameRef == 'gff':
|
|
|
38 -l gff
|
|
|
39 #elif $formatTypeRef.FormatInputFileNameRef == 'gff2':
|
|
|
40 -l gff2
|
|
|
41 #elif $formatTypeRef.FormatInputFileNameRef == 'gff3':
|
|
|
42 -l gff3
|
|
|
43 #elif $formatTypeRef.FormatInputFileNameRef == 'sam':
|
|
|
44 -l sam
|
|
|
45 #elif $formatTypeRef.FormatInputFileNameRef == 'gtf':
|
|
|
46 -l gtf
|
|
|
47 #end if
|
|
|
48
|
|
|
49 -o $outputFileGff
|
|
|
50
|
|
|
51 $simple
|
|
|
52 $adjusted
|
|
|
53
|
|
|
54 #if $optionSimplePara.simplePara == 'Yes':
|
|
|
55 -S $optionSimplePara.paraValue
|
|
|
56 #end if
|
|
|
57
|
|
|
58 #if $optionFixedSizeFactor.FSF == 'Yes':
|
|
|
59 -x $optionFixedSizeFactor.FSFValue
|
|
|
60 #end if
|
|
|
61
|
|
|
62 #if $optionFDR.FDR == 'Yes':
|
|
|
63 -d $optionFDR.FDRValue
|
|
|
64 #end if
|
|
|
65 </command>
|
|
|
66
|
|
|
67 <inputs>
|
|
|
68 <conditional name="formatType">
|
|
|
69 <param name="FormatInputFileName1" type="select" label="Input File Format 1">
|
|
|
70 <option value="bed">bed</option>
|
|
|
71 <option value="gff">gff</option>
|
|
|
72 <option value="gff2">gff2</option>
|
|
|
73 <option value="gff3">gff3</option>
|
|
|
74 <option value="sam">sam</option>
|
|
|
75 <option value="gtf">gtf</option>
|
|
|
76 </param>
|
|
|
77 <when value="bed">
|
|
|
78 <param name="inputFileName1" format="bed" type="data" label="Input File 1"/>
|
|
|
79 </when>
|
|
|
80 <when value="gff">
|
|
|
81 <param name="inputFileName1" format="gff" type="data" label="Input File 1"/>
|
|
|
82 </when>
|
|
|
83 <when value="gff2">
|
|
|
84 <param name="inputFileName1" format="gff2" type="data" label="Input File 1"/>
|
|
|
85 </when>
|
|
|
86 <when value="gff3">
|
|
|
87 <param name="inputFileName1" format="gff3" type="data" label="Input File 1"/>
|
|
|
88 </when>
|
|
|
89 <when value="sam">
|
|
|
90 <param name="inputFileName1" format="sam" type="data" label="Input File 1"/>
|
|
|
91 </when>
|
|
|
92 <when value="gtf">
|
|
|
93 <param name="inputFileName1" format="gtf" type="data" label="Input File 1"/>
|
|
|
94 </when>
|
|
|
95 </conditional>
|
|
|
96
|
|
|
97 <conditional name="formatType2">
|
|
|
98 <param name="FormatInputFileName2" type="select" label="Input File Format 2">
|
|
|
99 <option value="bed">bed</option>
|
|
|
100 <option value="gff">gff</option>
|
|
|
101 <option value="gff2">gff2</option>
|
|
|
102 <option value="gff3">gff3</option>
|
|
|
103 <option value="sam">sam</option>
|
|
|
104 <option value="gtf">gtf</option>
|
|
|
105 </param>
|
|
|
106 <when value="bed">
|
|
|
107 <param name="inputFileName2" format="bed" type="data" label="Input File 2"/>
|
|
|
108 </when>
|
|
|
109 <when value="gff">
|
|
|
110 <param name="inputFileName2" format="gff" type="data" label="Input File 2"/>
|
|
|
111 </when>
|
|
|
112 <when value="gff2">
|
|
|
113 <param name="inputFileName2" format="gff2" type="data" label="Input File 2"/>
|
|
|
114 </when>
|
|
|
115 <when value="gff3">
|
|
|
116 <param name="inputFileName2" format="gff3" type="data" label="Input File 2"/>
|
|
|
117 </when>
|
|
|
118 <when value="sam">
|
|
|
119 <param name="inputFileName2" format="sam" type="data" label="Input File 2"/>
|
|
|
120 </when>
|
|
|
121 <when value="gtf">
|
|
|
122 <param name="inputFileName2" format="gtf" type="data" label="Input File 2"/>
|
|
|
123 </when>
|
|
|
124 </conditional>
|
|
|
125
|
|
|
126 <conditional name="formatTypeRef">
|
|
|
127 <param name="FormatInputFileNameRef" type="select" label="Input Ref File Format ">
|
|
|
128 <option value="bed">bed</option>
|
|
|
129 <option value="gff">gff</option>
|
|
|
130 <option value="gff2">gff2</option>
|
|
|
131 <option value="gff3">gff3</option>
|
|
|
132 <option value="sam">sam</option>
|
|
|
133 <option value="gtf">gtf</option>
|
|
|
134 </param>
|
|
|
135 <when value="bed">
|
|
|
136 <param name="inputFileNameRef" format="bed" type="data" label="Input Ref File"/>
|
|
|
137 </when>
|
|
|
138 <when value="gff">
|
|
|
139 <param name="inputFileNameRef" format="gff" type="data" label="Input Ref File"/>
|
|
|
140 </when>
|
|
|
141 <when value="gff2">
|
|
|
142 <param name="inputFileNameRef" format="gff2" type="data" label="Input Ref File"/>
|
|
|
143 </when>
|
|
|
144 <when value="gff3">
|
|
|
145 <param name="inputFileNameRef" format="gff3" type="data" label="Input Ref File"/>
|
|
|
146 </when>
|
|
|
147 <when value="sam">
|
|
|
148 <param name="inputFileNameRef" format="sam" type="data" label="Input Ref File"/>
|
|
|
149 </when>
|
|
|
150 <when value="gtf">
|
|
|
151 <param name="inputFileNameRef" format="gtf" type="data" label="Input Ref File"/>
|
|
|
152 </when>
|
|
|
153 </conditional>
|
|
|
154
|
|
15
|
155 <param name="simple" type="boolean" truevalue="-s" falsevalue="" checked="false" label="Normalize using the number of reads in each condition"/>
|
|
|
156 <param name="adjusted" type="boolean" truevalue="-a" falsevalue="" checked="false" label="Normalize using the number of reads of interquartile expression region"/>
|
|
6
|
157
|
|
|
158 <conditional name="optionSimplePara">
|
|
|
159 <param name="simplePara" type="select" label="provide the number of reads" >
|
|
|
160 <option value="Yes">Yes</option>
|
|
|
161 <option value="No" selected="true">No</option>
|
|
|
162 </param>
|
|
|
163 <when value="Yes">
|
|
|
164 <param name="paraValue" type="text" value="None" label="provide the number of reads" />
|
|
|
165 </when>
|
|
|
166 <when value="No">
|
|
|
167 </when>
|
|
|
168 </conditional>
|
|
|
169
|
|
|
170 <conditional name="optionFixedSizeFactor">
|
|
15
|
171 <param name="FSF" type="select" label="Give the magnification factor for the normalization using fixed size sliding windows in reference regions (leave empty for no such normalization)">
|
|
6
|
172 <option value="Yes">Yes</option>
|
|
|
173 <option value="No" selected="true">No</option>
|
|
|
174 </param>
|
|
|
175 <when value="Yes">
|
|
|
176 <param name="FSFValue" type="integer" value="0" />
|
|
|
177 </when>
|
|
|
178 <when value="No">
|
|
|
179 </when>
|
|
|
180 </conditional>
|
|
|
181
|
|
|
182 <conditional name="optionFDR">
|
|
|
183 <param name="FDR" type="select" label="use FDR">
|
|
|
184 <option value="Yes">Yes</option>
|
|
|
185 <option value="No" selected="true">No</option>
|
|
|
186 </param>
|
|
|
187 <when value="Yes">
|
|
|
188 <param name="FDRValue" type="float" value="0.0"/>
|
|
|
189 </when>
|
|
|
190 <when value="No">
|
|
|
191 </when>
|
|
|
192 </conditional>
|
|
|
193
|
|
|
194 </inputs>
|
|
|
195
|
|
|
196 <outputs>
|
|
|
197 <data name="outputFileGff" format="gff3" label="[GetDifferentialExpression]out file"/>
|
|
|
198 </outputs>
|
|
|
199
|
|
|
200 <help>
|
|
15
|
201 This tool compares two sets of data and find the differential expression. One very important component of the tool is the reference set. Actually, to use the tool, you need the two input sets of data, of course, and the reference set. The reference set is a set of genomic coordinates and, for each interval, it will count the number of feature on each sample and compute the differential expression. For each reference interval, it will output the direction of the regulation (up or down, with respect to the first input set), and a *p*-value from a Fisher exact test.
|
|
|
202
|
|
|
203 This reference set seems boring. Why not computing the differential expression without this set? The answer is: the differential expression of what? I cannot guess it. Actually, you might want to compare the expression of genes, of small RNAs, of transposable elements, of anything... So the reference set can be a list of genes, and in this case, you can compute the differential expression of genes. But you can also compute many other things.
|
|
|
204
|
|
|
205 Suppose that you cluster the data of your two input samples (you can do it with the *clusterize* and the *mergeTranscriptLists* tools). You now have a list of all the regions which are transcribed in at least one of the input samples. This can be your reference set. This reference set is interesting since you can detect the differential expression of data which is outside any annotation.
|
|
|
206
|
|
|
207 Suppose now that you clusterize using a sliding window the two input samples (you can do it with the *clusterizeBySlidingWindows* and the *mergeSlidingWindowsClusters* tools). You can now select all the regions of a given size which contain at least one read in one of the two input samples (do it with *selectByTag* and the tag **nbElements**). Again, this can be an other interesting reference set.
|
|
|
208
|
|
|
209 In most cases, the sizes of the two input samples will be different, so you should probably normalize the data, which is an available option. The ---rather crude--- normalization increases the number of data in the least populated sample and decreases the number of data in the most populated sample to the average number of data.
|
|
6
|
210 </help>
|
|
|
211 </tool>
|
