Модели эволюции ДНК

Эта статья включает в себя список общих ссылок , но он остается в значительной степени непроверенным, поскольку в нем отсутствует достаточное количество соответствующих встроенных ссылок . Пожалуйста, помогите улучшить эту статью, добавив более точные цитаты. ( Ноябрь 2010 г. ) ( Узнайте, как и когда удалить этот шаблон сообщения )

Был предложен ряд различных марковских моделей эволюции последовательности ДНК . Эти модели замещения различаются параметрами, используемыми для описания скоростей, с которыми один нуклеотид заменяет другой в ходе эволюции. Эти модели часто используются в молекулярно-филогенетическом анализе . В частности, они используются при вычислении вероятности дерева (в байесовском подходе и подходе максимального правдоподобия к оценке дерева), и они используются для оценки эволюционного расстояния между последовательностями на основе наблюдаемых различий между последовательностями.

Введение [ править ]

Эти модели являются феноменологическим описанием эволюции ДНК как цепочки из четырех дискретных состояний. Эти марковские модели не описывают явно ни механизм мутации, ни действие естественного отбора. Скорее они описывают относительную скорость различных изменений. Например, мутационные смещения и очищающий отбор в пользу консервативных изменений, вероятно, оба ответственны за относительно высокую скорость переходов по сравнению с трансверсиями в развивающихся последовательностях. Однако модель Кимуры (K80), описанная ниже, пытается только уловить эффект обеих сил в параметре, который отражает относительную скорость переходов к трансверсиям.

Эволюционный анализ последовательностей проводится в самых разных временных масштабах. Таким образом, эти модели удобно выражать в терминах мгновенных скоростей изменения между различными состояниями ( матрицы Q ниже). Если нам дано начальное (наследственное) состояние в одной позиции, Q- матрица модели и длина ветви, выражающая ожидаемое количество изменений, которые должны произойти с момента предка, то мы можем вывести вероятность того, что дочерняя последовательность будет иметь каждый из четырех состояния. Математические детали этого преобразования от скорости матрицы к матрице вероятности описаны в математике моделей замещения раздела в модели замещениястраница. Выражая модели в терминах мгновенных скоростей изменений, мы можем избежать оценки большого количества параметров для каждой ветви филогенетического дерева (или каждого сравнения, если анализ включает множество парных сравнений последовательностей).

Модели, описанные на этой странице, описывают эволюцию отдельного сайта в рамках набора последовательностей. Они часто используются для анализа эволюции всего локуса , делая упрощающее предположение, что разные сайты развиваются независимо и одинаково распределены . Это предположение может быть оправдано, если можно предположить, что объекты развиваются нейтрально . Если первичный эффект естественного отбора на эволюцию последовательностей заключается в ограничении некоторых сайтов, то можно использовать модели межсайтовой гетерогенности по скорости. Этот подход позволяет оценить только одну матрицу относительных скоростей замещения и другой набор параметров, описывающих дисперсию общей скорости замещения по сайтам.

Эволюция ДНК как цепь Маркова с непрерывным временем [ править ]

Цепи Маркова с непрерывным временем [ править ]

Цепи Маркова с непрерывным временем имеют обычные переходные матрицы, которые, кроме того, параметризованы временем ,. В частности, если есть состояния, то матрица перехода${\displaystyle t}$${\displaystyle E_{1},E_{2},E_{3},E_{4}}$

${\displaystyle P(t)={\big (}P_{ij}(t){\big )}}$где каждая отдельная запись относится к вероятности того, что состояние изменится на состояние со временем .${\displaystyle P_{ij}(t)}$${\displaystyle E_{i}}$${\displaystyle E_{j}}$${\displaystyle t}$

Пример: мы хотели бы смоделировать процесс замены в последовательностях ДНК ( например, Джакса – Кантора , Кимуры и т. Д. ) В непрерывном режиме времени. Соответствующие матрицы переходов будут выглядеть так:

${\displaystyle P(t)={\begin{pmatrix}p_{\mathrm {AA} }(t)&p_{\mathrm {AG} }(t)&p_{\mathrm {AC} }(t)&p_{\mathrm {AT} }(t)\\p_{\mathrm {GA} }(t)&p_{\mathrm {GG} }(t)&p_{\mathrm {GC} }(t)&p_{\mathrm {GT} }(t)\\p_{\mathrm {CA} }(t)&p_{\mathrm {CG} }(t)&p_{\mathrm {CC} }(t)&p_{\mathrm {CT} }(t)\\p_{\mathrm {TA} }(t)&p_{\mathrm {TG} }(t)&p_{\mathrm {TC} }(t)&p_{\mathrm {TT} }(t)\end{pmatrix}}}$

