В популяционной генетике статистики частот аллелей обычно изучают при слабом отборе или нейтральности. Однако во многих биологических сценариях, таких как адаптация вирусов, бактерий или искусственный отбор, отбор силён относительно дрейфа. Рядом с этой работой лежат исследования, которые либо фокусируются на долговременных результатах при сильном отборе, либо используют одну детерминированную траекторию. Данная статья заполняет пробел: она рассматривает переходные (зависимые от времени) статистики в режиме сильного отбора, где стандартные приближения дают большие ошибки из-за игнорирования траекторий, ведущих к потере аллеля.
Автор работает в рамках диффузионного приближения модели Райта-Фишера. Вводится параметр R = 2Ne|s|, где Ne — эффективный размер популяции, s — коэффициент отбора. При R >> 1 отбор силён. Статистика частоты аллеля представляется как взвешенная сумма двух условных ожиданий: при условии фиксации и при условии потери. Затем каждое условное ожидание переписывается через модифицированный процесс с изменённым отбором. Далее применяется разложение по малому параметру ε = 1/R, где √ε — амплитуда шума (дрейфа). Приближение удерживает члены до порядка ε, что даёт вклад от двух детерминированных траекторий и флуктуаций вокруг них. Уравнения для детерминированных траекторий и флуктуационных поправок решаются численно.
Получены явные выражения для среднего и дисперсии частоты аллеля как функций времени. Среднее приближается суммой взвешенных детерминированных траекторий и поправки порядка 1/R. Дисперсия содержит член, возникающий из-за различия двух детерминированных траекторий, и член от флуктуаций. Численное сравнение с точными результатами модели Райта-Фишера показывает, что чисто детерминированное приближение может давать ошибку до 100% для дисперсии, особенно когда начальная частота мала. Включение флуктуационной поправки снижает ошибку до нескольких процентов. Например, при R=50 и начальной частоте y=0.001 ошибка среднего падает с 10.6% до 3.2%, а дисперсии — с 12.6% до 4.9%. Флуктуационный вклад особенно важен на ранних временах, когда детерминированная дисперсия почти нулевая.
Работа показывает, что даже при сильном отборе нельзя пренебрегать траекториями, ведущими к потере аллеля. В отличие от соседних исследований, которые используют одну детерминированную траекторию или рассматривают только долговременные пределы, здесь предложен систематический метод для переходных статистик. Это заполняет важный пробел: теперь можно корректно вычислять временную динамику частот аллелей в сильном режиме отбора, что было невозможно с помощью стандартных подходов. Метод также выявляет, что флуктуации не всегда являются малыми поправками: для дисперсии они могут доминировать на определённых интервалах времени.
Метод может быть расширен на случай множественных аллелей, где потребуется суммирование по фиксациям разных аллелей. Также возможно включение мутаций, что потребует иного подхода, поскольку стационарное распределение не описывается двумя детерминированными траекториями. Обобщение на изменяющийся размер популяции и переменный отбор также представляет интерес. Кроме того, можно исследовать поведение в пограничных слоях вблизи частот 0 и 1, где асимптотические разложения могут требовать отдельного рассмотрения.
Метод применим для анализа временных рядов частот аллелей в экспериментальной эволюции, вирусных и бактериальных популяциях, а также в программах искусственного отбора. Он может улучшить оценку параметров отбора из данных и помочь в предсказании судьбы полезных аллелей в малых популяциях, что важно для природоохранной генетики.
Следующим шагом может быть применение метода к реальным данным временных рядов аллельных частот для проверки его эффективности в оценке коэффициентов отбора. Также стоит исследовать расширение на случай множественных аллелей и включение мутаций.
В популяционной генетике остаётся открытой проблема точного описания переходной динамики частот аллелей при сильном отборе, особенно для дисперсии и высших моментов. Данная работа напрямую решает эту проблему для двухаллельного случая без мутаций, предлагая вычислительно эффективный метод. Она также поднимает вопрос о применимости асимптотических разложений вблизи границ частотного диапазона.
Даже при сильном отборе, когда детерминированная сила велика, дисперсия частоты аллеля может почти полностью определяться случайными флуктуациями в течение длительного времени. Например, при R=50 и начальной частоте y=0.001 детерминированный вклад в дисперсию практически нулевой до времени τ≈5, тогда как флуктуационный вклад точно описывает реальную дисперсию.