Компьютерра: Компьютерра PDA N160 (18.02.2012-24.02.2012)

В модели рассматривается популяция, состоящая из свободно скрещивающихся организмов. Как это и принято в СТЭ, фенотипы (и приспособленность организмов) определяются их генами. Рассмотрим два гена, A и B, каждый из который представлен двумя аллелями: A и a для гена A, B и b для гена B. Итак, в популяции могут существовать девять генотипов: AABB, AABb, AAbb, AaBB, AaBb, Aabb, aaBB, aaBb, aabb.

Слишком просто? И этого минимума хватит, чтобы комбинаторика генов оказалась не вполне тривиальной.

Итак, зададим численность популяции (обозначим её K – ёмкость среды – максимальное количество особей, которое может обитать в данных условиях). Зададим начальные доли генотипов (P AAbb; P AaBBи так далее). Что дальше?

По СТЭ, источник новых признаков – новые гены, возникающие в результате мутаций; не будем пока с этим спорить. Зададим в модели частоты, с которыми происходят переходы одних аллелей в другие: P a→A, P A→a, P b→Bи P B→b.

Осталось предусмотреть ту силу, которая должна менять соотношение аллелей в определённом направлении, – отбор. Для каждого генотипа укажем его приспособленность (fitness): F AAbb; F AaBBи прочие.

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

При вычислении состава гамет учитывается вероятность мутаций. Комбинации гамет определяют генотипы потомства. Шансы на выживание потомков зависят от их приспособленности.

Всё, цикл работы модели окончен. Осталось определить состав следующего поколения. Округление долей генотипов в популяции до единиц особей носит вероятностный характер. К примеру, величина 1,4 округлится до 2 с вероятностью 0,4, и до 1 - с вероятностью 0,6.

Осталось повторить описанный цикл много раз (модель построена для 500 поколений). Выведем на график самый интересный, с точки зрения СТЭ, показатель – динамику соотношения альтернативных аллелей. Компоненты модели я расположил на листе Excel так, чтобы в один экран попадали и график, и окошки для ввода входных значений (обозначения перечислены выше). Скрины этого экрана и иллюстрируют дальнейшее изложение.

Ну что, начнём.

Начальный состав популяции полностью состоит из особей aabb. Приспособленность всех генотипов одинакова, состав популяции постепенно меняется вследствие мутаций (по гену A чуть быстрее, чем по гену B, из-за разной частоты мутирования).

Если в начальном состоянии популяции по каждому гену представлен лишь один аллель, мутации постепенно повышают долю альтернативных аллеей. А что будет, если соотношение аллелей окажется равновесным?

Начальное соотношение аллелей равное и колеблется под влиянием случайностей (мутаций и смещения при дискретизации). Обратите внимание на изменение масштаба по оси ординат.

В отсутствие отбора частоты аллелей начнут "гулять" вокруг среднего значения. Кстати, иногда такие случайные смещения могут привести к потере одного из аллелей.

Пора "включать" отбор. Пусть особи, которые обладают "удачным" генотипом, имеют на 1 процент большие шансов на выживание, чем все прочие. Вам кажется, что такое малое преимущество не может сыграть свою роль в эволюции?

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

Если приспособленность особи определяется единственным геном, преимущества в 1% для носителей доминантного гена достаточно для довольно быстрого изменения соотношения аллелей. А как отразится на скорости эволюции рецессивность аллеля, которому способствует отбор?

Отбор в пользу рецессивного аллеля (из состояния равенства частот двух альтернатив) идёт ещё эффективнее!

Ускорит (при условии начального равенства частот аллелей). Смотрите: за время, рассмотренное в модели, аллель, которому благоприятствует отбор, почти вытеснил свою альтернативу!

В действительности и популяции, и виды отличаются не по одному, а по множеству генов. В нашей модели мы можем рассмотреть только одновременный отбор по двум из них. Как скажется на скорости эволюции вовлечение второго гена с ещё одной парой аллелей? По гену B тоже поддержим отбором рецессивный аллель (b).