Статистический анализ полученных данных показал, что изменения приспособленности соответствуют степенной модели с неограниченным (хотя и замедляющимся) ростом. Это значит, что приспособленность вовсе не собирается выходить на плато. Иными словами, бактерии, по всей видимости, не намерены останавливаться на достигнутом. Несмотря на долгую жизнь в одних и тех же условиях, у бактерий все еще время от времени закрепляются новые мутации, в данных условиях полезные (напомним, что в эволюционной биологии слово “полезный” обычно используется в смысле “повышающий приспособленность”).
рис. 1.1. Рост приспособленности бактерий за 50 000 поколений. Точками (первые 20 000 поколений) и кружочками (следующие 30 000 поколений) показаны усредненные значения по всем подопытным популяциям. По вертикальной оси – относительная приспособленность (скорость роста по сравнению с предковым штаммом). Пунктиром показаны предсказания моделей: гиперболической (приспособленность асимптотически стремится к предельному уровню; нижняя кривая) и степенной (неограниченный, хотя и замедляющийся рост приспособленности). Параметры моделей основаны на данных по первым 20 000 поколений. Видно, что гиперболическая модель хуже предсказывает динамику приспособленности за последующие 30 000 поколений, чем степенная. По рисунку из Wiser et al., 2013.
Удивительная неисчерпаемость эволюционных возможностей, впервые столь наглядно продемонстрированная в эксперименте, имеет огромное значение для понимания эволюции. Ведь можно было предположить (многие так и думали), что в постоянных условиях, тем более в монокультуре, возможности для адаптивной эволюции быстро исчерпаются и система придет в состояние стазиса. В таком случае объяснить продолжающуюся эволюцию жизни на нашей планете можно было бы только непостоянством среды (абиотической и биотической).
По-видимому, запаса потенциально полезных мутаций хватит не только на 50 000 поколений бактерий, но и на 50 000 поколений наблюдающих за ними исследователей. Как ни странно, идущий безостановочно эволюционный процесс не приведет к какому-то невероятному, запредельному росту приспособленности. Ученые рассчитали, что если приспособленность будет и впредь увеличиваться по той же траектории (описываемой степенной моделью), то пятидесятитысячный по счету преемник Ленски увидит бактерий, удваивающих свою численность каждые 23 минуты. Это высокая скорость размножения для бактерий, но не какая-то неслыханная (численность бактерий предкового штамма удваивается за 55 минут).
Откуда же берутся новые полезные мутации? Ведь за время эксперимента каждая возможная нуклеотидная замена уже наверняка успела произойти хотя бы у одной бактерии! Основных причин, по-видимому, две. Первая связана с эпистазом – своего рода взаимодействием между мутациями. Например, мутация Б становится полезной только после того, как бактерия уже обзавелась мутацией А. Здесь нужно помнить еще и о том, что многие полезные мутации имеют негативные побочные эффекты. Мутация А может улучшить работу одной молекулярной системы, немного повредив другой. Допустим, первый эффект сильнее второго, поэтому в суммарном зачете мутация А оказывается “полезной” и поддерживается отбором. Но теперь, когда в популяции закрепилась мутация А, полезной окажется любая компенсаторная мутация Б, сглаживающая негативный побочный эффект мутации А.
Вторая причина неисчерпаемости запаса полезных мутаций состоит в том, что мутации с очень слабым полезным эффектом требуют очень долгого времени для закрепления в популяции. Такая мутация, появившись, с большой вероятностью будет вскоре потеряна из-за генетического дрейфа – случайных колебаний частот аллелей. Чтобы отбор начал хоть немного “помогать” такой мутации распространяться, число ее носителей должно стать довольно большим, а до тех пор она будет находиться целиком под властью безжалостного, неразборчивого дрейфа (см. Исследование № 4). Расчеты показывают, что подавляющее большинство вновь возникающих полезных мутаций теряются. Чтобы данная слабополезная мутация в конце концов все-таки закрепилась, она должна многократно появиться у разных бактерий независимо. Ленски и его коллеги рассчитали, что мутация, повышающая приспособленность бактерий на одну миллионную, должна возникнуть примерно 250 000 раз, прежде чем она наконец распространится и зафиксируется в популяции. Если принять во внимание частоту мутирования подопытных бактерий (в среднем одна мутация на 10
нуклеотидов за поколение) и размер их популяций (эффективная численность каждой из них, с учетом ежедневных бутылочных горлышек, составляет примерно 33 000 000 клеток), получится, что такой мутации понадобится 100 000 000 поколений, чтобы “спастись” от дрейфа, и еще миллионы поколений, чтобы зафиксироваться (достичь стопроцентной частоты). Таким образом, накопление слабополезных мутаций в эксперименте Ленски вряд ли закончится в обозримом (и даже необозримом) будущем.
Причины замедления роста приспособленности связаны с тем, что каждая закрепившаяся полезная мутация делает последующие генетические усовершенствования в среднем менее полезными. Авторы называют этот эффект “эпистазом убывающей доходности” (по аналогии с экономическим законом убывающей доходности). Иными словами, сначала происходят наиболее радикальные адаптивные изменения, а затем идет все более тонкая настройка и оптимизация фенотипа.
Еще один фактор, способствующий замедлению роста приспособленности и тесно связанный с предыдущим, – клональная интерференция, то есть конкуренция между клонами бактерий с разными полезными мутациями. Напомним, что бактерии, участвующие в эксперименте, – бесполые, они не могут обмениваться генами. Поэтому, если у одной бактерии возникнет одна полезная мутация, а у другой – другая, эти мутации не смогут объединиться в одном геноме (как это произошло бы у нормальных микробов, способных к горизонтальному переносу генов, и уж тем более у организмов, размножающихся половым путем; см. Исследование № 7). Вместо этого потомки первой бактерии начнут не на жизнь, а на смерть (не на свою, конечно, а своей линии) конкурировать с потомками второй. Победит, естественно, тот клон, чья мутация окажется полезнее. Вторая, менее полезная мутация будет вытеснена и потеряна. Придется теперь ждать, пока она случайно появится снова у бактерии, уже имеющей первую мутацию. А закрепиться она сможет, только если у нее не окажется более удачливых конкурентов. Клональная интерференция во многом объясняет, почему на первых этапах эксперимента фиксировались в основном мутации с сильным полезным эффектом (в среднем первая закрепившаяся в каждой популяции мутация повышала приспособленность на 10 % – это очень много), а затем уже постепенно начинали фиксироваться все менее и менее полезные мутации.
Разумеется, у бактерий, способных к горизонтальному переносу генов, адаптация шла бы гораздо быстрее (см. Исследование № 7). Но и анализировать результаты было бы труднее, потому что к мутационной изменчивости добавилась бы комбинативная, а на эволюцию стали бы влиять такие факторы, как наследственные различия по склонности к приему и передаче генов, избирательность при выборе партнеров и т. д.
Еще один важный результат дало сравнение роста приспособленности у популяций, сохранивших исходный (низкий) темп мутагенеза, и тех, где закрепились аллели-мутаторы, несущие мутации, которые резко (в среднем на два порядка) повысили темп мутирования. На сегодняшний день “гипермутаторами” стали шесть из двенадцати популяций.
рис. 1.2. Повышение темпа мутагенеза ускоряет адаптацию. Черные точки – усредненные данные по шести популяциям, в которых темп мутагенеза остался на исходном низком уровне. Серые треугольники – усредненные данные по трем популяциям, в которых на ранних этапах эксперимента (за первые 20 000 поколений) закрепились аллели-мутаторы. Видно, что приспособленность у вторых росла быстрее, чем у первых. Еще три популяции, у которых мутаторы закрепились поздно, в данном случае не рассматривались. По рисунку из Wiser et al., 2013.
Оказалось, что у “гипермутаторов” адаптация протекала ускоренными темпами (рис. 1.2). Получается, повышенный темп мутагенеза пошел бактериям на пользу. Это противоречит распространенной идее о том, что в стабильных условиях организмам выгодно снизить темп мутирования до нуля. Ведь большинство вновь возникающих мутаций вредны, поэтому мутация, повышающая темп мутагенеза, в краткосрочной перспективе приносит больше вреда, чем пользы, и отбор, по идее, должен действовать против нее.
Почему же аллели-мутаторы все-таки распространяются? Дело в том, что в бесполой популяции они часто спасаются от отбраковки благодаря эффекту, который называют генетическим автостопом. Если какой-то гипермутабельной бактерии повезет и у нее появится редкая и очень полезная мутация, отбор начнет поддерживать последнюю. При этом вместе с полезной мутацией будет автоматически распространяться и сидящий в том же геноме аллель-мутатор. Он будет, подобно пассажиру-безбилетнику, ехать прицепом к хромосоме, распространение которой в генофонде популяции “оплачено” полезной мутацией. А все потому, что в бесполой популяции, где нет перетасовки генов и хромосомных участков между отдельными клетками, отбираться могут только целые геномы, но не отдельные гены.
В популяции, способной к генетическому обмену, ситуация будет иной. Связка мутатора с полезной мутацией не будет неразрывной, что позволит отбору работать с ними индивидуально. В итоге отбор сможет закрепить полезную мутацию и отбраковать породивший ее аллель-мутатор (см. Исследование № 7). Но это пока лишь теория: эволюционный эксперимент на организмах, способных к генетическому обмену, сопоставимый по масштабу с экспериментом Ленски, еще не поставлен.
Итак, мы увидели, что предела приспособленности, по всей видимости, нет. Об этом нам говорят экспериментальные данные. Этот вывод кажется контринтуитивным, потому что “высшую точку” приспособленности можно легко себе представить как нечто реальное и достижимое. Казалось бы, для заданных постоянных условий должно существовать некое оптимальное, наиболее приспособленное состояние – такое, что у организма, его достигшего, никакая мутация уже не сможет повысить приспособленность. Даже если для реальной эволюционирующей популяции этот оптимум недостижим, он продолжает оставаться для биологов удобной абстракцией, упрощением, помогающим решать определенные задачи. Эволюционисты привыкли представлять себе этот оптимум в виде горной вершины на воображаемом ландшафте. Но при этом важно понимать, какие пути ведут к высшей точке, или, образно выражаясь, какой рельеф имеют склоны горы. Упрощенное понимание естественного отбора рисует нам гладкие склоны и прямой путь к вершине. Но это, очевидно, не соответствует биологическим реалиям. Благодаря стремительному развитию науки, “путь к вершине приспособленности” постепенно перестает быть абстракцией и начинает поддаваться картированию. Об этом – следующая глава.
Исследование № 2
Трудный путь к совершенству по ландшафту приспособленности
Одним из способов представить себе, как идет эволюция, является картирование ландшафта приспособленности конкретных белков. В схематичном виде ландшафт приспособленности показан на илл. i (см. цветную вклейку). Конечно, в идеале хотелось бы просканировать все “пространство последовательностей” (все возможные аминокислотные последовательности белка), определив для каждой его точки (то есть для каждой последовательности) эффективность выполнения белком интересующей нас функции. Однако такая задача невыполнима, поскольку даже для маленькой белковой молекулы в 100 аминокислот число точек в пространстве последовательностей (то есть число всех возможных белков такого размера) на много порядков превышает число атомов во Вселенной (первое равно 20
? 10
, второе – в 10
раз меньше).
Поэтому лучшее, на что мы можем пока рассчитывать, – это картирование очень маленьких, специально подобранных, удобных для изучения участков пространства последовательностей. Даже такая задача невероятно трудоемка. Но все же сегодня, благодаря развитию биотехнологий, она уже выполнима. Об этом свидетельствует статья Анны Подгорной и Майкла Лауба из Массачусетского технологического института (США), опубликованная в 2015 году в журнале Science (Podgornaia, Laub, 2015). Авторы сосредоточили усилия на четырех ключевых аминокислотах, входящих в состав белка PhoQ все той же кишечной палочки (Escherichia coli). Белки, гомологичные PhoQ, есть и у других гамма-протеобактерий, таких как сальмонелла и чумная палочка. PhoQ – это рецептор, располагающийся на внутренней мембране бактерии. Он реагирует на изменения концентрации ионов Mg
снаружи от мембраны (в так называемом периплазматическом пространстве, которое у гамма-протеобактерий и других грамотрицательных бактерий находится между внутренней и наружной мембранами) и передает сигнал в цитоплазму посредством активации или инактивации другого белка, называемого PhoP. Белок PhoP, в свою очередь, включает и выключает гены, необходимые для адаптации бактерии к меняющимся условиям среды.
Ключевую роль в работе двухкомпонентной регуляторной системы PhoQ-PhoP играют четыре аминокислоты, входящие в состав PhoQ и образующие так называемый белок-белковый интерфейс. Это значит, что они необходимы для того, чтобы PhoQ мог безошибочно опознать своего партнера PhoP и правильно провзаимодействовать с ним.
В белке PhoQ кишечной палочки это следующие четыре аминокислоты: аланин в позиции 284, валин в позиции 285, серин в позиции 288, треонин в позиции 289 (сокращенно Ala284, Val285, Ser288, Thr289 или еще короче – AVST).
Но AVST – лишь одна из 160 000 (20
) возможных комбинаций четырех аминокислот. Ведь чисто теоретически на каждой из этих четырех позиций могла бы стоять любая из 20 аминокислот. Исследователи поставили перед собой амбициозную задачу откартировать все пространство последовательностей для данных четырех аминокислотных позиций. То есть для каждой из 160 000 комбинаций проверить, обеспечивает ли она нормальное взаимодействие PhoQ с PhoP, а значит, и правильную работу всей регуляторной системы.
Для этого были изготовлены плазмиды с геном phoQ, в котором триплеты нуклеотидов, кодирующие аминокислоты белок-белкового интерфейса AVST, были полностью рандомизированы, то есть заменены на случайные тройки нуклеотидов. Затем эти плазмиды вставили в клетки E. coli, из которых их собственный ген phoQ был заблаговременно удален. Кроме того, в геном подопытного штамма вставили ген желтого флуоресцирующего белка под управлением такой регуляторной последовательности, которая обеспечивает активацию этого гена в присутствии активной формы белка PhoP (чтобы клетки, в которых система PhoQ-PhoP сработала, можно было отличить по флуоресценции). В итоге получилась разношерстная популяция кишечных палочек, у каждой из которых в белке PhoQ на месте четырех аминокислот AVST находилась какая-то из 160 000 возможных аминокислотных комбинаций.
Теперь нужно было выбрать из этой популяции те клетки, у которых система PhoQ-PhoP по-прежнему исправно работала. Отбор проводился в два этапа. Сначала микробов выращивали при разных концентрациях ионов магния и отбирали тех, у которых интенсивность желтого свечения соответствовала таковой у контрольных бактерий с последовательностью AVST. Затем отобранных микробов помещали в среду, не содержащую ионов магния. Пережить такое “магниевое голодание” способны только бактерии с исправно работающей системой PhoQ-PhoP.
Геномы бактерий, выдержавших эти испытания, подвергли секвенированию, чтобы получить полный список аминокислотных последовательностей, обеспечивающих нормальную работу системы PhoQ-PhoP. Таких последовательностей оказалось неожиданно много: целых 1659. Ученые выборочно проверили некоторые из мутантных штаммов на способность конкурировать с контрольными кишечными палочками при разных концентрациях ионов магния. Эти эксперименты показали, что приспособленность мутантов осталась на том же уровне, что и у исходного штамма. Таким образом, по крайней мере некоторые (а скорее всего, большинство) из выявленных 1659 последовательностей не имеют серьезных дефектов по сравнению с исходной последовательностью AVST.
Это один из главных результатов работы. Он показывает, что “плато” ландшафта приспособленности, соответствующее функциональному белок-белковому интерфейсу, весьма обширно: оно занимает более 1 % всего пространства последовательностей (1659 точек из 160 000). Это значит, что, подобно генетическому коду, в котором 20 аминокислот и знак “стоп” кодируются 64 триплетами, “белок-белковый” код, обеспечивающий узнавание белками друг друга, тоже обладает высокой избыточностью. Есть много разных вариантов, способных работать одинаково, что, по идее, должно облегчать эволюцию белков. Чем обширнее плато ландшафта приспособленности, соответствующее данной функции, тем больше шансов, что с какой-нибудь его точки удастся перейти (или перепрыгнуть) на другое плато. Иными словами, повышается вероятность смены функции белка в ходе эволюции.
Полученные данные позволили изучить топографию плато и оценить его “эволюционную проходимость” (илл. i на цветной вклейке). Принципиальный вопрос, от ответа на который зависит ход эволюции, а конкретнее динамика перемещения эволюционирующей молекулы по ландшафту приспособленности, состоит в том, насколько сильно влияние четырех аминокислотных позиций друг на друга. Это влияние называется эпистазом. Если эпистаз отсутствует, то плато ландшафта приспособленности, соответствующее данной функции белка, представляет собой ровную, плоскую поверхность, по которой эволюционирующий белок может гулять абсолютно свободно. Есть некий набор “разрешенных” аминокислотных замен, и приобретать их можно в любом порядке. При сильном эпистазе плато становится похоже на лабиринт узких тропинок, разделенных пропастями. По такому плато передвигаться трудно, потому что аминокислотные замены, безвредные в одних контекстах, оказываются вредными в других. В результате порядок приобретения замен становится крайне важным: к каждой точке пространства последовательностей можно теперь добраться лишь строго определенными немногочисленными маршрутами. В таком случае многие работоспособные последовательности, скорее всего, до сих пор не найдены эволюцией просто потому, что у нее не хватило на это времени.
Анализ полученного списка из 1659 работоспособных последовательностей показал сильнейшее взаимное влияние четырех аминокислотных позиций, то есть очень сильный эпистаз. Это второй главный результат исследования. В наличии эпистаза можно убедиться разными способами. Например, можно сравнить друг с другом последовательности, отличающиеся от AVST одной, двумя, тремя и четырьмя аминокислотами. Исследователи обнаружили среди функциональных (“разрешенных”) последовательностей 16 одиночных мутантов, 100 двойных, 544 тройных и 998 четверных, то есть таких, в которых все четыре аминокислоты отличаются от исходных A, V, S и T. Если бы эпистаза не было, то разрешенные множественные мутации представляли бы собой просто-напросто произвольные комбинации одиночных. Однако ничего подобного не наблюдается. Каждая отдельная аминокислотная замена оказывается разрешенной в одних сочетаниях и запрещенной в других. Из-за этого кратчайший разрешенный эволюционный путь от одной последовательности к другой в большинстве случаев оказывается длиннее, чем кратчайший из потенциально возможных (рис. 2.1).
Рассмотрим, например, переход от AVST к разрешенной последовательности MLAI (рис. 2.1, вверху слева). Чтобы превратить AVST в MLAI, нужно заменить четыре аминокислоты. Соответственно, если считать аминокислотную замену единичным эволюционным событием, то кратчайший путь от AVST к MLAI будет состоять из четырех шагов. Разных путей такой длины потенциально существует 24, поскольку аминокислоты можно заменять в разном порядке. Однако из этих 24 кратчайших путей разрешенными оказались только два. Все остальные пути требуют пересечения “пропасти” на ландшафте приспособленности, то есть прохождения через нефункциональный промежуточный этап.
рис. 2.1. Кратчайшие разрешенные эволюционные пути от AVST к последовательностям MLAI и SIQC. Функциональные (разрешенные) последовательности обведены рамками. Контурными буквами обозначены аминокислоты, отсутствующие в начальной и конечной последовательностях. Два левых рисунка в качестве единичного эволюционного события предполагают замену аминокислоты в белке, на правом рисунке таким событием считается замена одного нуклеотида в ДНК. В последнем случае путь получается длиннее из-за структуры генетического кода: от каждой аминокислоты путем замены одного нуклеотида можно прийти не ко всем, а лишь к некоторым из остальных 19 аминокислот. По рисунку из Podgornaia, Laub, 2015.
Во многих случаях все кратчайшие пути оказываются запрещенными. Тогда маршрут из одной точки пространства последовательностей в другую (например, путь от AVST к SIQC, показанный на рис. 2.1 слева внизу) оказывается длиннее, чем кратчайший из потенциально возможных (кратчайший путь предполагает три промежуточных состояния, а на рисунке их четыре).
Ну а если учесть, что реальные единичные мутации обычно представляют собой замену одного нуклеотида в ДНК (а не триплета на любой другой триплет), то проходимость ландшафта окажется еще ниже. На рис. 2.1 справа показан кратчайший разрешенный путь от AVST к MLAI с учетом этого обстоятельства. Теперь он состоит не из четырех, а из целых десяти шагов. Это объясняется структурой генетического кода, которая позволяет путем единичных нуклеотидных замен переходить от данной аминокислоты только к некоторым, но не к любым другим аминокислотам. Например, фенилаланин кодируется триплетами UUU и UUC. Изменив в любом из этих триплетов только одну букву, можно заменить в белковой молекуле фенилаланин на лейцин, изолейцин, валин, серин, тирозин или цистеин. Для перехода к любой из оставшихся 13 аминокислот потребуется не одна, а две или три нуклеотидных замены. Таким образом, эпистаз и структура генетического кода совместными усилиями делают ландшафт приспособленности труднопро- ходимым.
В целом топография изученного плато ландшафта приспособленности оказалась довольно замысловатой. Это плато мало похоже на гладкое, легкопроходимое поле. Оно подразделяется на несколько областей, внутри каждой из которых эволюционирующий белок может передвигаться относительно свободно, однако переход в соседнюю область возможен лишь по немногочисленным тропкам. Некоторые разрешенные последовательности вообще оказываются недосягаемыми из обжитой реальными бактериями области пространства последовательностей. По-видимому, эволюция белок-белкового интерфейса действительно настолько затруднена эпистазом и структурой генетического кода, что эволюция просто не успела разведать многие потенциально доступные области пространства последовательностей.
Данное исследование кому-то может показаться слишком скучным и специальным. Но только не биологам, привыкшим считать “пространство последовательностей” и “ландшафт приспособленности” абстракциями, полезными для теоретических построений, но недосягаемыми для реального изучения. Больше всего впечатляет сам факт, что стало в принципе возможным прямое экспериментальное картирование отдельных областей пространства последовательностей. Хотя, конечно, речь пока идет лишь о четырех аминокислотных позициях, то есть о переборе 20
вариантов. Полное картирование пространства всех возможных белков навсегда останется неразрешимой задачей.
Исследование № 3
“Эволюция умнее, чем ты”: рождение экологического разнообразия
Эксперимент Ленски изначально был спланирован так, чтобы свести к минимуму все “осложняющие обстоятельства”: изменения среды, генетический обмен, экологические взаимодействия между организмами. Ученые хотели получить в чистом виде самый главный эволюционный процесс – адаптацию к среде на основе мутаций и отбора. Однако, как метко заметил биохимик Лесли Орджел, “эволюция умнее, чем ты”. Он имел в виду, что исследователям, утверждающим, будто эволюция на что-то не способна, скорее всего, просто не хватает воображения. Как выясняется, эволюция не боится сложностей и “в чистом виде” ничего не демонстрирует, порождая, вопреки чаяниям ученых, куда более замысловатые результаты, чем от нее ждут. В подопытных популяциях Ленски, существующих, казалось бы, в самых простых условиях, какие только можно придумать, стали сами собой зарождаться экологические взаимодействия, основанные на диверсификации (разделении) ниш. А это, в свою очередь, заставляет бактерий заново приспосабливаться к меняющейся биотической обстановке (Good et al., 2017).
На этот раз Ленски и его коллеги провели генетический анализ всей замороженной “ископаемой летописи” эксперимента, накопившейся за 60 000 бактериальных поколений и насчитывающей около 1440 проб (по 120 проб на каждую из двенадцати популяций). Для каждой пробы был проведен метагеномный анализ с 50-кратным покрытием. Это значит, что из пробы выделяли ДНК и секвенировали случайные фрагменты геномов до тех пор, пока каждый участок генома кишечной палочки не оказался “прочтен” в среднем 50 раз. Этого оказалось достаточно, чтобы идентифицировать все новые мутации, которые возникали в подопытных популяциях и достигали частоты не менее 10 % (то есть встречались как минимум у каждой десятой бактерии) хотя бы в двух пробах. Мутации, не получившие столь широкого распространения, не учитывались, потому что их трудно отличить от случайных ошибок секвенирования. В итоге получилась детальная реконструкция эволюционного процесса в двенадцати популяциях.
Выводы о том, что рост приспособленности замедлился, но не прекратился, подтвердились (см. Исследование № 1). Темп накопления новых мутаций остался высоким.
Главное же открытие состоит вот в чем. Динамика накопления мутаций не вписывается в простейшую модель, согласно которой эволюция монокультуры бесполых организмов в стабильных условиях сводится к последовательной фиксации отбором вновь возникающих полезных мутаций. Эта модель не может объяснить наблюдаемую картину даже с учетом таких осложняющих обстоятельств, как генетический автостоп и клональная интерференция, о которых мы говорили выше.
Оказалось, что многие мутации, достигнув некоторой частоты, вдруг перестают распространяться, то есть двигаться дальше в сторону фиксации (стопроцентной частоты встречаемости). А ведь именно таков должен быть естественный ход событий, если клон с данной мутацией имеет более высокую приспособленность, чем другие бактерии. Может быть, распространение мутации остановилось из-за того, что появился более приспособленный конкурент? Но тогда прежние чемпионы должны постепенно вытесняться из популяции и исчезать. Однако этого тоже не происходит. Частота мутации начинает колебаться около какого-то промежуточного значения. Эти колебания могут продолжаться десятки тысяч поколений. В чем же дело?
Метагеномные данные, полученные для каждой из 1440 проб, представляют собой множество отсеквенированных кусочков ДНК, принадлежащих разным бактериям. Поэтому нельзя сразу понять, какие мутации относятся к одному клону, а какие – к разным. Однако ученым удалось разобраться в этом, проанализировав согласованность изменений частот мутаций во времени (поскольку частоты мутаций, находящихся в одном и том же геноме, меняются синхронно). В итоге выяснилось, что по крайней мере в девяти из двенадцати подопытных популяций в течение длительного времени (свыше 10 000 поколений) имело место устойчивое сосуществование как минимум двух разных клад (эволюционных линий, ветвей). Внутри этих клад шли свои собственные эволюционные процессы, то есть появлялись и фиксировались различные мутации.
Это значит, что в большинстве популяций произошла диверсификация. Разные клады как-то поделили между собой экологические ниши и стали устойчиво сосуществовать, приспосабливаясь теперь уже не только к изначально заданным условиям среды, но и к специфическому и переменчивому биотическому окружению.