Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://kodomo.cmm.msu.ru/~lu.andreeva/ukorenenie.html
Дата изменения: Tue Mar 16 00:37:47 2010 Дата индексирования: Tue Oct 2 04:32:28 2012 Кодировка: Windows-1251 |
В предыдущем задании (задания 1 и 2) были построены филогенетические деревья
выбранных семи бактерий. Невозможно укоренить дерево, построенное методом максимальной бережливости, потому что
он не предполагает различия в расстояниях между листьями (не предполагает молекулярных часов),
и fneighbor с алгоритмом UPGMA, так как оно уже укоренено.
В предыдущем задании дерево было построено алгоритмом Neighbor-joining:
+----------ENO_NEIMA +-1 ! +-------ENO_RALPJ ! ! +ENO_HAEIN ! +--3 ! ! +ENO_PASMU 2-----------4 ! ! +----ENO_VIBCH ! +-5 ! +--ENO_YERPS ! +---------------ENO_BRAJA
Полученную формулу:
(((ENO_HAEIN:0.02338,ENO_PASMU:0.02066):0.04994,(ENO_VIBCH:0.07559,
ENO_YERPS:0.04895):0.02194):0.18703,(ENO_BRAJA:0.27318,(ENO_NEIMA:0.18782,
ENO_RALPJ:0.1305):0.03317):0.01138);
откроем программой TreeView и получим изображение:
г====3:ENO HAEIN г=========11 Ѓ L====4:ENO PASMU г=========================================10 Ѓ Ѓ г================5:ENO VIBCH Ѓ L===12 =13 L==========6:ENO YERPS Ѓ Ѓ г=============================================================7:ENO BRAJA L==8 Ѓ г==========================================1:ENO NEIMA L======9 L=============================2:ENO RALPJУкоренение прошло в ветвь 2-4 (см. дерево, построенное Neighbor-joining). Полученное дерево имеет такие же нетривиальные ветви, как и правильное дерево.
Реконструируем дерево по семейству белков энолаз, в качестве аутгруппы используем ENO_BACSU (Bacillus subtilis).
Fprotpars выдала следующий результат:
+--ENO_RALPJ +--------------7 ! +--ENO_NEIMA ! ! +--ENO_YERPS +--6 +-----5 ! ! ! +--ENO_VIBCH ! ! +--4 ! ! ! ! +--ENO_PASMU 1 +-----2 +-----3 ! ! +--ENO_HAEIN ! ! ! +-----------ENO_BACSU ! +--------------------ENO_BRAJAПосле обработки и укоренения дерева программой retree и выделения аутгруппы получили следующее:
г>>>>>>>>>>>>>>>>>>>>7:ENO BACSU Ѓ Ѓ г>>>>>8:ENO BRAJA ==9 г>>>>>>>>>>10 Ѓ Ѓ Ѓ г>>1:ENO RALPJ Ѓ Ѓ L>11 Ѓ Ѓ L>>2:ENO NEIMA L>12 Ѓ г>>3:ENO YERPS Ѓ г>>>>14 Ѓ Ѓ L>>4:ENO VIBCH L>>>>>>>13 Ѓ г>>5:ENO PASMU L>>>>15 L>>6:ENO HAEIN
С помощью программы fseqboot было создано 100 бутстрэп-реплик выравнивания белков выбранных бактерий,
с помощью fprotpars по полученным репликам были составлены деревья, из которых по принципу "расширенного большинства"
было создано единое дерево, совпавшее с правильным:
+------ENO YERPS +-54.2-| | +------ENO VIBCH +-100.0-| | | +------ENO HAEIN +-60.6-| +-100.0-| | | +------ENO PASMU +------| | | | +--------------------ENO BRAJA | | | +---------------------------ENO NEIMA | +----------------------------------ENO RALPJ4 ветви не получили поддержку: