stringtranslate.com

Вычислительная филогенетика

Вычислительная филогенетика , вывод филогении или филогенетический вывод фокусируется на вычислительных и оптимизационных алгоритмах , эвристике и подходах, используемых в филогенетическом анализе. Цель состоит в том, чтобы найти филогенетическое дерево , представляющее оптимальное эволюционное происхождение между набором генов , видов или таксонов . Максимальное правдоподобие , экономия , байесовский и минимальная эволюция являются типичными критериями оптимальности, используемыми для оценки того, насколько хорошо топология филогенетического дерева описывает данные последовательности. [1] [2] Обмен ближайшими соседями (NNI), обрезка и пересадка поддерева (SPR) и бисекция и повторное соединение дерева (TBR), известные как перестановки деревьев , являются детерминированными алгоритмами для поиска оптимального или наилучшего филогенетического дерева. Пространство и ландшафт поиска оптимального филогенетического дерева известны как пространство поиска филогении.

Критерий оптимальности максимального правдоподобия (также вероятности) — это процесс нахождения топологии дерева вместе с длинами его ветвей, которые обеспечивают наивысшую вероятность наблюдения данных последовательности, в то время как критерий оптимальности экономии — это наименьшее количество изменений состояний и эволюции, необходимых для филогенетического дерева для объяснения данных последовательности. [1] [2]

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

Многие формы молекулярной филогенетики тесно связаны с выравниванием последовательностей и широко используют его при построении и уточнении филогенетических деревьев, которые используются для классификации эволюционных связей между гомологичными генами, представленными в геномах расходящихся видов. Филогенетические деревья, построенные с помощью вычислительных методов, вряд ли будут идеально воспроизводить эволюционное дерево , представляющее исторические связи между анализируемыми видами. [ необходима цитата ] Историческое дерево видов также может отличаться от исторического дерева отдельного гомологичного гена, общего для этих видов.

Типы филогенетических деревьев и сетей

Филогенетические деревья, созданные вычислительной филогенетикой, могут быть как укорененными , так и неукорененными в зависимости от входных данных и используемого алгоритма. Укорененное дерево — это направленный граф , который явно идентифицирует последнего общего предка (MRCA), [ требуется ссылка ] обычно это введенная последовательность, которая не представлена ​​во входных данных. Меры генетического расстояния могут использоваться для построения дерева с входными последовательностями в качестве конечных узлов и их расстояниями от корня, пропорциональными их генетическому расстоянию от предполагаемого MRCA. Идентификация корня обычно требует включения во входные данные по крайней мере одной «внешней группы», о которой известно, что она лишь отдаленно связана с интересующими последовательностями.

Напротив, некорневые деревья отображают расстояния и отношения между входными последовательностями без предположений относительно их происхождения. Некорневое дерево всегда можно получить из корневого дерева, но корень обычно нельзя поместить на некорневое дерево без дополнительных данных о скоростях расхождения, таких как предположение о гипотезе молекулярных часов . [3]

Набор всех возможных филогенетических деревьев для заданной группы входных последовательностей можно концептуализировать как дискретно определенное многомерное «пространство деревьев», через которое можно проследить пути поиска с помощью алгоритмов оптимизации . Хотя подсчет общего числа деревьев для нетривиального числа входных последовательностей может быть затруднен из-за вариаций в определении топологии дерева, всегда верно, что для заданного числа входных последовательностей и выбора параметров существует больше корневых деревьев, чем некорневых. [2]

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

Кодирование символов и определение гомологии

Морфологический анализ

Основная проблема в морфологической филогенетике — это сборка матрицы, представляющей отображение каждого из сравниваемых таксонов с репрезентативными измерениями для каждой из фенотипических характеристик, используемых в качестве классификатора. Типы фенотипических данных, используемых для построения этой матрицы, зависят от сравниваемых таксонов; для отдельных видов они могут включать измерения среднего размера тела, длины или размеров отдельных костей или других физических особенностей, или даже поведенческих проявлений. Конечно, поскольку не все возможные фенотипические характеристики могут быть измерены и закодированы для анализа, выбор характеристик для измерения является основным неотъемлемым препятствием для метода. Решение о том, какие признаки использовать в качестве основы для матрицы, обязательно представляет собой гипотезу о том, какие признаки вида или более высокого таксона являются эволюционно значимыми. [4] Морфологические исследования могут быть озадачены примерами конвергентной эволюции фенотипов. [5] Основной проблемой при построении полезных классов является высокая вероятность межтаксонного перекрытия в распределении вариации фенотипа. Включение вымерших таксонов в морфологический анализ часто затруднено из-за отсутствия или неполноты ископаемых записей, но было показано, что это оказывает значительное влияние на полученные деревья; в одном исследовании только включение вымерших видов обезьян дало морфологически полученное дерево, которое соответствовало дереву, полученному на основе молекулярных данных. [6]

Некоторые фенотипические классификации, особенно те, которые используются при анализе очень разнообразных групп таксонов, являются дискретными и недвусмысленными; например, классификация организмов как обладающих или не обладающих хвостом, в большинстве случаев проста, как и подсчет таких признаков, как глаза или позвонки. Однако наиболее подходящее представление непрерывно изменяющихся фенотипических измерений является спорной проблемой без общего решения. Распространенный метод заключается в том, чтобы просто отсортировать интересующие измерения по двум или более классам, делая непрерывно наблюдаемые вариации дискретно классифицируемыми (например, все примеры с плечевыми костями длиннее заданного порогового значения, оцениваются как члены одного состояния, а все члены, плечевые кости которых короче порогового значения, оцениваются как члены второго состояния). Это приводит к легко манипулируемому набору данных , но подвергается критике за плохое представление основы для определений классов и за принесение в жертву информации по сравнению с методами, которые используют непрерывное взвешенное распределение измерений. [7]

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

Молекулярный анализ

Проблема кодирования символов сильно отличается в молекулярном анализе, поскольку символы в данных биологических последовательностей определяются немедленно и дискретно - отдельные нуклеотиды в последовательностях ДНК или РНК и отдельные аминокислоты в последовательностях белков . Однако определение гомологии может быть сложным из-за присущих трудностей множественного выравнивания последовательностей . Для заданного MSA с зазорами можно построить несколько укорененных филогенетических деревьев, которые различаются по своей интерпретации того, какие изменения являются « мутациями » по сравнению с предковыми признаками, и какие события являются мутациями вставки или мутациями делеции . Например, имея только попарное выравнивание с областью разрыва, невозможно определить, несет ли одна последовательность мутацию вставки или другая несет делецию. Проблема усугубляется в MSA с невыровненными и неперекрывающимися зазорами. На практике значительные области вычисленного выравнивания могут быть проигнорированы при построении филогенетического дерева, чтобы избежать интеграции шумных данных в расчет дерева.

