Практическое занятие Нахождение гомологов, выравнивание последовательностей, построение и статистическая оценка филогенетического дерева
1. Отобрать 15 -20 гомологичных последовательностей • По выбору, например: 16 S RNA 50 S RNA t. RNAs (для каждой аминокислоты своя t. RNA) ATP synthase (alpha, beta, gamma, delta, epsilon subintis) t. RNA synthases ((для каждой аминокислоты своя) DNA polymerase (alpha and beta subunits) RNA polymerase Ribosomal proteins (их ~50, так что следите за номером, буква L – large или S - small субъединицы рибосомы) – Любой ген, который вы хотите проанализировать – – – –
1. Отобрать 15 -20 гомологичных последовательностей • • Идете на сайт http: //www. ncbi. nlm. nih. gov/ Среди баз данных выбираете либо protein, либо nucleotide. Например:
Измените формат на FASTA Скопируйте в текстовый файл
Идем в BLAST
В зависимости от того, какие вы гены выбрали, кодирующие или РНК
Вставпоследова Удобно, чтобы результаты высвечивались в новом окне
Или даже 5000, 10000 Можно выбрать любую матрицу, запомните какую.
Прокрутите результаты вниз, до начала выравниваний. Начинайте отбирать последовательности снизу вверх, чтобы они были <=80% идентичности Щелкните на последовательность, это вас перенесет на локальное выравнивание.
Соберите 15 -20 последовательностей с условием, что их % identity < 80 Сохраните все по в один текстовый файл
2. Выровнять отобранные последовательности
Скачайте программу множественного выравнивания
Скачиваем Clustalx (графический вариант) – есть дистрибутивы для всех систем
Ваше выравнивание не должно быть таким хорошим, так как вы должны были выбрать более разошедшиеся последовательности.
2. Построить филогенетическое дерево методом ближайшего соседа (NJ) и выполнить bootstrap-анализ
Все в той же программе Clustal. X Формат phylip читается почти всеми программами отображения деревьев (Tree. Viewers) Change Label To Node here
Скачиваем какой-нибудь Tree Viewer
Tree. View под Windows
Открываете phb (bootstrapped) файл Fig. Tree или Tree. Viewer или в любой другой программе визуализации деревьев В Fig. Tree – появится сообщение Нажмите Ok Чтобы увидеть бутстрэп значения, надо в Branch Labels в разделе Display выбрать label
Bootstrap-значения должны быть в диапазоне от 1 до 1000 (мы так задали). Могли бы задать от 1 до 100.
Отчетность Распечатать или послать по email: 1. Текстовый файл с отобранными последовательностями. Указать использованные параметры BLAST (какая матрица была использована, сколько последовательностей было выбрано для поиска) 2. Скриншот выравнивания 3. Скриншот дерева, на котором видны бутстрэп-значения