5
|
1 <tool id="SAINT_preprocessing_v5" name="SAINT pre-processing">
|
|
2 <description></description>
|
|
3 <command interpreter="python">
|
|
4 #if (str($type) == 'Scaffold'):
|
|
5 SAINT_preprocessing.py $input $preybool $fasta_db $Inter_file $Prey_file
|
|
6 "
|
|
7 #for $ba in $bait
|
|
8 ${ba.bait1}
|
|
9 ${ba.assign}
|
|
10 ${ba.T_C}
|
|
11 #end for
|
|
12 "
|
|
13 $Bait_file \$INSTALL_RUN_PATH/ $bait_bool $bait_file_in
|
|
14 #elif (str($type) == 'MaxQuant'):
|
|
15 SAINT_preprocessing_mq_pep.py $input $preybool $fasta_db $Inter_file $Prey_file
|
|
16 "
|
|
17 #for $ba in $bait
|
|
18 ${ba.bait1}
|
|
19 ${ba.assign}
|
|
20 ${ba.T_C}
|
|
21 #end for
|
|
22 "
|
|
23 $Bait_file \$INSTALL_RUN_PATH/ $bait_bool $bait_file_in
|
55
|
24 #elif (str($type) == 'Peptideshaker'):
|
|
25 Protein_report_processing.py "$input" $Bait_file $use_metric $fasta_db $preybool
|
|
26 "
|
|
27 #for $ba in $bait
|
|
28 ${ba.bait1}
|
|
29 ${ba.assign}
|
|
30 ${ba.T_C}
|
|
31 #end for
|
|
32 "
|
|
33 \$INSTALL_RUN_PATH/ $bait_bool $bait_file_in $Inter_file $Prey_file
|
5
|
34 #end if
|
|
35 </command>
|
|
36 <requirements>
|
|
37 <requirement type="set_environment">INSTALL_RUN_PATH</requirement>
|
|
38 </requirements>
|
|
39 <inputs>
|
55
|
40 <conditional name="type_select">
|
|
41 <param type="select" name="type" label="MaxQuant or Scaffold">
|
|
42 <option value="MaxQuant">MaxQuant</option>
|
|
43 <option value="Scaffold">Scaffold</option>
|
|
44 <option value="Peptideshaker">Peptideshaker</option>
|
|
45 </param>
|
|
46 <when value="MaxQuant">
|
|
47 <param format="dat" name="input" type="data" label="MaxQuant peptides Output"/>
|
|
48 <param type="boolean" name="preybool" checked="true" label="Create Prey File"/>
|
|
49 <param type="data" name="fasta_db" format="fasta" label="Provide Uniprot Fasta database" optional="true"/>
|
|
50 <param name="bait_bool" type="boolean" checked="true" label="Are You Providing Your Own bait file?"/>
|
|
51 <param type="data" format="dat" name="bait_file_in" label="Bait File" optional="true"/>
|
|
52 <repeat name="bait" title="Bait Create">
|
|
53 <param name="bait1" type="text" size="100" label="Bait"/>
|
|
54 <param name="assign" type="text" size="100" label="Group Assignment"/>
|
|
55 <param name="T_C" type="boolean" checked="true" label="Is this a Control?"/>
|
|
56 </repeat>
|
59
|
57 </when>
|
55
|
58 <when value="Scaffold">
|
|
59 <param format="dat" name="input" type="data" label="Scaffold Output"/>
|
|
60 <param type="boolean" name="preybool" checked="true" label="Create Prey File"/>
|
|
61 <param type="data" name="fasta_db" format="fasta" label="Provide Uniprot Fasta database" optional="true"/>
|
|
62 <param name="bait_bool" type="boolean" checked="true" label="Are You Providing Your Own bait file?"/>
|
|
63 <param type="data" format="dat" name="bait_file_in" label="Bait File" optional="true"/>
|
|
64 <repeat name="bait" title="Bait Create">
|
|
65 <param name="bait1" type="text" size="100" label="Bait"/>
|
|
66 <param name="assign" type="text" size="100" label="Group Assignment"/>
|
|
67 <param name="T_C" type="boolean" checked="true" label="Is this a Control?"/>
|
|
68 </repeat>
|
|
69 </when>
|
|
70 <when value="Peptideshaker">
|
|
71 <param format="dat" name="input" type="data" label="Peptideshaker Output" multiple="true"/>
|
|
72 <param type="select" name="use_metric" label="Select Report File Value for Quantification">
|
|
73 <option value="Validated_Peptides">#Validated Peptides</option>
|
|
74 <option value="Peptides">#Peptides</option>
|
|
75 <option value="Unique">#Unique</option>
|
|
76 <option value="Validated_PSMs">#Validated PSMs</option>
|
|
77 <option value="PSMs">#PSMs</option>
|
|
78 </param>
|
|
79 <param type="boolean" name="preybool" checked="true" label="Create Prey File"/>
|
|
80 <param type="data" name="fasta_db" format="fasta" label="Provide Uniprot Fasta database" optional="true"/>
|
|
81 <param name="bait_bool" type="boolean" checked="true" label="Are You Providing Your Own bait file?"/>
|
|
82 <param type="data" format="dat" name="bait_file_in" label="Bait File" optional="true"/>
|
|
83 <repeat name="bait" title="Bait Create">
|
|
84 <param name="bait1" type="text" size="100" label="Bait"/>
|
|
85 <param name="assign" type="text" size="100" label="Group Assignment"/>
|
|
86 <param name="T_C" type="boolean" checked="true" label="Is this a Control?"/>
|
|
87 </repeat>
|
|
88 </when>
|
|
89 </conditional>
|
5
|
90 </inputs>
|
|
91 <outputs>
|
|
92 <data format="txt" name="Inter_file" label="Inter File"/>
|
|
93 <data format="txt" name="Prey_file" label="Prey File" />
|
|
94 <data format="txt" name="Bait_file" label="Bait File" />
|
|
95 </outputs>
|
|
96 <stdio>
|
52
|
97 <regex match="Error: bad bait"
|
|
98 source="stdout"
|
|
99 level="fatal"
|
|
100 description="Error: bad bait"/>
|
|
101 <regex match="IOError: [Errno 2] No such file or directory: './tukeys_output.txt'"
|
|
102 source="stderr"
|
|
103 level="fatal"
|
|
104 description="Error: Scaffold/MaxQuant mismatch. Check job settings."/>
|
|
105 <regex match="bait_temp_file = open(sys.argv[10], 'r')"
|
|
106 source="stderr"
|
|
107 level="fatal"
|
|
108 description="Error: Bait create settings mismatch. Check job settings."/>
|
5
|
109 <regex match="Error|error"
|
|
110 source="stdout"
|
|
111 level="fatal"
|
|
112 description="Unknown error"/>
|
|
113 <regex match="Error|error"
|
|
114 source="stderr"
|
|
115 level="fatal"
|
|
116 description="Unknown error"/>
|
|
117 </stdio>
|
|
118
|
|
119 <tests>
|
|
120 <test>
|
|
121 <param name="input" value="fa_gc_content_input.fa"/>
|
|
122 <output name="out_file1" file="fa_gc_content_output.txt"/>
|
|
123 </test>
|
|
124 </tests>
|
|
125 <help>
|
21
|
126 Pre-Processing
|
25
|
127 ^^^^^^^^^^^^^^
|
20
|
128
|
21
|
129 This tool will read in a Scaffold *Samples Report* file (tab-delimited
|
55
|
130 txt file), multiple Peptideshaker *Protein Report* files or a MaxQuant
|
|
131 *peptides.txt* file and process them to generate a *Bait File, Prey
|
|
132 File,* and *Inter File* for SAINTexpress analysis.
|
20
|
133
|
25
|
134 --------------
|
|
135
|
27
|
136 **1) MaxQuant or Scaffold**
|
20
|
137
|
21
|
138 APOSTL is able to recognize either a Scaffold *Samples Report* file
|
|
139 (tab-delimited txt file) or the *peptides.txt* file output in the
|
|
140 MaxQuant *txt* output folder. No modifications should be made to these
|
|
141 files. Please designate which one is being provided.
|
20
|
142
|
27
|
143 **2) Scaffold or MaxQuant File Input**
|
20
|
144
|
21
|
145 Select the corresponding *Samples Report* or *peptides.txt* here.
|
|
146
|
27
|
147 **3) Create Prey File**
|
5
|
148
|
20
|
149 Select whether or not you would like APOSTL to generate a prey file.
|
|
150
|
21
|
151 When making a prey file, APOSTL queries a user provided FASTA database
|
|
152 (see below) in order to extract protein amino acid lengths and gene
|
|
153 names. This may take several minutes depending on your computer and if
|
|
154 your Galaxy distribution is cluster enabled. Some users may want to run
|
|
155 SAINTexpress using the same data set while changing which baits are
|
|
156 considered test or control. It is useful to toggle **Create Prey File**
|
|
157 off in order to save time by circumventing this step as the same prey
|
|
158 file can be used for both SAINTexpress runs.
|
20
|
159
|
27
|
160 **4) Provide Uniprot FASTA database**
|
20
|
161
|
21
|
162 Please specify a FASTA file containing the protein sequence. You can
|
|
163 download the latest release of Uniprot's reviewed sequences (Swiss-Prot)
|
24
|
164 http://www.uniprot.org/downloads. For computing efficiency,
|
21
|
165 please limit your FASTA file to your species of interest. It is
|
|
166 generally advised to use the same FASTA database that was used for your
|
|
167 database searches during peptide and protein ID.
|
5
|
168
|
27
|
169 **5) Are You Providing Your Own Bait File?**
|
20
|
170
|
21
|
171 Users have two options for providing a *Bait* file. If you would like to
|
|
172 create your own *Bait* file, select Yes and designate the file below. If
|
|
173 you would like to use the **Bait Create** tool to create a *Bait* file
|
|
174 within Galaxy, please select No and proceed to the **Bait Create** tool.
|
20
|
175
|
27
|
176 **6) Bait File**
|
21
|
177
|
|
178 If you selected Yes above, please specify the *Bait* file you would like
|
|
179 to provide. If you selected No above, ignore this argument and proceed
|
|
180 to the **Bait Create** tool.
|
5
|
181
|
21
|
182 **Note:** Individual bait names must match to the sample names
|
|
183 within your Scaffold or MaxQuant output exactly and must contain no
|
|
184 whitespace characters (e.g. spaces) or dashes.
|
59
|
185
|
55
|
186 **Note:** When using Peptideshaker output the baits in the baitfile
|
|
187 need to be in the same order as the files selected.
|
5
|
188
|
27
|
189 **7) Bait Create**
|
21
|
190
|
|
191 Using the **Bait Create** tool, you can create your *bait.txt* file.
|
20
|
192
|
21
|
193 **Note:** Individual bait names must match to the sample names
|
|
194 within your Scaffold or MaxQuant output exactly and must contain no
|
|
195 whitespace characters (e.g. spaces) or dashes.
|
5
|
196
|
21
|
197 APOSTL uses this bait file to find the user's baits of interest within a
|
|
198 *Samples Report* or *peptides.txt* file and when preparing the
|
|
199 *inter.txt* file.
|
20
|
200
|
21
|
201 --------------
|
20
|
202
|
21
|
203 Once your parameters have been finalized, please press the Execute
|
|
204 button to start processing. This may take a few minutes. Once your
|
|
205 process has been completed, you will see your *Bait File, Prey File,*
|
|
206 and *Inter File* on the right hand side of your panel highlighted in
|
|
207 green.
|
5
|
208 </help>
|
|
209 </tool>
|