|
44
|
1
|
|
|
2 <tool id="IsoRead" name="IsoRead: miR and isomiR identification and classification" version="0.1.1">
|
|
|
3 <requirements>
|
|
|
4 <requirement type="package" version="1.7">fpdf</requirement>
|
|
|
5 <requirement type="package" version="0.8">logomaker</requirement>
|
|
|
6 <requirement type="package" version="3.7.4">python</requirement>
|
|
|
7 <requirement type="package" version="1.17.3">numpy</requirement>
|
|
|
8 <requirement type="package" version="3.1.2">matplotlib</requirement>
|
|
|
9 <requirement type="package" version="1.0.3">pandas</requirement>
|
|
|
10 </requirements>
|
|
|
11 <command>
|
|
|
12 #set controls=[]
|
|
|
13 #for $input in $control#
|
|
|
14 $controls.extend([str($input.element_identifier),str($input)])
|
|
|
15 #end for#
|
|
|
16 #set treateds=[]
|
|
|
17 #for $input in $treated#
|
|
|
18 $treateds.extend([str($input.element_identifier),str($input)])
|
|
|
19 #end for#
|
|
|
20 #if $mir_input.database == "1":
|
|
|
21 #if $f.fil == "1":
|
|
|
22 #set path=$mir_input.genome1.fields.path
|
|
|
23 python -W ignore $__tool_directory__/mirbase_opt.py -con $controls -tre $treateds -analysis $analysis -tool_dir $__tool_directory__ -gen "$path" -f "$mir_input.database" -percentage "-1" -counts "-1" -name1 "$fal1" -name2 "$fal2"
|
|
|
24 #end if
|
|
|
25 #if $f.fil == "2":
|
|
|
26 #set path=$mir_input.genome1.fields.path
|
|
|
27 python -W ignore $__tool_directory__/mirbase_opt.py -con $controls -tre $treateds -analysis $analysis -tool_dir $__tool_directory__ -gen "$path" -f "$mir_input.database" -percentage "$f.fil1" -counts "$f.fil2" -name1 "$fal1" -name2 "$fal2"
|
|
|
28 #end if
|
|
|
29 #else:
|
|
|
30 #if $f.fil == "1":
|
|
|
31 #set path=$mir_input.genome2.fields.value
|
|
|
32 python -W ignore $__tool_directory__/mirgene_with_new_names.py -con $controls -tre $treateds -analysis $analysis -tool_dir $__tool_directory__ -gen "$path" -f "$mir_input.database" -percentage "-1" -counts "-1" -name1 "$fal1" -name2 "$fal2"
|
|
|
33 #end if
|
|
|
34 #if $f.fil == "2":
|
|
|
35 #set path=$mir_input.genome2.fields.value
|
|
|
36 python -W ignore $__tool_directory__/mirgene_with_new_names.py -con $controls -tre $treateds -analysis $analysis -tool_dir $__tool_directory__ -gen "$path" -f "$mir_input.database" -percentage "$f.fil1" -counts "$f.fil2" -name1 "$fal1" -name2 "$fal2"
|
|
|
37 #end if
|
|
|
38 #end if
|
|
|
39
|
|
|
40 </command>
|
|
|
41 <inputs>
|
|
|
42 <param name="analysis" type="select" label="Discover miR with templated or/and non-templated isomiRs" help="Choose the category of miRNAs for detection">
|
|
|
43 <option value="1" selected="true">Detection of only templated miRNAs</option>
|
|
|
44 <option value="2">Detection of templated and non-templated miRNAs</option>
|
|
|
45 </param>
|
|
|
46 <conditional name="mir_input">
|
|
|
47 <param name="database" type="select" label="Reference miRNA DB" help="Choose which database prefer to be used.">
|
|
|
48 <option value="1" selected="true">MirBase</option>
|
|
|
49 <option value="2">MirGene</option>
|
|
|
50 </param>
|
|
|
51 <when value="1">
|
|
|
52 <param name="genome1" type="select" label="Reference miRNAs (organism)" help="If your genome coordinates of interest are not listed, contact the Galaxy team">
|
|
|
53 <options from_data_table="mirbase" />
|
|
|
54 </param>
|
|
|
55 </when>
|
|
|
56 <when value="2">
|
|
|
57 <param name="genome2" type="select" label="Reference miRNAs (organism)" help="If your genome coordinates of interest are not listed, contact the Galaxy team">
|
|
|
58 <options from_data_table="mirgene" />
|
|
|
59 </param>
|
|
|
60 </when>
|
|
|
61 </conditional>
|
|
|
62 <param name="fal1" type="text" value="FactorLevel" label="Specify a factor level, typical values could be 'tumor', 'normal', 'treated' or 'control'"/>
|
|
|
63 <param name="control" format="sam" type="data" multiple="True" label="Select BAM files of the factor level samples" />
|
|
|
64 <param name="fal2" type="text" value="FactorLevel" label="Specify a factor level, typical values could be 'tumor', 'normal', 'treated' or 'control'"/>
|
|
|
65 <param name="treated" format="sam" type="data" multiple="True" label="Select BAM files of the factor level samples" />
|
|
|
66 <conditional name="f">
|
|
|
67 <param name="fil" type="select" label="Filter low counts" help="Treat genes with very low expression as unexpressed and filter out">
|
|
|
68 <option value="1" selected="true">No</option>
|
|
|
69 <option value="2">Yes</option>
|
|
|
70 </param>
|
|
|
71 <when value="2">
|
|
|
72 <param name="fil1" type="integer" value="0" label="Minimum percentage of the samples" help="Filter out all genes that do not meet the Minimum counts in at least this many samples of every category"/>
|
|
|
73 <param name="fil2" type="integer" value="0" label="Minimum counts" help="Filter out all genes that do not meet this minimum count"/>
|
|
|
74 </when>
|
|
|
75 <when value="1">
|
|
|
76 </when>
|
|
|
77 </conditional>
|
|
|
78 <param name="db" type="boolean" checked="true" truevalue="1" falsevalue="0" label="Output Database files" />
|
|
|
79 <param name="cmatrix" type="boolean" checked="false" truevalue="1" falsevalue="0" label="Output Matrix files, one for each factor level" />
|
|
|
80 <param name="c_files" type="boolean" checked="true" truevalue="1" falsevalue="0" label="Output Count tables, one for each sample" />
|
|
|
81 </inputs>
|
|
|
82 <outputs>
|
|
|
83 <collection name="list_output1" type="list" label="Database ${fal1} templated" >
|
|
|
84 <discover_datasets pattern="__name__" format="tabular" directory="split1" />
|
|
|
85 <filter>db == 1 and (analysis == "1" or analysis == "2")</filter>
|
|
|
86 </collection>
|
|
|
87 <collection name="list_output2" type="list" label="Database ${fal2} templated" >
|
|
|
88 <discover_datasets pattern="__name__" format="tabular" directory="split2" />
|
|
|
89 <filter>db == 1 and (analysis == "1" or analysis == "2")</filter>
|
|
|
90 </collection>
|
|
|
91 <collection name="list_output3" type="list" label="Database ${fal1} non-templated" >
|
|
|
92 <discover_datasets pattern="__name__" format="tabular" directory="split3" />
|
|
|
93 <filter>db == 1 and analysis == "2"</filter>
|
|
|
94 </collection>
|
|
|
95 <collection name="list_output4" type="list" label="Database ${fal2} non-templated" >
|
|
|
96 <discover_datasets pattern="__name__" format="tabular" directory="split4" />
|
|
|
97 <filter>db == 1 and analysis == "2"</filter>
|
|
|
98 </collection>
|
|
|
99
|
|
|
100 <collection name="Counts" type="list" label="Count Matrices" >
|
|
|
101 <discover_datasets pattern="__name__" format="tabular" directory="Counts" />
|
|
|
102 <filter>cmatrix==1</filter>
|
|
|
103 </collection>
|
|
|
104 <collection name="list_output9" type="list" label="Count files ${fal1} for Differential Expression" >
|
|
|
105 <discover_datasets pattern="__name__" format="tabular" directory="Diff/temp_con" />
|
|
|
106 <filter>c_files==1 and (analysis == "1")</filter>
|
|
|
107 </collection>
|
|
|
108 <collection name="list_output10" type="list" label="Count files ${fal2} for Differential Expression" >
|
|
|
109 <discover_datasets pattern="__name__" format="tabular" directory="Diff/temp_tre" />
|
|
|
110 <filter>c_files==1 and (analysis == "1")</filter>
|
|
|
111 </collection>
|
|
|
112 <collection name="list_output11" type="list" label="Count files ${fal1} for Differential Expression" >
|
|
|
113 <discover_datasets pattern="__name__" format="tabular" directory="Diff/n_temp_con" />
|
|
|
114 <filter>c_files==1 and analysis == "2"</filter>
|
|
|
115 </collection>
|
|
|
116 <collection name="list_output12" type="list" label="Count files ${fal2} for Differential Expression" >
|
|
|
117 <discover_datasets pattern="__name__" format="tabular" directory="Diff/n_temp_tre" />
|
|
|
118 <filter>c_files==1 and analysis == "2"</filter>
|
|
|
119 </collection>
|
|
|
120
|
|
|
121 <data name="Results non templated treated1" format="pdf" label="IsomiR Profile Report" from_work_dir="$__tool_directory__/report1.pdf" />
|
|
|
122 </outputs>
|
|
|
123 <help><![CDATA[
|
|
|
124 **Inputs**
|
|
|
125
|
|
|
126 IsoRead takes SAM files from alignment tool and RefSeq miRNA sequences from miRbase or miRGene as input. It discovers and classifies miRs and their isoforms in two main types—templated isomiRs (_t) and non-templated isomiRs (_nont) across samples. The coordinates in the name extension show the sequence offset relative to the RefSeq sequence.
|
|
|
127
|
|
|
128 **Examples**
|
|
|
129
|
|
|
130 ============================= ========================
|
|
|
131 **RefSeq miRNA** **Sequence**
|
|
|
132 ----------------------------- ------------------------
|
|
|
133 hsa-let-7c-5p_chr21(+) TGAGGTAGTAGGTTGTATGGTT
|
|
|
134 **Templated miRNA** **Sequence**
|
|
|
135 hsa-let-7c-5p_chr21(+)_t_0_-1 TGAGGTAGTAGGTTGTATGGT
|
|
|
136 **Non Templated miRNA** **Sequence**
|
|
|
137 hsa-let-7c-5p_nont_0_+2_AT TGAGGTAGTAGGTTGTATGGTTAT
|
|
|
138 ============================= ========================
|
|
|
139
|
|
|
140 -----
|
|
|
141
|
|
|
142 **Outputs**
|
|
|
143
|
|
|
144 As output, the tool produces tabular files with all the detected miRs and isomiRs per category for every SAM file as a database for the user and two types of count matrices containing the copy number of each entry. One of the count matrix types is compatible with the DESeq2 and EdgeR tools so that the user can continue with differential expression analysis. The other types of count matrices generated per miR category (RefSeq miRs, template isomiRs, and non-template isomiRs) and per group are enriched with read sequences. Apart from the tabular files, the tool produces a PDF report file with several custom graphs and statistical information about the isomiR profiles in the user data.
|
|
|
145
|
|
|
146 **Examples**
|
|
|
147
|
|
|
148 Database file
|
|
|
149
|
|
|
150 .. image:: isoread.png
|
|
|
151
|
|
|
152 Count matrix enriched with read sequences
|
|
|
153
|
|
|
154 =============================== ====================== =========== ===========
|
|
|
155 **Name** **Sequence** **Sample1** **Sample2**
|
|
|
156 ------------------------------- ---------------------- ----------- -----------
|
|
|
157 hsa-miR-320a-3p_chr8(-) AAAAGCTGGGTTGAGAGGGCGA 260 609
|
|
|
158 hsa-miR-320a-3p_chr8(-)_t_+1_+1 AAAGCTGGGTTGAGAGGGCGAA 192 30
|
|
|
159 =============================== ====================== =========== ===========
|
|
|
160
|
|
|
161
|
|
|
162
|
|
|
163 ]]></help>
|
|
|
164 <citations>
|
|
|
165 <citation type="doi">10.3390/cancers13225663</citation>
|
|
|
166 </citations>
|
|
|
167 </tool>
|
|
|
168
|
