Занятие 3.

Укоренение в среднюю точку

Укорените в среднюю точку дерево, построенное при выполнении задания 7 предыдущего занятия методом neighbor-joining.

Почему это нельзя сделать с деревьями, построенными методом максимальной экономии?

Почему это не имеет смысла делать с деревом, построенным методом UPGMA?

Воспользуйтесь программой retree пакета PHYLIP. Для этого:

• Скопируйте файл retree.exe с диска P в свою рабочую директорию.
• Скопируйте файл со скобочной формулой, выданный программой neighbor, в файл с именем intree (именно так, без расширения!)
• Запустите retree (из командной строки Windows; например можно в Far manager навести курсор на exe-файл и нажать <Enter>)
• Постарайтесь сами разобраться, как работать с программой retree. Для этого внимательно читайте то, что появляется в окне. Нужное вам действие — укоренение в среднюю точку ("Midpoint root the tree").
Вставьте в отчёт изображение укоренённого дерева и опишите его (в какую ветвь произошло укоренение, можно ли это укоренение счиать правильным, может быть, ещё какие-то наблюдения).

Использование внешней группы

Реконструируйте программой fprotpars укоренённое дерево отобранных вами бактерий, используя то же семейство белков, что и в предыдущем задании, а в качестве внешней группы — белок того же семейства из кишечной палочки (Escherichia coli, ECOLI).

Подсказка: необходимо добавить к файлу с невыровненными последовательностями белков фирмикут последовательность белка из кишечной палочки. Это удобнее всего сделать на kodomo командой

seqret sw:xxxx_ecoli stdout >> sequences.fasta

(вместо "xxxx" надо подставить мнемонику функции белков, по которым вы строили деревья, а вместо "sequences.fasta" — имя файла с невыровненными последовательностями, заранее скопированного в рабочую директорию).

Затем надо построить выравнивание всех последовательностей, отредактировать имена (оставив только мнемонику видов) и результат подать на вход программе fprotpars. Полученное дерево лучше обработать программой retree, указав в качестве действия "select an Outgroup", а в качестве номера — тот, что программа retree присвоит листу ECOLI (смотрите внимательно на монитор!).

Вставьте в отчёт изображение (можно вместе с листом ECOLI; но можно, прежде чем получать изображение дерева, отредактировать скобочную формулу вручную, убрав внешнюю группу). Опишите результат (по образцу предыдущего пункта).

Бутстрэп

Проведите бутстрэп-анализ филогении своих белков, используя программу fprotpars. Этапы работы:
• Создайте 100 бутстрэп-реплик выравнивания белков программой fseqboot.
• Создайте по полученным репликам деревья программой fprotpars (для этого просто подайте выходной файл программы fseqboot на вход программе fprotpars).
• Создайте из полученных деревьев единое дерево по принципу "расширенного большинства" (extended majority rule tree). Для этого файл с деревьями, выданный программой fprotpars, подайте на вход программе fconsense.
Приведите в отчёте (как преформатированный текст) дерево из выходного файла программы fconsense (не скобочную формулу, а дерево "нарисованное" ASCII-символами, с числами, означающими поддержку ветвей). Опишите свои наблюдения (в частности, улучшилась ли реконструкция филогении по сравнению с результатом fprotpars на исходном выравнивании?) Проанализируйте в выходном файле ветви, не получившие большинства (описаны в виде последовательностей точек и звёздочек). Есть ли среди них верные и какова их поддержка?