Hi all!
I assembled two different genomes and wanted to see how similar they are on both nucleotide and protein levels so I aligned their nucleotide and translated nucleotide sequences. Here are the results I obtained:
Nucleotide identity=90% Protein identity=57%
How would one make sense of this high nucleotide yet low protein identity result? I have been doing a lot of reading and it seems that if the species are close its better to use the DNA sequence to compare, and I believe these two species should be fairly close. However, I am still confused as to why the values would differ so much.
Thanks for your input!
•
written
by
el97004 •
1 answer
Hi,
Little changes on nucleotide level can lead to drastic changes on protein level. In a worst case scenario, you might introduce a frame shift with a mutation in a gene's 5' region which lead to a totally different products. You'll have in such a case nearly 100%identity on nucleotide level but nearly none for the protein.
Depending on your species you have more or less "junk DNA" intergenic region, introns, etc. These non-coding regions can increase the overall nucleotide identity, but not that of the proteins.
Cheers,
Michael
Log in to answer this question.
-
ANI for inter-species analysis and genome tree phylogeny
written by blur
-
Check Conservation of a Given Protein Domain
written by DVA
-
In Jalview, how can I colour just the nucleotides that are different between sequences
written by SaltedPork
-
primer desing for qPCR on de novo assembly output
written by pablo61991 •
-
Phylogeny Analysis Data
written by l.souza
-
ortholog finding with unsequenced species
written by Mehmet
-
How to remove transcripts that have poor alignment scores in exonerate analysis
written by Ginsea Chen
-
CORE genes detection in two closely related species
written by Mehmet
-
Multiple Alignment: Protein Or Nucleotide Sequence?
written by tomas.iesmantas •
-
Est Sequence Database: Which Frame To Choose ?
written by Woa
There are lots of reasons for this, and all else being equal this is to be expected.
You need to clarify whether these are DNA sequences of genes or the whole genome etc.
Sorry I should have clarified. Whole genomes!
It doesn't make any sense to translate the whole genome, and consequently even less to align/compare them.
Exactly!! Only translate and compare protein-coding regions. For non-coding regions, DNA similarity can be high but when ERRONEOUSLY translated, the "protein" sequences could be from different frames and therefore very low similarity. Again, only translate and compare protein-coding regions.