Wikiomics:Sequence motifs

From OpenWetWare

(Difference between revisions)
Jump to: navigation, search
m (Multiple algorithms / metaservers)
m (Multiple algorithms / metaservers)
Line 211: Line 211:
==Multiple algorithms / metaservers==
==Multiple algorithms / metaservers==
-
* [http://www.comp.nus.edu.sg/~bioinfo/MotifVoter/ MotifVoter] [http://precedings.nature.com/documents/1251/version/1/files/npre20071251-1.pdf PDF] 10! algorithms from National University of Singapore. Web server, result by email.  
+
* [http://compbio.ddns.comp.nus.edu.sg/~edward/MotifVoter2/ MotifVoter] [http://precedings.nature.com/documents/1251/version/1/files/npre20071251-1.pdf PDF] 11! algorithms from National University of Singapore. Web server, result by email.  

Revision as of 20:53, 2 March 2010

Sequence motifs - for more general background see Wikipedias Sequence_motif.


It is important to use few complementary programs as well as allow predictions of several motifs per program. There is threshold when it comes to number of sequences used, above which there is no improvement in sensitivity.


Contents

Motif finding programs

There are obvious trade offs between speed and accuracy/sensitivity when running these programs. Few rules based on review by Hu et al.:

  • avoid using longer sequencess than necessary. It does increase noise signal and increases running time.
  • for programs like MEME there is a plateu of motifs found after about 10 input sequences, so one can randomly select 10 sequences out of a larger set and

then check for occurence of detected motifs in the remaining ones.

  • scores between different programs are incompatible

Tested set

AlignACE

Gibbs sampling algorithm (stochastic), requires multiple runs to get all top hits.

  • Web server: [1] Gibbs sampling algorithm
  • command line C++, requires compilation
  • to run it:
./AlignACE -i GAL.seq > test.ace
  • Accessory programs for motif comparisons and motif finding

CompareACE [2] Compares set of found motifs to a database of TFs (in yeast on the web)


MEME

Web: [3]

  • Command line
meme -nmotifs numberof_motifs_2_find -mod oops  -protein -sequence your_fasta_file -outfile your_fasta_file.meme

Where 'model' could be oops/zoops/anr.

oops: One Occurrence Per Sequence

zoops: Zero or One Occurrence Per Sequence

anr: Any Number of Repetitions


Weeder

./weederlauncher.out your_promoters_file.fasta  MM large A M S t15 

Where MM stands for Mus musculus, HS Homo Sapiens atc. See ./FreqFiles/ directory for more.

The output files: your_promoters_file.mix your_promoters_file.wee (output as text) your_promoters_file.html (output summary in HTML)


Bioprospector

  • Web server: [4]

accepts only FastA files with single sequence line

>sequence1 name
GGTGACGAC sequence1 as ONE LINE
>sequence2 name
GTAGCCTCATG sequence2 as ONE LINE
  • fixed widith of motifs (default 10, range: 4-50)
  • binaries for Linux, Sun and Cygwin.
  • creata background file:
./genomebg.linux -i your_background -o your_background.genomeBG
  • run BioProspector
./BioProspector.linux -n 200 -d 1 -r 30 -i target_fasta_seqs -f your_background.genomeBG -o outputbiop1

Improbizer

Web takes into account position of the motifs Command line utility called ameme

 ./ameme good=your_test_set.fasta bad=set_of_random_promoter_sequences.fasta numMotifs=10 mo tifOutpu=ameme.out outputLogo background=m2


NestedMica

(java + C program, requires compilation ) [5]

  • Calculate background:
makemosaicbg -seqs your_input_sequences.fasta -mosaicClasses 1 -mosaicOrder 1 -out your_sequences_background.sbg
  • Calculate motifs:
motiffinder -seqs your_input_sequences.fasta -backgroundModel your_sequences_background.sbg -numMotifs 2
  • View the motifs:
motifviewer motifs.xms

Output in xms (an XML variant) format.

  • MotifExplorer (Java motif viewer compatible with NestedMica): [6]

as for Jan 22nd 2007 works on MacOS. Problems on Windows and Ubuntu Linux.

SOMBRERO

Self Organising Maps [7]

  • create background model (takec >15min on 35Mbytes fasta file/2.26GHz Pentium) :
perl ./BackExtract.pl -seq your_background_sequences.fasta
  • Create SOM with motif lenght 8, 10, 12, 14 &16 nucleotides:
./SOMBRERO -t your_target_sequences.fasta -b out.back -lm 8 16 -out target_sombrero.out
  • View the SOM (requires installation of Tkperl)
perl ./SOMBREROView.pl  target_sombrero.ou

MDScan

  • uses restricted fasta format!
>sequence1 name
ATGGTGACGAC sequence1 as ONE LINE
>sequence2 name
GTAGCCTCATG sequence2 as ONE LINE
  • requires sequence background file
./genomebg -i inputSequenceFile -o outputDistributionFile
  • running it:
./MDscan -i inseq -w 15 -f yeast_all.bg -t 10 -c 80 -r 10 -n 5 -g 1 

find motif of width 15 from sequence file inseq, use yeast_all.bg as the background distribution. Find candidate motifs from top 10 sequences, and refine 5 iterations from the top 0 sequences. Report the final best 5 motifs to stdout, and do not print out progress messages on the way.

YMF (not working)

  • C++ program, requires compilation
  • running it:
./stats stats.config 800 6 ../ymftables/yeast -sort ../examples/abf/abf1 ../examples/abf/cha1

segfaults.

DIPS

Obtainable from author. Requires compilation

  • example ( ~50 Kbp takes 70mins)
dmotif -positive positive.fna -negative negative.fna -len 9 -bkg fly_background.fna -niter 5 -nmotif 1 > dmotifoutput 2> dmotiflog

Change '-nmotif 1' to a '-nmotif 10' if you want to get top ten motifs.

Novel algorithms

more sensitive than MEME according to authors


  • The AMADEUS Motif Discovery Tool (whole platform)

http://acgt.cs.tau.ac.il/amadeus/

  • GAME (java, genetic algorithm) [8] (authors page)

paper http://bioinformatics.oxfordjournals.org/cgi/content/abstract/22/13/1577

  • MotifCut (maximum density subgraphs) [9]

Paper: http://bioinformatics.oxfordjournals.org/cgi/reprint/22/14/e150 motif lenght: fixed, between 6,31

  • GibbsST [10] Not working yet

Paper(PDF): [11]

  • THEME (ChIP-chip only??)

http://bioinformatics.oxfordjournals.org/cgi/content/abstract/22/4/423

Ortholog sequences

Paper: [13]

Paper: [15]

Paper: [17]

Multiple algorithms / metaservers

  • MotifVoter PDF 11! algorithms from National University of Singapore. Web server, result by email.


  • Credo [18] (broken link)

Visualisation of AlignACE, DIALIGN, FootPrinter, MEME and MotifSampler results. Paper: [19]

  • BEST

Download: [20] (Qt on linux) BEST includes four commonly used motif-finding programs: AlignACE, BioProspector, CONSENSUS and MEME, as well as the optimization program BioOptimizer. Paper: [21]

Combining/comparing motifs

  • MatCompare
    • set of user provided PWMs: [22]
    • single PWM against TRANSFAC or JASPAR [23]
    • simple correlation vlue between each apir of motifs
  • CompareACE [24] supporting program for AlignAce

Handles 12 various output formats from a wide range of motif finding programs. Compares these motifs to known TFBS from JASPAR and other, also user-defined databases of motifs. Converts the output to intermediate format accessible on the server under "X motifs loaded".

APIs

Python package with interfaces to i.e. MEME, AlignACE, Co-Bind, and FootPrinter Paper [26]

Python/C++ package with Interfaces to MDscan, AlignACE, and MEME. Paper: [27]

Databases

  • ORegAnno database as a dynamic collection of literature-curated regulatory regions, transcription factor binding sites and regulatory mutations (polymorphisms and haplotypes).

Paper: [28]


To Read

See also: http://www.connotea.org/user/darked89/tag/motif


External courses/tutorials

  • EMBL 2001 [31]
  • Harward (AlignAce) [32]

stuff to incorporate

  • Melina II Metaserver (uses four out of five programs: CONSENSUS, MEME, Gibbs Sampler, MDScan and Weeder) by Kenta Nakai @Tokyo University. Paper

Credits

Template:Credits


Personal tools