Эффект притяжения длинных ветвей – серьезная помеха для систематики. Он проявляется везде, где много конвергенций и реверсий. К сожалению, эту проблему нельзя решить простым увеличением объема рассматриваемого текста. Наоборот, чем больше текст, тем выше вероятность обнаружения случайных совпадений. Про такие деревья говорят, что они лежат в “зоне Фельзенстайна” (звучит устрашающе!), названной в честь американского биолога Джо Фельзенстайна. Увы, ДНК особенно подвержена эффекту притяжения длинных ветвей. Основная причина в том, что в ДНК всего четыре “буквы”. Поскольку большинство изменений затрагивают всего одну “букву”, случайные мутации с высокой вероятностью могут привести к совпадениям. Так возникает притяжение длинных ветвей. Очевидно, что для таких случаев нужна альтернатива методу парсимонии. Она существует – это метод правдоподобия. В последнее время он используется все чаще.
Оценка правдоподобия требует больше вычислительных мощностей, чем метод парсимонии, поскольку здесь мы учитываем длину ветвей. Таким образом, приходится иметь дело с еще большим количеством деревьев: вдобавок к рассмотрению возможных схем ветвления мы должны учитывать возможные длины ветвей. Геркулесов труд! Поэтому, несмотря на упрощенные методы вычисления, компьютеры пока могут подвергнуть анализу небольшое количество видов.
Термин “правдоподобие” здесь имеет вполне точное значение. Возьмем древо определенной формы (с учетом длины ветвей). Из всех возможных эволюционных траекторий, посредством которых может сформироваться филогенетическое древо данной формы, всего несколько могут привести к тому тексту, который мы сейчас видим. "Правдоподобие” данного древа – это ничтожно малая вероятность получения реально существующих текстов, а не каких-нибудь текстов, которые могут появиться на таком древе. Величина правдоподобия для древа очень мала, однако это не мешает сравнить одну малую величину с другой, чтобы выбрать нужную.
Неукорененное филогенетическое древо первых 250 строк 24 списков «Кентерберийских рассказов». Здесь представлен набор списков, изученный в рамках проекта «Кентерберийские рассказы». Сокращения соответствуют тем, что использованы в проекте. Схема построена методом парсимонии, на каждой ветви указаны индексы бутстреп-поддержки. Для четырех списков, которые обсуждаются нами, указаны их полные названия.
Выбирать "лучшее” древо методом правдоподобия можно по-разному. Самый простой способ – искать наиболее правдоподобное древо. Это метод максимального правдоподобия. Однако то, что это наиболее правдоподобное древо, вовсе не означает, что другие деревья не окажутся почти столь же правдоподобными. Совсем недавно было предложено не искать одно самое правдоподобное древо, а рассматривать все возможные. При этом степень "доверия” к древу должна зависеть от его правдоподобия. Этот подход представляет собой альтернативу методу правдоподобия и известен как байесовский метод. Если схема ветвления подтверждается большим количеством правдоподобных деревьев, мы заключаем, что эта схема с высокой вероятностью верна. Конечно, как и в методе максимального правдоподобия, мы не можем проверить все деревья. Но существуют способы упрощения вычислений, и они довольно неплохо работают.
Степень нашего доверия древу, которое мы в итоге выберем, зависит от того, насколько мы уверены в правильности каждого разветвления. Поэтому возле точек ветвления часто указывают степень “уверенности” в них. При использовании байесовского метода правдоподобие точек ветвления вычисляется автоматически, однако для других методов, таких как парсимония или максимальное правдоподобие, необходимы альтернативные способы подсчета. Чаще всего используют метод бутстрепа: многократно обсчитываются выборки данных, и оценки сравниваются с результатами для всего древа. Так мы можем понять, насколько древо устойчиво к ошибкам. Чем больше индекс бутстреп-поддержки, тем надежнее точка ветвления. Правда, точно интерпретировать полученные индексы бывает непросто. По сходному алгоритму работают методы “складного ножа” (jackknife) и “поддержки Бремера”. Все они служат для оценки достоверности точек ветвления.
Прежде чем оставить литературу, рассмотрим итоговое древо, построенное для первых 250 строк в 24 манускриптах Чосера. Это филограмма, на которой информативна не только схема ветвления, но и длина ветвей. На схеме видно, какие списки почти идентичны, а какие сильно отличаются от остальных. Эта филограмма неукорененная, то есть не указывает на то, какой из 24 манускриптов ближе всех к “оригиналу”.
Вернемся к гиббонам. Принцип парсимонии предполагает существование четырех групп. Ниже приведена укорененная диаграмма, основанная на морфологических признаках. Здесь виды рода Hylobates (настоящие гиббоны) группируются вместе, как и виды рода Nomascus. Обе группы поддерживаются высокими индексами бутстреп-поддержки (указаны над ветвями). Однако в нескольких местах порядок ветвлений не определен. Хотя Hylobates и Bunopithecus вроде бы формируют группу, индекс бутстреп-поддержки (63) представляется неубедительным для тех, кто умеет читать подобные руны. Морфологических признаков для построения древа недостаточно.
Укорененная кладограмма гиббонов, построенная на основе морфологии. Geissmann [100].
По этой причине Кристиан Роос и Томас Гайсман, ученые из Германии, обратились к молекулярной генетике, а именно к участку митохондриальной ДНК, который называют контрольным регионом. Взяв ДНК шести гиббонов, они расшифровали последовательности, выровняли их и провели анализ с помощью методов поиска ближайшего соседа, парсимонии и максимального правдоподобия. Самый убедительный результат был получен с помощью метода максимального правдоподобия, который лучше других методов справляется с эффектом притяжения длинных ветвей. Итоговое древо, где показаны отношения между четырьмя группами, приведено здесь. Значения бутстреп-поддержки на этом древе вполне убедительны. Так что, на мой взгляд, это то, что нам нужно.
Кладограмма гиббонов, построенная с помощью метода максимального правдоподобия по данным ДНК. Roos and Geissmann [246].
Видообразование у гиббонов произошло сравнительно недавно. Однако если изучать все более удаленные виды, которые будут разделены все более длинными ветвями, в конце концов даже изощренные методы Байеса и максимального правдоподобия откажутся нам служить. В определенный момент недопустимо большая доля сходств окажется случайной. Когда это происходит с ДНК, говорят, что наступило насыщение. И тогда ни один метод не поможет реконструировать схему родственных отношений: действие времени заглушает “филогенетический” сигнал. Особенно остро этот вопрос стоит в отношении нейтральных мутаций ДНК. Давление естественного отбора не позволяет генам сбиваться с пути, удерживая их в узком диапазоне. В некоторых случаях самые важные функциональные гены могут оставаться практически неизменными сотни миллионов лет. Однако для псевдогена, с которого никогда ничего не считывается, таких промежутков времени достаточно для безнадежно сильного насыщения. В таких случаях нам приходится искать другие данные. Одна из самых перспективных идей – использование редких геномных мутаций, о которых я упоминал. Эти изменения затрагивают значительные участки ДНК, а не одну “букву”. Поскольку такие перестройки редки и, как правило, уникальны, проблема случайного сходства не возникает. Эти мутации могут выявлять неожиданные родственные связи. Мы убедимся в этом, когда к толпе пилигримов присоединятся гиппопотамы. (Вот увидите, они расскажут удивительные вещи!)