Страница 8 из 22
Так бывает всегда. Если аллель нейтрален, его частота будет «случайно блуждать» между нулем и единицей до тех пор, пока не упрется либо в верхний, либо в нижний предел. Рано или поздно нейтральный аллель либо зафиксируется (достигнет частоты 1), либо элиминируется — исчезнет из генофонда. Третьего не дано. Таким образом, блуждания хоть и случайны, их исход предсказуем. Мы точно знаем, что дело кончится либо фиксацией, либо элиминацией. Чем больше популяция, тем в среднем дольше придется ждать исхода[17], но он все равно неизбежен.
Здесь самое время сказать, что случайные, не зависящие от приспособленности изменения частот аллелей называются генетическим дрейфом. Под властью дрейфа находятся все нейтральные мутации (а в маленьких популяциях отчасти также и вредные с полезными, но об этом позже).
Можно ли вычислить вероятность того, что нейтральный аллель в конце концов зафиксируется, а не элиминируется? Да, это просто. В нашем примере аллели А1 и А2 исходно находились в равном положении, потому что у обоих была частота 0,5. Очевидно, что в такой ситуации шансы на фиксацию у них должны быть одинаковыми и равными 0,5. В половине случаев А1 зафиксируется, А2 элиминируется, в половине — наоборот.
Ну а если конкурирующих нейтральных аллелей не два, а, скажем, четыре и начальные частоты у них тоже одинаковые (0,25)? В этом случае дело кончится фиксацией одного из аллелей и элиминацией трех остальных, причем для каждого аллеля вероятность фиксации равна 0,25. Таким образом, очевидно, что вероятность фиксации нейтральной мутации в будущем просто-напросто равна ее частоте в данный момент: Pfix = q.
Если вы, дорогие читатели, еще не устали от этой примитивной математики, то позвольте познакомить вас еще с двумя простыми, интересными и полезными формулами.
Сколько нейтральных мутаций будет фиксироваться в популяции в каждом поколении? (Имеются в виду мутации уже не в одном и том же, а в разных генах.) Если мы сумеем это вычислить, то получим прекраснейший инструмент — молекулярные часы. Тогда мы сможем по количеству нейтральных генетических различий определить, когда жил последний общий предок сравниваемых организмов.
Вывод этой формулы — подлинный шедевр «биологической математики». Судите сами. Определим сначала, сколько новых мутаций появляется в популяции в каждом поколении. Будем считать для простоты, что подавляющее большинство мутаций нейтральны (это недалеко от истины). Ответ очевиден: U × N, где U — темп мутагенеза (среднее число новых мутаций у каждой новорожденной особи), N — численность популяции. Определить U можно, просто сравнивая геномы детей и родителей.
Теперь нужно понять, какая часть из этих U × N только что появившихся мутаций в итоге зафиксируется. Это и будет искомая величина — темп фиксации нейтральных мутаций в генофонде популяции за одно поколение. Обозначим ее буквой V.
Мы уже знаем, что вероятность фиксации мутации равна ее частоте: Pfix = 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]. Быстро меняющиеся участки ДНК используются для датировки недавних событий, а те, что меняются медленно, хороши для датировки событий глубокой древности.
Дрейф и отбор: кто кого?
Генетический дрейф царствует над нейтральными мутациями (аллелями), отбор — над полезными и вредными. Отбор, повышающий частоту полезных мутаций, называют положительным. Отбор, отбраковывающий вредные мутации, — отрицательным, или очищающим.
Впрочем, все так хорошо и просто только в больших популяциях. В маленьких ситуация сложнее, потому что отбор и дрейф — две главные движущие силы эволюции — начинают конкурировать друг с другом за контроль над слабовредными и слабополезными мутациями.
Посмотрим, как это происходит. Воспользуемся для этого той же моделью с мюмзиками. Изменим только характер мутации, которая привела к появлению аллеля А2. До сих пор мы считали, что мутация была нейтральной. Поскольку она была нейтральной, обладатели аллелей А1 и А2 имели одинаковую плодовитость. Давайте теперь предположим, что мутация была полезной, что она повысила плодовитость мюмзиков на 5 %. Это можно смоделировать так: пусть мюмзики с генотипом А1 рожают по 20 детенышей, а мюмзики с генотипом А2 — по 21. Начальная частота аллеля А2 пусть будет по-прежнему равна 0,5. Только теперь мы рассмотрим популяции с разной численностью (N).
На рисунке [см. ниже] показано, как будет меняться частота А2 в том случае, если его носители имеют пятипроцентное адаптивное преимущество. Мы видим, что в большой популяции (N = 5000) частота А2 неуклонно растет, приближаясь к единице. Так работает отбор в идеальных для него условиях, т. е. в больших популяциях, где влияние дрейфа на полезные и вредные аллели пренебрежимо мало. Форма у кривой довольно правильная, и это наводит на мысль, что ее можно описать какой-нибудь математической формулой. Это действительно так, но выводить формулу мы не будем, чтобы не утомить читателей (а любители математики могут сделать это самостоятельно)[19].
Перед нами, между прочим, важнейший эволюционный процесс — аллельное замещение, т. е. вытеснение более приспособленным аллелем менее приспособленного. Процесс идет не слишком быстро. В большой популяции, например состоящей из миллиона особей, для того чтобы зафиксировалась новая полезная мутация, дающая 5-процентное адаптивное преимущество, требуется около 560 поколений. А ведь 5 % — это серьезное преимущество. Такие мутации — редкость. Ждать, пока зафиксируется мутация с преимуществом в 1 %, придется уже 2800 поколений! Тем не менее в большой популяции отбор «чувствует» даже самую незначительную разницу в приспособленности. Это обеспечивает эффективную (пусть и медленную) фиксацию слабополезных мутаций и элиминацию слабовредных[20].
17
Не вдаваясь в сложные вычисления, заметим лишь, что число поколений, которое (в среднем) должно пройти до фиксации или элиминации нейтральной мутации, имеет примерно тот же порядок величин, что и численность популяции. Если в популяции несколько тысяч особей, то и ждать придется несколько тысяч поколений, в миллионной популяции речь будет идти о миллионах поколений. Поэтому на судьбу нейтральных мутаций сильное влияние оказывают «бутылочные горлышки» — периоды сокращения численности популяции, во время которых нейтральный аллель может потеряться или зафиксироваться особенно быстро.
18
Потому что разные участки ландшафта приспособленности имеют разную «проходимость»: одни похожи на ровные плато (там нейтральная эволюция идет быстро), другие — на лабиринт узких тропинок над пропастями (там нейтральные изменения накапливаются медленно). Подробнее об этом см. ниже.
19
Ответ такой: q = (Fx × k)/(1 + Fx × k), где q — частота аллеля А2 по прошествии X поколений, F — относительная приспособленность аллеля А2 по сравнению с конкурирующим аллелем А1 (в нашем случае F = 21/20 = 1,05), k = q0/(1 − q0), где q0 — частота аллеля А2 в начальный момент времени. В бесконечно большой популяции динамика q должна в точности соответствовать этой формуле. В маленьких популяциях соответствие будет неточным из-за генетического дрейфа.
20
Медленность элиминации слабовредных мутаций может пойти на пользу популяции. Ведь некоторые мутации, вредные здесь и сейчас, могут оказаться полезными в перспективе. Главное, чтобы это потенциальное преимущество успело реализоваться до того, как отбор элиминирует мутацию. Мы рассмотрим такие коллизии в следующих главах.