• A
  • A
  • A
  • АБB
  • АБB
  • АБB
  • А
  • А
  • А
  • А
  • А
Обычная версия сайта

Московский институт электроники
и математики им. А.Н. Тихонова

"Вычислительные среды": доклад "О распределении длины трактов адаптивной интрогрессии"

2 октября 2019 г. прошел очередной семинар «Суперкомпьютерное моделирование в науке и инженерии, или Вычислительные среды». Владимир Львович Щур, к.ф.-м.н., доцент департамента прикладной математики МИЭМ, выступил с докладом "О распределении длины трактов  адаптивной интрогрессии" ("On the distribution of tract lengths during adaptive introgression").

Поскольку изучаемый вопрос находится в области геномики и многие слушатели мало знакомы с этой областью, начало доклада было посвящено разбору основных терминов. Тракты, о которых идет речь, это участок генома, т.е. некоторая последовательность нуклеотидов, имеющая иное происхождение, чем остальная часть генома. Адаптивная интрогрессия – это сценарий, по которому популяция приспосабливается к окружающей среде. Модель адаптивной интрогрессии заключается в следующем. Имеется две популяции, одна из которых вырабатывает адаптацию, а вторая – является донором для первой. Во второй популяции есть аллель (вариант гена), которая улучшает приспособленность носителя к окружающей среде, а в первой популяции такой аллели нет. При интрогрессии часть популяции-донора «пульсово», то есть одномоментно, попадает в первую популяцию. Задача исследования заключается в том, чтобы определить вероятность того, что тракт генома, содержащий адаптационную аллель, пришел из популяции-донора, а также предсказать длину этого тракта.

Известно, что ожидаемая траектория частоты выбранной аллели может быть приближена логистической функцией w(t)=1/(1+exp(-st/2)). При помощи компьютерного моделирования было показано, что распределение длины трактов можно вычислять вдоль ожидаемой траектории частоты аллели, таким образом пренебрегая стохастическим поведением. Так, относительная ошибка при оценке математического ожидания длины тракта составляет не более 5%. Одним из параметров моделирования являлась сила естественного отбора, поскольку естественный отбор сильно влияет на распределение длин путей. Положительный более сильный отбор будет увеличивать ожидаемую длину интрогрессировавшего тракта. Однако, возможно, неинтуитивным образом, при условии фиксированного значения частоты аллели в наблюдаемой популяции, более сильный отбор на самом деле приведет к более короткой ожидаемой длине тракта.

Поскольку наличие отбора влияет на распределение длин трактов, это явление может иметь значение при датировании интрогрессии. В конце доклада был разобраны пример сценария интрогрессии гена EPAS1 от денисовцев в предков тибетцев, позволившего выработать адаптацию к высокогорному существованию. Было показано, что оценка времени в предположении нейтрального сценария (без естественного отбора) вместо адаптивного приводит к ошибке в 20% при датировании времени интрогрессии.