Non-Synonmous Snps
1
0
Entering edit mode
13.1 years ago
Elena ▴ 250

Hi Friends,

I have a doubt. For finding out the non-synonymous SNPs between various organisms, I have translated the nucleotide sequence to protein sequence, found out the mismatches in the protein sequences between various organisms, and multiplied the the mismatch locations in the protein sequences by three to get the nucleotide position.

Now my doubt is if the protein sequences, which I have taken, should be translated from the "multiple alignment of nucleotide sequences" or from "the normal nucleotide sequences and then aligned the multiple protein sequences to check the mismatches"?

snp non mutation multiple • 3.2k views
ADD COMMENT
1
Entering edit mode

I think you are not really using the completely correct definition of SNP. A SNP is connected to a population and the frequency of alleles in those, then when at least two alleles exist and one of them is found above a certain frequency threshold in the population this is called a SNP. Did you compare two organisms published genome sequences for mismatches? That is a slightly different topic.

ADD REPLY
0
Entering edit mode

you are right Michael. I've edited again my answer in order to reflect the possibility that Elena may not be talking about SNPs but about sequence mismatches.

ADD REPLY
3
Entering edit mode
13.1 years ago

I see here a slight misconception of the term "synonymous SNP", so I think this should be clarified in advance. the strict meaning of this term is of a SNP that changes a base on the DNA sequence which does not affect the amino acid translation, so the protein stays the being same (definition #1). but I think you are using it in a slightly different way, trying to define a SNP that generates the same kind of amino acid change on all your organisms (definition #2). I will try to address both cases.

If you are looking for non-synonymous SNPs in their strict meaning (definition #1) you need to translate each DNA sequence into its coded protein considering the SNPs' alleles, and looking at the amino acid changes get the SNPs that affected such protein translation. this should be done for each sequence separately, since the strict definition doesn't depend on anything else.

But if you are looking for SNPs that create different protein sequences on different organisms then you will need to perform the alignments you mention. choosing whether to align DNA or protein sequence shouldn't matter, as long as your translating script takes into account the results of the alignment, but if you are asking for advice I see aligning protein sequences and check for mismatches on them a little bit easier. but of course the SNP position wouldn't necessarily have to be amino acid position multiplied by 3, as it could be placed on any of the 3 positions of the triplet, so you should store that information too in order to reverse position such SNP.

EDIT: I tried to address 2 different ideas I thought you could have in mind, but I was always considering that you were talking about SNPs. after reading Michael's comment on your question I realize that you are maybe not intending to deal with SNPs at all but with sequence mismatches. if that is the case then I would suggest you to translate your DNA sequences into proteins, align both DNA and protein sequences (DNA with DNA and protein with protein, of course), and check whether DNA mismatches correspond to protein mismatches. you may do the position check by evaluating the following expression:

if ( $proteinPosition == int ( ($dnaPosition-1) / 3 ) + 1 )
ADD COMMENT
1
Entering edit mode

isn't non-sysnonymous = DNA different AA different, synoymous = DNA different, AA same? Your second sentence says the opposite.

ADD REPLY
0
Entering edit mode

whoops! I was indeed defining "synonymous SNP", instead of "non-synonymous SNP". thanks for the correction, I've just edited the text.

ADD REPLY
0
Entering edit mode

@jorge abd @michael thanks for the information

ADD REPLY
0
Entering edit mode

@jorge abd @michael thanks for the information

ADD REPLY

Login before adding your answer.

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