где верхний левый и нижний правый блоки 2 × 2 соответствуют вероятностям перехода, а верхний правый и нижний левый блоки 2 × 2 соответствуют вероятностям трансверсии .

Предположение: если в какой-то момент цепь Маркова находится в состоянии , то вероятность того, что в какой- то момент она будет в состоянии, зависит только от , и . Это позволяет нам записать эту вероятность как .${\displaystyle t_{0}}$${\displaystyle E_{i}}$${\displaystyle t_{0}+t}$${\displaystyle E_{j}}$${\displaystyle i}$${\displaystyle j}$${\displaystyle t}$${\displaystyle p_{ij}(t)}$

Теорема. Матрицы переходов с непрерывным временем удовлетворяют:

${\displaystyle P(t+\tau )=P(t)P(\tau )}$

Примечание: здесь возможна путаница между двумя значениями слова « переход» . (i) В контексте цепей Маркова переход - это общий термин для перехода между двумя состояниями. (ii) В контексте нуклеотидных изменений в последовательностях ДНК переход - это особый термин для обмена между двумя пуринами (A ↔ G) или двумя пиримидинами (C ↔ T) (дополнительные сведения см. в статье о переходах в генетике ). Напротив, обмен между одним пурином и одним пиримидином называется трансверсией .

Вывод динамики замещения [ править ]

Рассмотрим последовательность ДНК фиксированной длины m, эволюционирующую во времени за счет замены оснований. Предположим, что процессы, за которыми следуют m узлов, являются марковскими независимыми, одинаково распределенными и что процесс постоянен во времени. Для конкретного сайта пусть

${\displaystyle {\mathcal {E}}=\{A,\,G,\,C,\,T\}}$

- набор возможных состояний сайта, и

${\displaystyle \mathbf {p} (t)=(p_{A}(t),\,p_{G}(t),\,p_{C}(t),\,p_{T}(t))}$

их соответствующие вероятности во времени . Для двух различных , позвольте быть скоростью перехода из состояния в состояние . Аналогичным образом , для любого , пусть суммарной скорости изменения от В${\displaystyle t}$${\displaystyle x,y\in {\mathcal {E}}}$${\displaystyle \mu _{xy}\ }$${\displaystyle x}$${\displaystyle y}$${\displaystyle x}$${\displaystyle x}$

${\displaystyle \mu _{x}=\sum _{y\neq x}\mu _{xy}\,.}$

Изменения в распределении вероятностей для малых приращений времени даются выражением${\displaystyle p_{A}(t)}$${\displaystyle \Delta t}$

${\displaystyle p_{A}(t+\Delta t)=p_{A}(t)-p_{A}(t)\mu _{A}\Delta t+\sum _{x\neq A}p_{x}(t)\mu _{xA}\Delta t\,.}$

Другими словами (на частотном языке), частота «s» в определенный момент времени равна частоте «s» во времени минус частота потерянных «s» плюс частота вновь созданных «s».${\displaystyle A}$${\displaystyle t+\Delta t}$${\displaystyle t}$ ${\displaystyle A}$ ${\displaystyle A}$

Аналогично для вероятностей , и . Эти уравнения можно компактно записать как${\displaystyle p_{G}(t)}$${\displaystyle p_{C}(t)}$${\displaystyle p_{T}(t)}$

${\displaystyle \mathbf {p} (t+\Delta t)=\mathbf {p} (t)+\mathbf {p} (t)Q\Delta t\,,}$

где

${\displaystyle Q={\begin{pmatrix}-\mu _{A}&\mu _{AG}&\mu _{AC}&\mu _{AT}\\\mu _{GA}&-\mu _{G}&\mu _{GC}&\mu _{GT}\\\mu _{CA}&\mu _{CG}&-\mu _{C}&\mu _{CT}\\\mu _{TA}&\mu _{TG}&\mu _{TC}&-\mu _{T}\end{pmatrix}}}$

называется матрицей ставок . Обратите внимание, что по определению сумма записей в каждой строке равна нулю. Это следует из того${\displaystyle Q}$

