Анатолий Клёсов - Кому мешает ДНК-генеалогия?

Но это еще не всё, и погрешность часто еще значительно больше. Дело в том, что снип-мутации происходят не через равные промежутки времени, а неупорядоченно. Если при этом рассматриваются цепочки снипов протяженностью десятки и сотни тысяч лет, то это аналогично рассмотрению ситуации с сотнями бросков монеты, и погрешности из-за неупорядоченности нивелируются. Но для относительно недавних снипов метод расчета очень неточный. Например, если снип образовался 750 лет назад, то в цепочке нисходящих снипов их всего 750/144 = 5 снипов, и погрешность определения времени образования снипа, согласно формулам математической статистики, не меньше ±46 % (с 68 %-ной достоверностью) или ±92 % (с 95%ной достоверностью). В итоге погрешности у большинства датировок по снипам в таблице значительно более высокие, чем у датировок по мутациям в гаплотипах, которые рассматривает ДНК-генеалогия.

Есть еще серьезная проблема – в расчетах компании YFull часто игнорируют фактическое количество снипов, и ставят ту датировку, которая им представлется более правильной, «по понятиям». При этом они всегда постулируют, что при разветвлении снипов они образовались точно в одно и то же время, хотя это события совершенно независимые. По аналогии, если у отца двое сыновей, то их возраст далеко не всегда одинаковый, один мог родиться, когда отцу было 18 лет, другой – когда отцу было 70 лет. Такой постулат, «о равенстве возраста сыновей» часто приводит к ситуациям, когда к одному снипу от вышестоящего ведет 5 мутаций (в среднем 720 лет), к другому 20 мутаций (в среднем 2880 лет), а YFull записывает их «возраст» как одинаковый, либо 720, либо 2880 лет, выбор датировки фактически произвольный. В итоге получаются несуразные датировки.

Примеры – в переходе R1a-Y35 > YP7278 зафиксировано 9 промежуточных снипов (это приведено в таблице на сайте YFull), что соответствует 1300 лет разницы, однако датировка для обоих снипов, родительского и нисходящего, приведена как одинаковая, 4300 лет назад. В переходе CTS3402 > CTS2613 зафиксировано 17 промежуточных снипов (2450 лет разницы), однако датировка для обоих снипов в списке YFull – одинаковые 4300 лет назад, вместо 1850 лет назад для нисходящего, как должно быть (и как дает более точная датировка по гаплотипам). И так в таблицах YFull повторяется десятки и сотни раз.

Приходится на этом столь подробно останавливаться, потому что в последнее время усиленно распространяется легенда, которую безграмотно пропагандирует Балановский, что расчеты по снипам – значительно более точные, чем любыми другими методами. Историки могут принять эту легенду за чистую монету, и опять повторится ситуация последних 15 лет, когда практически все расчеты популяционными генетиками, опубликованные в академических журналах, были кардинально искажены, с завышениями до 300–400 %. Для расчетов на времена удаленные, как правило, более 5-10 тысяч лет назад, и тем более на 100–200 тысяч лет назад, когда число снип-мутаций исчисляется многими сотнями, расчеты по снипам действительно могут оказаться полезными, но вводимые постулаты, как равенство датировок для «параллельных» снипов, опять сводит эту пользу к нулю, если не к отрицательным величинам.

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

Но, как мы видим, попгенетики опять в своем амплуа – Балановский уже провозгласил, что метод расчетов по снипам самый точный, и ненавистные ему расчетные подходы ДНК-генеалогии можно отставить. Это опять было провозглашено без какого-либо исследования, и единственная статья, в которой Балановский этот метод использовал, дала неверные результаты. Дело в том, что Балановский в своей манере выхватил из многих вариантов лишь одну скорость мутации – 122 года на снип, без малейшего ее обоснования или проверки, и опять получил неверные датировки (Underhill et al, 2015, 32 автора, Eur. J. Human Genetics).

Теперь обратимся к статье Балановского с 22 соавторами (PLOS ONE, 2015), которую он представляет в своей «дискуссии» наряду с образцовыми, в которой, как объявляет, определил константу скорости мутации, равную 0.78Х10 -9мутаций на нуклеотид в год. Надо же, и точность примечательная, с надежностью до второго знака, то есть почти 1 %, что, конечно, не бывает и быть не может, но он так объявляет. И на диаграмме в своей «дискуссии» он своё значение помещает прямо в яблочко, остальные светила зарубежной попгенетики гуляют со своими значениями вокруг, а «наша – лучше всех». Зная уровень Балановского как ниже среднего, я позволил себе усомниться, и данные пересчитал. И что получилось? Балановский взял данные по девяти казахам-аргынам, общий предок которых (по легендам) упомянут в хрониках в 1405 году, набросил 50 лет, выбрал из упомянутых девяти казахов шесть, привел их средний год рождения, и записал, что от общего предка до их современных шестерых потомков прошло 606 лет – замечаете, как работают попгенетики? Они возятся с мелочами, развивая якобы бурную деятельность, чтобы на «настоящую науку» было похоже, но пролетают по крупному там, в чем безграмотны. В итоге получается пшик. Так и здесь, после всей этой возни с единицами лет Балановский написал, что среднее число снип-мутаций у шестерых выбранных казахов от их общего предка равно 4.67.

Когда я вижу, что среднее число снип-мутаций указывают с точностью до третьего знака, у меня, как у известного персонажа Пелевина, рука тянется к пистолету. Те, кто так «подают данные», либо жулики, либо безграмотные. Либо оба в одном флаконе. Смотрим на данные. Вот – числа мутаций от «общего предка» до каждого из шести современных потомков: 5, 5, 4, 7, 5, 2. Любой, имеющий элементарное арифметическое образование должен знать, что среднего из этого набора с точностью до третьего знака (то есть фактически до 0.1 %) быть просто не может. Среднее здесь равно 4.67±1.11 без округления, с натяжкой – 4.7±1.1, а на самом деле 5±1. Но пусть пока без округления, тогда средняя скорость мутаций равна не 0.78×10 -9, как объявляет Балановский, а (0.77±0.18) × 10 -9мутаций на нуклеотид в год, или между 0.59 × 10 -9и 0.95 × 10 -9, то есть в том же интервале, где это определения были еще несколько лет назад, то есть между 0.6 и 1.0 миллиардных долей. И нужно ему было считать годы рождения казахов-аргынов с точностью до одного года? Мог спокойно им всем по сто лет поставить, и ничего бы в его расчетах не изменилось. Ничего нового Балановский со своим коллективом в уточнение константы скорости мутации не внес. Просто еще раз продемонстрировал вопиющую математическую безграмотность. И, конечно, полное отсутствие научной школы, что, впрочем, давно было ясно.