This is a test version of Biostars. For the public version, visit https://www.biostars.org.
A Pwm With Gapped Alignments In Biopython

I'm trying to generate a Position-Weighted Matrix (PWM) in Biopython from Clustalw multiple sequence alignments. I get a "Wrong Alphabet" error every time I do it with gapped alignments. From reading the documentation, I think I need to utilize the Gapped Alphabet to deal with the '-' character in gapped alignments. But when I do this, it still doesn't resolve the error. Does anyone see the problem with this code, or have a better way to generate a PWM from gapped Clustal alignments?

from Bio.Alphabet import Gapped
alignment = AlignIO.read("filename.clustalw", "clustal", alphabet=Gapped)
m = Motif.Motif()
for a in alignment:
    m.add_instance(a.seq)
m.pwm()
biopython multiple clustalw
• 5,873 views • link
updated by Brad Chapman Brad Chapman • written by RossCampbell RossCampbell

Here's an example how to do it with with unambiguous_dna:

m = motifs.create(list, alphabet=Gapped(IUPAC.unambiguous_dna))

• 1 views • link
written by jykzel • jykzel

2 answers

Paul-Michael's answer is exactly right and you were nearly there with Gapped. Just pass the same alphabet to AlignIO and Motif; here is a working example with a protein alignment:

from Bio import AlignIO, Motif
from Bio.Alphabet import IUPAC, Gapped
alphabet = Gapped(IUPAC.protein)
alignment = AlignIO.read("cw02.aln", "clustal", alphabet=alphabet)
m = Motif.Motif(alphabet)
for a in alignment:
    m.add_instance(a.seq)
print m.pwm()[0]

With a DNA alignment you'd want IUPAC.unambiguous_dna.

• 1 views • link
written by Brad Chapman Brad Chapman

Just quickly: the proximate problem here is that the alphabets of the Motif and and aligned sequence being added have to be the same (see line 44 of _Motif.py). So you could create the motif with the same alphabet, or perhaps try to set it to None, which with avoid the alphabet comparison.

Biopython Alphabets are a tricky subject - I regularly have to pore over the source code to understand how they work, only to forget it all by the next time something odd happens.

• 1 views • link
written by Agapow Agapow

Log in to answer this question.

More posts like this