Методы дистанционных матриц

Методы филогенетического анализа на основе матрицы расстояний явно полагаются на меру «генетического расстояния» между классифицируемыми последовательностями, и поэтому они требуют MSA в качестве входных данных. Расстояние часто определяется как доля несовпадений в выровненных позициях, при этом пробелы либо игнорируются, либо считаются несовпадениями. [3] Методы расстояний пытаются построить матрицу «все-все» из набора запросов последовательностей, описывающую расстояние между каждой парой последовательностей. Из этого строится филогенетическое дерево, которое помещает близкородственные последовательности под один и тот же внутренний узел и чьи длины ветвей близко воспроизводят наблюдаемые расстояния между последовательностями. Методы матриц расстояний могут создавать как корневые, так и некорневые деревья, в зависимости от алгоритма, используемого для их вычисления. Они часто используются в качестве основы для прогрессивных и итеративных типов множественных выравниваний последовательностей . Основным недостатком методов матриц расстояний является их неспособность эффективно использовать информацию о локальных регионах с высокой вариабельностью, которые появляются в нескольких поддеревьях. [2]

UPGMA и WPGMA

Методы UPGMA ( метод невзвешенных парных групп со средним арифметическим ) и WPGMA ( метод взвешенных парных групп со средним арифметическим ) создают корневые деревья и требуют предположения о постоянной скорости, то есть они предполагают ультраметрическое дерево, в котором расстояния от корня до каждой верхушки ветви равны. [9]

Присоединение соседей

Методы объединения соседей применяют общие методы кластерного анализа к анализу последовательностей, используя генетическое расстояние в качестве метрики кластеризации. Простой метод объединения соседей создает некорневые деревья, но он не предполагает постоянную скорость эволюции (т. е. молекулярные часы ) по всем линиям. [10]

Метод Фитча–Марголиаша

Метод Фитча-Марголиаша использует метод взвешенных наименьших квадратов для кластеризации на основе генетического расстояния. [11] Близкородственным последовательностям придается больший вес в процессе построения дерева для коррекции возросшей неточности измерения расстояний между отдаленно родственными последовательностями. Расстояния, используемые в качестве входных данных для алгоритма, должны быть нормализованы для предотвращения больших артефактов при вычислении взаимосвязей между близкородственными и отдаленно родственными группами. Расстояния, вычисляемые этим методом, должны быть линейными ; критерий линейности для расстояний требует, чтобы ожидаемые значения длин ветвей для двух отдельных ветвей были равны ожидаемому значению суммы двух расстояний ветвей — свойство, которое применяется к биологическим последовательностям только тогда, когда они были скорректированы на возможность обратных мутаций на отдельных участках. Эта коррекция выполняется с помощью матрицы подстановки, такой как полученная из модели эволюции ДНК Джукса-Кантора . Коррекция расстояния необходима на практике только тогда, когда скорости эволюции между ветвями различаются. [2] Другая модификация алгоритма может быть полезна, особенно в случае концентрированных расстояний (см. явление концентрации меры и проклятие размерности ): было показано, что эта модификация, описанная в [12], повышает эффективность алгоритма и его надежность.

Критерий наименьших квадратов, применяемый к этим расстояниям, более точен, но менее эффективен, чем методы объединения соседей. Дополнительное улучшение, которое корректирует корреляции между расстояниями, возникающими из многих тесно связанных последовательностей в наборе данных, также может быть применено при увеличении вычислительных затрат. Нахождение оптимального дерева наименьших квадратов с любым поправочным коэффициентом является NP-полным [ 13], поэтому эвристические методы поиска, подобные тем, которые используются в анализе максимальной экономии, применяются к поиску в пространстве деревьев.

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

Независимая информация о взаимосвязи между последовательностями или группами может быть использована для сокращения пространства поиска дерева и корневых некорневых деревьев. Стандартное использование методов матрицы расстояний подразумевает включение по крайней мере одной последовательности внешней группы , которая, как известно, имеет лишь отдаленное родство с интересующими последовательностями в наборе запросов. [3] Такое использование можно рассматривать как тип экспериментального контроля . Если внешняя группа была выбрана правильно, она будет иметь гораздо большее генетическое расстояние и, следовательно, большую длину ветви, чем любая другая последовательность, и она будет отображаться вблизи корня укорененного дерева. Выбор подходящей внешней группы требует выбора последовательности, которая умеренно связана с интересующими последовательностями; слишком близкое родство сводит на нет цель внешней группы, а слишком отдаленное добавляет шум в анализ. [3] Также следует проявлять осторожность, чтобы избегать ситуаций, в которых виды, из которых были взяты последовательности, имеют отдаленное родство, но ген, кодируемый последовательностями, высококонсервативен между линиями. Горизонтальный перенос генов , особенно между в противном случае расходящимися бактериями , также может затруднить использование внешней группы.

Максимальная бережливость

Максимальная экономия (MP) — это метод определения потенциального филогенетического дерева, которое требует наименьшего общего числа эволюционных событий для объяснения наблюдаемых данных о последовательностях. Некоторые способы оценки деревьев также включают «стоимость», связанную с определенными типами эволюционных событий, и пытаются найти дерево с наименьшей общей стоимостью. Это полезный подход в случаях, когда не все возможные типы событий одинаково вероятны — например, когда известно, что определенные нуклеотиды или аминокислоты более изменчивы, чем другие.

Самый наивный способ определения наиболее экономного дерева — это простое перечисление — рассмотрение каждого возможного дерева по порядку и поиск дерева с наименьшим счетом. Однако это возможно только для относительно небольшого числа последовательностей или видов, поскольку известно, что проблема определения наиболее экономного дерева является NP-трудной ; [2] следовательно, был разработан ряд эвристических методов поиска для оптимизации, чтобы найти высокоэкономное дерево, если не лучшее в наборе. Большинство таких методов включают механизм минимизации в стиле наискорейшего спуска , работающий на критерии перестановки деревьев .

Ветви и границы

