How can we make this plot? (sequence logo)
4
2
Entering edit mode
6.4 years ago
zx8754 10k

The plot is from this publication, see page 2 the bottom of the plot with ACTG letters in different colours and sizes.

  1. Could anyone clarify in layman's terms what it represents?
  2. More importantly, how do we plot it using R or any other tools?

Here is the plot, if you can't access the paper:

Found related post at SO: Plotting a “sequence logo” using ggplot2?

R plot sequence-logo logo seqlogo • 5.6k views
ADD COMMENT
9
Entering edit mode
6.4 years ago

That type of plot is called Sequence Logo. The X axis represent the nucleotide positions of a sequence, and the Y the entropy of each possible nucleotide. The size of each letter indicates how frequently it is found at that position: for example position 8 is almost always an A, while position 5 can be a T or C.

You can create sequence logos from a web interface using a popular application called WebLogo, or using R packages such as SequenceLogo.

ADD COMMENT
0
Entering edit mode

Thank you, "sequence logo" was the term I was missing :)

ADD REPLY
4
Entering edit mode
6.4 years ago
Chirag Nepal ★ 2.3k

The plot shows the comparison of the reference human sequence with N-number of samples (patients). The samples turn to have variation of sequences in that region. Now from these samples, you can compare the probability how often a given nucleotide is mismatched/mutated in sapmles. If the letter is big it means the sequence is invariant.

You can use weblogo either in their webpage or in R to make that figure.

Or if you have Chip-seq peaks, you can do motif analysis and make those plots using HOMER tools.

ADD COMMENT
0
Entering edit mode

X axis 1 to 10, is it SNP position, why 1 to 10 and not genomic positions?

ADD REPLY
3
Entering edit mode
6.4 years ago
Zhilong Jia ★ 2.0k

Here is motiflogo, which can make a SNP-specific motif logo representation. See figure, please.

 

ADD COMMENT
0
Entering edit mode

You might want to add this as an answer to this post, too. Plotting a “sequence logo” using ggplot2?

ADD REPLY
1
Entering edit mode

Thanks. Done.

ADD REPLY
2
Entering edit mode
6.4 years ago
benaneely ▴ 70

Maybe not as helpful for simply plotting, but in proteomics we often use IceLogo.

ADD COMMENT

Login before adding your answer.

Traffic: 2389 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