Statistical alignment of retropseudogenes and their functional paralogs.
We describe a model for the sequence evolution of a processed pseudogene and its paralog from a common protein-coding ancestor. The model accounts for substitutions, insertions, and deletions and combines nucleotide- and codon-level mutation models. We give a dynamic programming method for calculating the likelihood of homology between two sequences in the model and describe the accompanying alignment algorithm. We also describe how ancestral codons can be computed when the same gene produced multiple pseudogene homologs. We apply our methods to the evolution of human cytochrome c.
Csurös
Miklós
M
Department of Computer Science and Operations Research, Université de Montréal, succursale Centre-Ville, Montréal, Québec H3C 3J7, Canada. csuros@iro.umontreal.ca
Miklós
István
I
D001482:Base Composition; D003062:Codon; D045304:Cytochromes c / Q000235:genetics*; D019143:Evolution, Molecular; D008957:Models, Genetic; D011544:Pseudogenes; D012689:Sequence Homology, Nucleic Acid
D003062:Codon; D045304:Cytochromes c