Алгоритм ветвей и границ — это общий метод, используемый для повышения эффективности поиска почти оптимальных решений NP-трудных задач, впервые примененный в филогенетике в начале 1980-х годов. [14] Ветви и границы особенно хорошо подходят для построения филогенетического дерева, поскольку он по своей сути требует разделения проблемы на древовидную структуру , поскольку подразделяет проблемное пространство на более мелкие области. Как следует из его названия, он требует в качестве входных данных как правило ветвления (в случае филогенетики — добавление следующего вида или последовательности к дереву), так и границу (правило, которое исключает определенные области пространства поиска из рассмотрения, тем самым предполагая, что оптимальное решение не может занимать эту область). Определение хорошей границы — самый сложный аспект применения алгоритма в филогенетике. Простой способ определения границы — это максимальное количество предполагаемых эволюционных изменений, разрешенных для одного дерева. Набор критериев, известный как правила Жарких [15], серьезно ограничивает пространство поиска, определяя характеристики, общие для всех кандидатов на «наиболее экономные» деревья. Два самых основных правила требуют устранения всех, кроме одной избыточной последовательности (для случаев, когда множественные наблюдения дали идентичные данные) и устранения характерных участков, в которых два или более состояний не встречаются по крайней мере у двух видов. В идеальных условиях эти правила и связанный с ними алгоритм полностью определяли бы дерево.

Алгоритм Санкоффа-Мореля-Седергрена

Алгоритм Санкоффа-Мореля-Седергрена был одним из первых опубликованных методов для одновременного создания MSA и филогенетического дерева для последовательностей нуклеотидов. [16] Метод использует расчет максимальной экономии в сочетании с функцией оценки, которая штрафует пробелы и несовпадения, тем самым отдавая предпочтение дереву, которое вводит минимальное количество таких событий (альтернативная точка зрения гласит, что деревья, которым следует отдавать предпочтение, - это те, которые максимизируют количество сходства последовательностей, которое можно интерпретировать как гомологию, точка зрения, которая может привести к различным оптимальным деревьям [17] ). Вмененные последовательности во внутренних узлах дерева оцениваются и суммируются по всем узлам в каждом возможном дереве. Сумма дерева с наименьшей оценкой обеспечивает как оптимальное дерево, так и оптимальное MSA с учетом функции оценки. Поскольку метод является очень вычислительно интенсивным, приближенный метод, в котором начальные предположения для внутренних выравниваний уточняются по одному узлу за раз. Как полная, так и приближенная версия на практике вычисляются с помощью динамического программирования. [2]

MALIGN и POY

Более поздние методы филогенетического дерева/MSA используют эвристики для выделения деревьев с высокими оценками, но не обязательно оптимальных. Метод MALIGN использует технику максимальной экономии для вычисления множественного выравнивания путем максимизации оценки кладограммы , а его компаньон POY использует итеративный метод, который объединяет оптимизацию филогенетического дерева с улучшениями в соответствующей MSA. [18] Однако использование этих методов при построении эволюционных гипотез подвергалось критике как предвзятое из-за преднамеренного построения деревьев, отражающих минимальные эволюционные события. [19] Это, в свою очередь, было опровергнуто мнением, что такие методы следует рассматривать как эвристические подходы для поиска деревьев, которые максимизируют количество сходства последовательностей, которое можно интерпретировать как гомологию. [17] [20]

Максимальная вероятность

Метод максимального правдоподобия использует стандартные статистические методы для вывода распределений вероятностей , чтобы назначить вероятности конкретным возможным филогенетическим деревьям. Метод требует модели замещения для оценки вероятности конкретных мутаций ; грубо говоря, дерево, которому требуется больше мутаций во внутренних узлах для объяснения наблюдаемой филогении, будет оцениваться как имеющее более низкую вероятность. Это в целом похоже на метод максимальной экономии, но максимальное правдоподобие обеспечивает дополнительную статистическую гибкость, допуская различные скорости эволюции как по линиям, так и по участкам. Фактически, метод требует, чтобы эволюция на разных участках и вдоль разных линий была статистически независимой . Таким образом, максимальное правдоподобие хорошо подходит для анализа отдаленно связанных последовательностей, но считается, что его вычислительно трудно вычислить из-за его NP-трудности. [21]

Алгоритм «отсечения», вариант динамического программирования , часто используется для сокращения пространства поиска путем эффективного вычисления вероятности поддеревьев. [2] Метод вычисляет вероятность для каждого сайта «линейным» образом, начиная с узла, единственными потомками которого являются листья (то есть кончики дерева), и работая в обратном направлении к «нижнему» узлу во вложенных наборах. Однако деревья, полученные с помощью метода, укореняются только в том случае, если модель подстановки необратима, что обычно не относится к биологическим системам. Поиск дерева максимального правдоподобия также включает компонент оптимизации длины ветвей, который трудно улучшить алгоритмически; часто используются общие глобальные инструменты оптимизации, такие как метод Ньютона–Рафсона .

Некоторые инструменты, которые используют максимальное правдоподобие для выведения филогенетических деревьев из данных о частоте аллелей вариантов (VAF), включают AncesTree и CITUP. [22] [23]

Байесовский вывод

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

Реализации байесовских методов обычно используют алгоритмы выборки Монте-Карло с цепями Маркова , хотя выбор набора перемещений варьируется; выборки, используемые в байесовской филогенетике, включают циклическую перестановку листовых узлов предлагаемого дерева на каждом шаге [24] и обмен поддеревьями потомков случайного внутреннего узла между двумя связанными деревьями. [25] Использование байесовских методов в филогенетике было спорным, в основном из-за неполной спецификации выбора набора перемещений, критерия приемлемости и предварительного распределения в опубликованных работах. [2] Байесовские методы, как правило, считаются превосходящими методы, основанные на экономии; они могут быть более склонны к притяжению длинных ветвей, чем методы максимального правдоподобия, [26] хотя они лучше способны учитывать отсутствующие данные. [27]

В то время как методы правдоподобия находят дерево, которое максимизирует вероятность данных, байесовский подход восстанавливает дерево, которое представляет наиболее вероятные клады, опираясь на апостериорное распределение. Однако оценки апостериорной вероятности кладов (измерение их «поддержки») могут быть довольно далеки от цели, особенно в кладах, которые не являются подавляюще вероятными. Таким образом, были предложены другие методы для оценки апостериорной вероятности. [28]

Некоторые инструменты, которые используют байесовский вывод для выведения филогенетических деревьев из данных о частоте аллелей вариантов (VAF), включают Canopy, EXACT и PhyloWGS. [29] [30] [31]

Выбор модели

