|
88
|
1 <tool id="QDNAseq-regioning" name="QDNAseq-regioning" version="0.0.1" force_history_refresh="True">
|
|
|
2
|
|
|
3 <requirements>
|
|
|
4
|
|
|
5 <!-- R 3.1.0 dependency will be used instead when available, now default R is used, see command -->
|
|
|
6 <!-- <requirement type="package" version="3.1.0">R</requirement> -->
|
|
|
7 <!-- <requirement type="package" version="1.2.2">qdnaseq</requirement> -->
|
|
|
8
|
|
|
9 </requirements>
|
|
|
10
|
|
|
11 <description>Perform regioning on QDNAseq data</description>
|
|
|
12
|
|
|
13 <!-- change to /full/path/to/Rscript if required (eg /ccagc/lib/R/R-3.1.0/bin/Rscript) -->
|
|
|
14 <command interpreter="Rscript">
|
|
|
15 QDNAseq-regioning.R
|
|
|
16 $cfg
|
|
|
17 </command>
|
|
|
18
|
|
|
19 <version_command interpreter="Rscript">QDNAseq-version.R</version_command>
|
|
|
20
|
|
|
21 <stdio>
|
|
|
22 <!-- Anything higher than 0 means the R script didnt finish (correctly) -->
|
|
|
23 <!-- Because different R packages deal with err/warn differently unable to waterproof this -->
|
|
|
24 <exit_code range="1:" level="fatal" description="R script finished too early, check log" />
|
|
|
25 </stdio>
|
|
|
26
|
|
|
27 <inputs>
|
|
|
28
|
|
|
29 <!-- ==================== -->
|
|
|
30 <!-- General inputs -->
|
|
|
31 <!-- ==================== -->
|
|
|
32
|
|
|
33 <!-- Job name: must contain non-whitespace chars -->
|
|
|
34 <param name="output_name" type="text" optional="false" label="Analysis/ouput name" help="Supply a name for the outputs to remind you what they contain" value="TEST">
|
|
|
35 <!-- <validator type="empty_field" /> -->
|
|
|
36 <validator type="regex" message="No whitespace characters allowed">^[^\s\\]*$</validator>
|
|
|
37 </param>
|
|
|
38
|
|
|
39 <!-- ==================== -->
|
|
|
40 <!-- Input RDS -->
|
|
|
41 <!-- ==================== -->
|
|
|
42 <param name="rds_file_path" type="data" optional="False" format="rds" label="Input RDS file" help="RDS file should contain a QDNAseq R object (output of QDNAseq tool)" />
|
|
|
43
|
|
|
44 <param name="av_error" size="10" type="float" value="0.00001" optional="False" label="Average Error" help="By default low so that most if not all possible breakpoints in the data are kept, output will then contain smaller regions." />
|
|
|
45
|
|
|
46 </inputs>
|
|
|
47
|
|
|
48 <!-- ==================== -->
|
|
|
49 <!-- Config file to pass params to R script -->
|
|
|
50 <!-- ==================== -->
|
|
|
51 <configfiles>
|
|
|
52 <configfile name="cfg">
|
|
|
53 ## Desc: this file is sourced in QDNAseq-region.R wrapper script
|
|
|
54 ## as means to pass all galaxy params to R
|
|
|
55
|
|
|
56 "${output_name}" -> outputName
|
|
|
57 "${output_file}" -> outputFile
|
|
|
58 "${rds_file_path}" -> rdsFilePath
|
|
|
59 as.double( "${av_error}" ) -> avError
|
|
|
60
|
|
|
61 </configfile>
|
|
|
62 </configfiles>
|
|
|
63
|
|
|
64 <!-- ==================== -->
|
|
|
65 <!-- One text file as output -->
|
|
|
66 <!-- ==================== -->
|
|
|
67 <outputs>
|
|
|
68 <data format="tabular" name="output_file" label="QDNAseq: ${output_name} regions file" />
|
|
|
69 </outputs>
|
|
|
70
|
|
|
71 <help>
|
|
|
72
|
|
|
73 **Introduction**
|
|
|
74
|
|
|
75 This tool is a wrapper for the "CGHregions" function of the R Bioconductor package CGHregions_ and works on objects of the R bioconductor package QDNAseq_
|
|
|
76
|
|
|
77 -----
|
|
|
78
|
|
|
79 **What it does**
|
|
|
80
|
|
|
81 **Input:** The input for this tool is a QDNAseq R object in RDS (R data structure) format, an (optional) output file of the QDNAseq galaxy tool. **Output:** Running this regioning tool overlaps the called called segments from all samples present in the input object. The output contains the merged segments with their call values for each sample.
|
|
|
82
|
|
|
83 -----
|
|
|
84
|
|
|
85 **Output examples**
|
|
|
86
|
|
|
87 *Example TSV output:*
|
|
|
88
|
|
|
89 ::
|
|
|
90
|
|
|
91 Chromosome Start End Nclone AveDist sample1 sample2
|
|
|
92 1 6000001 7000000 42 0 1 -1
|
|
|
93 1 7000001 8000000 11 0 0 -1
|
|
|
94
|
|
|
95 -----
|
|
|
96
|
|
|
97 .. class:: warningmark
|
|
|
98
|
|
|
99 As there is no R 3.1.0 galaxy-package yet (a requirement for QDNAseq) that works with all requirements, the **dependencies** need to be installed by hand and available to the user under which galaxy runs: R (>= 3.1.0) and bioconductor package QDNAseq (>= 1.2.2). In case the path to this R installation is not "R", also the wrapper xml must be updated to include the correct path during installation of this tool.
|
|
|
100
|
|
|
101 -----
|
|
|
102
|
|
|
103 **Citation**
|
|
|
104
|
|
|
105 For the underlying QDNAseq R package please cite:
|
|
|
106 Scheinin I, Sie D, Bengtsson H, van de Wiel MA, Olshen AB, van Thuijl HF, van Essen HF, Eijk PP, Rustenburg F, Meijer GA, Reijneveld JC, Wesseling P, Pinkel D, Albertson DG and Ylstra B (2014). “DNA copy number analysis of fresh and formalin-fixed specimens by shallow whole-genome sequencing with identification and exclusion of problematic regions in the genome assembly.” Genome Research. doi:10.1101/gr.175141.114.
|
|
|
107
|
|
|
108 See also the bioconductor documentation of QDNAseq_ and CGHregions_
|
|
|
109
|
|
|
110 .. _CGHregions: http://http://www.bioconductor.org/packages/release/bioc/html/CGHregions.html
|
|
|
111 .. _QDNAseq: http://www.bioconductor.org/packages/release/bioc/html/QDNAseq.html
|
|
|
112
|
|
|
113 .. image:: LGG150_copynumber.png
|
|
|
114 .. image:: LGG150_copynumberSegmented.png
|
|
|
115
|
|
|
116 </help>
|
|
|
117
|
|
|
118 </tool>
|
