Кому мешает ДНК-генеалогия? Ложь, инсинуации, и русофобия в современной российской науке
на парах отец-сын, потому что подавляющее большинство из них вообще не дадут мутаций даже на 2000 парах, в лучшем случае некоторые дадут одну мутацию, из которых константу скорости не вычислить, погрешность таких расчетов составит плюс-минус 100 % даже при доверительном интервале в 68 % (одна сигма). Поразительно, что популяционные генетики используют такие «скорости мутаций», расчитанные из одной мутации, или даже без единой мутации в парах отец-сын (приведенные, например, в работах[69], [70]), даже не задумываясь, что смысла в этом нет практически никакого. Естественно, получают при этом невоспроизводимые «времена жизни общих предков», и объявляют, что по мутациям в гаплотипах считать нельзя (D. Pontikos[71]; Busby et al[72]).
Самые «быстрые» маркеры мутируют в среднем от частот (1/0.0073) один раз в 137 условных поколений, то есть примерно раз в 3400 лет, или один раз на 137 рождений мальчиков, до примерно (1/0.004) один раз в 250 условных поколений. Давайте посмотрим, как это выглядит на практическом примере. Воспользуемся тем же набором из 3466 гаплотипов гаплогруппы R1b-L21. В маркере DYS472, самом «медленном», аллели образуют следующий набор:
7 – 1 (то есть аллель 7 встречается в 3466 маркерах DYS472 всего один раз)
8 – 3461 раз
9 – 4 раза
Таким образом, мы наблюдаем всего пять мутаций на 3466 маркеров DYS472 за 152 условных поколений (расчеты см. выше), прошедших со времени жизни общего предка этих гаплотипов. Это дает константу скорости мутации, равную примерно 5/3466/152 = 0.95 × 10-5 мутаций на условное поколение. Если давать более строгую формулировку, то при пяти мутациях надо написать (0.95±0.43) × 10-5 мутаций на условное поколение, или, что более правильно, (1.0±0.4) × 10-5 мутаций на условное поколение. Погрешности здесь рассчитываются по обычным правилам статистики, как квадрат обратной величины квадратного корня из числа мутаций. Здесь мы пренебрегли поправкой на возвратные мутации, потому что она в данном случае практически ничего не меняет, поправка составит менее одной сотой доли процента.
Но этот расчет проведен только на одной серии гаплотипов. В серии из 976 гаплотипов гаплогруппы R1a в маркере DYS472 прошла всего одна мутация, что при 154 условных поколениях до общего предка дает константу скорости мутации 1/976/154 = 0.67 × 10-5 на условное поколение. В целом, при подобном рассмотрении серий гаплотипов из 24 разных субкладов средняя величина для константы скорости мутации для DYS472 оказалась равной 0.8 × 10-5 мутаций на условное поколение.
Заметим, что это – самая «медленная скорость мутации из всех 111 маркеров. Для остальных задача решается еще проще, там мутаций больше.
Приведем еще несколько примеров «медленных» маркеров, расчеты по которым самые сложные. По парам отец-сын они вообще не определяются. Например, маркер DYS455. В серии из 3466 аллелей этого маркера наблюдаем следующее распределение:
9 – 6 (то есть аллель 9 встречается в 3466 маркерах DYS455 всего 6 раз)
10 – 28
11 – 3409
12 – 23
Считая, что все мутации одношаговые, получаем 63 мутации на 3466 аллелей DYS455 за те же 152 условных поколения, прошедших со времени жизни общего предка этих гаплотипов. Это дает константу скорости мутации, равную 63/3466/152 = 0.00012±0.00002 мутаций на маркер DYS455 за условное поколение. Поправка на возвратные мутации здесь составляет менее одного процента, поэтому вводить ее бессмысленно. В таблице выше дана константа