Mostra la font per a Models d'evolució de l'ADN

S'ha proposat un nombre de '''models d'evolució de l'ADN''' de [[cadena de Màrkov|Màrkov]]. Aquests [[model de substitució|models de substitució]] difereixen pel que fa als paràmetres emprats per a descriure les taxes de [[mutació]], és a dir en les que els [[nucleòtid]]s reemplacen a altres amb el pas de les [[generació|generacions]]. Aquests models s'utilitzen freqüentment en anàlisis de [[filogènia molecular]]. En particular, s'utilitzen en el càlcul de la versemblança d'un arbre, en aproximacions de l'estimació dels arbres [[inferència Bayesiana|bayesians]] i de [[màxima versemblança]] (''Maximum Likehood''), i s'utilitzen per a estimar la distància evolutiva entre seqüències a partir de les diferències observades entre aquestes.

Aquests models són descripcions fenomenològiques de l'evolució de l'ADN com una cadena de quatre estats discrets. Aquests models de Màrkov no representen explícitament el mecanisme de mutació ni l'acció de la [[selecció natural]]. Més aviat descriuen les taxes relatives dels diferents canvis. Per exemple, biaixos de mutació i [[selecció purificadora]] afavoreixen els cansi conservatius i probablement ambdós són responsables de la taxa relativament elevada de [[transició genètica|transicions]] en comparació amb les [[tranversió|transversions]] en les seqüències que evolucionen. Tanmateix, el model de [[Kimura]] (K80) que es descriu més avall senzillament intenta capturar l'efecte de les dues forces en un paràmetre que reflecteix la taxa relativa de transicions i transversions.

Les anàlisis evolutives de seqüències es duen a terme en una àmplia varietat d'escales temporals. Per això és convenient expressar aquests models en termes de taxes instantànies de canvis entre diferents estats (les matrius Q de sota). Si donem un estat inicial (ancestral) a la primera posició, la matriu del model Q i la longitud de la branca expressant el nombre esperat de canvis que s'han donat des de l'ancestre, aleshores podem derivar la probabilitat de la seqüència descendent tenint cadascun dels quatre estadis. Els detalls matemàtics d'aquesta transformació des de taxa-matriu a matriu de probabilitats en descriuen a la secció de [[Model matemàtic|models matemàtics]] de substitució de la pàgina de [[models de substitució]]. Expressant en termes de taxes instantànies de canvi podem evitar estimar grans nombres de paràmetres per a cada branca en un [[arbre filogenètic]] (o cada comparació si l'anàlisi inclou moltes comparacions de seqüències per parelles).

== Evolució de l'ADN com una cadena de Màrkov contínua en el temps ==
=== Cadenes de Màrkov contínues en el temps ===
Les cadenes de Màrkov contínues en el temps tenen les típiques matrius de transició que són, a més a més, parametritzades pel temps, <math> t\ </math>. Específicament, si <math> E_1,\ldots,E_4\ </math> són els estats, aleshores la matriu de transició

: <math> P(t) = \big(P_{ij}(t)\big)</math> en la que cada entrada, <math> P_{ij}(t)\ </math> fa referència a la probabilitat que l'estat <math> E_i\ </math> canviï a l'estat <math> E_j\ </math> en el temps <math> t\ </math>.

'''Exemple:''' Ens agradaria obtenir el model del procés de substitució en seqüències d'ADN (p. ex.: [[Jukes–Cantor]], Kimura, etc.) en una moda contínua en el temps. Les matrius de transició corresponents seran com les que segueixen:

: <math>
P(t) = \begin{pmatrix} p_{AA}(t) & p_{GA}(t) & p_{CA}(t) & p_{TA}(t) \\
 p_{AG}(t) & p_{GG}(t) & p_{CG}(t) & p_{TG}(t) \\
 p_{AC}(t) & p_{GC}(t) & p_{CC}(t) & p_{TC}(t) \\
 p_{AT}(t) & p_{GT}(t) & p_{CT}(t) & p_{TT}(t)
 \end{pmatrix} </math>

en les que els blocs de dalt a l'esquerra i a baix a la dreta 2 &times;&nbsp;2 corresponen a ''probabilitats de transició'' i les els blocs de dalt a la dreta i a baix a l'esquerra 2 &times;&nbsp;2 corresponen a les ''probabilitats de transversió''.

'''Supòsit:''' Si en algun moment <math> t_0\ </math>, la cadena de Màrkov està a l'estat <math> E_i\ </math>, aleshores la probabilitat que al temps <math> t_0+t\ </math>, sigui un estat <math>E_j\ </math> depèn només de <math>i\ </math>, <math>j\ </math> i <math>t\ </math>. Aleshores això ens permet escriure la probabilitat com a <math>p_{ij}(t)\ </math>.

'''Teorema:''' Les matrius de transició contínues en el temps satisfan:

::<math>P(t+\tau) = P(t)P(\tau)\ </math>

== Models més comuns d'evolució de l'ADN ==
=== Model JC69 (Jukes i Cantor, 1969)<ref>{{ref-llibre |autor=Jukes TH and Cantor CR |any=1969 |títol=Evolution of Protein Molecules |pàgines=21–132 |editorial=Academic Press |lloc=Nova York}}</ref> ===
JC69 és el model de substitució més simple. Hi ha diferents supòsits. Assumeix freqüències de bases i [[taxa de mutació|taxes de mutació]] són iguals (<math>\pi_1 = \pi_2 = \pi_3 = \pi_4 = {1\over4}</math>). L'únic paràmetre d'aquest model és per tant <math>\mu</math>, la taxa general de substitució. Com s'ha mencionat prèviament, aquesta variable esdevé constant quan normalitzem la taxa mitjana a 1.

: <math>Q = \begin{pmatrix} {*} & {\mu\over 4} & {\mu\over 4} & {\mu\over 4} \\ {\mu\over 4} & {*} & {\mu\over 4}& {\mu\over 4}\\ {\mu\over 4}& {\mu\over 4}& {*} & {\mu\over 4}\\ {\mu\over 4}& {\mu\over 4}& {\mu\over 4}& {*} \end{pmatrix}</math>

: <math> P= \begin{pmatrix} {{1\over4} + {3\over4}e^{-t\mu}} & {{1\over4} - {1\over4}e^{-t\mu}} & {{1\over4} - {1\over4}e^{-t\mu}} & {{1\over4} - {1\over4}e^{-t\mu}} \\\\ {{1\over4} - {1\over4}e^{-t\mu}} & {{1\over4} + {3\over4}e^{-t\mu}} & {{1\over4} - {1\over4}e^{-t\mu}} & {{1\over4} - {1\over4}e^{-t\mu}} \\\\ {{1\over4} - {1\over4}e^{-t\mu}} & {{1\over4} - {1\over4}e^{-t\mu}} & {{1\over4} + {3\over4}e^{-t\mu}} & {{1\over4} - {1\over4}e^{-t\mu}} \\\\ {{1\over4} - {1\over4}e^{-t\mu}} & {{1\over4} - {1\over4}e^{-t\mu}} & {{1\over4} - {1\over4}e^{-t\mu}} & {{1\over4} + {3\over4}e^{-t\mu}} \end{pmatrix}</math>

Quan la longitud de branca, <math>\nu</math>, és mesurada en el nombre esperat de canvis per lloc:
: <math>P_{ij}(\nu) = \left\{
\begin{array}{cc}
{1\over4} + {3\over4}e^{-4\nu/3} & \mbox{ if } i = j \\
{1\over4} - {1\over4}e^{-4\nu/3} & \mbox{ if } i \neq j 
\end{array}
\right.</math>

La estimació Jukes-Cantor de la [[distància genètica]] (en termes de nombre de canvis esperats) entre dues seqüències ve donada per
:<math>\hat{d}=-{3\over4} \ln({1-{4\over3}p})</math>
en que <math>p</math> és la proporció de lloc que difereixen entre dues seqüències. La <math>p</math> en aquesta fórmula està freqüentment referida com a distància <math>p</math>. És un estadístic suficient per a la correcció calculada de la distància Jukes-Cantor, però no és suficient per al càlcul de la distància evolutiva sota models més complexes que segueixen (la <math>p</math> emprada en les següents fórmules no és idèntica a la "distància <math>p</math>").

== Referències ==
{{Referències}}

== Enllaços externs ==
* [http://scit.us/projects/dawg DAWG: DNA Assembly With Gaps] {{Webarchive|url=https://web.archive.org/web/20160815184338/http://scit.us/projects/dawg/ |date=2016-08-15}} — free software for simulating sequence evolution
* {{cite journal |author=Gu X, Li W |year=1992 |title=Higher rates of amino acid substitution in rodents than in man |journal=Molecular Phylogenetics and Evolution |volume=1 |pages=211–214 |doi=10.1016/1055-7903(92)90017-B |pmid=1342937}}
* {{cite journal |author=Li W-H, Ellsworth DL, Krushkal J, Chang BH-J, Hewett-Emmett D |year=1996 |title=Rates of nucleotide substitution in primates and rodents and the generation-time effect hypothesis |journal=Molecular Phylogenetics and Evolution |volume=5 |pages=182–187 |doi=10.1006/mpev.1996.0012 |pmid=8673286}}

{{ORDENA:Models D'Evolucio De L'ADN}}
[[Categoria:Filogènia]]
[[Categoria:Bioinformàtica]]