Эволюционный кейс: на океаническом атолле изоляция популяции моллюсков длилась несколько тысяч лет; какие признаки и механизмы (генетический дрейф, естественный отбор, миграция) вы ожидаете увидеть при сравнении с материковой популяцией и как проверить это генетически
Кратко — чего ожидать и как это проверить генетически.
1) Каких биологических признаков можно ожидать на атолле посравнениюсматериковойпопуляциейпо сравнению с материковой популяциейпосравнениюсматериковойпопуляцией
Нейтральные изменения, вызванные дрейфом/эффектом основателя: снижение генетического разнообразия меньшегетерозигот,меньшеаллельногобогатстваменьше гетерозигот, меньше аллельного богатстваменьшегетерозигот,меньшеаллельногобогатства, наличие «частных» аллелей или их потеря, усиление инбридинга. Могут появиться фиксированные или почти фиксированные аллели.Адаптивные изменения под локальную среду: изменение формы/толщины раковины волноваянагрузка,абразияволновая нагрузка, абразияволноваянагрузка,абразия, окраски маскировка,УФмаскировка, УФмаскировка,УФ, размеров тела и скорости роста, сроков размножения, соле/термостойкости, поведения, устойчивости к хищникам и патогенам.Морфологические/физиологические сдвиги могут быть быстрыми несколькодесятков—сотенпоколенийнесколько десятков — сотен поколенийнесколькодесятков—сотенпоколений если сильный отбор.При периодическом миграционном обмене реже—посравнениюсматерикомреже — по сравнению с материкомреже—посравнениюсматериком ожидают смешанные признаки: в основном изоляция + низкий уровень поступления аллелей, возможно адмиксия последних лет.
2) Механизмы и какие генетические сигналы они дают
Генетический дрейф
Сигналы: пониженная гетерозиготность Hobs,HexpHobs, HexpHobs,Hexp, уменьшение аллельного богатства, высокие значения FST между популяциями, случайные изменения частот аллелей, длинные участки гомозиготности ROHROHROH при долгой изоляции/малом Ne.Тесты/метрики: Hobs/Hexp, allelic richness, private alleles count, LD-based Ne, тесты на их потери bottlenecktests—M−ratio,heterozygosity−excessbottleneck tests — M-ratio, heterozygosity-excessbottlenecktests—M−ratio,heterozygosity−excess, демографические модели fastsimcoal2,dadifastsimcoal2, dadifastsimcoal2,dadi чтобы оценить величину снижения Ne и длительность изоляции.
Естественный отбор локальнаяадаптациялокальная адаптациялокальнаяадаптация
Сигналы: локализованные участки генома с высокой дивергенцией FST−выбросыFST-выбросыFST−выбросы, пониженная локальная нуклеотидная дивергенция вокруг положительно отобранных аллелей свипысвипысвипы, отличие по частотам функциональных/кодирующих аллелей, ассоциация аллелей с экологическими переменными или адаптивными фенотипами.Тесты/методы: FST outlier BayeScan,FDIST,pcadaptBayeScan, FDIST, pcadaptBayeScan,FDIST,pcadapt, environmental association LFMM,Bayenv2LFMM, Bayenv2LFMM,Bayenv2, селекционные сканы XP−EHH,iHS,XP−CLRдляхаплотипныхданныхXP-EHH, iHS, XP-CLR для хаплотипных данныхXP−EHH,iHS,XP−CLRдляхаплотипныхданных, McDonald–Kreitman/πN/πS для белковых генов, комбинированные подходы FST+разрежениедиверсииFST + разрежение диверсииFST+разрежениедиверсии. Рекоммендуется также провести общие эксперименты reciprocaltransplant,commongardenreciprocal transplant, common gardenreciprocaltransplant,commongarden чтобы подтвердить адаптивность фенотипов.
Сигналы: снижение FST, наличие смешанных генотипов адмиксияадмиксияадмиксия, участки генома с разной историей локальныеаутлайерыадмиксиилокальные аутлайеры адмиксиилокальныеаутлайерыадмиксии, ассиметричный поток аллелей от материка к атоллу или наоборот.Тесты/методы: STRUCTURE/ADMIXTURE/PCA/DAPC для присвоения индивидов; Treemix и D-statistics ABBA−BABAABBA-BABAABBA−BABA — для детекции адмиксии; MIGRATE-n, Migrate, BayesAss, G-PhoCS для оценки скорости миграции и направления; EEMS — визуализация потоков; демографическое моделирование fastsimcoal2,dadifastsimcoal2, dadifastsimcoal2,dadi с параметризованными миграционными потоками.
3) Практический план исследования пошаговопошаговопошагово
Дизайн выборки
Количество популяций: минимум одна атолльская и несколько материковых реплик чтобыотличитьлокальнуюдифференциациюотширокойструктурычтобы отличить локальную дифференциацию от широкой структурычтобыотличитьлокальнуюдифференциациюотширокойструктуры.Количество индивидов: >=20–30 индивидуумов на популяцию желательно30+дляпопуляционнойгенетикижелательно 30+ для популяционной генетикижелательно30+дляпопуляционнойгенетики; больше нужно для точных оценок Ne и сканов на отбор.Если возможно: временные серии/исторические образцы или несколько атоллов для репликации.
Молекулярные маркеры
Для общей картины: RAD‑seq / GBS → тысячи—десятки тысяч SNP хорошийкомпромиссцена/покрытиехороший компромисс цена/покрытиехорошийкомпромиссцена/покрытие.Для глубокой истории/отборов: WGS еслибюджетпозволяетесли бюджет позволяетеслибюджетпозволяет.Дополнительно: mtDNA/микросателлиты для сравнения с предыдущими работами и материнской истории.Функциональные данные: секвенирование транскриптомов RNA‑seqRNA‑seqRNA‑seq для выявления кандидатов по экспрессии.
Анализы
Общая популяционная статистика: Hobs/Hexp, allelic richness, π, FIS, private alleles.Структура популяций: PCA, ADMIXTURE/STRUCTURE, DAPC, pairwise FST.Демография: LD‑based Ne NeEstimatorNeEstimatorNeEstimator, fastsimcoal2/dadi/∂a∂i/SMC++ для оценки времени дивергенции, величины Ne и скорости миграции; тесты на бутылочные горлышки.Скан на отбор: FST outliers, environmental association LFMM/Bayenv2LFMM/Bayenv2LFMM/Bayenv2, XP-EHH/iHS/XP-CLR при наличии WGS/фазированных данных.Гаплотипные/функциональные анализы для кандидатов: локальная дивергенция, снижение π вокруг кандидата, проверка аннотаций генов.
Экспериментальная проверка адаптации
Common garden / reciprocal transplant / лабораторные испытания на соли/температуру/волновую нагрузку.Ассоциации генотип–фенотип GWASпридостаточномsamplesizeGWAS при достаточном sample sizeGWASпридостаточномsamplesize.
4) Интерпретация: как отделить дрейф от отбора
Дрейф влияет на весь геном случайно, даёт распространённое снижение diversity и случайные FST; отбор даёт локальные, систематические сигналы узкие«пики»FST,снижениеdiversityвозлекандидатныхлокусов,согласованныеэкологическиеассоциацииузкие «пики» FST, снижение diversity возле кандидатных локусов, согласованные экологические ассоциацииузкие«пики»FST,снижениеdiversityвозлекандидатныхлокусов,согласованныеэкологическиеассоциации.Используйте нейтральные маркеры/геномное распределение: если дивергенция равномерна по геному — в основном дрейф; если разбросанная — есть кандидаты отбора.Моделирование/ABC: смоделируйте сценарии чистыйдрейф,дрейф+низкаямиграция,дрейф+отбор,периодическаямиграциячистый дрейф, дрейф+низкая миграция, дрейф+отбор, периодическая миграциячистыйдрейф,дрейф+низкаямиграция,дрейф+отбор,периодическаямиграция и сравните с наблюдаемыми данными.
5) Практические замечания и подводные камни
Малые выборки и малое число маркеров дают слабую мощность для сканов на отбор и оценки миграции.mtDNA может давать искажённую картину материнскаялиния,селективныеэффектыматеринская линия, селективные эффектыматеринскаялиния,селективныеэффекты.Повторяющиеся адаптации еслиестьдругиеостроваесли есть другие островаеслиестьдругиеострова помогают подтверждать локальную адаптацию.Гибридизация/человек могут вносить недавно привнесённые аллели — учитывайте исторические данные.
6) Примерные ожидания через «несколько тысяч лет»
Если Ne атолла низок и миграция редка: сильный дрейф → заметно пониженное разнообразие, высокие FST, возможно фиксация некоторых нейтральных аллелей; если сильный локальный отбор — обнаружите несколько геномных регионов под сильным отбором, ассоциированных с локальными условиями.Если миграция периодична редкиегостиредкие гостиредкиегости: общая структура сохранится, но будут следы адмиксии и смешанные генотипы у некоторых индивидуумов; Ne будет выше, diversity — ближе к материковому.
Если хотите, могу:
Предложить конкретный протокол выборки и минимальные цифры Nинд.,покрытиеN инд., покрытиеNинд.,покрытие под RAD‑seq или WGS для вашего бюджета;Составить список программ/команд для основных анализов примерныйрабочийпайплайнпримерный рабочий пайплайнпримерныйрабочийпайплайн.
Кратко — чего ожидать и как это проверить генетически.
1) Каких биологических признаков можно ожидать на атолле посравнениюсматериковойпопуляциейпо сравнению с материковой популяциейпосравнениюсматериковойпопуляцией
Нейтральные изменения, вызванные дрейфом/эффектом основателя: снижение генетического разнообразия меньшегетерозигот,меньшеаллельногобогатстваменьше гетерозигот, меньше аллельного богатстваменьшегетерозигот,меньшеаллельногобогатства, наличие «частных» аллелей или их потеря, усиление инбридинга. Могут появиться фиксированные или почти фиксированные аллели.Адаптивные изменения под локальную среду: изменение формы/толщины раковины волноваянагрузка,абразияволновая нагрузка, абразияволноваянагрузка,абразия, окраски маскировка,УФмаскировка, УФмаскировка,УФ, размеров тела и скорости роста, сроков размножения, соле/термостойкости, поведения, устойчивости к хищникам и патогенам.Морфологические/физиологические сдвиги могут быть быстрыми несколькодесятков—сотенпоколенийнесколько десятков — сотен поколенийнесколькодесятков—сотенпоколений если сильный отбор.При периодическом миграционном обмене реже—посравнениюсматерикомреже — по сравнению с материкомреже—посравнениюсматериком ожидают смешанные признаки: в основном изоляция + низкий уровень поступления аллелей, возможно адмиксия последних лет.2) Механизмы и какие генетические сигналы они дают
Генетический дрейф
Сигналы: пониженная гетерозиготность Hobs,HexpHobs, HexpHobs,Hexp, уменьшение аллельного богатства, высокие значения FST между популяциями, случайные изменения частот аллелей, длинные участки гомозиготности ROHROHROH при долгой изоляции/малом Ne.Тесты/метрики: Hobs/Hexp, allelic richness, private alleles count, LD-based Ne, тесты на их потери bottlenecktests—M−ratio,heterozygosity−excessbottleneck tests — M-ratio, heterozygosity-excessbottlenecktests—M−ratio,heterozygosity−excess, демографические модели fastsimcoal2,dadifastsimcoal2, dadifastsimcoal2,dadi чтобы оценить величину снижения Ne и длительность изоляции.Естественный отбор локальнаяадаптациялокальная адаптациялокальнаяадаптация
Сигналы: локализованные участки генома с высокой дивергенцией FST−выбросыFST-выбросыFST−выбросы, пониженная локальная нуклеотидная дивергенция вокруг положительно отобранных аллелей свипысвипысвипы, отличие по частотам функциональных/кодирующих аллелей, ассоциация аллелей с экологическими переменными или адаптивными фенотипами.Тесты/методы: FST outlier BayeScan,FDIST,pcadaptBayeScan, FDIST, pcadaptBayeScan,FDIST,pcadapt, environmental association LFMM,Bayenv2LFMM, Bayenv2LFMM,Bayenv2, селекционные сканы XP−EHH,iHS,XP−CLRдляхаплотипныхданныхXP-EHH, iHS, XP-CLR для хаплотипных данныхXP−EHH,iHS,XP−CLRдляхаплотипныхданных, McDonald–Kreitman/πN/πS для белковых генов, комбинированные подходы FST+разрежениедиверсииFST + разрежение диверсииFST+разрежениедиверсии. Рекоммендуется также провести общие эксперименты reciprocaltransplant,commongardenreciprocal transplant, common gardenreciprocaltransplant,commongarden чтобы подтвердить адаптивность фенотипов.Миграция генетическийпотокгенетический потокгенетическийпоток
Сигналы: снижение FST, наличие смешанных генотипов адмиксияадмиксияадмиксия, участки генома с разной историей локальныеаутлайерыадмиксиилокальные аутлайеры адмиксиилокальныеаутлайерыадмиксии, ассиметричный поток аллелей от материка к атоллу или наоборот.Тесты/методы: STRUCTURE/ADMIXTURE/PCA/DAPC для присвоения индивидов; Treemix и D-statistics ABBA−BABAABBA-BABAABBA−BABA — для детекции адмиксии; MIGRATE-n, Migrate, BayesAss, G-PhoCS для оценки скорости миграции и направления; EEMS — визуализация потоков; демографическое моделирование fastsimcoal2,dadifastsimcoal2, dadifastsimcoal2,dadi с параметризованными миграционными потоками.3) Практический план исследования пошаговопошаговопошагово
Дизайн выборки
Количество популяций: минимум одна атолльская и несколько материковых реплик чтобыотличитьлокальнуюдифференциациюотширокойструктурычтобы отличить локальную дифференциацию от широкой структурычтобыотличитьлокальнуюдифференциациюотширокойструктуры.Количество индивидов: >=20–30 индивидуумов на популяцию желательно30+дляпопуляционнойгенетикижелательно 30+ для популяционной генетикижелательно30+дляпопуляционнойгенетики; больше нужно для точных оценок Ne и сканов на отбор.Если возможно: временные серии/исторические образцы или несколько атоллов для репликации.Молекулярные маркеры
Для общей картины: RAD‑seq / GBS → тысячи—десятки тысяч SNP хорошийкомпромиссцена/покрытиехороший компромисс цена/покрытиехорошийкомпромиссцена/покрытие.Для глубокой истории/отборов: WGS еслибюджетпозволяетесли бюджет позволяетеслибюджетпозволяет.Дополнительно: mtDNA/микросателлиты для сравнения с предыдущими работами и материнской истории.Функциональные данные: секвенирование транскриптомов RNA‑seqRNA‑seqRNA‑seq для выявления кандидатов по экспрессии.Анализы
Общая популяционная статистика: Hobs/Hexp, allelic richness, π, FIS, private alleles.Структура популяций: PCA, ADMIXTURE/STRUCTURE, DAPC, pairwise FST.Демография: LD‑based Ne NeEstimatorNeEstimatorNeEstimator, fastsimcoal2/dadi/∂a∂i/SMC++ для оценки времени дивергенции, величины Ne и скорости миграции; тесты на бутылочные горлышки.Скан на отбор: FST outliers, environmental association LFMM/Bayenv2LFMM/Bayenv2LFMM/Bayenv2, XP-EHH/iHS/XP-CLR при наличии WGS/фазированных данных.Гаплотипные/функциональные анализы для кандидатов: локальная дивергенция, снижение π вокруг кандидата, проверка аннотаций генов.Экспериментальная проверка адаптации
Common garden / reciprocal transplant / лабораторные испытания на соли/температуру/волновую нагрузку.Ассоциации генотип–фенотип GWASпридостаточномsamplesizeGWAS при достаточном sample sizeGWASпридостаточномsamplesize.4) Интерпретация: как отделить дрейф от отбора
Дрейф влияет на весь геном случайно, даёт распространённое снижение diversity и случайные FST; отбор даёт локальные, систематические сигналы узкие«пики»FST,снижениеdiversityвозлекандидатныхлокусов,согласованныеэкологическиеассоциацииузкие «пики» FST, снижение diversity возле кандидатных локусов, согласованные экологические ассоциацииузкие«пики»FST,снижениеdiversityвозлекандидатныхлокусов,согласованныеэкологическиеассоциации.Используйте нейтральные маркеры/геномное распределение: если дивергенция равномерна по геному — в основном дрейф; если разбросанная — есть кандидаты отбора.Моделирование/ABC: смоделируйте сценарии чистыйдрейф,дрейф+низкаямиграция,дрейф+отбор,периодическаямиграциячистый дрейф, дрейф+низкая миграция, дрейф+отбор, периодическая миграциячистыйдрейф,дрейф+низкаямиграция,дрейф+отбор,периодическаямиграция и сравните с наблюдаемыми данными.5) Практические замечания и подводные камни
Малые выборки и малое число маркеров дают слабую мощность для сканов на отбор и оценки миграции.mtDNA может давать искажённую картину материнскаялиния,селективныеэффектыматеринская линия, селективные эффектыматеринскаялиния,селективныеэффекты.Повторяющиеся адаптации еслиестьдругиеостроваесли есть другие островаеслиестьдругиеострова помогают подтверждать локальную адаптацию.Гибридизация/человек могут вносить недавно привнесённые аллели — учитывайте исторические данные.6) Примерные ожидания через «несколько тысяч лет»
Если Ne атолла низок и миграция редка: сильный дрейф → заметно пониженное разнообразие, высокие FST, возможно фиксация некоторых нейтральных аллелей; если сильный локальный отбор — обнаружите несколько геномных регионов под сильным отбором, ассоциированных с локальными условиями.Если миграция периодична редкиегостиредкие гостиредкиегости: общая структура сохранится, но будут следы адмиксии и смешанные генотипы у некоторых индивидуумов; Ne будет выше, diversity — ближе к материковому.Если хотите, могу:
Предложить конкретный протокол выборки и минимальные цифры Nинд.,покрытиеN инд., покрытиеNинд.,покрытие под RAD‑seq или WGS для вашего бюджета;Составить список программ/команд для основных анализов примерныйрабочийпайплайнпримерный рабочий пайплайнпримерныйрабочийпайплайн.