Представьте популяцию, в которой частота аллеля резко изменилась за несколько поколений без очевидного отбора; какие демографические и генетические процессы (дрейф, миграция, бутылочное горлышко, целенаправленное размножение) могут объяснить наблюдения и как их отличить эмпирически
Коротко: возможные процессы — генетический дрейф, миграция (внесение аллеля извне), бутылочное горлышко (резкое сокращение NeN_eNe) и целенаправленное размножение/отбор. Как отличить — по масштабу (локус vs геном), темпу изменений и характерным геномным сигналам. Ниже — признаки каждого процесса и эмпирические тесты. 1) Генетический дрейф - Механизм: случайные флуктуации частот в конечной популяции; вероятность фиксации равна начальной частоте. - Ожидаемая дисперсия смены частоты за одну генерацию: Var(Δp)=p(1−p)2Ne\mathrm{Var}(\Delta p)=\dfrac{p(1-p)}{2N_e}Var(Δp)=2Nep(1−p). - Через ttt поколений: Var(pt)≈p0(1−p0)(1−(1−12Ne)t)\mathrm{Var}(p_t)\approx p_0(1-p_0)\left(1-\left(1-\dfrac{1}{2N_e}\right)^t\right)Var(pt)≈p0(1−p0)(1−(1−2Ne1)t). - Признаки: изменения случайны по всему геному, одинаковая амплитуда для многих негарантированных локусов; нет сильных локус-специфичных паттернов LD или длинных гаплотипов. - Тесты: оценить эффективный размер NeN_eNe (темповой метод по временным сэмплам), сравнить наблюдаемую дисперсию Δp\Delta pΔp с ожидаемой при дрейфе; симуляции Wright–Fisher. 2) Миграция (генетический поток) - Модель простейшая за одну генерацию: p′=(1−m)p+mpmp'=(1-m)p + m p_mp′=(1−m)p+mpm, где mmm — доля мигрантов, pmp_mpm — частота у источника. - Признаки: сдвиг частот направленный и коррелирован между локусами в сторону частот источника; появление блоков аутосомного происхождения (адмикция), новые аллели/гаплотипы с характерной длиной в зависимости от времени миграции. - Тесты: PCA/ADMIXTURE/treemix для обнаружения компонент происхождения; локусные FST между до/после и со «источником»; анализ длин интродуцированных гаплотипов (local ancestry); асимметрия смен частот между локусами (согласованность с pmp_mpm). 3) Бутылочное горлышко (резкое падение NeN_eNe) - Механизм: внезапное уменьшение численности повышает дрейф и теряет редкие аллели. - Признаки: геномно-широкое уменьшение диверситета (гетерозиготности), уменьшение числа полиморфных сайтов, повышенное LD и длинные участки гомозиготности (ROH), усиление вариативности частот между локусами. - Формула убыли гетерозиготности: Ht+1=Ht(1−12Ne)H_{t+1}=H_t\left(1-\dfrac{1}{2N_e}\right)Ht+1=Ht(1−2Ne1). - Тесты: сравнить SSS (число segregating sites), π \piπ и HHH до/после; LD-декей; распределение длин ROH; SFS — потеря редких вариантов; демографические методы (SMC++, stairway plot, dadi) для восстановления истории Ne(t)N_e(t)Ne(t). 4) Целенаправленное размножение / искусственный отбор - Механизм: люди/скрытый отбор целенаправленно увеличивают частоту конкретных аллелей; действует как сильный отбор. - Признаки: быстрый, локус-специфичный сдвиг частоты у функциональных/связанных локусов; сигнатуры отжима (selective sweep): падение диверситета вокруг локуса, длинные гаплотипы высокого частотного класса, сильный локусный FST. - Тесты: сканы на отбор (iHS, XP-EHH, CLR, PBS, Fay & Wu, Tajima's D), поиск длинных однородных гаплотипов, проверка родословных/помещённых данных о селекции, связь с фенотипами и известной селекционной программой. Как эмпирически отличать в практической последовательности 1. Соберите широкомасштабные геномные данные (много локусов), желательно временные сэмплы. 2. Оцените геномно-широкие метрики: π \piπ, HHH, число segregating sites, LD, ROH. Глобальные изменения → бутылочное горлышко/дрейф; локус-специфичные → миграция или отбор. 3. Оцените темповую дисперсию частот и NeN_eNe методом временных сэмплов; если наблюдаемая дисперсия >> ожидания при предполагаемом NeN_eNe, подумать об миграции или отборе. 4. Выполните PCA/ADMIXTURE и локальную анценстри-аналитику — сильная адмикция укажет на миграцию. 5. Сделайте сканы на отбор и изучите гаплотипную структуру вокруг интересующих локусов — сигналы sweep говорят за целенаправленный отбор/размножение. 6. Моделирование/ABC: смоделируйте альтернативные сценарии (дрейф, бутылочное горло, миграция, отбор) и сравните с данными по множеству статистик (SFS, LD, ROH, FST, temporal variance). Краткие эвристики - Если изменение быстро и локус-специфично + сильные гаплотипные сигналы → целенаправленное размножение/отбор. - Если изменение быстро и совпадает с появлением чужих гаплотипов/компонент → миграция. - Если снижение диверситета, рост LD и ROH по всему геному → бутылочное горлышко (усиленный дрейф). - Если изменения случайны по геному и соответствуют ожидаемой дисперсии при малом NeN_eNe → дрейф. Если хотите, могу предложить конкретный аналитический план и список команд/программ (e.g. PLINK, ADMIXTURE, PCAngsd, SMC++, dadi, selscan) для ваших данных.
1) Генетический дрейф
- Механизм: случайные флуктуации частот в конечной популяции; вероятность фиксации равна начальной частоте.
- Ожидаемая дисперсия смены частоты за одну генерацию: Var(Δp)=p(1−p)2Ne\mathrm{Var}(\Delta p)=\dfrac{p(1-p)}{2N_e}Var(Δp)=2Ne p(1−p) .
- Через ttt поколений: Var(pt)≈p0(1−p0)(1−(1−12Ne)t)\mathrm{Var}(p_t)\approx p_0(1-p_0)\left(1-\left(1-\dfrac{1}{2N_e}\right)^t\right)Var(pt )≈p0 (1−p0 )(1−(1−2Ne 1 )t).
- Признаки: изменения случайны по всему геному, одинаковая амплитуда для многих негарантированных локусов; нет сильных локус-специфичных паттернов LD или длинных гаплотипов.
- Тесты: оценить эффективный размер NeN_eNe (темповой метод по временным сэмплам), сравнить наблюдаемую дисперсию Δp\Delta pΔp с ожидаемой при дрейфе; симуляции Wright–Fisher.
2) Миграция (генетический поток)
- Модель простейшая за одну генерацию: p′=(1−m)p+mpmp'=(1-m)p + m p_mp′=(1−m)p+mpm , где mmm — доля мигрантов, pmp_mpm — частота у источника.
- Признаки: сдвиг частот направленный и коррелирован между локусами в сторону частот источника; появление блоков аутосомного происхождения (адмикция), новые аллели/гаплотипы с характерной длиной в зависимости от времени миграции.
- Тесты: PCA/ADMIXTURE/treemix для обнаружения компонент происхождения; локусные FST между до/после и со «источником»; анализ длин интродуцированных гаплотипов (local ancestry); асимметрия смен частот между локусами (согласованность с pmp_mpm ).
3) Бутылочное горлышко (резкое падение NeN_eNe )
- Механизм: внезапное уменьшение численности повышает дрейф и теряет редкие аллели.
- Признаки: геномно-широкое уменьшение диверситета (гетерозиготности), уменьшение числа полиморфных сайтов, повышенное LD и длинные участки гомозиготности (ROH), усиление вариативности частот между локусами.
- Формула убыли гетерозиготности: Ht+1=Ht(1−12Ne)H_{t+1}=H_t\left(1-\dfrac{1}{2N_e}\right)Ht+1 =Ht (1−2Ne 1 ).
- Тесты: сравнить SSS (число segregating sites), π \piπ и HHH до/после; LD-декей; распределение длин ROH; SFS — потеря редких вариантов; демографические методы (SMC++, stairway plot, dadi) для восстановления истории Ne(t)N_e(t)Ne (t).
4) Целенаправленное размножение / искусственный отбор
- Механизм: люди/скрытый отбор целенаправленно увеличивают частоту конкретных аллелей; действует как сильный отбор.
- Признаки: быстрый, локус-специфичный сдвиг частоты у функциональных/связанных локусов; сигнатуры отжима (selective sweep): падение диверситета вокруг локуса, длинные гаплотипы высокого частотного класса, сильный локусный FST.
- Тесты: сканы на отбор (iHS, XP-EHH, CLR, PBS, Fay & Wu, Tajima's D), поиск длинных однородных гаплотипов, проверка родословных/помещённых данных о селекции, связь с фенотипами и известной селекционной программой.
Как эмпирически отличать в практической последовательности
1. Соберите широкомасштабные геномные данные (много локусов), желательно временные сэмплы.
2. Оцените геномно-широкие метрики: π \piπ, HHH, число segregating sites, LD, ROH. Глобальные изменения → бутылочное горлышко/дрейф; локус-специфичные → миграция или отбор.
3. Оцените темповую дисперсию частот и NeN_eNe методом временных сэмплов; если наблюдаемая дисперсия >> ожидания при предполагаемом NeN_eNe , подумать об миграции или отборе.
4. Выполните PCA/ADMIXTURE и локальную анценстри-аналитику — сильная адмикция укажет на миграцию.
5. Сделайте сканы на отбор и изучите гаплотипную структуру вокруг интересующих локусов — сигналы sweep говорят за целенаправленный отбор/размножение.
6. Моделирование/ABC: смоделируйте альтернативные сценарии (дрейф, бутылочное горло, миграция, отбор) и сравните с данными по множеству статистик (SFS, LD, ROH, FST, temporal variance).
Краткие эвристики
- Если изменение быстро и локус-специфично + сильные гаплотипные сигналы → целенаправленное размножение/отбор.
- Если изменение быстро и совпадает с появлением чужих гаплотипов/компонент → миграция.
- Если снижение диверситета, рост LD и ROH по всему геному → бутылочное горлышко (усиленный дрейф).
- Если изменения случайны по геному и соответствуют ожидаемой дисперсии при малом NeN_eNe → дрейф.
Если хотите, могу предложить конкретный аналитический план и список команд/программ (e.g. PLINK, ADMIXTURE, PCAngsd, SMC++, dadi, selscan) для ваших данных.