${\displaystyle \mathbf {p} '(t)=\mathbf {p} (t)Q\,.}$

Для стационарного процесса , не зависящего от времени t , это дифференциальное уравнение может быть решено. Первый,${\displaystyle Q}$

${\displaystyle P(t)=\exp(tQ),}$

где обозначает экспоненту матрицы . Как результат,${\displaystyle \exp(tQ)}$${\displaystyle tQ}$

${\displaystyle \mathbf {p} (t)=\mathbf {p} (0)P(t)=\mathbf {p} (0)\exp(tQ)\,.}$

Эргодичность [ править ]

Если цепь Маркова неприводима , т. Е. Если всегда можно перейти из состояния в состояние (возможно, за несколько шагов), то она также эргодична . В результате оно имеет уникальное стационарное распределение , которое соответствует доле времени, проведенного в состоянии после того, как цепь Маркова проработала бесконечное количество времени. В эволюции ДНК, в предположении общего процесса для каждого сайта, стационарные частоты соответствуют равновесным базовым композициям. В самом деле, обратите внимание, что, поскольку стационарное распределение удовлетворяет , мы видим, что когда текущее распределение является стационарным, мы имеем${\displaystyle x}$${\displaystyle y}$ ${\displaystyle {\boldsymbol {\pi }}=\{\pi _{x},\,x\in {\mathcal {E}}\}}$${\displaystyle \pi _{x}}$${\displaystyle x}$${\displaystyle \pi _{A},\,\pi _{G},\,\pi _{C},\,\pi _{T}}$${\displaystyle {\boldsymbol {\pi }}}$${\displaystyle {\boldsymbol {\pi }}Q=0}$${\displaystyle \mathbf {p} (t)}$${\displaystyle {\boldsymbol {\pi }}}$

${\displaystyle {\mathbf {p} '(t)=\mathbf {p} (t)Q={\boldsymbol {\pi }}}Q=0\,.}$

Другими словами, частоты не меняются.${\displaystyle p_{A}(t),\,p_{G}(t),\,p_{C}(t),\,p_{T}(t)}$

Обратимость времени [ править ]

Определение : Стационарный марковский процесс обратим во времени, если (в установившемся состоянии) величина изменения от состояния до равна величине изменения от до (хотя эти два состояния могут возникать с разными частотами). Это означает, что:${\displaystyle x\ }$${\displaystyle y\ }$${\displaystyle y\ }$${\displaystyle x\ }$

${\displaystyle \pi _{x}\mu _{xy}=\pi _{y}\mu _{yx}\ }$

Не все стационарные процессы обратимы, однако наиболее часто используемые модели эволюции ДНК предполагают обратимость во времени, что считается разумным предположением.

В предположении обратимости времени let легко увидеть, что:${\displaystyle s_{xy}=\mu _{xy}/\pi _{y}\ }$

${\displaystyle s_{xy}=s_{yx}\ }$

Определение Симметричный член называется взаимозаменяемостью между состояниями и . Другими словами, это доля частоты состояния, которая является результатом переходов из состояния в состояние .${\displaystyle s_{xy}\ }$${\displaystyle x\ }$${\displaystyle y\ }$${\displaystyle s_{xy}\ }$${\displaystyle x\ }$${\displaystyle y\ }$${\displaystyle x\ }$

Следствие 12 недиагональных элементов матрицы скорости (обратите внимание, что недиагональные элементы определяют диагональные элементы, так как строки суммы равны нулю) могут быть полностью определены 9 числами; это: 6 членов взаимозаменяемости и 3 стационарные частоты (так как сумма стационарных частот равна 1).${\displaystyle Q\ }$${\displaystyle Q\ }$${\displaystyle \pi _{x}\ }$

Масштабирование длин ветвей [ править ]

Сравнивая существующие последовательности, можно определить степень расхождения последовательностей. Это необработанное измерение расхождения предоставляет информацию о количестве изменений, произошедших на пути, разделяющем последовательности. Простой подсчет различий (The расстояние Хэмминга ) между последовательностями часто недооценивают количество замещения из - за многократных ударов (см homoplasy). Пытаться оценить точное количество произошедших изменений сложно и обычно не требуется. Вместо этого длина ветвей (и длина пути) в филогенетическом анализе обычно выражается в ожидаемом количестве изменений на сайт. Длина пути - это произведение продолжительности пути во времени и средней скорости замен. Хотя их продукт можно оценить, скорость и время невозможно определить по расхождению последовательностей.

Описания матриц ставок на этой странице точно отражают относительную величину различных замен, но эти матрицы ставок не масштабируются таким образом, что длина ветви, равная 1, дает одно ожидаемое изменение. Это масштабирование может быть выполнено путем умножения каждого элемента матрицы на один и тот же коэффициент или просто путем масштабирования длин ветвей. Если мы используем β для обозначения коэффициента масштабирования и ν для обозначения длины ветви, измеренной в ожидаемом количестве замен на сайт, тогда βν используется в формулах вероятности перехода ниже вместо μ t . Обратите внимание, что ν - это параметр, который оценивается по данным, и называется длиной ветви, в то время как β - это просто число, которое можно вычислить из матрицы скорости (это не отдельный бесплатный параметр).

Значение β можно найти, установив ожидаемую скорость потока состояний равной 1. Диагональные элементы матрицы скорости ( Q- матрицы) представляют -1 кратную скорость выхода из каждого состояния. Для моделей с обратимым во времени мы знаем частоты равновесного состояния (это просто значение параметра π _i для состояния i ). Таким образом, мы можем найти ожидаемую скорость изменения, вычислив сумму потоков из каждого состояния, взвешенную по доле сайтов, которые, как ожидается, будут в этом классе. Установка β равной величине, обратной этой сумме, гарантирует, что масштабируемый процесс будет иметь ожидаемый поток 1:

${\displaystyle \beta =1/\left(-\sum _{i}\pi _{i}\mu _{ii}\right)}$

Например, в Jukes-Cantor коэффициент масштабирования будет 4 / (3μ), потому что скорость выхода из каждого состояния составляет 3μ / 4 .

Наиболее распространенные модели эволюции ДНК [ править ]

Модель JC69 (Джукс и Кантор, 1969) [ править ]

JC69, модель Джукса и Кантора 1969 года ^[1], является простейшей моделью замещения . Есть несколько предположений. Он предполагает одинаковые базовые частоты и одинаковые частоты мутаций . Таким образом, единственным параметром этой модели является общий коэффициент замещения. Как упоминалось ранее, эта переменная становится постоянной, когда мы нормализуем среднюю скорость до 1.${\displaystyle \left(\pi _{A}=\pi _{G}=\pi _{C}=\pi _{T}={1 \over 4}\right)}$${\displaystyle \mu }$

${\displaystyle 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}}}$

Вероятность перехода от начального состояния к конечному как функция длины ветви ( ) для JC69. Красная кривая: состояния нуклеотидов и разные. Синяя кривая: начальное и конечное состояния одинаковы. По прошествии длительного времени вероятности стремятся к частотам равновесия нуклеотидов (0,25: пунктирная линия).${\displaystyle P_{ij}}$${\displaystyle i}$${\displaystyle j}$${\displaystyle \nu }$${\displaystyle i}$${\displaystyle j}$

${\displaystyle P={\begin{pmatrix}{{1 \over 4}+{3 \over 4}e^{-t\mu }}&{{1 \over 4}-{1 \over 4}e^{-t\mu }}&{{1 \over 4}-{1 \over 4}e^{-t\mu }}&{{1 \over 4}-{1 \over 4}e^{-t\mu }}\\\\{{1 \over 4}-{1 \over 4}e^{-t\mu }}&{{1 \over 4}+{3 \over 4}e^{-t\mu }}&{{1 \over 4}-{1 \over 4}e^{-t\mu }}&{{1 \over 4}-{1 \over 4}e^{-t\mu }}\\\\{{1 \over 4}-{1 \over 4}e^{-t\mu }}&{{1 \over 4}-{1 \over 4}e^{-t\mu }}&{{1 \over 4}+{3 \over 4}e^{-t\mu }}&{{1 \over 4}-{1 \over 4}e^{-t\mu }}\\\\{{1 \over 4}-{1 \over 4}e^{-t\mu }}&{{1 \over 4}-{1 \over 4}e^{-t\mu }}&{{1 \over 4}-{1 \over 4}e^{-t\mu }}&{{1 \over 4}+{3 \over 4}e^{-t\mu }}\end{pmatrix}}}$

Когда длина ветви, измеряется в ожидаемом количестве изменений на сайт, тогда:${\displaystyle \nu }$

${\displaystyle P_{ij}(\nu )=\left\{{\begin{array}{cc}{1 \over 4}+{3 \over 4}e^{-4\nu /3}&{\mbox{ if }}i=j\\{1 \over 4}-{1 \over 4}e^{-4\nu /3}&{\mbox{ if }}i\neq j\end{array}}\right.}$

Стоит отметить, что то , что обозначает сумму любого столбца (или строки) матрицы, умноженную на время и, таким образом, означает ожидаемое количество замен во времени (длительность ветвления) для каждого конкретного сайта (на сайт), когда скорость замены равна .${\displaystyle \nu ={3 \over 4}t\mu =({\mu \over 4}+{\mu \over 4}+{\mu \over 4})t}$${\displaystyle Q}$${\displaystyle t}$${\displaystyle \mu }$

Учитывая долю сайтов, которые различаются между двумя последовательностями, оценка Джукса-Кантора эволюционного расстояния (с точки зрения ожидаемого числа изменений) между двумя последовательностями дается выражением ${\displaystyle p}$

${\displaystyle {\hat {d}}=-{3 \over 4}\ln({1-{4 \over 3}p})={\hat {\nu }}}$

В этой формуле часто называют -дистанцией. Это достаточный статистический показатель для расчета поправки на расстояние Джукса-Кантора, но его недостаточно для расчета эволюционного расстояния в более сложных моделях, которые следуют ниже (также обратите внимание, что используемые в последующих формулах не идентичны « -дистанции») .${\displaystyle p}$${\displaystyle p}$${\displaystyle p}$${\displaystyle p}$

Модель K80 (Kimura 1980) [ править ]

K80, модель Кимуры 1980 ^[2], часто называемая двухпараметрической моделью Кимуры (или моделью K2P ), различает переходы ( то есть от пурина к пурину или , то есть от пиримидина к пиримидину) и трансверсии (от пурина к пиримидин или наоборот). В первоначальном описании модели Кимурой α и β использовались для обозначения скоростей этих типов замен, но теперь более распространено устанавливать скорость трансверсий на 1 и использовать κ для обозначения отношения скорости перехода / трансверсии (как делается ниже). Модель K80 предполагает, что все базы одинаковы ( ).${\displaystyle A\leftrightarrow G}$${\displaystyle C\leftrightarrow T}$${\displaystyle \pi _{A}=\pi _{G}=\pi _{C}=\pi _{T}=0.25}$

Скорость матрица , столбцы , , и , соответственно.${\displaystyle Q={\begin{pmatrix}{*}&{\kappa }&{1}&{1}\\{\kappa }&{*}&{1}&{1}\\{1}&{1}&{*}&{\kappa }\\{1}&{1}&{\kappa }&{*}\end{pmatrix}}}$${\displaystyle A}$${\displaystyle G}$${\displaystyle C}$${\displaystyle T}$

Двухпараметрическое расстояние Кимуры определяется как:

${\displaystyle K=-{1 \over 2}\ln((1-2p-q){\sqrt {1-2q}})}$

где p - доля сайтов, показывающих переходные различия, а q - доля сайтов, показывающих трансверсионные различия.

Модель K81 (Kimura 1981) [ править ]

K81, то Кимура 1 981 модель, ^[3] часто называют три модели параметра Кимуров (модель K3P) или Кимура три типа замещения (K3ST) моделью, имеют различные тарифы для переходов и два различных типа трансверсий . Два типа трансверсии - это те, которые сохраняют слабые / сильные свойства нуклеотидов (т. Е. И , обозначаются символом ^[3] ), и те, которые сохраняют амино / кето-свойства нуклеотидов (то есть, и , обозначаются символом ^{[3 ]} ). Модель K81 предполагает, что все равновесные базовые частоты равны (т. Е. ).${\displaystyle A\leftrightarrow T}$${\displaystyle C\leftrightarrow G}$${\displaystyle \gamma }$ ${\displaystyle A\leftrightarrow C}$${\displaystyle G\leftrightarrow T}$${\displaystyle \beta }$ ${\displaystyle \pi _{A}=\pi _{G}=\pi _{C}=\pi _{T}=0.25}$

Скорость матрица , столбцы , , и , соответственно.${\displaystyle Q={\begin{pmatrix}{*}&{\alpha }&{\beta }&{\gamma }\\{\alpha }&{*}&{\gamma }&{\beta }\\{\beta }&{\gamma }&{*}&{\alpha }\\{\gamma }&{\beta }&{\alpha }&{*}\end{pmatrix}}}$${\displaystyle A}$${\displaystyle G}$${\displaystyle C}$${\displaystyle T}$

Модель K81 используется гораздо реже, чем модель K80 (K2P) для оценки расстояния, и она редко является наиболее подходящей моделью в филогенетике максимального правдоподобия. Несмотря на эти факты, модель K81 продолжала изучаться в контексте математической филогенетики. ^{[4] [5] [6]} Одним из важных свойств является способность выполнять преобразование Адамара, предполагая, что шаблоны сайтов были сгенерированы на дереве с нуклеотидами, эволюционирующими в рамках модели K81. ^{[7] [8] [9]}

При использовании в контексте филогенетики преобразование Адамара обеспечивает элегантное и полностью обратимое средство для вычисления ожидаемых частот паттернов сайтов с учетом набора длин ветвей (или наоборот). В отличие от многих вычислений максимального правдоподобия, относительные значения для , и могут варьироваться в зависимости от ветви, и преобразование Адамара может даже предоставить свидетельство того, что данные не соответствуют дереву. Преобразование Адамара также можно комбинировать с широким спектром методов для учета неоднородности скорости между сайтами ^[10], используя непрерывные распределения, а не дискретные приближения, обычно используемые в филогенетике максимального правдоподобия ^[11]${\displaystyle \alpha }$${\displaystyle \beta }$${\displaystyle \gamma }$(хотя нужно пожертвовать обратимостью преобразования Адамара, чтобы использовать определенные распределения неоднородности скорости между сайтами ^[10] ).

Модель F81 (Felsenstein 1981) [ править ]

F81, модель Фельзенштейна 1981 г. ^[12], является расширением модели JC69, в которой базовые частоты могут изменяться от 0,25 ( )${\displaystyle \pi _{A}\neq \pi _{G}\neq \pi _{C}\neq \pi _{T}}$

Матрица оценок:

${\displaystyle Q={\begin{pmatrix}{*}&{\pi _{G}}&{\pi _{C}}&{\pi _{T}}\\{\pi _{A}}&{*}&{\pi _{C}}&{\pi _{T}}\\{\pi _{A}}&{\pi _{G}}&{*}&{\pi _{T}}\\{\pi _{A}}&{\pi _{G}}&{\pi _{C}}&{*}\end{pmatrix}}}$

Когда длина ветви, ν, измеряется в ожидаемом количестве изменений на сайт, тогда:

${\displaystyle \beta =1/(1-\pi _{A}^{2}-\pi _{C}^{2}-\pi _{G}^{2}-\pi _{T}^{2})}$

${\displaystyle P_{ij}(\nu )=\left\{{\begin{array}{cc}e^{-\beta \nu }+\pi _{j}\left(1-e^{-\beta \nu }\right)&{\mbox{ if }}i=j\\\pi _{j}\left(1-e^{-\beta \nu }\right)&{\mbox{ if }}i\neq j\end{array}}\right.}$

Модель HKY85 (Hasegawa, Kishino и Yano 1985) [ править ]

HKY85, модель Hasegawa, Kishino и Yano 1985, ^[13] может рассматриваться как объединение расширений, сделанных в моделях Kimura80 и Felsenstein81. А именно, он различает скорость переходов и трансверсий (с использованием параметра κ) и допускает неравные базовые частоты ( ). [Фельзенштейн описал аналогичную (но не эквивалентную) модель в 1984 году, используя другую параметризацию; ^[14] эта последняя модель упоминается как модель F84. ^[15] ]${\displaystyle \pi _{A}\neq \pi _{G}\neq \pi _{C}\neq \pi _{T}}$

Матрица ставок ${\displaystyle Q={\begin{pmatrix}{*}&{\kappa \pi _{G}}&{\pi _{C}}&{\pi _{T}}\\{\kappa \pi _{A}}&{*}&{\pi _{C}}&{\pi _{T}}\\{\pi _{A}}&{\pi _{G}}&{*}&{\kappa \pi _{T}}\\{\pi _{A}}&{\pi _{G}}&{\kappa \pi _{C}}&{*}\end{pmatrix}}}$

Если мы выразим длину ветви ν через ожидаемое количество изменений на сайте, тогда:

${\displaystyle \beta ={\frac {1}{2(\pi _{A}+\pi _{G})(\pi _{C}+\pi _{T})+2\kappa [(\pi _{A}\pi _{G})+(\pi _{C}\pi _{T})]}}}$

${\displaystyle P_{AA}(\nu ,\kappa ,\pi )=\left[\pi _{A}\left(\pi _{A}+\pi _{G}+(\pi _{C}+\pi _{T})e^{-\beta \nu }\right)+\pi _{G}e^{-(1+(\pi _{A}+\pi _{G})(\kappa -1.0))\beta \nu }\right]/(\pi _{A}+\pi _{G})}$

${\displaystyle P_{AC}(\nu ,\kappa ,\pi )=\pi _{C}\left(1.0-e^{-\beta \nu }\right)}$

${\displaystyle P_{AG}(\nu ,\kappa ,\pi )=\left[\pi _{G}\left(\pi _{A}+\pi _{G}+(\pi _{C}+\pi _{T})e^{-\beta \nu }\right)-\pi _{G}e^{-(1+(\pi _{A}+\pi _{G})(\kappa -1.0))\beta \nu }\right]/\left(\pi _{A}+\pi _{G}\right)}$

${\displaystyle P_{AT}(\nu ,\kappa ,\pi )=\pi _{T}\left(1.0-e^{-\beta \nu }\right)}$

а формула для других комбинаций состояний может быть получена путем подстановки соответствующих базовых частот.

Модель Т92 (Тамура 1992) [ править ]

T92, модель Тамуры 1992 ^[16], представляет собой математический метод, разработанный для оценки количества нуклеотидных замен на сайт между двумя последовательностями ДНК путем расширения двухпараметрического метода Кимуры (1980) на случай, когда существует систематическая ошибка содержания G + C. . Этот метод будет полезен при наличии сильных смещений трансверсии перехода и содержания G + C, как в случае митохондриальной ДНК дрозофилы . ^[16]

T92 включает в себя один составной параметр базовой частоты (также отмечен )${\displaystyle \theta \in (0,1)}$${\displaystyle \pi _{GC}}$${\displaystyle =\pi _{G}+\pi _{C}=1-(\pi _{A}+\pi _{T})}$

Поскольку T92 перекликается с вторым правилом четности Чаргаффа - спаривающиеся нуклеотиды имеют одинаковую частоту на одной цепи ДНК, G и C с одной стороны, и A и T с другой стороны - из этого следует, что четыре базовые частоты могут быть выражены как функция${\displaystyle \pi _{GC}}$

${\displaystyle \pi _{G}=\pi _{C}={\pi _{GC} \over 2}}$ и ${\displaystyle \pi _{A}=\pi _{T}={(1-\pi _{GC}) \over 2}}$

Матрица ставок ${\displaystyle Q={\begin{pmatrix}{*}&{\kappa \pi _{GC}/2}&{\pi _{GC}/2}&{(1-\pi _{GC})/2}\\{\kappa (1-\pi _{GC})/2}&{*}&{\pi _{GC}/2}&{(1-\pi _{GC})/2}\\{(1-\pi _{GC})/2}&{\pi _{GC}/2}&{*}&{\kappa (1-\pi _{GC})/2}\\{(1-\pi _{GC})/2}&{\pi _{GC}/2}&{\kappa \pi _{GC}/2}&{*}\end{pmatrix}}}$

Эволюционное расстояние между двумя последовательностями ДНК согласно этой модели определяется выражением

${\displaystyle d=-h\ln(1-{p \over h}-q)-{1 \over 2}(1-h)\ln(1-2q)}$

где и - содержание G + C ( ).${\displaystyle h=2\theta (1-\theta )}$${\displaystyle \theta }$${\displaystyle \pi _{GC}=\pi _{G}+\pi _{C}}$

Модель TN93 (Тамура и Ней 1993) [ править ]

TN93, модель Тамуры и Нея 1993, ^[17] различает два разных типа перехода ; ie ( ) может иметь скорость, отличную от ( ). Предполагается, что все трансверсии происходят с одинаковой скоростью, но эта скорость может отличаться от обеих скоростей переходов.${\displaystyle A\leftrightarrow G}$${\displaystyle C\leftrightarrow T}$

TN93 также допускает неравные базовые частоты ( ).${\displaystyle \pi _{A}\neq \pi _{G}\neq \pi _{C}\neq \pi _{T}}$

Матрица ставок ${\displaystyle Q={\begin{pmatrix}{*}&{\kappa _{1}\pi _{G}}&{\pi _{C}}&{\pi _{T}}\\{\kappa _{1}\pi _{A}}&{*}&{\pi _{C}}&{\pi _{T}}\\{\pi _{A}}&{\pi _{G}}&{*}&{\kappa _{2}\pi _{T}}\\{\pi _{A}}&{\pi _{G}}&{\kappa _{2}\pi _{C}}&{*}\end{pmatrix}}}$

Модель GTR (Tavaré 1986) [ править ]

ОТО, Обобщенная обратимая во времени модель Таваре 1986, ^[18], является наиболее общей нейтральной, независимой, конечной моделью с обратимой во времени возможной. Впервые он был описан в общей форме Саймоном Таваре в 1986 году ^[18].

Параметры ОТО состоят из вектора равновесной базовой частоты , задающего частоту, с которой каждое основание встречается в каждом узле, и матрицы скоростей${\displaystyle \Pi =(\pi _{A},\pi _{G},\pi _{C},\pi _{T})}$

${\displaystyle Q={\begin{pmatrix}{-(\alpha \pi _{G}+\beta \pi _{C}+\gamma \pi _{T})}&{\alpha \pi _{G}}&{\beta \pi _{C}}&{\gamma \pi _{T}}\\{\alpha \pi _{A}}&{-(\alpha \pi _{A}+\delta \pi _{C}+\epsilon \pi _{T})}&{\delta \pi _{C}}&{\epsilon \pi _{T}}\\{\beta \pi _{A}}&{\delta \pi _{G}}&{-(\beta \pi _{A}+\delta \pi _{G}+\eta \pi _{T})}&{\eta \pi _{T}}\\{\gamma \pi _{A}}&{\epsilon \pi _{G}}&{\eta \pi _{C}}&{-(\gamma \pi _{A}+\epsilon \pi _{G}+\eta \pi _{C})}\end{pmatrix}}}$

Где

${\displaystyle {\begin{aligned}\alpha =r(A\rightarrow G)=r(G\rightarrow A)\\\beta =r(A\rightarrow C)=r(C\rightarrow A)\\\gamma =r(A\rightarrow T)=r(T\rightarrow A)\\\delta =r(G\rightarrow C)=r(C\rightarrow G)\\\epsilon =r(G\rightarrow T)=r(T\rightarrow G)\\\eta =r(C\rightarrow T)=r(T\rightarrow C)\end{aligned}}}$

параметры скорости перехода.

Следовательно, ОТО (для четырех символов, как это часто бывает в филогенетике) требует 6 параметров скорости замещения, а также 4 параметра равновесной базовой частоты. Однако обычно это устраняется до 9 параметров плюс общее количество замен в единицу времени. При измерении времени в подстановках ( = 1) остается всего 8 свободных параметров.${\displaystyle \mu }$${\displaystyle \mu }$

В общем, чтобы вычислить количество параметров, нужно подсчитать количество записей над диагональю в матрице, то есть для n значений признаков на сайт , а затем добавить n для равновесных базовых частот и вычесть 1, потому что это фиксировано. Один получает${\displaystyle {{n^{2}-n} \over 2}}$${\displaystyle \mu }$

${\displaystyle {{n^{2}-n} \over 2}+n-1={1 \over 2}n^{2}+{1 \over 2}n-1.}$

Например, для аминокислотной последовательности (существует 20 «стандартных» аминокислот, из которых состоят белки ), можно найти 209 параметров. Однако при изучении кодирующих областей генома чаще работают с моделью замены кодонов (кодон состоит из трех оснований и кодирует одну аминокислоту в белке). Есть кодоны, но скорости для переходов между кодонов , которые отличаются более чем одним основанием принимается равным нулю. Значит, есть параметры.${\displaystyle 4^{3}=64}$${\displaystyle {{20\times 19\times 3} \over 2}+64-1=633}$

См. Также [ править ]

• Молекулярная эволюция
• Молекулярные часы
• UPGMA

Ссылки [ править ]

Дальнейшее чтение [ править ]

Внешние ссылки [ править ]

• DAWG: DNA Assembly With Gaps - бесплатное программное обеспечение для моделирования эволюции последовательностей