Методы молекулярной филогенетики опираются на определенную модель замещения , которая кодирует гипотезу об относительных скоростях мутаций на различных участках вдоль изучаемых последовательностей гена или аминокислот. В простейшем случае модели замещения направлены на корректировку различий в скоростях переходов и трансверсий в последовательностях нуклеотидов. Использование моделей замещения обусловлено тем фактом, что генетическое расстояние между двумя последовательностями увеличивается линейно только в течение короткого времени после того, как две последовательности расходятся друг от друга (в качестве альтернативы расстояние становится линейным только незадолго до слияния ). Чем больше времени проходит после расхождения, тем больше вероятность того, что две мутации произойдут в одном и том же нуклеотидном участке. Таким образом, простые расчеты генетического расстояния будут недооценивать число событий мутаций, которые произошли в эволюционной истории. Степень этого недоучета увеличивается с увеличением времени с момента расхождения, что может привести к явлению притяжения длинных ветвей или неправильному назначению двух отдаленно связанных, но конвергентно эволюционирующих последовательностей как тесно связанных. [32] Метод максимальной экономии особенно восприимчив к этой проблеме из-за его явного поиска дерева, представляющего минимальное количество различных эволюционных событий. [2]

Типы моделей

Все модели замен назначают набор весов каждому возможному изменению состояния, представленному в последовательности. Наиболее распространенные типы моделей неявно обратимы, поскольку они назначают тот же вес, например, мутации нуклеотида G>C, что и мутации C>G. Самая простая возможная модель, модель Джукса-Кантора , назначает равную вероятность каждому возможному изменению состояния для данного нуклеотидного основания. Скорость изменения между любыми двумя различными нуклеотидами будет составлять одну треть от общей скорости замены. [2] Более продвинутые модели различают переходы и трансверсии . Самая общая возможная обратимая во времени модель, называемая моделью GTR, имеет шесть параметров скорости мутации. Еще более обобщенная модель, известная как общая 12-параметрическая модель, нарушает обратимость во времени ценой значительного дополнительного усложнения расчета генетических расстояний, которые согласованы между несколькими линиями. [2] Одна из возможных вариаций на эту тему корректирует скорости таким образом, что общее содержание GC — важная мера стабильности двойной спирали ДНК — меняется со временем. [33]

Модели также могут допускать вариацию скоростей в зависимости от позиций во входной последовательности. Наиболее очевидный пример такой вариации следует из расположения нуклеотидов в генах, кодирующих белок, в трехосновных кодонах . Если местоположение открытой рамки считывания (ORF) известно, скорости мутаций можно скорректировать для позиции заданного сайта в кодоне, поскольку известно, что спаривание оснований с колебаниями может допускать более высокие скорости мутаций в третьем нуклеотиде заданного кодона, не влияя на значение кодона в генетическом коде . [32] Менее основанный на гипотезах пример, который не полагается на идентификацию ORF, просто назначает каждому сайту скорость, случайно выбранную из предопределенного распределения, часто гамма-распределения или логнормального распределения . [2] Наконец, более консервативная оценка вариаций скоростей, известная как метод ковариона, допускает автокоррелированные вариации скоростей, так что скорость мутаций заданного сайта коррелирует между сайтами и линиями. [34]

Выбор лучшей модели

Выбор подходящей модели имеет решающее значение для проведения хорошего филогенетического анализа, как потому, что недостаточно параметризованные или чрезмерно ограничительные модели могут давать аномальное поведение, когда нарушаются их основные предположения, так и потому, что чрезмерно сложные или чрезмерно параметризованные модели требуют больших вычислительных затрат, а параметры могут быть переобучены. [32] Наиболее распространенным методом выбора модели является тест отношения правдоподобия (LRT), который дает оценку правдоподобия, которую можно интерпретировать как меру « хорошего соответствия » между моделью и входными данными. [32] Однако следует проявлять осторожность при использовании этих результатов, поскольку более сложная модель с большим количеством параметров всегда будет иметь более высокое правдоподобие, чем упрощенная версия той же модели, что может привести к наивному выбору моделей, которые являются слишком сложными. [2] По этой причине компьютерные программы выбора модели выберут самую простую модель, которая не будет значительно хуже, чем более сложные модели замещения. Существенным недостатком LRT является необходимость проведения серии попарных сравнений между моделями; Было показано, что порядок, в котором сравниваются модели, оказывает большое влияние на выбор той, которая в конечном итоге будет выбрана. [35]

Альтернативным методом выбора модели является информационный критерий Акаике (AIC), формально представляющий собой оценку расхождения Кульбака–Лейблера между истинной моделью и проверяемой моделью. Его можно интерпретировать как оценку правдоподобия с поправочным коэффициентом для штрафования сверхпараметризованных моделей. [32] AIC рассчитывается для отдельной модели, а не для пары, поэтому он не зависит от порядка оценки моделей. Связанная альтернатива, информационный критерий Байеса (BIC), имеет схожую базовую интерпретацию, но сильнее штрафует сложные модели. [32] Определение наиболее подходящей модели для реконструкции филогении является фундаментальным шагом во многих эволюционных исследованиях. Однако различные критерии выбора модели приводят к спорам о том, какой критерий предпочтительнее. Недавно было показано, что когда желаемым результатом являются топологии и реконструкция предковой последовательности, выбор одного критерия вместо другого не имеет решающего значения. Вместо этого использование самой сложной модели замены нуклеотидов, GTR+I+G, приводит к аналогичным результатам для вывода топологии дерева и предковых последовательностей. [36]

Полный пошаговый протокол построения филогенетических деревьев, включая сборку смежных последовательностей ДНК/аминокислот, множественное выравнивание последовательностей, модельный тест (тестирование наиболее подходящих моделей замен) и реконструкцию филогении с использованием максимального правдоподобия и байесовского вывода, доступен на сайте Protocol Exchange [37]

Нетрадиционный способ оценки филогенетического дерева — сравнить его с результатом кластеризации. Можно использовать технику многомерного масштабирования, так называемое интерполяционное соединение, чтобы уменьшить размерность и визуализировать результат кластеризации для последовательностей в 3D, а затем сопоставить филогенетическое дерево с результатом кластеризации. Лучшее дерево обычно имеет более высокую корреляцию с результатом кластеризации. [38]

Оценка поддержки деревьев

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

Узловая поддержка

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

Дерево консенсуса

Многие методы оценки поддержки узлов включают рассмотрение множественных филогений. Дерево консенсуса суммирует узлы, которые являются общими для набора деревьев. [39] В *строгом консенсусе* показаны только узлы, найденные в каждом дереве, а остальные свернуты в неразрешенную политомию . Менее консервативные методы, такие как дерево *консенсуса большинства*, рассматривают узлы, которые поддерживаются заданным процентом рассматриваемых деревьев (например, не менее 50%).

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

Бутстрэппинг и складной нож

В статистике бутстрап — это метод вывода изменчивости данных с неизвестным распределением с использованием псевдорепликаций исходных данных. Например, если задан набор из 100 точек данных, псевдореплика — это набор данных того же размера (100 точек), случайно выбранный из исходных данных с заменой. То есть каждая исходная точка данных может быть представлена ​​в псевдореплике более одного раза или вообще не представлена. Статистическая поддержка включает оценку того, имеют ли исходные данные схожие свойства с большим набором псевдорепликаций.

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

