|
0
|
1 <tool id="pyclone" name="PyClone" version="0.13.0.01762d9">
|
|
|
2 <description>
|
|
|
3 PyClone is a tool for inferring the cellular prevalence of point mutations from deeply sequenced data
|
|
|
4 </description>
|
|
|
5
|
|
|
6 <requirements>
|
|
|
7 <requirement type="set_environment">MORINLAB_SCRIPTS</requirement>
|
|
|
8 <requirement type="package" version="14.2">morinlab_scripts</requirement>
|
|
|
9 <requirement type="package" version="0.13.0.01762d9">pyclone</requirement>
|
|
|
10 <requirement type="package" version="1.0">pyclone_dependencies</requirement>
|
|
|
11 <!-- <requirement type="package" version="3.11">pyyaml</requirement>
|
|
|
12 <requirement type="package" version="0.2.3">pydp</requirement>
|
|
|
13 <requirement type="package" version="0.7.1">seaborn</requirement> -->
|
|
|
14 </requirements>
|
|
|
15
|
|
|
16 <command>
|
|
|
17 <![CDATA[
|
|
|
18 PyClone run_analysis_pipeline
|
|
|
19 --in_files
|
|
|
20 #for $file in $infiles:
|
|
|
21 '$file'
|
|
|
22 #end for
|
|
|
23 --working_dir ./pyclone_results
|
|
|
24 #if $adv_opts.adv_opts_selector == "advanced" and $adv_opts.TUMOUR_CONTENTS != "":
|
|
|
25 --tumour_contents "$adv_opts.TUMOUR_CONTENTS"
|
|
|
26 #end if
|
|
|
27 #if $adv_opts.adv_opts_selector == "advanced" and $adv_opts.SAMPLES != "":
|
|
|
28 --samples $adv_opts.SAMPLES
|
|
|
29 #end if
|
|
|
30 #if $adv_opts.adv_opts_selector == "advanced":
|
|
|
31 --density $adv_opts.DENSITY
|
|
|
32 --num_iters $adv_opts.NUM_ITERS
|
|
|
33 --prior $adv_opts.PRIOR
|
|
|
34 --burnin $adv_opts.BURNIN
|
|
|
35 --thin $adv_opts.THIN
|
|
|
36 --mesh_size $adv_opts.MESH_SIZE
|
|
|
37 --min_cluster_size $adv_opts.MIN_CLUSTER_SIZE
|
|
|
38 #end if
|
|
|
39 ;
|
|
|
40
|
|
|
41 gs -dBATCH -dNOPAUSE -q -sDEVICE=pdfwrite -sOutputFile=lociplots ./pyclone_results/plots/loci/parallel_coordinates.pdf ./pyclone_results/plots/loci/scatter.pdf ./pyclone_results/plots/loci/vaf_parallel_coordinates.pdf ./pyclone_results/plots/loci/vaf_scatter.pdf;
|
|
|
42
|
|
|
43 gs -dBATCH -dNOPAUSE -q -sDEVICE=pdfwrite -sOutputFile=clusterplots ./pyclone_results/plots/cluster/parallel_coordinates.pdf ./pyclone_results/plots/cluster/scatter.pdf ./pyclone_results/plots/cluster/density.pdf;
|
|
|
44
|
|
|
45 cat ./clusterplots > $cluster_plots;
|
|
|
46 cat ./lociplots > $loci_plots;
|
|
|
47 cat ./pyclone_results/plots/loci/density.pdf > $loci_density_plot;
|
|
|
48
|
|
|
49 cat ./pyclone_results/tables/cluster.tsv > $cluster_table;
|
|
|
50 cat ./pyclone_results/tables/loci.tsv > $loci_table;
|
|
|
51
|
|
|
52 #if $plot_custom.plot_options == "multisample":
|
|
|
53 Rscript \$MORINLAB_SCRIPTS/plot_pyclone_output/plot_ccfs.R
|
|
|
54
|
|
|
55 $plot_custom.patient
|
|
|
56
|
|
|
57 #if $adv_opts.adv_opts_selector == "advanced" and $adv_opts.SAMPLES != "":
|
|
|
58 --samples $adv_opts.SAMPLES
|
|
|
59 #end if
|
|
|
60
|
|
|
61 ./pyclone_results/tables/loci.tsv
|
|
|
62
|
|
|
63 $plot_custom.genes
|
|
|
64
|
|
|
65 #if $plot_custom.maf_files:
|
|
|
66 --mafs
|
|
|
67 #for $maf in $plot_custom.maf_files:
|
|
|
68 '$maf'
|
|
|
69 #end for
|
|
|
70 #end if
|
|
|
71
|
|
|
72 #if $plot_custom.effects
|
|
|
73 --effects $plot_custom.effects
|
|
|
74 #end if
|
|
|
75 ;
|
|
|
76 #end if
|
|
|
77 ]]>
|
|
|
78
|
|
|
79 </command>
|
|
|
80
|
|
|
81 <stdio>
|
|
|
82 <exit_code range="1:" />
|
|
|
83 <exit_code range=":-1" />
|
|
|
84 <regex match="Error:" />
|
|
|
85 <regex match="Exception:" />
|
|
|
86 </stdio>
|
|
|
87
|
|
|
88 <inputs>
|
|
|
89 <param type="data" name="infiles" format="tabular" multiple="true" label="input_file" help="PyClone input file with mutation and copy number info, as generated by generate_pyclone_input_tool" />
|
|
|
90
|
|
|
91
|
|
|
92 <conditional name="adv_opts">
|
|
|
93 <param name="adv_opts_selector" type="select" label="Advanced Options">
|
|
|
94 <option value="basic" selected="True">Hide Advanced Options</option>
|
|
|
95 <option value="advanced">Show Advanced Options</option>
|
|
|
96 </param>
|
|
|
97
|
|
|
98 <when value="basic" />
|
|
|
99
|
|
|
100 <when value="advanced">
|
|
|
101 <param type="text" name="TUMOUR_CONTENTS" value="" label="Space delimited list of tumour contents. Should match the order of --in_files (the list of inputs in mutations_file). If not given tumour content is assumed to 1.0 in all samples."/>
|
|
|
102
|
|
|
103 <param type="text" name="SAMPLES" value="" label="Space delimited list of sample names. Should be in the same order as --in_files (the list of inputs in mutations_file). If not set sample name will be inferred from file names and ordering in plots will be arbitrary."/>
|
|
|
104
|
|
|
105 <param type="text" name="DENSITY" value="pyclone_beta_binomial" label="Emission density for the model. Default is pyclone_beta_binomial."/>
|
|
|
106
|
|
|
107 <param type="integer" name="NUM_ITERS" value="10000" label="Emission density for the model. Default is pyclone_beta_binomial."/>
|
|
|
108
|
|
|
109 <param type="select" name="PRIOR" value="major_copy_number" label="Method used to set the possible genotypes.">
|
|
|
110 <option value="major_copy_number">major_copy_number: Considers all possible genotypes with up to the major copy number of B alleles.</option>
|
|
|
111 <option value="parental_copy_number">parental_copy_number: Considers all possible genotypes compatible with the predicted parental copy number.</option>
|
|
|
112 <option value="total_copy_number">total_copy_number: Use this if information about major and minor copy number is not available.</option>
|
|
|
113 </param>
|
|
|
114
|
|
|
115 <param type="integer" name="BURNIN" value="0" label="Number of samples to discard as burning for the MCMC chain. Default is 0."/>
|
|
|
116
|
|
|
117 <param type="integer" name="THIN" value="1" label="Number of samples to thin MCMC trace. For example if thin=10 every tenth sample after burning will be used for inference. Default is 1."/>
|
|
|
118
|
|
|
119 <param type="integer" name="MESH_SIZE" value="101" label="Number of points to use for approximating the cluster posteriors. Default is 101."/>
|
|
|
120
|
|
|
121 <param type="integer" name="MIN_CLUSTER_SIZE" value="5" label="Clusters with fewer mutations than this value will not be plotted."/>
|
|
|
122
|
|
|
123 </when>
|
|
|
124
|
|
|
125 </conditional>
|
|
|
126
|
|
|
127 <conditional name="plot_custom">
|
|
|
128 <param name="plot_options" type="select" label="Custom Plotting Options" help="If running PyClone with more than one sample, use this menu to customize CCF plots.">
|
|
|
129 <option value="single">Hide options (Running PyClone with a single sample)</option>
|
|
|
130 <option value="multisample">Show options (Running PyClone with two or more samples)</option>
|
|
|
131 </param>
|
|
|
132
|
|
|
133 <when value="single" />
|
|
|
134
|
|
|
135 <when value="multisample">
|
|
|
136 <param type="text" name="patient" option="false" size="20" label="Patient ID"/>
|
|
|
137
|
|
|
138 <param type="data" name="genes" optional="false" format="text" label="Label mutations in these genes (provide a file with one gene per line)"/>
|
|
|
139
|
|
|
140 <param type="data" name="maf_files" optional="true" format="tabular" multiple="true" label="maf_file" help="MAF file for each sample, in same order as the samples were input" />
|
|
|
141
|
|
|
142 <param type="select" name="effects" optional="true" multiple="true" label="Select effect criteria to apply to mutations to label" display="checkboxes">
|
|
|
143 <option value="Splice_Site">Splice_Site</option>
|
|
|
144 <option value="Nonsense_Mutation">Nonsense_Mutation</option>
|
|
|
145 <option value="Frame_Shift_Del">Frame_Shift_Del</option>
|
|
|
146 <option value="Frame_Shift_Ins">Frame_Shift_Ins</option>
|
|
|
147 <option value="Nonstop_Mutation">Nonstop_Mutation</option>
|
|
|
148 <option value="Translation_Start_Site">Translation_Start_Site</option>
|
|
|
149 <option value="In_Frame_Ins">In_Frame_Ins</option>
|
|
|
150 <option value="In_Frame_Del">In_Frame_Del</option>
|
|
|
151 <option value="Missense_Mutation">Missense_Mutation</option>
|
|
|
152 <option value="Intron">Intron</option>
|
|
|
153 <option value="Splice_Region">Splice_Region</option>
|
|
|
154 <option value="Silent">Silent</option>
|
|
|
155 <option value="RNA">RNA</option>
|
|
|
156 <option value="5'UTR">5'UTR</option>
|
|
|
157 <option value="3'UTR">3'UTR</option>
|
|
|
158 <option value="IGR">Intergenic/regulatory region</option>
|
|
|
159 <option value="5'Flank">5'Flank</option>
|
|
|
160 <option value="3'Flank">3'Flank</option>
|
|
|
161 </param>
|
|
|
162
|
|
|
163 </when>
|
|
|
164
|
|
|
165 </conditional>
|
|
|
166
|
|
|
167
|
|
|
168 </inputs>
|
|
|
169
|
|
|
170 <outputs>
|
|
|
171 <data name="cluster_plots" label="${tool.name} on ${on_string} - Cluster Plots" format="pdf">
|
|
|
172 <filter>output_files == "both" or output_files == "plots"</filter>
|
|
|
173 </data>
|
|
|
174
|
|
|
175 <data name="loci_plots" label="${tool.name} on ${on_string} - Loci Plots" format="pdf">
|
|
|
176 <filter>output_files == "both" or output_files == "plots"</filter>
|
|
|
177 </data>
|
|
|
178
|
|
|
179 <data name="loci_density_plot" label="${tool.name} on ${on_string} - Loci Density Plots" format="pdf">
|
|
|
180 <filter>output_files == "both" or output_files == "plots"</filter>
|
|
|
181 </data>
|
|
|
182
|
|
|
183 <data name="cluster_table" label="${tool.name} on ${on_string} - Cluster.tsv" format="tsv">
|
|
|
184 <filter>output_files == "both" or output_files == "tables"</filter>
|
|
|
185 </data>
|
|
|
186
|
|
|
187 <data name="loci_table" label="${tool.name} on ${on_string} - Loci.tsv" format="tsv">
|
|
|
188 <filter>output_files == "both" or output_files == "tables"</filter>
|
|
|
189 </data>
|
|
|
190
|
|
|
191 <data format="png" name="ccf_plot">
|
|
|
192 <discover_datasets pattern="__designation__" ext="png" visible="true" />
|
|
|
193 </data>
|
|
|
194
|
|
|
195 </outputs>
|
|
|
196
|
|
|
197 <help>
|
|
|
198 <![CDATA[
|
|
|
199
|
|
|
200 .. class:: infomark
|
|
|
201
|
|
|
202 **What this tool does**
|
|
|
203
|
|
|
204 PyClone is statistical model and software tool designed to infer the prevalence of point mutations in heterogeneous cancer samples. The input data for PyClone consists of a set read counts from a deep sequencing experiment, the copy number of the genomic region containing the mutation and an estimate of tumour content.
|
|
|
205
|
|
|
206 .. _PyClone: https://bitbucket.org/aroth85/pyclone/wiki/Home
|
|
|
207
|
|
|
208 -----
|
|
|
209
|
|
|
210 .. class:: infomark
|
|
|
211
|
|
|
212 **Input**
|
|
|
213
|
|
|
214 This tool take a tab seperated value file with six required fields. This file can also be generated using the 'Generate PyClone input' tool which takes a takes a MAF file and a segments.txt file produced by sequenza.
|
|
|
215
|
|
|
216 ** The required fields of the header in this input file are**
|
|
|
217
|
|
|
218 1. mutation_id - A unique ID to identify the mutation. Good names are thing such a the genomic co-ordinates of the mutation i.e. chr22:12345. Gene names are not good IDs because one gene may have multiple mutations, in which case the ID is not unique and PyClone will fail to run or worse give unexpected results. If you want to include the gene name I suggest adding the genomic coordinates i.e. TP53_chr17:753342.
|
|
|
219
|
|
|
220 2. ref_counts - The number of reads covering the mutation which contain the reference (genome) allele.
|
|
|
221
|
|
|
222 3. var_counts - The number of reads covering the mutation which contain the variant allele.
|
|
|
223
|
|
|
224 4. normal_cn - The copy number of the cells in the normal population. For autosomal chromosomes this will be 2 and for sex chromosomes it could be either 1 or 2. For species besides human other values are possible.
|
|
|
225
|
|
|
226 5. minor_cn - The minor copy number of the cancer cells. Usually this value will be predicted from WGSS or array data.
|
|
|
227
|
|
|
228 6. major_cn - The major copy number of the cancer cells. Usually this value will be predicted from WGSS or array data.
|
|
|
229
|
|
|
230 -----
|
|
|
231
|
|
|
232 .. class:: infomark
|
|
|
233
|
|
|
234 **Cite**
|
|
|
235
|
|
|
236 Roth et al. PyClone: statistical inference of clonal population structure in cancer PMID: 24633410
|
|
|
237
|
|
|
238 ]]>
|
|
|
239 </help>
|
|
|
240
|
|
|
241 </tool> |
