changeset 1:e851d9021494 draft

Uploaded
author jbrayet
date Tue, 22 Sep 2015 05:14:36 -0400
parents a0a25763b618
children dce9495ac542
files peakMotifs_wrapper.xml
diffstat 1 files changed, 239 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/peakMotifs_wrapper.xml	Tue Sep 22 05:14:36 2015 -0400
@@ -0,0 +1,239 @@
+<!--Peak Motifs - developed by Jocelyn Brayet <jocelyn.brayet@curie.fr>
+Copyright (C) 2015  Institut Curie
+
+This program is free software: you can redistribute it and/or modify
+it under the terms of the GNU General Public License as published by
+the Free Software Foundation, either version 3 of the License, or
+(at your option) any later version.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+GNU General Public License for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program.  If not, see <http://www.gnu.org/licenses/>.-->
+<tool id="rsat_peak_motifs" name="RSAT Peak Motifs" version="1.0">
+  <description>Create Logo images from peaks</description>
+  <command interpreter="python">
+    peakMotifs_wrapper.py -test $fastaFile -outGalaxy $htmlFile
+    #if $control_conditional.controle_file == "yes":
+      -control $control_conditional.controlFastaFile
+    #end if
+    -server $serverChoice.server
+    -motif_db $serverChoice.database
+    #if $advanced_parameters.adv_param == "show":
+      #if $advanced_parameters.restrict.restrictPeaks == "yes":
+        -top_peaks $advanced_parameters.restrict.top_peaks
+        -max_seq_length $advanced_parameters.restrict.max_seq_length
+      #end if
+      #if $advanced_parameters.motif.motifDiscovery == "yes":
+        #if $advanced_parameters.motif.algorithms.value: 
+          -disco $advanced_parameters.motif.algorithms
+        #end if
+        -max_motif_number $advanced_parameters.motif.max_motif_number
+        #if $advanced_parameters.motif.oligo_lengths.oligomer_lengths == "yes":
+          -min_length $advanced_parameters.motif.oligo_lengths.min_length
+          -max_length $advanced_parameters.motif.oligo_lengths.max_length
+        #end if
+      #end if
+      #if $advanced_parameters.bedFile.resultBedFile == "yes":
+        -source galaxy
+        -outGalaxy2 $bedFile
+      #end if
+    #end if
+  </command>
+  <inputs>
+    <param name="fastaFile" type="data" format="fasta" label="Peak sequences (FASTA file)" />
+    <conditional name="control_conditional" >
+      <param name="controle_file" type="select" label="Control file" help="" >
+        <option value="no" selected="true">No</option>
+        <option value="yes">Yes</option>
+      </param>
+      <when value="no" />
+      <when value="yes">
+        <param name="controlFastaFile" type="data" format="fasta" label="Control peak sequences (FASTA file)" />
+      </when>
+    </conditional>
+    <conditional name="serverChoice">
+       <param name="server" type="select" label="RSAT server (if Time out error, change server)" >
+         <option value="fr_ens" selected="true">RSAT Protists</option>
+         <option value="fr_mrs" >RSAT Fungis</option>
+         <option value="fr_ro" >RSAT Metazoa</option>
+         <option value="fr_mrs_2" >RSAT Teaching</option>
+         <option value="es" >RSAT Plants</option>
+         <option value="mx" >RSAT Prokaryotes</option>
+       </param>
+       <when value="fr_ens">
+         <param name="database" type="select" label="Select a motif database" help="If your motif database of interest is not listed, contact the RSAT team">
+           <options from_file="RSAT_motif_databases_fr_ens.loc">
+             <column name="name" index="0" />
+             <column name="value" index="0" />
+           </options>
+         </param>
+       </when>
+       <when value="fr_mrs">
+         <param name="database" type="select" label="Select an motif database" help="If your motif database of interest is not listed, contact the RSAT team">
+           <options from_file="RSAT_motif_databases_fr_mrs.loc">
+             <column name="name" index="0" />
+             <column name="value" index="0" />
+           </options>
+         </param>
+       </when>
+       <when value="fr_mrs_2">
+         <param name="database" type="select" label="Select an motif database" help="If your motif database of interest is not listed, contact the RSAT team">
+           <options from_file="RSAT_motif_databases_fr_mrs_2.loc">
+             <column name="name" index="0" />
+             <column name="value" index="0" />
+           </options>
+         </param>
+       </when>
+       <when value="fr_ro">
+         <param name="database" type="select" label="Select an motif database" help="If your motif database of interest is not listed, contact the RSAT team">
+           <options from_file="RSAT_motif_databases_fr_ro.loc">
+             <column name="name" index="0" />
+             <column name="value" index="0" />
+          </options>
+        </param>
+       </when>
+       <when value="mx">
+         <param name="database" type="select" label="Select an motif database" help="If your motif database of interest is not listed, contact the RSAT team">
+           <options from_file="RSAT_motif_databases_mx.loc">
+             <column name="name" index="0" />
+             <column name="value" index="0" />
+           </options>
+         </param>
+       </when>
+       <when value="es">
+         <param name="database" type="select" label="Select an motif database" help="If your motif database of interest is not listed, contact the RSAT team">
+           <options from_file="RSAT_motif_databases_es.loc">
+             <column name="name" index="0" />
+             <column name="value" index="0" />
+           </options>
+         </param>
+       </when>
+    </conditional>
+    <conditional name="advanced_parameters" >
+      <param name="adv_param" type="select" label="Advanced Parameters" help="" >
+        <option value="hide" selected="true">Hide</option>
+        <option value="show">Show</option>
+      </param>
+      <when value="hide" />
+      <when value="show">
+        <conditional name="restrict" >
+          <param name="restrictPeaks" type="select" label="Reduce peak sequences" help="" >
+            <option value="no" selected="true">No</option>
+            <option value="yes">Yes</option>
+          </param>
+          <when value="no" />
+          <when value="yes">
+            <param name="top_peaks" type="integer" value="50" label="Number of top sequences to retain"/>
+            <param name="max_seq_length" type="integer" value="500" label="Cut peak sequences (+/- x bp on each side of peak centers)"/>
+          </when>
+        </conditional>
+        <conditional name="motif" >
+          <param name="motifDiscovery" type="select" label="Change motif discovery parameters" help="" >
+            <option value="no" selected="true">No</option>
+            <option value="yes">Yes</option>
+          </param>
+          <when value="no" />
+          <when value="yes">
+            <param name="algorithms" type="select" display="checkboxes" multiple="True" label="Algorithm(s) to use">
+              <option value="oligos" selected="True">Discover over-represented words : oligo-analysis</option>
+              <option value="positions" selected="True">Discover words with a positional bias : position-analysis</option>
+              <option value="local_words">Discover words with local over-representation : local-word-analysis</option>
+              <option value="dyads">Discover over-represented spaced word pairs : dyad-analysis</option>
+              <option value="merged_words">Merge all word lengths for assembly into matrix</option>
+              <validator type="no_options" message="You must select at least one algorithm." />
+            </param>
+            <param name="max_motif_number" type="integer" value="3" label="Number of motifs per algorithm"/>
+            <conditional name="oligo_lengths" >
+              <param name="oligomer_lengths" type="select" label="Change oligomer lengths (oligo-analysis, position-analysis and local-word-analysis)" help="" >
+                <option value="no" selected="true">No</option>
+                <option value="yes">Yes</option>
+              </param>
+              <when value="no" />
+              <when value="yes">
+                <param name="min_length" type="integer" value="6" label="Minimal oligomer length"/>
+                <param name="max_length" type="integer" value="8" label="Maximal oligomer length"/>
+              </when>
+            </conditional>
+          </when>
+        </conditional>
+        <conditional name="bedFile" >
+          <param name="resultBedFile" type="select" label="Generate predicted binding sites (BED file)" help="" >
+            <option value="no" selected="true">No</option>
+            <option value="yes">Yes</option>
+          </param>
+        </conditional>
+      </when>
+    </conditional>
+  </inputs>
+  <outputs>
+    <data format="html" name="htmlFile" label="HMTL RSAT output - PeakMotifs results" />
+    <data format="bed" name="bedFile" label="BED RSAT output - PeakMotifs results" >
+      <filter>((advanced_parameters['adv_param'] == 'show') and (advanced_parameters['bedFile']['resultBedFile'] == 'yes'))</filter>
+    </data>
+  </outputs>
+  <help>
+
+**What it does**
+
+.. class:: infomark
+
+This tool is a pipeline for discovering motifs in massive ChIP-seq peak sequences.
+
+------
+
+*Example with FASTA file (1000 sequences)*:: 
+
+ >mm9_chr3_121848111_121848740_+
+ tggggtgggttccaggacagccaggactgtcacacaaagaaaccttgtctcaaaaaaaca
+ aaaCAAAAAACAAAACAAAACCAAGCAAGCAAGCAAACATGGGCTTAAATctggatacag
+ tggcctttatttctagttccagtgactgggagactgaaacaagagagtcacttgagtaca
+ ggagtgcaaggctagcttgagcaatatagtaagactatctcaaaaTGTGAATTtagatca
+ acagaattgacatcaagaaaaatactgatatcactcaaagcaatctacagattcaacaca
+ atctccatcaacatgacaatgacttccatcaGCATGACAATGACTCCATCAACATGCCAA
+ TGGGCCCCATCAACATAACAATGACCCCTATCATCATGACAATGATCCCCATCAACATGA
+ CAATGACCTCCATCAACATGACAATTACTCCTGTCAACATGCCAATtgttggggttcaga
+ agtcaccctgcaaaccacaagaacactaatctcagtcaagcagggatggtttactgaacg
+ tatatccaaagactgagtgaccaagggaacagctcagactctagagctgaaagctagctg
+ tgcgctggacatttctcggggccaactta
+ >mm9_chr14_86795691_86796311_+
+ CTCAAGGAGGATCCAGAAGTTGGCAGTTTCTGAGGCGAGTCCCATATTCCTCCCCTAAGG
+ GGTCAGGATTTTTCAGGTCTGGGCTCTTCTTGTTCTTTTGACAGCGACATTAATAATTGT
+ ACCAGCTCTCCCCTGGCAGGGCCGCACCACAGAGTAAAGCCTGGAGTAGGAGCTGTGCCC
+ AGCGCAATAATACCAGTTAAATAAGTACGTTCATTACCTCCCAACAGTCAAGGAGTTTAA
+ AATCCGTCAATTCACCCCACATGAGGGAGATTATGTGATTTACATGTTAAAGTGCCCCTG
+ TGGTTTGATTTGCATAGCAAAGACTTTGGGGGCACAGAAACAAAGCATCCGCATGCATGA
+ CAAGAGGACTATTAGCATAGAGAGAGCAGGTTTTCCGACAGCCCAGCCTGGCAAACAATG
+ CTGCACCTTCGCTGCTCGCTGGAGTTTATAGGATTTGACAGTTTCTCATCTGAGGGAGGA
+ GAAGGTTAGGGAGTTGGGGTGGAATGAGGTTCGCTAGATTGGCTGTCTGCCTTAAGCACA
+ ATAATTTGTCTTTCTCTGAAGACCTCCCCTCCTCTTATCACCTTTATCGTTTCTTTCTGA
+ TGTTCATTCAAAGACGTCTT
+
+------
+
+*Partial results :*
+
+.. image:: ${static_path}/images/result_peak_motifs_1.png
+  :height: 350
+  :width: 600
+
+.. image:: ${static_path}/images/result_peak_motifs_2.png
+  :height: 250
+  :width: 600
+
+  </help>
+  <citations>
+    <citation type="bibtex">@article{Medina-Rivera22042015,
+      author = {Medina-Rivera, Alejandra and Defrance, Matthieu and Sand, Olivier and Herrmann, Carl and Castro-Mondragon, Jaime A. and Delerce, Jeremy and Jaeger, Sébastien and Blanchet, Christophe and Vincens, Pierre and Caron, Christophe and Staines, Daniel M. and Contreras-Moreira, Bruno and Artufel, Marie and Charbonnier-Khamvongsa, Lucie and Hernandez, Céline and Thieffry, Denis and Thomas-Chollier, Morgane and van Helden, Jacques}, 
+      title = {RSAT 2015: Regulatory Sequence Analysis Tools},
+      year = {2015},
+      doi = {10.1093/nar/gkv362}, 
+      URL = {http://nar.oxfordjournals.org/content/early/2015/04/21/nar.gkv362.abstract}, 
+      eprint = {http://nar.oxfordjournals.org/content/early/2015/04/21/nar.gkv362.full.pdf+html}, 
+      journal = {Nucleic Acids Research} 
+      }</citation>
+  </citations>
+</tool>