Статистическая строгость теста bootstrap была эмпирически оценена с использованием вирусных популяций с известной эволюционной историей [41] , обнаружив, что 70% поддержки bootstrap соответствует 95% вероятности того, что клада существует. Однако это было проверено в идеальных условиях (например, отсутствие изменений в эволюционных скоростях, симметричные филогении). На практике значения выше 70% обычно поддерживаются и оставляются исследователю или читателю для оценки достоверности. Узлы с поддержкой ниже 70% обычно считаются неразрешенными.

Jackknife в филогенетике — это похожая процедура, за исключением того, что столбцы матрицы выбираются без замены. Псевдореплики генерируются путем случайной подвыборки данных — например, «10% jackknife» будет включать случайную выборку 10% матрицы много раз для оценки узловой поддержки.

Апостериорная вероятность

Реконструкция филогений с использованием байесовского вывода генерирует апостериорное распределение высоковероятных деревьев с учетом данных и эволюционной модели, а не одно «лучшее» дерево. Деревья в апостериорном распределении обычно имеют много различных топологий. Когда входные данные представляют собой данные о частоте аллелей вариантов (VAF), инструмент EXACT может точно вычислять вероятности деревьев для небольших биологически значимых размеров деревьев, путем исчерпывающего поиска по всему пространству деревьев. [29]

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

Ожидается, что статистическая поддержка узла в байесовском выводе будет отражать вероятность того, что клад действительно существует с учетом данных и эволюционной модели. [42] Поэтому порог для принятия узла как поддерживаемого, как правило, выше, чем для бутстреппинга.

Методы подсчета шагов

Поддержка Бремера подсчитывает количество дополнительных шагов, необходимых для противоречия кладу.

Недостатки

У каждой из этих мер есть свои слабости. Например, клады меньшего или большего размера, как правило, привлекают большие значения поддержки, чем клады среднего размера, просто из-за количества таксонов в них. [43]

Поддержка Bootstrap может обеспечить высокие оценки поддержки узлов из-за шума в данных, а не из-за истинного существования клады. [44]

Ограничения и обходные пути

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

Гомоплазия

Некоторые признаки с большей вероятностью будут развиваться конвергентно , чем другие; логически, таким признакам следует придавать меньший вес при реконструкции дерева. [45] Веса в форме модели эволюции можно вывести из наборов молекулярных данных, чтобы можно было использовать методы максимального правдоподобия или байесовские методы для их анализа. Для молекулярных последовательностей эта проблема усугубляется, когда изучаемые таксоны существенно расходятся. По мере того, как время после расхождения двух таксонов увеличивается, увеличивается и вероятность множественных замен на одном и том же сайте или обратных мутаций, все из которых приводят к гомоплазиям. Для морфологических данных, к сожалению, единственным объективным способом определения конвергенции является построение дерева — несколько круговой метод. Тем не менее, взвешивание гомоплазированных признаков [ как? ] действительно приводит к более поддерживаемым деревьям. [45] Дальнейшее уточнение может быть достигнуто путем взвешивания изменений в одном направлении выше, чем изменений в другом; Например, наличие грудных крыльев почти гарантирует размещение среди крыловидных насекомых, поскольку, хотя крылья часто теряются вторично, нет никаких доказательств того, что они были приобретены более одного раза. [46]

Горизонтальный перенос генов

В целом, организмы могут наследовать гены двумя способами: вертикальный перенос генов и горизонтальный перенос генов . Вертикальный перенос генов — это передача генов от родителя к потомству, а горизонтальный (также называемый латеральным) перенос генов происходит, когда гены перескакивают между неродственными организмами, что является обычным явлением, особенно у прокариот ; хорошим примером этого является приобретенная устойчивость к антибиотикам в результате обмена генами между различными бактериями, что приводит к появлению видов бактерий с множественной лекарственной устойчивостью. Также были хорошо документированы случаи горизонтального переноса генов между эукариотами .

Горизонтальный перенос генов усложнил определение филогении организмов, и были зарегистрированы несоответствия в филогении среди определенных групп организмов в зависимости от генов, используемых для построения эволюционных деревьев. Единственный способ определить, какие гены были приобретены вертикально, а какие горизонтально, — это экономно предположить, что наибольший набор генов, которые были унаследованы вместе, был унаследован вертикально; это требует анализа большого количества генов.

Гибриды, видообразование, интрогрессии и неполная сортировка линий

Основное предположение, лежащее в основе математической модели кладистики, заключается в ситуации, когда виды аккуратно разделяются бифуркационным образом. Хотя такое предположение может иметь место в большем масштабе (за исключением горизонтального переноса генов, см. выше), видообразование часто гораздо менее упорядочено. Исследования с момента введения кладистического метода показали, что гибридное видообразование , когда-то считавшееся редким, на самом деле довольно распространено, особенно у растений. [47] [48] Также распространено парафилетическое видообразование , что делает предположение о бифуркационном шаблоне непригодным, что приводит к филогенетическим сетям , а не деревьям. [49] [50] Интрогрессия также может перемещать гены между в остальном различными видами, а иногда даже родами, [51] усложняя филогенетический анализ, основанный на генах. [52] Это явление может способствовать «неполной сортировке по родословной» и считается распространенным явлением в ряде групп. При анализе на уровне видов с этим можно справиться с помощью более крупной выборки или лучшего анализа всего генома. [53] Часто эту проблему можно обойти, ограничив анализ меньшим количеством не тесно связанных между собой образцов.

Выборка таксона

Благодаря развитию передовых методов секвенирования в молекулярной биологии стало возможным собирать большие объемы данных (ДНК или аминокислотные последовательности) для выведения филогенетических гипотез. Например, нередко можно найти исследования с матрицами признаков, основанными на целых митохондриальных геномах (~16 000 нуклеотидов у многих животных). Однако моделирование показало, что важнее увеличить количество таксонов в матрице, чем увеличить количество признаков, поскольку чем больше таксонов, тем точнее и надежнее полученное филогенетическое дерево. [54] [55] Это может быть отчасти связано с разрывом длинных ветвей .

Филогенетический сигнал

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

Непрерывные символы

Морфологические признаки, которые являются выборкой континуума, могут содержать филогенетический сигнал, но их трудно кодировать как дискретные признаки. Было использовано несколько методов, одним из которых является кодирование пробелов, и существуют вариации кодирования пробелов. [57] В исходной форме кодирования пробелов: [57]

