Полиморфизм генов FADS и история формирования коренного населения Сибири

Полный текст

Кластер FADS-генов (FADS1, FADS2 и FADS3), кодирующих десатуразы жирных кислот, представляет большой интерес в плане изучения генетических механизмов адаптации к условиям природной среды и особенностям питания различных народов [1‒6]. В области расположения FADS-генов (11q12-13.1) находятся два гаплоблока: первый определяется вариантами полиморфизма в гене FADS1 и первой части гена FADS2, второй – вариантами полиморфизма оставшейся последовательности гена FADS2 [1]. Установлено также, что первый гаплоблок представлен двумя основными гаплотипами – A и D, которые кодируют ферменты соответственно с пониженной и повышенной десатуразной активностью.

Гаплотип A довольно редко встречается в популяциях Африки (примерно 1‒3%) и Южной Азии (примерно 8%), а также Европы (12‒16%) (по сведениям базы данных dbSNP (https://www.ncbi.nlm.nih.gov/snp/)). Однако у населения Восточной Азии и Америки частота гаплотипа A достигает высоких частот (более 50%). У гренландских эскимосов наблюдается почти полная фиксация этого гаплотипа в популяции – 98% по данным работы [2]. Предполагается, что распространение менее активного в плане десатурации жирных кислот гаплотипа A началось еще со времен заселения Евразии и Америки верхнепалеолитическими людьми, потреблявшими богатую липидами и белками пищу [7–9]. Однако позже в неолите, по мере появления технологий сельского хозяйства более распространенными в некоторых регионах мира (Европе, Южной Азии) стали D-гаплотипы, позволяющие с большей скоростью синтезировать полиненасыщенные жирные кислоты (ПНЖК) из липидов растений [1, 8, 9]. Между тем по палеогеномным данным известно, что у донеолитических европейцев был распространен в основном гаплотип A [10]. Таким образом, частота гаплотипа D в Европе увеличилась от менее 10% 10 тыс. лет тому назад до 60‒75% в настоящее время [7]. Предполагается, что столь резкое изменение частоты этого гаплотипа у европейцев связано с действием положительного отбора [8‒10]. Более того, недавно установлено, что действие отбора на локусы FADS-генов у европейцев продолжается до сих пор – согласно результатам анализа данных UK Biobank, частота гаплотипа D у британцев европейского происхождения продолжает расти на 0.009% в год [11]. Одной из причин этого может быть связь с репродуктивным успехом, поскольку рост частоты гаплотипа D ассоциируется с увеличением коэффициента рождаемости.

Результаты популяционно-генетических исследований показали, что с наиболее высокими частотами гаплотип А распространен среди коренного населения Сибири и Америки [4‒6, 12‒14]. Причем в сибирских популяциях частота гаплотипа A увеличивается в направлении с юга на север [5, 14]. У американских индейцев частота этого гаплотипа составляет более 90% как на севере, так и на юге Америки [4]. Предполагается, что высокая распространенность гаплотипа A у америндов связана с отбором наиболее оптимальных вариантов FADS-полиморфизма еще у населения Берингии, предкового по отношению к американским индейцам, в связи с необходимостью адаптироваться к холоду и ограниченным пищевым ресурсам, предоставленным Арктикой [2, 4]. Другая гипотеза, основанная на палеогеномных данных, предполагает, что отбор гаплотипа A произошел еще раньше – во время выхода первых групп человека из Африки и расселения на территориях Евразии [8]. Поэтому к моменту формирования берингийской популяции (20–15 тыс. лет тому назад) вполне вероятно, что альтернативы гаплотипу A просто не было, однако отбор мог происходить на уровне его субгаплотипов.

В популяционно-генетическом исследовании гренландских эскимосов был выявлен сильный сигнал отбора в области FADS-генов, а также оказалось, что у эскимосов с частотой 98% распространен единственный гаплотип, сформированный вариантами полиморфизма FADS-генов, связанными с рядом метаболических и антропометрических фенотипов (например, ростом, распределением жира по телу, композицией жирных кислот в мембранах эритроцитов) [2]. Этот гаплотип представлен вариантами полиморфизма, маркирующими гаплотип A (rs174570-T гена FADS2), а также вариантами полиморфизма, выходящими за пределы области этого гаплоблока (rs74771917-T гена FADS2 и rs7115739-T гена FADS3). Поскольку информация о полиморфизме этих локусов в популяциях коренного населения Сибири очень фрагментарна, в настоящей работе нами исследован полиморфизм локусов rs174570, rs74771917 и rs7115739 в популяциях северной и южной частей Сибири, а также получены сведения о распространении гаплотипов для указанных локусов FADS-генов у современного и древнего населения Сибири.

МАТЕРИАЛЫ И МЕТОДЫ

Исследованы выборки коренного населения Сибири: чукчи (N = 14) из Чукотского автономного округа (г. Анадырь), коряки (N = 96) и эвены (N = 96) из Северо-Эвенского района Магаданской области, эвенки (N = 56) из различных поселков Эвенкийского района Красноярского края, буряты (N = 123) из различных районов Бурятии, алтайцы (N = 125) из различных районов Республики Алтай, тувинцы (N = 84) из различных районов Республики Тыва.

Выделение и очистку геномной ДНК проводили как описано нами ранее [6]. Генотипирование локусов rs174570, rs74771917 и rs7115739 выполняли с помощью секвенирования ДНК по Сэнгеру с использованием набора для циклического секвенирования ДНК BigDye Terminator (Applied Biosystems v 3.1) и генетического анализатора ABI Prism 3500xL (Applied Biosystems). Участок ДНК, включающий в свой состав полиморфный локус rs174570 (ген FADS2, позиция 61597212 хромосомы 11), амплифицировали с помощью пары олигонуклеотидных праймеров FA1 (5ꞌ-TCCTGCCAAGAACACCCTTC-3ꞌ) и FA2 (5ꞌ-CACCACCTGTCCATTGTCCA -3ꞌ). Для генотипирования локуса rs74771917 (ген FADS2, позиция 61627960 хромосомы 11) использовали праймеры FA5 (5ꞌ-AGAACCGCCCATGACAACAA-3ꞌ) и FA6 (5ꞌ-GTTTCGCCCATGTTGCTCAG-3ꞌ), для локуса rs7115739 (ген FADS3, позиция 61641717 хромосомы 11) применяли праймеры FA7 (5ꞌ-AACACCCCTGTATCCCCCAT-3ꞌ) и FA8 (5ꞌ-CGCCTTTGTACACGCTGTTC-3ꞌ). Нуклеотидные последовательности праймеров подобраны с помощью программы Primer3 [15]. Нумерация нуклеотидов приводится согласно референтной последовательности генома человека GRCh37.p13 (hg19).

Для анализа нуклеотидных последовательностей участков, включающих локусы rs174570, rs74771917 и rs7115739, использовали пакет программ MEGA5 [16]. Распределение аллелей, генотипов и гаплотипов, гетерозиготность и генетическую дифференциацию популяций (в значениях F_st) исследовали с помощью пакета программ Arlequin 3.5 [17]. Отклонение от равновесия Харди‒Вайнберга в популяциях и различия в частоте аллелей и генотипов оценивали с помощью точного теста Фишера. Частоты гаплотипов определяли с помощью EM-алгоритма. Неравновесие по сцеплению между парами полиморфных локусов оценивали с помощью коэффициентов D’ и r² (пакет программ Arlequin 3.5).

Для анализа распределения генотипов и аллелей локуса rs174570 использовали опубликованные ранее результаты мультилокусного генотипирования в популяциях эскимосов, якутов и шорцев [18]. Частоты аллелей для локусов rs174570, rs74771917 и rs7115739 в различных популяциях мира определяли с помощью баз генетических данных dbSNP (https://www.ncbi.nlm.nih.gov/snp/) и ALFRED (https://alfred.med.yale.edu/alfred/). Сведения о полиморфизме этих локусов у древних индивидов взяты из базы данных Allen Ancient DNA Resource (https://reich.hms.harvard.edu/). Значения генных потоков (долю привнесенных аллелей (M) и интенсивность генного потока (m)) рассчитаны по формулам (10.3) и (10.4) работы [19].

Для визуализации распределения частоты аллеля rs174570-T в географическом пространстве использовали пакет программ Surfer (Golden Software, LLC) (www.goldensoftware.com) и геоинформационную систему QGIS (https://qgis.org/).

РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

Полиморфизм локуса rs174570 гена FADS2 в современных популяциях Сибири

Результаты исследования полиморфизма локуса rs174570 в популяциях коренного населения Сибири показали, что частота аллеля rs174570-T и генотипа TT почти в два раза выше в популяциях Северо-Восточной Сибири в сравнении с югом Сибири (табл. 1). У эскимосов, чукчей и коряков частота варианта rs174570-T достигает почти 90%, что сопоставимо с аналогичными данными для эскимосов Аляски, Канады и Гренландии [2, 12, 20]. Для распределений генотипов локуса rs174570 в исследованных популяциях Сибири отклонений от равновесия Харди‒Вайнберга не обнаружено (табл. 1).

 

Таблица 1. Частоты генотипов и аллелей локуса rs174570 гена FADS2 в популяциях Сибири

Популяция (N)	Генотипы			Аллели		He	P
	TT	CT	CC	T	C
Эскимосы (19)1	0.789	0.211	0	0.895	0.105	0.194	1.0
Чукчи (14)2	0.714	0.214	0.072	0.821	0.179	0.304	0.35
Коряки (96)2	0.771	0.208	0.021	0.875	0.125	0.22	0.64
Эвены (96)2	0.646	0.302	0.052	0.797	0.203	0.325	0.53
Эвенки (56)2	0.625	0.286	0.089	0.768	0.232	0.36	0.14
Якуты (22)1	0.546	0.409	0.045	0.75	0.25	0.384	1.0
Буряты (123)2	0.22	0.536	0.244	0.488	0.512	0.502	0.47
Тувинцы (84)2	0.226	0.512	0.262	0.482	0.518	0.502	1.0
Алтайцы (125)2	0.216	0.504	0.28	0.468	0.532	0.5	1.0
Шорцы (24)1	0.083	0.584	0.333	0.375	0.625	0.479	0.39

Примечание. N – размер выборки, He – ожидаемая гетерозиготность, P – статистическая значимость отклонения от равновесия Харди‒Вайнберга (достоверно при P < 0.05). 1 – по данным работы [18]; 2 – настоящая работа.

 

Анализ межпопуляционной генетической дифференциации (в виде F_st-значений) показал, что по распределению аллелей локуса rs174570 популяции северной части Сибири (эскимосы, чукчи, коряки, эвены, эвенки, якуты) статистически значимо отличаются от популяций Южной Сибири (буряты, алтайцы, шорцы, тувинцы) (табл. 1 доп. материалов). Результаты пространственной интерполяции значений частоты аллеля rs174570-T в популяциях Северной Евразии также показали, что частота этого варианта гена FADS2 постепенно повышается к северу, достигая максимальных значений у коренного населения Чукотки (рис. 1; табл. 2 доп. материалов).

 

Таблица 2. Частота варианта rs174570-T гена FADS2 у древнего населения Сибири

Регион (N)	Координаты	Возраст (тыс. лет)	Частота варианта rs174570-T
Северо-Восточная Сибирь (26)	60.0°‒66.2° с.ш., 132.0° в.д.–169.8° з.д.	0.6‒9.8	1.0
Дальний Восток (Приморье) (11)	42.8°‒44.5° с.ш., 131.0°‒135.0° в.д.	5.6‒7.7	1.0
Юго-Восточная Сибирь (55)	50.0°‒57.5° с.ш., 102.0°‒115.0° в.д.	0.6‒8.8	0.873
Юго-Западная Сибирь (65)	48.8°‒56.2° с.ш., 83.9°‒95.8° в.д.	1.1‒5.1	0.554

 

Рис. 1. Распространение варианта rs174570-T у населения Сибири. Белыми кружками на карте отмечены опорные точки (исследованные популяции). Сбоку показана шкала интервалов значений частоты варианта полиморфизма (в %).

 

Полиморфизм локуса rs174570 гена FADS2 в древних популяциях Сибири

Анализ результатов генотипирования локуса rs174570, представленных в базе данных Allen Ancient DNA Resource, показывает, что в восточной части Евразии вариант rs174570-T издревле был самым частым (табл. 3 доп. материалов). Вариант rs174570-C зарегистрирован только в верхнепалеолитических образцах из Янской стоянки (~32 тыс. лет), которые, как установлено ранее [21], проявляли генетическое сходство с древним населением Западной Евразии, у которого уже 38‒33 тыс. лет тому назад отмечались C-аллели по локусу rs174570 (стоянки Костёнки, Сунгирь). У древних жителей Северо-Восточной Сибири и Приморья в широком диапазоне времени присутствовал исключительно вариант rs174570-T, а на юго-востоке Сибири (территория примерно соответствует этническому ареалу бурят) частота этого варианта составляла 87.3%, причем первые случаи появления аллеля C наблюдаются примерно 5.5 тыс. лет назад (табл. 2). У древнего населения Юго-Западной Сибири (в участке, примерно соответствующем этническим ареалам современных алтайцев и тувинцев) частота варианта rs174570-T составляла 55.4% (табл. 2). Таким образом, результаты показывают, что в Восточной Сибири в древности преобладал вариант rs174570-T, а появление C-аллелей у населения Байкальского региона, по всей видимости, связано с проникновением на восток Сибири, начиная с эпохи энеолита, племен афанасьевской культуры, у которых частота варианта rs174570-С достигала 40% (по данным для афанасьевцев Алтая [22]).

Между тем на северо-востоке Сибири появление С-аллелей регистрируется поздно (не раньше 600 лет назад) и связано, по всей видимости, с периодом освоения Северо-Восточной Сибири восточноевропейским по происхождению населением, начиная с XVII в. и особенно интенсивно в советское время [23]. По данным базы ALFRED для восточных европейцев (как и русских, в частности) характерна высокая частота варианта rs174570-С – в среднем, 75%. Таким образом, основываясь на данных о распространенности этого варианта в современных и древних популяциях Северо-Восточной Сибири (у эскимосов, чукчей, коряков и их предков), а также у пришлого восточноевропейского населения, принявшего участие в метисации с аборигенами Крайнего Севера, доля варианта rs174570-С у последних оценивается в 16.5%, а интенсивность генного потока составляет 1.5% на поколение. Вполне вероятно, что интенсивность генного потока еще выше, если считать, что исходная точка смешения коренных народов Северо-Восточной Сибири с пришлыми группами приходится лишь на XX в. [24]. Тогда она увеличивается до 4.4% на поколение. Для современного коренного населения Байкальского и Алтае-Саянского регионов оценить поток генов еще сложнее, поскольку приток С-аллелей осуществлялся в широком диапазоне времени и из разных популяционных источников. Необходимо отметить очень высокий прирост частоты варианта rs174570-С в Байкальском регионе – от менее 15% в промежутке от 8.8 до 0.6 тыс. лет назад до более 50% у современных бурят (табл. 1, 2).

Распределение гаплотипов по локусам rs174570, rs74771917 и rs7115739 в популяциях Сибири

Анализ полиморфизма локусов rs74771917 и rs7115739 показал отсутствие значимых различий в распределении аллелей в исследованных выборках чукчей, коряков, эвенов и бурят (табл. 4, 5 доп. материалов). Частота производных аллелей rs74771917-T и rs7115739-T в сибирских популяциях (20‒40%) находится на уровне частот, характерных для населения Восточной Азии (в среднем 25‒29% по данным dbSNP). Поэтому более информативным для целей исследования представляется анализ распределения гаплотипов по локусам rs174570, rs74771917 и rs7115739, так как ранее было обнаружено, например, что у гренландских эскимосов частота гаплотипа TTT (как комбинации аллельных вариантов для указанных выше локусов) составляет почти 100% [2]. Анализ неравновесия по сцеплению между локусами rs174570, rs74771917 и rs7115739 в исследованных нами популяциях Сибири показал высокий уровень сцепления между вариантами полиморфизма, хотя наиболее значимые величины коэффициентов D’ и r² наблюдались во всех популяциях для локусов rs74771917 и rs7115739 (табл. 6 доп. материалов).

В табл. 3 показано распределение гаплотипов по локусам rs174570, rs74771917 и rs7115739 у чукчей, коряков, эвенов и бурят. По распределению гаплотипов буряты статистически значимо отличаются от всех остальных, более северных выборок. У бурят более чем в два раза выше частота гаплотипа CCG и ниже частота гаплотипов TCG и TTT. У чукчей частота гаплотипа TTT достигает почти 40%, однако в сравнении с таковой для гренландских эскимосов (98%) она оказывается довольно умеренной, а у коряков частота гаплотипа TTT даже ниже, чем у эвенов. В табл. 4 показано распределение гаплотипов по локусам rs174570, rs74771917 и rs7115739 у более древнего населения Сибири и прилегающих территорий Евразии и Америки. Как видно, в прошлом гаплотип TTT был самым частым (более 60%) как на северо-востоке Азии, так и в Америке, в то время как на юго-востоке Сибири (в Байкальском регионе) его частота почти не изменилась со временем. При сравнении современного и древнего населения Сибири заметно, что на северо-востоке Сибири уменьшилась частота гаплотипа TTT и выросла частота гаплотипа TCG, а также произошел приток аллеля rs174570-C в виде гаплотипа CCG (рис. 1,а доп. материалов. В Байкальском регионе в последние почти 2 тыс. лет наблюдались высокий прирост частоты гаплотипа CCG и снижение частоты гаплотипа TCG (рис. 1,б доп. материалов). Между тем у американских индейцев распространенность гаплотипа TTT осталась на высоком уровне (65‒75%); в древних и современных популяциях сохраняется спектр гаплотипов, связанных с аллелем rs174570-T, хотя у современных индейских народов также наблюдается небольшой приток C-ассоциированных гаплотипов (рис. 2 доп. материалов).

 

Таблица 3. Частота гаплотипов по локусам rs174570, rs74771917 и rs7115739 в современных популяциях Сибири

Гаплотипы	Чукчи (N = 14)	Коряки (N = 96)	Эвены (N = 93)	Буряты (N = 112)
TCG	0.429	0.641	0.484	0.293
CCG	0.178	0.125	0.199	0.465
TTT	0.393	0.234	0.263	0.142
TCT	0	0	0.043	0.015
CCT	0	0	0.0055	0.048
TTG	0	0	0.0055	0.015
CTT	0	0	0	0.022

 

Таблица 4. Частота гаплотипов по локусам rs174570, rs74771917 и rs7115739 у древнего населения Сибири и Америки

Гаплотип	Юго-Восточная Сибирь, 1.9‒8.8 тыс. лет назад (N = 39)	Северо-Восточная Сибирь, 0.6‒9.8 тыс. лет назад (N = 22)	Америка, 0.3‒3.9 тыс. лет назад (N = 14)
TTT	18.0	63.7	64.3
TTG	12.8	13.6	7.1
TCG	48.7	13.6	14.3
TCT	5.1	9.1	14.3
CTG	7.7	0	0
CCT	2.6	0	0
CCG	5.1	0	0

 

Поскольку гаплотип TTT с высокой частотой был распространен в прошлом среди населения северо-восточной части Сибири и Америки, интересным представляется вопрос об истоках этого гаплотипа. Анализ данных о распределении гаплотипов по локусам rs174570, rs74771917 и rs7115739 в наиболее древних антропологических материалах (в диапазоне от 10 до 45 тыс. лет) показывает, что единичные случаи появления гаплотипа TTT наблюдаются у верхнепалеолитических европейцев (34 и 45 тыс. лет назад) (табл. 7 доп. материалов). Однако восточноазиатским источником этого гаплотипа, вероятнее всего, является Приамурье, где гаплотип TTT присутствовал еще 14‒19 тыс. лет тому назад (табл. 7 доп. материалов). Более того, результаты палеогеномных исследований показали, что древнее население Приамурья стало, по всей видимости, основным генетическим источником в процессах заселения северо-востока Азии и формирования населения Берингии [25].

Таким образом, результаты проведенного исследования свидетельствуют о том, что в современных популяциях коренного населения Сибири наблюдается увеличение частоты варианта rs174570-T, маркирующего гаплотип A, в направлении с юга на север Сибири – от около 50% у бурят, алтайцев и тувинцев до более 80% у эскимосов, чукчей и коряков. Однако у древнего населения восточной части Сибири частота варианта rs174570-T была еще выше (в диапазоне от более 80 до 100%), хотя на юго-западе Сибири частота этого варианта почти не изменилась на протяжении 5 тыс. лет. Таким образом, основные изменения генофондов населения Восточной Сибири обусловлены взаимодействием с популяциями, для которых характерны высокие частоты варианта rs174570-C. На северо-востоке Сибири это пришлое восточноевропейское по происхождению население, которое примерно с XVII в. начали осваивать Северную Сибирь, но основной период межэтнических контактов приходится, по всей видимости, на советское время [23, 24]. В результате частота варианта rs174570-C у современных эскимосов, чукчей и коряков достигает в среднем 16.5%. Насколько это полезное приобретение и способствует ли оно увеличению репродуктивного успеха (как предполагается в случае европейцев [11]), пока неизвестно. Примечательно, что у бурят также наблюдаются повышенные в сравнении с древним населением Байкальского региона частоты варианта rs174570-C, однако в этом случае, по всей видимости, прирост частоты этого аллеля связан с более широким кругом контактов с индо-европейскими народами, начиная с эпохи энеолита и бронзового века.

Интересны также результаты анализа распределения гаплотипов, сформированных вариантами полиморфизма локусов rs174570, rs74771917 и rs7115739 у современного и древнего населения Восточной Азии и Америки. Как на северо-востоке Сибири, так и в Байкальском регионе приток варианта rs174570-C в популяции связан лишь с одним гаплотипом CCG. Гаплотип TTT, который рассматривался как адаптивный вариант полиморфизма у гренландских эскимосов [2], имеет, согласно палеогеномным данным, верхнепалеолитическое происхождение. Присутствие этого гаплотипа примерно 14‒19 тыс. лет тому назад у населения Приамурья согласуется с гипотезой о том, что древние жители этого региона принимали участие в формировании населения Берингии, давшего начало американским индейцам [25].

Авторы заявляют, что все процедуры, выполненные в исследовании с участием людей, соответствуют этическим стандартам институционального и/или национального комитета по исследовательской этике и Хельсинкской декларации 1964 г. и ее последующим изменениям или сопоставимым нормам этики. От каждого из включенных в исследование участников было получено информированное добровольное согласие.

Авторы заявляют, что у них нет конфликта интересов.

Исследование выполнено за счет гранта Российского научного фонда № 22-24-00264, https://rscf.ru/project/22-24-00264/.

Об авторах

Б. А. Малярчук

Институт биологических проблем Севера Дальневосточного отделения Российской академии наук

Автор, ответственный за переписку.
Email: malyarchuk@ibpn.ru
Россия, Магадан, 685000

М. В. Деренко

Институт биологических проблем Севера Дальневосточного отделения Российской академии наук

Email: malyarchuk@ibpn.ru
Россия, Магадан, 685000

Г. А. Денисова

Институт биологических проблем Севера Дальневосточного отделения Российской академии наук

Email: malyarchuk@ibpn.ru
Россия, Магадан, 685000

А. Н. Литвинов

Институт биологических проблем Севера Дальневосточного отделения Российской академии наук

Email: malyarchuk@ibpn.ru
Россия, Магадан, 685000

И. К. Дамбуева

Институт биологических проблем Севера Дальневосточного отделения Российской академии наук

Email: malyarchuk@ibpn.ru
Россия, Магадан, 685000

Список литературы

Дополнительные файлы