|
0
|
1 <tool id="rbc_rnabob" name="RNABOB" version="2.2.1.0">
|
|
|
2 <description>Fast Pattern searching for RNA secondary structures</description>
|
|
|
3 <requirements>
|
|
|
4 <requirement type="package" version="2.2.1">rnabob</requirement>
|
|
|
5 </requirements>
|
|
|
6 <version_command>echo "2.2.1"</version_command>
|
|
|
7 <command>
|
|
|
8 <![CDATA[
|
|
1
|
9 rnabob
|
|
|
10 -q
|
|
|
11 $fancy
|
|
|
12 $compStrands
|
|
|
13 $skipOverlapping
|
|
|
14 $descriptorFile
|
|
|
15 $sequenceFile > $stdout
|
|
0
|
16 ]]>
|
|
|
17 </command>
|
|
|
18 <stdio>
|
|
|
19 <exit_code range="1:" level="fatal" description="Error occurred. Please check Tool Standard Error" />
|
|
|
20 <exit_code range=":-1" level="fatal" description="Error occurred. Please check Tool Standard Error" />
|
|
|
21 </stdio>
|
|
|
22 <inputs>
|
|
|
23 <param name="descriptorFile" type="data" format="txt" multiple="false" label="Motif Descriptor File" help="This file contains the description of the motif for which to search"/>
|
|
1
|
24 <param name="sequenceFile" type="data" format="fasta" multiple="false" label="Sequence File" help="This file specifies the sequence in which the motif will be searched"/>
|
|
|
25 <param name="compStrands" type="boolean" truevalue="-c" falsevalue="" checked="false" label="Also search on complementary strands" help="-c : Search both strands of the supplied sequence"/>
|
|
|
26 <param name="skipOverlapping" type="boolean" truevalue="-s" falsevalue="" checked="false" label="Skip overlapping matches" help="-s : This is a workaround to avoid a problem in the DNABANK, overlapping matches will be ignored"/>
|
|
|
27 <param name="fancy" type="boolean" checked="false" truevalue="-F" falsevalue="" label="Show Alignments" help="Display full alignments to pattern"/>
|
|
0
|
28 </inputs>
|
|
|
29 <outputs>
|
|
|
30 <data format="txt" name="stdout" label="${tool.name} on ${on_string}" />
|
|
|
31 </outputs>
|
|
|
32 <tests>
|
|
|
33 <test>
|
|
|
34 <param name="descriptorFile" value="r17.des" />
|
|
|
35 <param name="sequenceFile" value="F22B7.fa" />
|
|
|
36 <param name="compStrands" value="True" />
|
|
|
37 <param name="skipOverlapping" value="False" />
|
|
|
38 <param name="fancy" value="False" />
|
|
|
39 <output name="stdout" file="r17.bob" />
|
|
|
40 </test>
|
|
|
41 <test>
|
|
|
42 <param name="descriptorFile" value="trna.des" />
|
|
|
43 <param name="sequenceFile" value="F22B7.fa" />
|
|
|
44 <param name="compStrands" value="True" />
|
|
|
45 <param name="skipOverlapping" value="False" />
|
|
|
46 <param name="fancy" value="False" />
|
|
|
47 <output name="stdout" file="trna.bob" />
|
|
|
48 </test>
|
|
|
49 </tests>
|
|
|
50 <help>
|
|
1
|
51 <![CDATA[
|
|
0
|
52 **What RNABOB does**
|
|
|
53
|
|
|
54 RNABOB allows searching a sequence database for RNA structural motifs.
|
|
|
55 The probe motif is specified in a *descriptor* file,
|
|
|
56 which describes its primary sequence, secondary structure, and tertiary constraints.
|
|
|
57 The source in its original packaging can be found at http://selab.janelia.org/software/#rnabob.
|
|
|
58
|
|
|
59 -----
|
|
|
60
|
|
|
61 **Sequence database format**
|
|
|
62
|
|
1
|
63 RNABOB is currently restricted to reading sequence files in FASTA format.
|
|
0
|
64 The command line version of RNABOB can also read sequence files in GCG, EMBL, GenBank and other formats.
|
|
|
65
|
|
|
66 -----
|
|
|
67
|
|
|
68 **Descriptor file syntax**
|
|
|
69
|
|
1
|
70 The descriptor file syntax is fairly powerful, and allows a great deal of freedom for specifying
|
|
0
|
71 RNA motifs. The syntax is therefore a bit complicated.
|
|
|
72
|
|
|
73 The descriptor file has two parts: a **topology** description and an **explicit** description.
|
|
|
74
|
|
1
|
75 The first non-blank, non-comment line of the file is the topology description. It defines the
|
|
|
76 order of occurrence of a series of single-stranded, double-stranded and related elements. Each
|
|
|
77 element must be given a unique name (a number, typically) and must be prefixed with '**s**',
|
|
|
78 '**h**', or '**r**', indicating single-strand, helical, or a relational element. Helical and
|
|
0
|
79 relational elements are paired to other elements, which are suffixed by a prime, **\'**.
|
|
|
80
|
|
|
81 For example::
|
|
|
82
|
|
1
|
83 \
|
|
|
84 h1 s1 h1'
|
|
0
|
85
|
|
1
|
86 describes a hairpin loop structure with a simple helix and single-stranded loop. If the helix
|
|
0
|
87 always contained a non-canonical base pair at one position, the topology coud be described as::
|
|
|
88
|
|
1
|
89 \
|
|
|
90 h1 r1 h2 s1 h2' r1' h1'
|
|
0
|
91
|
|
1
|
92 where r1,r1' indicate a correlation, where the sequence r1 constrains the sequence of r1'.
|
|
0
|
93 (Helices are a special case of this.)
|
|
|
94
|
|
1
|
95 The remaining non-comment, non-blank lines are explicit descriptions of each element in turn. Each
|
|
|
96 line contains 3 or 4 fields, separated by tabs or blank space. The first field is the name of the
|
|
|
97 element, from the topology description. The second field is the number of mismatches allowed in
|
|
0
|
98 this element. The third field is the primary sequence constraint to apply to this element.
|
|
|
99
|
|
1
|
100 Helices and relational element pairs are specified on a single line rather than two. Mismatches
|
|
|
101 and primary sequence constraints are given as pairs, separated by a colon '**:**'. The left side
|
|
|
102 is the constraint applied to the upstream element, and the right side is applied to the downstream
|
|
0
|
103 elements.
|
|
|
104
|
|
1
|
105 The primary sequence constraint is given as a sequence of nucleotides. Any IUPAC single-letter
|
|
|
106 code is recognized, including N if the position can have any base identity. Allowed length
|
|
|
107 variations are specified with asterisks ``'*'``, where each ``*`` will allow either 0 or 1 N at
|
|
0
|
108 that position.
|
|
|
109
|
|
|
110 For example::
|
|
|
111
|
|
1
|
112 \
|
|
|
113 GGAGG******NNNAUG
|
|
0
|
114
|
|
1
|
115 specifies a GGAGG Shine/Dalgarno site and an AUG initiation codon, separated by a spacer of 3 to 9
|
|
0
|
116 nucleotides of any sequence.
|
|
|
117
|
|
|
118 An alternative syntax can be used for very long gaps::
|
|
|
119
|
|
1
|
120 \
|
|
|
121 GGAGG[10]NNNAUG is the same as GGAGG**********NNNAUG
|
|
0
|
122
|
|
1
|
123 Be careful defining variable length helices and relational elements; if the number and type (gap
|
|
|
124 or identity) of position do not match on left and right sides, the program will refuse to accept
|
|
0
|
125 the descriptor.
|
|
|
126
|
|
1
|
127 Relational elements have an additional field which specifies a "transformation matrix" of four
|
|
|
128 nucleotides, specifying the rule for making the ``r'`` pattern from the ``r`` sequence in order
|
|
|
129 ``A-C-G-T``. For example, the transformation matrix for a simple helix is ``TGCA``; if you allow
|
|
|
130 ``G-U`` pairs, it is ``TGYR``. RNABOB allows ``G-U`` pairing by default and uses the ``TGYR``
|
|
0
|
131 matrix for helical elements.
|
|
|
132
|
|
|
133 For example, the explicit description of our hairpin might be:
|
|
|
134
|
|
|
135 ::
|
|
|
136
|
|
1
|
137 \
|
|
|
138 h1 0:0 NNN:NNN
|
|
|
139 r1 0:0 R:N GNAN
|
|
|
140 h2 0:0 **NC:GN**
|
|
|
141 s1 0 UUCG
|
|
0
|
142
|
|
1
|
143 This describes a stem of 6 to 8 base pairs, in which the 4th pair from the bottom of the stem must
|
|
|
144 be a non-canonical GA pair. Note that, in general, the left side of the primary constraint for
|
|
|
145 helices and relational elements is redundant, and should be given as all N. In some cases it is
|
|
|
146 convenient to constrain the right side to require a particular base pair (GU, for instance) at one
|
|
0
|
147 position.
|
|
|
148
|
|
1
|
149 A note on mismatches: The split format for helices and relational elements works like this. The
|
|
|
150 number on the left constrains the primary sequence match of the left side of the primary
|
|
|
151 constraint. The number on the right constrains the match of the right side of the primary
|
|
|
152 constraint, *after* that side has been constructed according to the sequence on the left. In other
|
|
|
153 words, the number on the left constrains the mismatches in primary sequence only, while the number
|
|
0
|
154 on the right will constrain the number of mispaired positions in the helix.
|
|
|
155
|
|
1
|
156 Finally: any line that begins with a pound sign '#' is a comment line, and will not be interpreted
|
|
0
|
157 by the pattern compiler.
|
|
|
158
|
|
|
159 **Options**
|
|
|
160
|
|
|
161 The behavior of RNABOB can be modified by use of the following options:
|
|
|
162
|
|
1
|
163 *Complement*: Selecting this option will cause RNABOB to search for the pattern also on the
|
|
0
|
164 complementary strands.
|
|
|
165
|
|
1
|
166 *Skip*: This is a workaround to avoid a problem in the DNABANK. There are some sequences in the
|
|
|
167 database which have long stretches of ambiguous sequence (N's). Descriptors with no primary
|
|
|
168 sequence constraints will match these garbage sequences at many, many positions, and generate huge
|
|
|
169 outputs. This option toggles a search strategy that skips forward a pattern-length rather than a
|
|
|
170 single base when a match is found, thus printing out only a single match when overlapping matches
|
|
|
171 are found.
|
|
0
|
172
|
|
|
173 **Examples**
|
|
|
174
|
|
1
|
175 The following example descriptors included in the source distribution
|
|
0
|
176 (http://selab.janelia.org/software/rnabob/rnabob.tar.gz):
|
|
|
177
|
|
1
|
178 - trna.des - a general descriptor of a tRNA structure
|
|
|
179 - r17.des - descriptor of the consensus binding site for the r17 phage coat protein
|
|
|
180 - pseudoknot.des - description of a simple pseudoknotted structure
|
|
0
|
181
|
|
1
|
182 An example cosmid ``F22B7.fa`` from the *C. elegans* genome sequencing project is also provided
|
|
0
|
183 for running these descriptors against.
|
|
|
184
|
|
|
185 ::
|
|
|
186
|
|
1
|
187 \
|
|
|
188 # trna.des
|
|
|
189 #
|
|
|
190 # Generalized descriptor of a tRNA cloverleaf. Doesn't
|
|
|
191 # find them all though.
|
|
|
192 #
|
|
0
|
193
|
|
1
|
194 h1 s1 h2 s2 h2' s3 h3 s4 h3' s5 h4 s6 h4' h1' s8
|
|
0
|
195
|
|
1
|
196 h1 0:2 NNNNNNN:NNNNNNN
|
|
|
197 h2 0:1 *NNN:NNN*
|
|
|
198 h3 0:1 NNNNN:NNNNN
|
|
|
199 h4 0:1 NNNNN:NNNNN
|
|
|
200 s1 0 TN
|
|
|
201 s2 0 NNNN**********
|
|
|
202 s3 0 N
|
|
|
203 s4 0 NNNNNN*
|
|
|
204 s5 0 NN********************
|
|
|
205 s6 0 TTC****
|
|
|
206 s8 0 NCCA
|
|
0
|
207
|
|
1
|
208 Running RNABOB with ``trna.des`` against ``F22B7.fa`` searches the top strand of the cosmid for
|
|
|
209 the above motif. ``trna.des`` hits twice, once on each strand. (F22B7 has several other tRNA genes
|
|
0
|
210 in it which the pattern fails to detect - this is *not* a pattern to use for tRNA genefinding!).
|
|
1
|
211 ]]>
|
|
|
212 </help>
|
|
0
|
213 <citations>
|
|
1
|
214 <citation type="doi">10.1093/bioinformatics/6.4.325</citation>
|
|
|
215 <citation type="bibtex">@UNPUBLISHED{rnabob,
|
|
0
|
216 author = {Eddy S.R},
|
|
|
217 title = {RNABOB: a program to search for RNA secondary structure motifs in sequence databases},
|
|
|
218 note = {}}</citation>
|
|
|
219 </citations>
|
|
|
220 </tool>
|