Групповые средние значения для символа сначала упорядочиваются по размеру. Вычисляется объединенное внутригрупповое стандартное отклонение... и сравниваются различия между соседними средними значениями... относительно этого стандартного отклонения. Любая пара соседних средних значений считается различной и получает различные целочисленные оценки... если средние значения разделены "пробелом", большим, чем внутригрупповое стандартное отклонение... умноженное на некоторую произвольную константу.

Если в анализ добавляется больше таксонов, промежутки между таксонами могут стать настолько малыми, что вся информация будет потеряна. Обобщенное кодирование промежутков решает эту проблему, сравнивая отдельные пары таксонов, а не рассматривая один набор, содержащий все таксоны. [57]

Отсутствующие данные

В общем, чем больше данных доступно при построении дерева, тем точнее и надежнее будет полученное дерево. Отсутствующие данные не более вредны, чем просто меньше данных, хотя влияние сильнее всего, когда большая часть отсутствующих данных находится в небольшом количестве таксонов. Концентрация отсутствующих данных по небольшому количеству признаков дает более надежное дерево. [58]

Роль окаменелостей

Поскольку многие признаки включают в себя эмбриологические или мягкотканные или молекулярные признаки, которые (в лучшем случае) почти никогда не окаменевают, а интерпретация ископаемых более неоднозначна, чем интерпретация живых таксонов , вымершие таксоны почти всегда имеют более высокие доли недостающих данных, чем живые. Однако, несмотря на эти ограничения, включение ископаемых бесценно, поскольку они могут предоставить информацию в редких областях деревьев, разбивая длинные ветви и ограничивая промежуточные состояния признаков; таким образом, ископаемые таксоны вносят такой же вклад в разрешение дерева, как и современные таксоны. [59] Ископаемые также могут ограничивать возраст родословных и, таким образом, демонстрировать, насколько дерево соответствует стратиграфической записи; [1] стратокладистика включает информацию о возрасте в матрицы данных для филогенетического анализа.

Смотрите также

