How to compare two sequence logo?
2
0
Entering edit mode
2.4 years ago
Andrew Liu ▴ 10

Hi, I want to compare motifs between two sets of sequences by drawing their sequence logo, should I align them together or separately?

sequence-logo motif • 1.6k views
ADD COMMENT
1
Entering edit mode

Can you expand on your question more? It's unclear what data you have and what you want to do with it.

ADD REPLY
0
Entering edit mode

should I align them together

If you want to compare then you will need at least two things to do so. So finding motifs separately would be the first thing to do. Then you may want to follow-up by doing a single alignment (if the results indicate that the two groups may share a motif) to confirm.

ADD REPLY
0
Entering edit mode
2.4 years ago
Mensur Dlakic ★ 28k

We don't know whether you are talking about protein or DNA/RNA motifs, so my answer will be general.

If you are talking about finding motifs in a common sense use of that term, and especially if those two groups of sequences are from different proteins or promoters, you may not be able to align them at all using global alignments of both groups. The whole point of motifs is that they are short stretches of sequence in what can otherwise be long and unrelated sequences.

Assuming you have enough signal to detect motifs individually in both groups of sequences, and that you can find all the motifs that are present (there may be more than one), it should be obvious if both groups share the same motif even if you do the motif finding separately. Combining the two groups before motif finding will likely give you the same result, but only if proper motif finding is done. Without knowing exactly what you have and how (un)related the two sequence groups are, there is no guarantee that they can be globally aligned in such a way where a common motif will be obvious.

ADD COMMENT
0
Entering edit mode
2.4 years ago
  1. Create two PWMs for each set separately,
  2. Use difflogo to generate sequence logo for each set, followed by sequence logo comparison and finally, export comparison figure (diff logo)

https://github.com/mgledi/DiffLogo and r package: https://bioconductor.org/packages/release/bioc/html/DiffLogo.html

Image from https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4650857/: diff_logo.

ADD COMMENT

Login before adding your answer.

Traffic: 2506 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6