Так бывает всегда. Если аллель нейтрален, его частота будет «случайно блуждать» между нулем и единицей до тех пор, пока не упрется либо в верхний, либо в нижний предел. Рано или поздно нейтральный аллель либо зафиксируется (достигнет частоты 1), либо элиминируется – исчезнет из генофонда. Третьего не дано. Таким образом, блуждания хоть и случайны, их исход предсказуем. Мы точно знаем, что дело кончится либо фиксацией, либо элиминацией. Чем больше популяция, тем в среднем дольше придется ждать исхода[17], но он все равно неизбежен.

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

Можно ли вычислить вероятность того, что нейтральный аллель в конце концов зафиксируется, а не элиминируется? Да, это просто. В нашем примере аллели А₁ и А₂ исходно находились в равном положении, потому что у обоих была частота 0,5. Очевидно, что в такой ситуации шансы на фиксацию у них должны быть одинаковыми и равными 0,5. В половине случаев А₁ зафиксируется, А₂ элиминируется, в половине – наоборот.

Ну а если конкурирующих нейтральных аллелей не два, а, скажем, четыре и начальные частоты у них тоже одинаковые (0,25)? В этом случае дело кончится фиксацией одного из аллелей и элиминацией трех остальных, причем для каждого аллеля вероятность фиксации равна 0,25. Таким образом, очевидно, что вероятность фиксации нейтральной мутации в будущем просто-напросто равна ее частоте в данный момент: P_fix = q.

Если вы, дорогие читатели, еще не устали от этой примитивной математики, то позвольте познакомить вас еще с двумя простыми, интересными и полезными формулами.

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

Вывод этой формулы – подлинный шедевр «биологической математики». Судите сами. Определим сначала, сколько новых мутаций появляется в популяции в каждом поколении. Будем считать для простоты, что подавляющее большинство мутаций нейтральны (это недалеко от истины). Ответ очевиден: U × N, где U – темп мутагенеза (среднее число новых мутаций у каждой новорожденной особи), N – численность популяции. Определить U можно, просто сравнивая геномы детей и родителей.

Теперь нужно понять, какая часть из этих U × N только что появившихся мутаций в итоге зафиксируется. Это и будет искомая величина – темп фиксации нейтральных мутаций в генофонде популяции за одно поколение. Обозначим ее буквой V.

Мы уже знаем, что вероятность фиксации мутации равна ее частоте: P_fix = q. Чему же равна частота только что появившейся мутации? Это совсем просто. Поскольку мутация только что появилась, она пока есть только у одной особи. Следовательно, ее частота равна 1/N. Вот, собственно, и все. Умножаем количество новых мутаций на вероятность фиксации каждой из них (т. е. на долю мутаций, которые в итоге зафиксируются) и получаем ответ: V = U × N × 1/N. Поразительно! Численность популяции волшебным образом сокращается, и величина N уходит из уравнения. Мы приходим к выводу, что темп фиксации нейтральных мутаций не зависит от численности и равен просто-напросто скорости мутагенеза: V = U. Не знаю, как вы, дорогие читатели, а мы, биологи, приходим в восторг от такой красоты.

Нам остался один шаг до молекулярных часов. После того как два вида, произошедшие от общего предка, разделились, в их генофондах независимо накапливаются нейтральные мутации. За время t (измеряемое в поколениях) первый вид накопит V × t = U × t нейтральных мутаций, и второй вид накопит столько же. Совместными усилиями они накопят 2U × t нейтральных отличий друг от друга. Зная темп мутагенеза U и подсчитав число различий между геномами сравниваемых видов (обозначим его буквой D), определяем время жизни их последнего общего предка: t = D/2U. Это и есть знаменитые молекулярные часы.

Если бы темп мутагенеза и скорость смены поколений были одинаковыми у всех живых существ, все было бы совсем просто. Но они, конечно, неодинаковы, поэтому в вычисления приходится вносить поправки. А еще есть проблема насыщения: в какой-то момент – обычно по прошествии десятков, а чаще сотен миллионов лет – гены разделившихся когда-то видов оказываются настолько «переполнены» нейтральными различиями, что величина D перестает расти, хотя нейтральные мутации продолжают фиксироваться. Ниже мы увидим пример исследования, показавшего, как уровень сходства между «случайно блуждающими» последовательностями приблизился к минимально возможному и дальше снижаться уже не мог.

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

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

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

Посмотрим, как это происходит. Воспользуемся для этого той же моделью с мюмзиками. Изменим только характер мутации, которая привела к появлению аллеля А₂. До сих пор мы считали, что мутация была нейтральной. Поскольку она была нейтральной, обладатели аллелей А₁ и А₂ имели одинаковую плодовитость. Давайте теперь предположим, что мутация была полезной, что она повысила плодовитость мюмзиков на 5 %. Это можно смоделировать так: пусть мюмзики с генотипом А₁ рожают по 20 детенышей, а мюмзики с генотипом А₂ – по 21. Начальная частота аллеля А₂ пусть будет по-прежнему равна 0,5. Только теперь мы рассмотрим популяции с разной численностью (N).

На рисунке (см. ниже) показано, как будет меняться частота А₂ в том случае, если его носители имеют пятипроцентное адаптивное преимущество. Мы видим, что в большой популяции (N = 5000) частота А₂ неуклонно растет, приближаясь к единице. Так работает отбор в идеальных для него условиях, т. е. в больших популяциях, где влияние дрейфа на полезные и вредные аллели пренебрежимо мало. Форма у кривой довольно правильная, и это наводит на мысль, что ее можно описать какой-нибудь математической формулой. Это действительно так, но выводить формулу мы не будем, чтобы не утомить читателей (а любители математики могут сделать это самостоятельно)[19].

Перед нами, между прочим, важнейший эволюционный процесс – аллельное замещение, т. е. вытеснение более приспособленным аллелем менее приспособленного. Процесс идет не слишком быстро. В большой популяции, например состоящей из миллиона особей, для того чтобы зафиксировалась новая полезная мутация, дающая 5-процентное адаптивное преимущество, требуется около 560 поколений. А ведь 5 % – это серьезное преимущество. Такие мутации – редкость. Ждать, пока зафиксируется мутация с преимуществом в 1 %, придется уже 2800 поколений! Тем не менее в большой популяции отбор «чувствует» даже самую незначительную разницу в приспособленности. Это обеспечивает эффективную (пусть и медленную) фиксацию слабополезных мутаций и элиминацию слабовредных[20].