Ссылки

  1. ^ abc Халафванд, Тайлер (2015). «Поиск структуры в пространстве поиска филогении». Университет Далхаузи .
  2. ^ abcdefghijklmno Felsenstein J (2004). Выводы о филогениях . Сандерленд, Массачусетс: Sinauer Associates. ISBN 978-0-87893-177-4.
  3. ^ abcd Mount DM (2004). Биоинформатика: Анализ последовательностей и генома (2-е изд.). Cold Spring Harbor, Нью-Йорк: Cold Spring Harbor Laboratory Press. ISBN 978-0-87969-712-9.
  4. ^ Swiderski DL, Zelditch ML, Fink WL (сентябрь 1998 г.). «Почему морфометрия не является чем-то особенным: кодирование количественных данных для филогенетического анализа». Systematic Biology . 47 (3): 508–19. JSTOR  2585256. PMID  12066691.
  5. ^ Gaubert P, Wozencraft WC, Cordeiro-Estrela P, Veron G (декабрь 2005 г.). «Мозаики конвергенций и шума в морфологических филогениях: что есть в виверровых-подобных плотоядных?». Systematic Biology . 54 (6): 865–94. doi :10.1080/10635150500232769. PMID  16282167.
  6. ^ Strait DS, Grine FE (декабрь 2004 г.). «Вывод филогении гоминоидов и ранних гоминид с использованием краниодентальных признаков: роль ископаемых таксонов». Журнал эволюции человека . 47 (6): 399–452. Bibcode : 2004JHumE..47..399S. doi : 10.1016/j.jhevol.2004.08.008. PMID  15566946.
  7. ^ Wiens JJ (2001). «Анализ признаков в морфологической филогенетике: проблемы и решения». Systematic Biology . 50 (5): 689–99. doi :10.1080/106351501753328811. PMID  12116939.
  8. ^ Jenner RA (2001). «Билатерийская филогения и некритическая переработка морфологических наборов данных». Systematic Biology . 50 (5): 730–42. doi : 10.1080/106351501753328857 . PMID  12116943.
  9. ^ Сокал Р., Миченер К. (1958). «Статистический метод оценки систематических связей». Научный бюллетень Канзасского университета . 38 : 1409–1438.
  10. ^ Saitou N, Nei M (июль 1987). «Метод присоединения соседей: новый метод реконструкции филогенетических деревьев». Молекулярная биология и эволюция . 4 (4): 406–25. doi : 10.1093/oxfordjournals.molbev.a040454 . PMID  3447015.
  11. ^ Fitch WM , Margoliash E (январь 1967). «Построение филогенетических деревьев». Science . 155 (3760): 279–84. Bibcode : 1967Sci...155..279F. doi : 10.1126/science.155.3760.279. PMID  5334057.
  12. ^ Lespinats S, Grando D, Maréchal E, Hakimi MA, Tenaillon O, Bastien O (2011). «Как алгоритм Фитча-Марголиаша может выиграть от многомерного масштабирования». Evolutionary Bioinformatics Online . 7 : 61–85. doi :10.4137/EBO.S7048. PMC 3118699. PMID  21697992 . 
  13. ^ Day WH (1987). «Вычислительная сложность вывода филогений из матриц различий». Бюллетень математической биологии . 49 (4): 461–7. doi :10.1007/BF02458863. PMID  3664032. S2CID  189885258.
  14. ^ Хенди МД, Пенни Д (1982). «Алгоритмы ветвей и границ для определения минимальных эволюционных деревьев». Математические биологические науки . 59 (2): 277–290. doi :10.1016/0025-5564(82)90027-X.
  15. ^ Ратнер В.А., Жарких А.А., Колчанов Н, Родин С, Соловьев С, Антонов А.С. (1995). Молекулярная эволюция . Серия Биоматематика. Том. 24. Нью-Йорк: Спрингер-Верлаг. ISBN 978-3-662-12530-4.
  16. ^ Sankoff D, Morel C, Cedergren RJ (октябрь 1973 г.). «Эволюция 5S РНК и неслучайность замены оснований». Nature . 245 (147): 232–4. doi :10.1038/newbio245232a0. PMID  4201431.
  17. ^ ab De Laet J (2005). «Экономия и проблема неприменимости в данных последовательностей». В Albert VA (ред.). Экономия, филогения и геномика . Oxford University Press. стр. 81–116. ISBN 978-0-19-856493-5.
  18. ^ Wheeler WC, Gladstein DS (1994). «MALIGN: программа выравнивания множественных последовательностей нуклеиновых кислот». Журнал наследственности . 85 (5): 417–418. doi :10.1093/oxfordjournals.jhered.a111492.
  19. ^ Simmons MP (июнь 2004 г.). «Независимость выравнивания и поиска по дереву». Молекулярная филогенетика и эволюция . 31 (3): 874–9. Bibcode :2004MolPE..31..874S. doi :10.1016/j.ympev.2003.10.008. PMID  15120385.
  20. ^ De Laet J (2015). «Анализ экономии невыровненных данных последовательностей: максимизация гомологии и минимизация гомоплазии, а не минимизация операционально определенной общей стоимости или минимизация равновзвешенных преобразований». Cladistics . 31 (5): 550–567. doi :10.1111/cla.12098. PMID  34772278. S2CID  221582410.
  21. ^ Chor B, Tuller T (июнь 2005 г.). «Максимальное правдоподобие эволюционных деревьев: твердость и приближение». Биоинформатика . 21 (Приложение 1): i97–106. doi : 10.1093/bioinformatics/bti1027 . PMID  15961504.
  22. ^ El-Kebir M, Oesper L, Acheson-Field H, Raphael BJ (июнь 2015 г.). «Реконструкция клональных деревьев и состава опухоли на основе данных секвенирования нескольких образцов». Биоинформатика . 31 (12): i62-70. doi :10.1093/bioinformatics/btv261. PMC 4542783. PMID  26072510 . 
  23. ^ Malikic S, McPherson AW, Donmez N, Sahinalp CS (май 2015 г.). «Вывод о клональности в образцах множественных опухолей с использованием филогении». Биоинформатика . 31 (9): 1349–56. doi : 10.1093/bioinformatics/btv003 . PMID  25568283.
  24. ^ Mau B, Newton MA (1997). «Филогенетический вывод для двоичных данных на дендрограммах с использованием цепей Маркова Монте-Карло». Журнал вычислительной и графической статистики . 6 (1): 122–131. doi :10.2307/1390728. JSTOR  1390728.
  25. ^ Yang Z, Rannala B (июль 1997 г.). «Байесовский филогенетический вывод с использованием последовательностей ДНК: метод Монте-Карло с использованием цепей Маркова». Молекулярная биология и эволюция . 14 (7): 717–24. doi : 10.1093/oxfordjournals.molbev.a025811 . PMID  9214744.
  26. ^ Kolaczkowski B, Thornton JW (декабрь 2009 г.). Delport W (ред.). «Смещение притяжения длинных ветвей и непоследовательность в байесовской филогенетике». PLOS ONE . 4 (12): e7891. Bibcode : 2009PLoSO ...4.7891K. doi : 10.1371/journal.pone.0007891 . PMC 2785476. PMID  20011052. 
  27. ^ Симмонс МП (2012). «Вводящие в заблуждение результаты филогенетических анализов, основанных на правдоподобии, при наличии отсутствующих данных». Кладистика . 28 (2): 208–222. doi : 10.1111/j.1096-0031.2011.00375.x . PMID  34872185. S2CID  53123024.
  28. ^ Ларджет Б. (июль 2013 г.). «Оценка апостериорных вероятностей деревьев с использованием условных распределений вероятностей клад». Systematic Biology . 62 (4): 501–11. doi :10.1093/sysbio/syt014. PMC 3676676 . PMID  23479066. 
  29. ↑ Ab Рэй ​​С., Цзя Б., Сафави С., ван Опийнен Т., Исберг Р., Рош Дж., Бенто Дж. (22 августа 2019 г.). «Точный вывод в рамках идеальной модели филогении». arXiv : 1908.08623 . Бибкод : 2019arXiv190808623R. {{cite journal}}: Цитировать журнал требует |journal=( помощь )
  30. ^ Jiang Y, Qiu Y, Minn AJ, Zhang NR (сентябрь 2016 г.). «Оценка гетерогенности внутри опухоли и отслеживание продольной и пространственной клональной эволюционной истории с помощью секвенирования следующего поколения». Труды Национальной академии наук Соединенных Штатов Америки . 113 (37): E5528-37. Bibcode : 2016PNAS..113E5528J. doi : 10.1073/pnas.1522203113 . PMC 5027458. PMID  27573852 . 
  31. ^ Deshwar AG, Vembu S, Yung CK, Jang GH, Stein L, Morris Q (февраль 2015 г.). "PhyloWGS: реконструкция субклонального состава и эволюции с помощью секвенирования всего генома опухолей". Genome Biology . 16 (1): 35. doi : 10.1186/s13059-015-0602-8 . PMC 4359439 . PMID  25786235. 
  32. ^ abcdef Салливан Дж., Джойс П. (2005). «Выбор модели в филогенетике». Annual Review of Ecology, Evolution, and Systematics . 36 (1): 445–466. doi :10.1146/annurev.ecolsys.36.102003.152633. PMC 3144157. PMID  20671039 . 
  33. ^ Galtier N, Gouy M (июль 1998 г.). «Вывод шаблона и процесса: реализация неоднородной модели эволюции последовательности ДНК с максимальным правдоподобием для филогенетического анализа». Молекулярная биология и эволюция . 15 (7): 871–9. doi : 10.1093/oxfordjournals.molbev.a025991 . PMID  9656487.
  34. ^ Fitch WM, Markowitz E (октябрь 1970 г.). «Улучшенный метод определения изменчивости кодонов в гене и его применение к скорости фиксации мутаций в эволюции». Biochemical Genetics . 4 (5): 579–93. doi :10.1007/bf00486096. PMID  5489762. S2CID  26638948.
  35. ^ Pol D (декабрь 2004 г.). «Эмпирические проблемы иерархического теста отношения правдоподобия для выбора модели». Systematic Biology . 53 (6): 949–62. doi : 10.1080/10635150490888868 . PMID  15764562.
  36. ^ Abadi S, Azouri D, Pupko T, Mayrose I (февраль 2019 г.). «Выбор модели может не быть обязательным шагом для реконструкции филогении». Nature Communications . 10 (1): 934. Bibcode :2019NatCo..10..934A. doi :10.1038/s41467-019-08822-w. PMC 6389923 . PMID  30804347. 
  37. ^ Баст Ф. (2013). «Поиск сходства последовательностей, множественное выравнивание последовательностей, выбор модели, матрица расстояний и реконструкция филогении». Протокол обмена . doi : 10.1038/protex.2013.065 .
  38. ^ Ruan Y, House GL, Ekanayake S, Schütte U, Bever JD, Tang H, Fox G (26 мая 2014 г.). «Интеграция кластеризации и многомерного масштабирования для определения филогенетических деревьев как сферических филограмм, визуализированных в 3 измерениях». 2014 14-й Международный симпозиум IEEE/ACM по кластерным, облачным и грид-вычислениям . IEEE. стр. 720–729. doi :10.1109/CCGrid.2014.126. ISBN 978-1-4799-2784-5. S2CID  9581901.
  39. ^ Baum DA, Smith SD (2013). Древовидное мышление: введение в филогенетическую биологию. Roberts. стр. 442. ISBN 978-1-936221-16-5.
  40. ^ Felsenstein J (июль 1985 г.). «Пределы уверенности в филогениях: подход с использованием Bootstrap». Эволюция; Международный журнал органической эволюции . 39 (4): 783–791. doi :10.2307/2408678. JSTOR  2408678. PMID  28561359.
  41. ^ Хиллис Д. М., Булл Дж. Дж. (1993). «Эмпирический тест бутстреппинга как метода оценки достоверности филогенетического анализа». Systematic Biology . 42 (2): 182–192. doi :10.1093/sysbio/42.2.182. ISSN  1063-5157.
  42. ^ Huelsenbeck J, Rannala B (декабрь 2004 г.). «Частотные свойства байесовских апостериорных вероятностей филогенетических деревьев в простых и сложных моделях замещения». Systematic Biology . 53 (6): 904–13. doi : 10.1080/10635150490522629 . PMID  15764559.
  43. ^ Chemisquy MA, Prevosti FJ (2013). «Оценка эффекта размера клады в альтернативных мерах поддержки ветвей». Журнал зоологической систематики и эволюционных исследований . 51 (4): 260–273. doi : 10.1111/jzs.12024 . hdl : 11336/4144 .
  44. ^ Филлипс М.Дж., Дельсюк Ф., Пенни Д. (июль 2004 г.). «Филогения в масштабе генома и обнаружение систематических смещений» (PDF) . Молекулярная биология и эволюция . 21 (7): 1455–8. doi : 10.1093/molbev/msh137 . PMID  15084674.
  45. ^ ab Goloboff PA, Carpenter JM, Arias JS, Esquivel DR (2008). «Взвешивание против гомоплазии улучшает филогенетический анализ морфологических наборов данных». Cladistics . 24 (5): 758–773. doi : 10.1111/j.1096-0031.2008.00209.x . hdl : 11336/82003 . S2CID  913161.
  46. ^ Goloboff PA (1997). «Оптимизация с самовзвешиванием: поиск по дереву и реконструкции состояний символов при подразумеваемых затратах на преобразование». Cladistics . 13 (3): 225–245. doi :10.1111/j.1096-0031.1997.tb00317.x. PMID  34911233. S2CID  196595734.
  47. ^ Арнольд МЛ (1996). Естественная гибридизация и эволюция . Нью-Йорк: Oxford University Press. стр. 232. ISBN 978-0-19-509975-1.
  48. ^ Вендел Дж. Ф., Дойл Дж. Дж. (1998). «Секвенирование ДНК». В Soltis DE, Soltis PS , Doyle JJ (ред.). Молекулярная систематика растений II . Бостон: Клювер. стр. 265–296. ISBN 978-0-19-535668-7.
  49. ^ Funk DJ, Omland KE (2003). «Парафилия и полифилия на уровне видов: частота, причины и последствия с учетом митохондриальной ДНК животных». Annual Review of Ecology, Evolution, and Systematics . 34 : 397–423. doi : 10.1146/annurev.ecolsys.34.011802.132421. S2CID  33951905.
  50. ^ "Генеалогия жизни (GoLife)". Национальный научный фонд . Получено 5 мая 2015 г. Программа GoLife основана на программе AToL, приспосабливаясь к сложности моделей диверсификации на протяжении всей истории жизни. Наши текущие знания о таких процессах, как гибридизация, эндосимбиоз и латеральный перенос генов, ясно показывают, что эволюционная история жизни на Земле не может быть точно изображена — для каждой ветви дерева — как единое, типологическое, разветвляющееся дерево.
  51. ^ Kutschera VE, Bidon T, Hailer F, Rodi J, Fain SR, Janke A (2014). «Медведи в лесу генных деревьев: филогенетический вывод осложнен неполной сортировкой линий и потоком генов». Молекулярная биология и эволюция . 31 (8): 2004–2017. doi :10.1093/molbev/msu186. PMC 4104321. PMID  24903145 . 
  52. ^ Qu Y, Zhang R, Quan Q, Song G, Li SH, Lei F (декабрь 2012 г.). «Неполная сортировка линий или вторичная примесь: распутывание исторического расхождения с недавним потоком генов у попугайчика-барсука (Paradoxornis webbianus)». Молекулярная экология . 21 (24): 6117–33. Bibcode : 2012MolEc..21.6117Q. doi : 10.1111/mec.12080. PMID  23095021. S2CID  22635918.
  53. ^ Pollard DA, Iyer VN, Moses AM, Eisen MB (октябрь 2006 г.). «Широко распространенное несоответствие генных деревьев видовому дереву у Drosophila: доказательства неполной сортировки линий». PLOS Genetics . 2 (10): e173. doi : 10.1371/journal.pgen.0020173 . PMC 1626107 . PMID  17132051. 
  54. ^ Zwickl DJ, Hillis DM (август 2002 г.). «Увеличение выборки таксонов значительно снижает филогенетическую ошибку». Systematic Biology . 51 (4): 588–98. doi : 10.1080/10635150290102339 . PMID  12228001.
  55. ^ Wiens JJ (февраль 2006 г.). «Отсутствующие данные и дизайн филогенетических анализов». Журнал биомедицинской информатики . 39 (1): 34–42. doi : 10.1016/j.jbi.2005.04.001 . PMID  15922672.
  56. ^ Blomberg SP, Garland T, Ives AR (апрель 2003 г.). «Тестирование филогенетического сигнала в сравнительных данных: поведенческие черты более лабильны». Эволюция; Международный журнал органической эволюции . 57 (4): 717–45. doi :10.1111/j.0014-3820.2003.tb00285.x. PMID  12778543. S2CID  221735844.
  57. ^ abc Archie JW (1985). «Методы кодирования переменных морфологических признаков для численного таксономического анализа». Systematic Zoology . 34 (3): 326–345. doi :10.2307/2413151. JSTOR  2413151.
  58. ^ Prevosti FJ, Chemisquy MA (2009). «Влияние отсутствующих данных на реальные морфологические филогении: влияние количества и распределения отсутствующих записей». Cladistics . 26 (3): 326–339. doi : 10.1111/j.1096-0031.2009.00289.x. hdl : 11336/69010 . PMID  34875786. S2CID  86850694.
  59. ^ Cobbett A, Wilkinson M, Wills MA (октябрь 2007 г.). «Ископаемые оказывают такое же сильное воздействие, как и живые таксоны в экономичном анализе морфологии». Systematic Biology . 56 (5): 753–66. doi : 10.1080/10635150701627296 . PMID  17886145.

Дальнейшее чтение

Внешние ссылки