Чихуахуа оказались отчасти койотами: в ДНК нашли следы древней гибридизации

Миниатюрная собачка весом менее трех килограммов, дрожащая в дорожной сумке, давно стала символом декоративного питомца. Но репутация чихуахуа как хрупкого аксессуара разбивается о реальность: эти собаки способны без колебаний атаковать овчарку, бросаться на незнакомцев и яростно защищать территорию. Кинологи десятилетиями списывали это на «синдром маленькой собаки» — комплекс поведенческих проблем, вызванный тем, что владельцы прощают карликовым питомцам агрессию, за которую крупного пса немедленно изолировали бы. Биологическая реальность оказалась сложнее.

Исследовательская группа из Стэнфордского университета опубликовала в журнале Nature Communications масштабную работу по деконволюции предкового состава геномов с помощью системы машинного обучения Gnomix. Анализируя геномы сотен пород собак и диких псовых, ученые зафиксировали в ДНК современной чихуахуа неопровержимые следы интрогрессии — включения генетического материала дикого североамериканского койота (Canis latrans). Это открытие меняет взгляд на эволюцию древнейших пород Нового Света: чихуахуа — не продукт европейской декоративной селекции, утративший рабочие качества, а миниатюрный хищник с подтвержденной примесью дикой крови.

Чтобы понять, как гены дикого зверя попали в комнатную собаку, нужно обратиться к истории коренных популяций Америки. Собаки заселили континент вместе с предками коренных народов, перейдя по Берингову перешейку тысячи лет назад. Предком чихуахуа считается течичи — небольшая собака, которую разводили тольтеки и ацтеки. В хрониках испанских конкистадоров XVI века упоминаются свидетельства того, что коренные народы практиковали преднамеренное скрещивание домашних собак с койотами. Индейские селекционеры стремились объединить послушание животного с адаптивными качествами дикого зверя: устойчивостью к болезням, обостренными органами чувств, выносливостью и выраженной территориальной агрессией.

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

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

Чтобы доказать факт гибридизации, генетикам требовался инструмент локального анализа предков. До настоящего времени алгоритмы локального анализа хромосом обладали тремя фатальными недостатками:

• Они требовали колоссальных вычислительных мощностей и обрабатывали полногеномные последовательности днями и неделями.
• Они полагались на жесткие математические допущения о времени смешения популяций, которые невозможно проверить для древних событий.
• Точность существующих программ резко падала, если гибридизация произошла глубже 20–30 поколений назад: короткие предковые фрагменты алгоритмы принимали за ошибки секвенирования.

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

На первом этапе хромосома разделяется на равные отрезки — окна. Для полных геномов размер окна составил 1000 однонуклеотидных полиморфизмов (SNP). Генетические последовательности не инвариантны к сдвигу: одна и та же цепочка нуклеотидов в начале хромосомы и в ее конце может указывать на разное происхождение. Поэтому инженеры обучили отдельную базовую модель машинного обучения для каждого конкретного окна. Чтобы классификатор не терял данные на стыках, применили концепцию расширенного контекста: каждое окно считывается моделью вместе с дополнительными участками с обеих сторон. Для полных геномов лидером стала логистическая регрессия — за счет скорости и способности учитывать взаимосвязи между соседними мутациями. Для низкоплотных данных создан специальный математический аппарат — строковое ядро со снижением ковариации (CovRSK).

На втором этапе подключается модуль-сглаживатель. Он скользит поверх полученных вероятностей, объединяя информацию из 75 соседних окон. Исследователи отказались от традиционных формул популяционной динамики и обучили модуль на базе алгоритма экстремального градиентного бустинга (XGBoost). Модель напрямую обучилась на миллионах симулированных примеров распознавать физические закономерности биологического кроссинговера.

Геном диплоиден: он состоит из материнской и отцовской нитей ДНК. Секвенаторы не умеют разделять считываемые мутации по родительским хромосомам. Стандартные программы сборки фаз регулярно совершают ошибки переключения. Вместе с Gnomix стэнфордцы создали алгоритм Gnofix. Он оценивает уверенность модели в присвоении предкового статуса сегменту. Если на хромосоме возникает неестественная чересполосица, Gnofix моделирует локальные перестановки участков между материнской и отцовской цепочками и оставляет ту конфигурацию, которая максимизирует предсказательную уверенность сглаживателя. Там, где специализированный модуль фазирования предыдущего отраслевого стандарта (программы RFMix) работал 37 минут, алгоритм Gnofix выполнил задачу за 7 секунд.

В результате ученые получили инструмент, способный без потери точности (выше 90%) находить предковые фрагменты даже на глубине в 100 поколений с момента гибридизации — то есть более чем на 300 лет в прошлое для собак и на 3000 лет для человека.

Исследователи применили систему Gnomix к геномам 531 собаки различных пород и диких псовых, используя обучающую панель из 191 высококачественного эталонного генома. Когда алгоритм построил филогенетическое дерево на основе среднего числа попарных нуклеотидных различий, сегменты чихуахуа, помеченные моделью как «койот», образовали четкий совместный кластер именно с дикими североамериканскими койотами, а не с волками или другими группами собак. Алгоритм локализовал эти участки на хромосомах чихуахуа, зафиксировав устойчивую интрогрессию. Это означает, что предки современных чихуахуа действительно скрещи

По материалам: iXBT