CC BY
35
DOI: 10.32521/2074-8132.2019.3.067-076
Агджоян А.Т12), Богунов Ю.В.2Д), Богунова А.А.2), Каменщикова Е.Н.3), Запорожченко В.В.2Д), Пылёв В.Ю.2), Короткова ЕА34), Утриван СА2Д), Схаляхо РА1Д5)
Кошель С.М.6), Балановский О.П.и4), Балановская Е.В.24
1ФГБУН«Институт общей генетики им. НИ. Вавилова РАН», 119991, ул. Губкина, д. 3, Москва, Россия; 2 ФГБНУ«Медико-генетический научный центр», 115522, ул. Москворечье, д. 1, Москва, Россия; 3) Амурский гуманитарно-педагогический государственный университет, кафедра педагогики профессионального образования, 681000, ул. Кирова, д. 17, корп. 2, Комсомольск-на-Амуре, Россия;
4>АНО «Биобанк Северной Евразии», 115201, ул. Котляковская, д. 3, Москва, Россия; 5)Научный парк Санкт-Петербургского государственного университета, Санкт-Петербург, 198504;
6) МГУ имени МВ. Ломоносова, географический факультет, кафедра картографии
и геоинформатики, 119991, ГСП-1, Москва, Россия
МОЗАИКА ГЕНОФОНДА ЭВЕНКОВ: ЗАБАЙКАЛЬСКИЙ И АМУРСКИЙ СЕГМЕНТЫ
Материалы и методы. Гзнофонд двух групп восточных эвенков - забайкальских (N=50) и амурских (N=81) - изучен по подробной панели 60 SNP маркеров У-хромосомы. Анализ частот выявленных гаплогрупп проведен методами многомерной статистики и геногеографии. Образцы венозной крови собраны в ходе экспедиционных обследований 2010-2015 гг. у мужчин, не связанных родственными отношениями, все предки которых на протяжении не менее трех поколений по мужской линии относили себя к данной популяции эвенков и родились в районе обследования. Сбор образцов сопровождался добровольным письменным информированным согласием под контролем Этической комиссии ФГБНУ «МГНЦ».
Результаты и обсуждение. Вопреки значительной географической удаленности ареалов забайкальских и амурских эвенков, сходство их генофондов по У-хромосоме выражено ярко: более 2/3 их генофондов составляют гаплогруппы С2-М48(хБК1066), Ы-М2118 и /^1а-М198. В многомерном генетическом пространстве Сибири забайкальские и амурские эвенки вместе с эвенами Охотского побережья и Камчатки вошли в кластер Приамурья, и при этом значительно удалены от популяций бурят, монголов, якутов и народов Южной Сибири - тувинцев, алтайцев и хакасов. Карты генетических расстояний и от забайкальских, и от амурских эвенков демонстрируют их связь с генофондами Приамурья и Камчатки и значительную отдаленность от генофондов Сибири, Восточной и Центральной Азии. Включение в анализ данных литературных источников позволило выявить различия генофондов западных и восточных эвенков, тяготеющих к разным регионам Сибири: западные эвенки - к якутам и отчасти к народам Южной Сибири, восточные - к популяциям Дальнего Востока.
Заключение. Не выявлены генетические следы ни ассимиляции бурятами эвенкийских родовых групп, ни включения бурятских родовых групп в состав эвенков. Полученные результаты позволяют выдвинуть гипотезу сохранения эвенками Забайкалья и Приамурья «восточнотунгус-ского» генетического компонента, распространенного между Байкалом и Охотским побережьем и включенного в систему генофондов Дальнего Востока либо в составе более обширного «приамурского» компонента, либо наряду с ним.
Ключевые слова: популяционная генетика человека; генофонд; Дальний Восток; эвенки; У-хромосома; гаплогруппы
Введение
Эвенки являются не самым многочисленным народом Сибири (в 2010 г. насчитывалось лишь 38 тыс. чел. на территории РФ [Итоги переписи ... 2010]), но определенно - одним из самых широко расселенных: их современный ареал охватывает территории Красноярского края, республик Якутия (Саха) и Бурятия, Забайкальского и Хабаровского краев, Амурской области, а также северные области Монголии и Китая. На всем своем этническом ареале эвенки проживают в близком соседстве со многими другими народами Сибири: якутами, долганами, юкагирами, кетами и селькупами - вдоль северной полосы ареала, с популяциями Южной Сибири (алтайцами, хакасами, тувинцами, бурятами, монголами) и Приамурья (ульчи, нанайцы, нивхи) - вдоль южной. Ареал эвенков условно подразделяют на западную и восточную части, которые пересекаются по линии р. Лена - оз. Байкал. Исследователи этногенеза эвенков отмечают, что за последние несколько веков прослеживались контакты с кетами и предками энцев, группами тувинцев и бурят, якутами и юкагирами. Одним из следов таких контактов являются аналогичные названия родов у эвенков и соседствующих с ними народов [Этническая история ... 1982].
Широта расселения эвенков не могла не отразиться на разнообразии хозяйственных типов (на севере - оленеводческий, на юге - скотоводческий, также встречались и группы эвенков, «единственным домашним животным которых была собака» [Этническая история ... 1982]) или особенностях языка (в эвенкийском языке - три наречия с большим числом говоров: северное, южное и восточное [Афанасьева, 2010]). Изменчивость внешних признаков у эвенков позволила физическим антропологам включить их в тип «североазиатских монголоидов» [История и культура дальневосточных эвенков ... 2010]. Полиморфизм ряда классических маркеров у западных и восточных эвенков выявляет небольшие (по частотам иммунологических маркеров - групп крови ABO, MN, Lewis) или разительные (по частотам аллелей генов Km и Gm иммуноглобулинов) различия в частотах аллелей, зачастую приближая генофонд западной группы к популяциям Средней и Южной Сибири, а восточных эвенков - к популяциям северо-востока Азии [Генофонд и ге-ногеография народонаселения. 2000, т.1].
Митохондриальный генофонд эвенков изучается на протяжении последних 25 лет, начиная с работы [Torroni et al., 1993], как в сравнении с генофондами других народов Сибири и Восточной Азии [Kong et al., 2003; Pakendorf et al., 2006;
Ракеле е! а1., 2007; йегепко е! а1., 2007], так и в связи с вопросом заселения Америки ^апкоуэкауа е! а1., 2005; Татт е! а1., 2007]. Разнообразие гап-лотипов популяций эвенков соответствует географической широте проживания конкретной субпопуляции: например, эвенки Тунгуски накопили разнообразие, сходное с таковым в популяциях якутов, и оно заметно ниже, чем разнообразие более крупных популяций, проживающих южнее (например, бурят и монголов) ^апкоуэкауа е! а1., 2005; Ракеп^И7 е! а1., 2006]. Примечательно взаимное сходство всех изученных групп эвенков. Наиболее близкими к ним по частотам гаплогрупп оказались якутские популяции [Ракеп^И7 е! а1., 2006], с которыми эвенки разделяют характерную черту - преобладание гаплогрупп С4а1, С4а2, С5а, й4, й5, которые в совокупности представляют более половины генофонда (в некоторых субпопуляциях - более 70%). Доля западноевразийских гаплогрупп невелика, и представлены они в основном в самой западной части ареала. Невелика и доля линий, характерных для бассейна Амура, например гаплогруппы У1 (2%). В следовых количествах обнаруживаются «палеоазиатские» гаплогруппы й2 и й3 (4%).
По данным другой гаплоидной системы - маркерам У-хромосомы - изучены популяции эвенков (N=57) и эвенов (N=24) из республики Якутия (Саха) в работе [Ре^гоуа е! а1., 2013]. В генофонде изученных групп тунгусов гаплогруппа М-М178 (по классификации 2013 года - Же) охватывает более половины разнообразия эвенков и 2/5 генофонда эвенов, вторая по частоте - гаплогруп-па С2-М48 (по классификации 2013 года - С3с) -достигает четверти разнообразия эвенков и 2/5 всех линий у эвенов. Опубликованная в работе [Ре^гоуа е! а1., 2013] панель маркеров актуальна на 2013 год, однако не позволяет провести сравнение с изученными в данной статье популяциями на современном уровне.
В работе [Ре^гоуа е! а1., 2013] по результатам анализа широкогеномных панелей маркеров аутосомного генома эвенки на графике главных компонент приближаются к якутам, юкагирам и долганам, в меньшей степени - к нганасанам и популяциям Южной Сибири - хакасам, тувинцам, эвены приближаются к группам Дальнего Востока - орочонам и нивхам. Стоит отметить, что положение тунгусских популяций на графике главных компонент в исследовании [Ре^гоуа е! а1., 2013], похоже, отражает географию сбора образцов для анализа ДНК эвенов и эвенков.
Полные геномы эвенков и эвенов опубликованы в работах [МаШек е! а1., 2016; Радап1 е! а1., 2016]. В статье [МаШек е! а1., 2016] на древе родства популяций эвены объединяются в кластер с ульчами,
орочонами и нанайцами Китая, что отражает дальневосточный круг генетических связей. В статье [Радат е! а1., 2016] при анализе полных геномов по вкладу различных региональных компонентов у эвенков (популяции из Красноярского края и республики Якутия (Саха) и эвенов (из Якутии и Магаданской области) преобладает компонент, выраженный в популяциях Средней (якуты) и Южной Сибири (алтайцы, тувинцы, буряты), Центральной Азии (монголы). Значительно в меньшей степени у эвенков и эвенов представлен компонент, характерный для коряков, чукчей и эскимосов (у эвенов Магаданской области - четверть генофонда, у эвенков и эвенов Якутии - десятая часть генофонда).
Иными словами, по разным системам генетических маркеров изучены преимущественно эвенки западной части ареала, и по сути - это несколько популяций из Красноярского края и Якутии, образцы из которых проанализированы в ходе нескольких исследований по У-хромосомным и/или аутосомным (широкогеномные панели маркеров или полногеномное секвенирование) маркерам.
В данной работе рассматривается генофонд двух популяций восточных эвенков - забайкальских и амурских, которые являются практически неизученными по генетическим маркерам как гаплоидных систем (мтДНК и У-хромосомы), так и аутосомного генома. Не претендуя на разрешение широких и давних проблем (например, происхождения и прародины эвенков), данная работа представляет данные о генофонде двух восточных сегментов обширной эвенкийской общности, являющейся важным звеном при реконструкции истории народонаселения Сибири и его связей с соседними регионами.
Материалы и методы
Популяции эвенков были обследованы в Забайкальском крае (пять районов - Каларский, Тун-гокоченский, Кыринский, Акшинский и Баунтовс-кий, суммарно N=50) и в Хабаровском крае (Тугу-ро-Чумиканский и Аяно-Майский районы, N=81). Суммарно материалом для изучения генофонда эвенков послужили образцы крови (N=131) мужчин, включенные в Биобанк народонаселения Северной Евразии [Биобанк Северной Евразии] и собранные согласно правилам формирования популяционных выборок Биобанка, изложенным в работе [Балановская с соавт., 2016]. В выборку включались только мужчины, предки которых на протяжении не менее чем трех поколений (до деда включительно) относили себя или к забайкальским, или к амурским эвенкам.
Выделение ДНК из венозной крови проведено методом фенол-хлороформной экстракции с использованием протеиназы К. Генотипирование 60 SNP-маркеров Y-хромосомы для образцов ДНК эвенков проведено на приборе 7900HT Real-Time PCR System (Applied Biosystems) с использованием TaqMan зондов. Наименования гаплогрупп даны согласно [ISOGG, 2019]. Популяции для сравнения получены из базы данных «Y-base», разработанной под руководством О.П. Балановс-кого [БД Y-base]. Генетические расстояния Нея рассчитаны по панели 30 гаплогрупп Y-хромосо-мы (C2-M217(xM407,F3918,F3791 ,M48), C2-M407, C2-F3918, C2-F3791x(F5481), C2-F5481, C2-M48(xSK1066), C2-SK1066, E-M35.1, G-M285, G-P15, I-M253, I-P37.2, J-M267, J-M172, L-M11, N-P43, N-B211, N-M2118, N-VL29, N-Z1936, N-F4205, N-B202, N-B479, O-M122, Q-M242, R1a-M198, R1b-M73, R1b-M269, R2-M124, T-L206) в программе DJgenetic [Balanovsky et al., 2008] и визуализированы методом многомерного шкалирования в программе Statistica [Statistica Software, 2005].
Картографический анализ проведен в программе GeneGeo [Balanovsky et al., 2011; Кошель, 2012]. Карты (модели) распространения 30 гаплогрупп Y-хромосомы (C2-M217x(M407, F3918, F3791, M48), C2-M407, C2-F3918, C2-F3791x(F5481), C2-F5481, C2-M48(xSK1066), C2-SK1066, E-M35, G-M285, G-P15, I-M253, I-P37.2, J-M267, J-M172, L-M11, N-P43, N-B211, N-M2118, N-VL29, N-Z1936, N-F4205, N-B202, N-B479, O-M122, Q-M242, R-M198, R-M73, R-M269, R-M124, T-L206) были построены методом средневзвешенной интерполяции с радиусом влияния 1700 км и значением степени весовой функции 3 [Русский генофонд ... 2007, с. 343-345; Кошель, 2012] по данным о частотах гаплогрупп Y-хромосомы из базы данных «Y-base», разработанной под руководством О.П. Ба-лановского [БД Y-base]. Были созданы две карты генетических расстояний Нея - для обеих групп эвенков, и алгоритм построения каждой карты состоял из двух этапов. На первом этапе для каждой из 30 гаплогрупп была построена карта генетических расстояний от средней частоты этой гап-логруппы в популяции эвенков до интерполированных значений частоты этой гаплогруппы в каждой точке карты. На втором этапе было найдено среднее значение генетических расстояний по всем гаплогруппам для каждой точки карты, и таким образом создана обобщенная карта генетических расстояний от эвенков Забайкалья или Приамурья.
Рисунок 1. Частоты гаплогрупп У-хромосомы в изученных популяциях эвенков Figure 1. The frequencies of Y-chromosomal haplogroups in the studied groups of Evenks
Результаты
Спектр гаплогрупп Y-хромосомы в изученных популяциях
Генофонды забайкальских и амурских эвенков оказались удивительно похожи (рис. 1) как по низкому генетическому разнообразию (более 90% генофонда каждой группы представлено лишь пятью гаплогруппами У-хромосомы), так и по вкладу основных гаплогрупп (С2-М48(хБК1066) и М-М2118, которые охватывают две трети генофонда каждой группы эвенков и представлены в сходных соотношениях).
Доминирующая у забайкальских и амурских эвенков гаплогруппа С2-М48 распространена в популяциях Северо-Восточной и Центральной Азии (рис. 1, врезка-карта), появилась около 13 тыс. лет назад и к настоящему времени включает в себя кластер линий с датировками разной хронологической глубины. В генофонде изученных популяций представлен, вероятно, кластер линий внутри гаплогруппы С2-М48, за исключением как минимум одной - линии С2-БК1066 (обнаружена с наибольшей частотой в Центральной Азии в популяциях Западной Монголии, реже - у казахов, узбеков; возраст составляет 1000±300 лет).
Вторая гаплогруппа, определяющая сходство изученных популяций эвенков - М-М2118, со-
ставляет в генофонде забайкальских эвенков 14%, а амурских - 9%. Эта гаплогруппа появилась около 4-5 тыс. лет назад и среди современного народонаселения наиболее характерна для популяций Средней Сибири, а пика частоты достигает у якутов [Шытае е! а1., 2016].
Различия в спектре гаплогрупп У-хромосомы у изученных популяций (рис. 1) связаны со вкладом двух линий: гаплогрупп N-8479 и Ма-М198. В генофонде амурских эвенков гаплогруппа N-8479 составляет 20% (у забайкальских - отсутствует), её возраст около 4-5 тыс. лет, а современная география охватывает преимущественно Дальний Восток [Шытае е! а1., 2016]. У обеих групп эвенков встречена гаплогруппа Ма-М198: у амурских эвенков - в 2 раза чаще (10%), чем у забайкальских (4%). Важно отметить, что гаплогруппа Ма-М198 является одной из наиболее широко распространенных в Евразии и включает в себя кластер линий разной хронологической глубины и с разной географией [Underhill е! а1., 2014]. Также среди эвенков Забайкалья обнаружено небольшое число индивидов - носителей гаплогрупп 11-М253 и 12-Р37.2 (рис. 1), что скорее всего является следом метисации с населением из Восточной Европы и Урало-Поволжья, где распространены данные линии У-хромосомы.
Рисунок 2. Положение эвенков Забайкалья и Приамурья в генетическом пространстве Сибири (график многомерного шкалирования, показатели стресса=0.15, алиенации=0.19) Figure 2. The position of the Evenks from Transbaikalia and from Amur Region in the genetic space of Siberia
(multidimensional scaling graph, stress = 0.15, alienation = 0.19) Примечания. Для возможности сравнения по широкой панели 30 гаплогрупп Y-хромосомы в анализ вовлечены только популяции, изученные нашим научным коллективом. В правой части рисунка приведены генетические расстояния Нея от усредненной популяции эвенков (забайкальских и амурских) до каждой из групп сравнения.
Notes. For the comparison on the wide panel of 30 Y chromosomal haplogroups only populations studied by the research team are involved in the analysis. The right side of the figure shows Nei's genetic distances from the average Evenk population (Transbaikalian and Amur) to each of the comparison groups.
Положение эвенков Забайкалья и Приамурья в многомерном генетическом пространстве
Генетические соотношения забайкальских и амурских эвенков в контексте популяций Сибири, учитывающие весь спектр гаплогрупп У-хромосо-мы, представлены на рис. 2. Отметим, что в одной из предыдущих работ [Агджоян с соавт., 2019] результат аналогичного анализа в масштабе Северной Азии (26 популяций) по узкой панели маркеров У-хромосомы (16 гаплогрупп высокого иерархического уровня, «стволов») выявил сходство генофондов эвенков Забайкалья и Приамурья между собой и с популяцией эвенов охотского побережья и Камчатки. В текущей статье многомерный анализ (рис. 2) проведен по иному принципу: меньшее число популяций (п=21) изучены по значительно более широкой панели маркеров (30 гаплогрупп У-хромосомы, включая 6 линий внутри гаплогруп-пы С2-М217 и 7 локальных линий внутри гапло-группы 1Ч-М178).
В генетическом пространстве популяций Сибири (рис. 2) выделяются три кластера: кластер А
объединяет изученные группы эвенков с популяциями Приамурья и Камчатки, кластер Б включает популяции Южной Сибири - северных и южных алтайцев, хакасов сагайцев, тувинцев и тоджинцев; кластер В - группы бурят Бурятии, Иркутской области и Забайкалья. Отдалённо от указанных кластеров на графике расположились популяции якутов, чукчей Камчатки, монголов, тофалар и хакасов-качинцев.
Внутри кластера А генетические связи между эвенками и другими популяциями неодинаковы (рис. 2, график и врезка). Наиболее близки друг к другу обе группы эвенков (Ь=0.02), к ним максимально тяготеют эвены (Ь=0.04). Промежуточная -в масштабе данного графика (рис. 2) - степень генетического сходства (0.14<Ь<0.29) с изученными эвенками прослеживается для ульчей, коряков, нивхов, негидальцев. Генетическое расстояние от эвенков до популяции нанайцев (которая отражает усредненное положение четырёх групп нанайцев - низовых, верховых, озерных и горинских)
Рисунок 3. Карты генетических расстояний от эвенков Забайкалья (рисунок А) и Приамурья (рисунок Б) Figure 3. Maps of genetic distances from the Evenks from Transbaikalia (Figure A) and from Amur Region (Figure B) Примечания. Ареалы популяций, генетически наиболее близких к эвенам, отражены желто-зелеными тонами и расположены в Приамурье. Карта построена по частотам 30 гаплогрупп У-хромосомы. Черными точками на карте обозначены популяции сравнения (N=103 популяции), сиреневой звездой - положение изученных популяций эвенков. Коэффициент корреляции Пирсона между картами на рисунках А и Б составляет r=0.95. Коэффициенты корреляции карт на рисунках А и Б с картой генетических расстояний от эвенов из работы [Агджоян с соавт., 2019] составляют r=0.87 и r=0.93 соответственно.
Notes. Areas of populations that are genetically closest to the Evenks are shown in yellow-green tones and are located in the Amur region. The map is based on the frequencies of 30 Y-chromosomal haplogroups. Black dots on the map indicate comparison populations (N = 103 populations), the purple star indicates the position of the studied Evenk population. The Pearson correlation coefficient between the maps in Figures A and B is r = 0.95. The correlation coefficients of the maps in Figures A and B with the genetic distance map for Evens from the paper [Agdzhoyan et al., 2019. In Russian] are r = 0.87 and r = 0.93.
составляет d=0.70, что тем не менее в 2 раза меньше расстояние от эвенков до монголов и в 3 раза меньше аналогичного расстояния до бурят Забайкалья ^=2.77) (рис. 2, врезка). Отметим, что генетическое расстояние между ближайшими соседями, столетиями проживающими в одном крае -эвенками и бурятами Забайкалья ^=2.55) - в 2,5 раза больше аналогичного показателя между эвенками Забайкалья и нанайцами ^=1.03), разделенными расстоянием более 1000 км хребтами гор и долинами рек .
От популяций сравнения за пределами кластера А эвенки Забайкалья и Приамурья расположены на значительном генетическом расстоянии (в интервале 1.46^<3.62), при этом значительно удалены как от северных групп (якуты d=1.51, чукчи d=1.85), так и от проживающих южнее (монголы d=1.46, южные алтайцы d=1.94).
Дальневосточный паттерн эвенков на картах генетических расстояний
Карты генетических расстояний от эвенков Забайкалья (рис. 3А) и Приамурья (рис. 3Б) созданы по такой же широкой панели маркеров, как и график многомерного шкалирования (рис. 2). При этом каждая карта отражает круг связей изученного генофонда с окружающими популяциями Сибири, что дополняет картину, представленную на графике многомерного шкалирования (рис. 2). На рис. 3 наглядно отражено, насколько похожи между собой обе карты: максимальное сходство прослеживается между восточными тунгусскими популяциями - забайкальскими и амурскими эвенками, а также с эвенами Охотского побережья и Камчатки ^<0.12), более слабое сходство -с генофондами Приамурья и севера Камчатки (0.14^<0.4, ульчи, нивхи, нанайцы, негидальцы, удэгейцы, коряки севера Камчатки), а также с популяциями правобережья Амура. Значительная генетическая отдаленность от эвенков Забайкалья
и Приамурья ^>0.6) прослеживается для всех остальных популяций сравнения из Сибири и восточной части Центральной Азии.
Обсуждение
Сходство генофондов изученных популяций эвенков как между собой и с эвенами Охотского побережья и Камчатки, так и с широким кругом популяций Дальнего Востока связано, в первую очередь, с высокой долей гаплогруппы С2-М48(хБК1066), которая отмечена с частотой 10-40% среди ульчей, нивхов, нанайцев, коряков и чукчей Камчатки. Значительная доля одного признака может означать последствия дрейфа генов, но это представляется более вероятным событием не для нескольких одновременно, а для одной, но - предковой для многих популяции. С другой стороны, внутри гаплогруппы С2-М48(хБК1066) может скрываться комплекс более молодых, чем 4 тыс. лет, линий (субвариантов). Соответственно, в разных популяциях Дальнего Востока может быть представлен как комплекс разных субвариантов в разных долях, так и накопление «своего» локального субварианта в каждой популяции или же накопление у всех одной и той же линии: ответ на эти вопросы может представить только дальнейший анализ филогеографии гаплогруппы С2-М48(хБК1066).
Другой общий для забайкальских и амурских эвенков признак - гаплогруппа Ы-М2118 - отражает их связь с популяциями Средней Сибири, что может указывать либо на контакты предков эвенков с якутами, либо на другие источники заимствования изученными группами «среднесибирского» субстрата генофонда.
Истоки третьего общего для изученных групп эвенков признака - гаплогруппы К1э-М198 (рис. 1) -могут как восходить к палеолиту, так и отражать след ассимиляции мигрантов из европейской части материка за последние несколько столетий. Прояснить этот вопрос может только анализ филогеографии гаплогруппы К1э-М198: если у эвенков будут обнаружены ветви с азиатской географией и прародиной - это довод в пользу первого предположения, если с европейской - в пользу второго.
Вклад гаплогруппы 11-М253 в генофонд эвенков Забайкалья указывает на следы ассимиляции переселенцев из европейской части материка: появилась гаплогруппа 11-М253 в Европе в неолите и среди современного народонаселения Северной Евразии практически не встречается восточнее Уральских гор.
Паттерн включения трех тунгусских популяций - забайкальских и амурских эвенков, эвенов Охотского побережья и Камчатки - в дальневос-
точный кластер на графике многомерного шкалирования (рис. 2, кластер А) также представлен и в географическом пространстве, на картах генетических расстояний (рис. 3). Довольно неожиданным оказалось столь близкое сходство эвенков Забайкалья с популяциями Дальнего Востока и при этом столь явная их отдаленность от генофондов разных групп бурят (и забайкальских, и прибайкальских, при этом значительно различающихся между собой по спектру гаплогрупп У-хромосомы). Иными словами, не выявлены генетические следы ни ассимиляции бурятами эвенкийских родовых групп, ни включения бурятских родовых групп в состав эвенков: иначе бы значения генетических расстояний между бурятами и эвенками Забайкалья были бы значительно ниже.
С другой стороны, полученный результат может указывать на более широкую (от Охотского моря до Байкала) географию «приамурского» генетического компонента либо же на выделение в его составе (или наряду с ним) «восточнотунгус-ского» компонента, наиболее характерного для восточных эвенков и эвенов: по причине неизученности других популяций восточных эвенков и эвенов по современной широкой панели маркеров география этого компонента пока отражается неполностью. Однако в пользу второго предположения свидетельствует некоторый «разрыв» между значениями генетических расстояний в пределах кластера А на графике многомерного шкалирования (рис. 2, график и врезка): среднее генетическое расстояние от усредненной популяции эвенков для группы «эвенки Забайкалья, Приамурья и эвены» составляет 0.03, в группе «ульчи, коряки, нивхи, негидальцы» - 0.22 (рис. 2, врезка). Иными словами, три популяции тунгусов среди заданного круга популяций наиболее близки к генофондам Приамурья и Камчатки. Но сохранится ли такая закономерность и в дальнейшем, при изучении большего числа популяций эвенков и эвенов в пределах восточной части их ареала? Можно предполагать, что при подробной изученности десятка популяций «восточных» эвенков и эвенов на карте или графике могут отразиться два пересекающихся по территории компонента - «амурский» и «восточнотунгусский», включенные в систему генофондов Дальнего Востока.
Согласно опубликованным данным по гапло-группам У-хромосомы, широко- и полногеномным маркерам [Ре^гоуа е! а1., 2013; Радап е! а1., 2016], генофонд западных эвенков (популяции из Якутии, Красноярского края, Забайкалья) тяготеет не к Дальнему Востоку, а к соседним популяциям Средней Сибири, т.е. представляет иной сегмент эвенкийского генофонда, чем описанный в данной работе «восточный» сегмент. Похожая ситуация выявлена и для эвенов: популяции из Магаданской
области и Якутии [Fedorova et al., 2013; Mallick et al., 2016; Pagani et al., 2016] занимают промежуточное положение между кластерами Средней Сибири и Северо-Восточной Азии, при этом тяготея к Приамурью (к нивхам в работе [Fedorova et al., 2013]), что также отличается от генетического портрета эвенов Охотского побережья и Камчатки [Агджоян с соавт., 2019], сходного с популяциями восточных эвенков и народов Амура. Результат, сходный с показанным в данной работе, отражен на графике главных компонент по частотам гаплогрупп мтДНК в работе [Derenko et al., 2007], где западные эвенки (представлены популяцией из Красноярского края) объединяются в кластер с алтайцами, тувинцами и якутами, а восточные эвенки (популяция из Забайкалья) располагаются в отдалении от этого кластера и от западных эвенков. Удаленность восточных эвенков на упомянутом графике из [Derenko et al., 2007] связана с меньшей долей второй главной компоненты, которая убывает в направлении запад-восток и, соответственно, может быть связана со вкладом «приамурского» компонента генофонда.
Из сопоставления с литературными данными вытекает заключение о различиях между генофондами эвенков на западном и восточном «полюсах» их ареала. Но можно ли прогнозировать, какими окажутся генофонды эвенков между этими «полюсами»? Если следовать подразделению всех эвенков на три антропологических типа («северозападный» - в центральной части Средней Сибири, «южный» с примесью «центральноазиатского компонента» - вдоль южных берегов озера Байкал, «восточный» - эвенки Забайкалья и Приамурья) [История и культура дальневосточных эвенков ... 2010, с. 25], то наряду с «западным» и «восточным» следует ожидать как минимум ещё один «полюс» - «южный». Принимая же во внимание масштаб ареалов каждого антропологического типа эвенков и наличие у тунгусов родовых подразделений, можно ожидать и внутреннюю подразделен-ность на кластеры. Преобладание среди современных эвенков сельского населения (3/4 всех эвенков, или 28 тыс. чел. из 38 тыс.чел. на 2010 год [Итоги переписи ... 2010]) оставляет небольшую надежду на возможность изучить многообразие популяций на всем ареале расселения этой этнической группы, и, возможно, реконструировать не только фрагменты «генетической мозаики» тунгусов, но и важные звенья в древней истории народонаселения Сибири.
Благодарности
Авторы выражают благодарность администрациям Тугуро-Чумиканского и Аяно-Майского районов Хабаровского края и Министерству здравоохранения Хабаровского края, а также администрациям Калар-ского, Тунгокоченского, Кыринского, Акшинского и Баунтовского районов Забайкальского края и Министерству здравоохранения Забайкальского края, при организационной поддержке которых проведено экспедиционное обследование генофондов амурских и забайкальских эвенков.
Авторы выражают глубокую благодарность: к.филос.н., доценту кафедры философии Забайкальского государственного университета Андрею Григорьевичу Букину, Антонине Андреевне Кузнецовой и Алексею Геннадьевичу Балаганскому за активное участие в экспедиционном обследовании и этапах пробоподготовки биологических образцов эвенков Забайкальского края.
Исследование выполнено при финансовой поддержке РНФ в рамках научного проекта № 17-1401345 (в части изучения генофонда эвенков Приамурья) и проекта РФФИ № 18-09-00487 (в части картографического анализа генофонда эвенков Забайкалья). Другие части исследования выполнены в рамках тем Государственного задания Минобрнауки России для ФГБНУ «МГНЦ» и Государственного задания Минобрнауки России для ИОГен РАН.
Библиография
Агджоян А.Т., Богунов Ю.В., Богунова A.A., Каменщикова E.H., Кагазежева Ж.А. с соавт. Генетический портрет охотских и камчатских эвенов // Вестник Московского университета. Серия 23. Антропология, 2019. № 2. С. 116-125. Афанасьева E.Ô. Фонология и фонетика эвенкийского языка: учеб. пособие. Улан-Удэ: Изд-во Бурятского гос. университета, 2010. 116 с.
БалановскаяЕ.В., Балановский О.П. Русский генофонд на Русской равнине. М.: Луч, 2007. 416 с.
БалановскаяЕ.В., Жабагин М.К., Агджоян А.Т., Чухряева М.И., Маркина Н.В. с соавт. Популяционные биобанки: принципы организации и перспективы применения в геногеографии и персонализированной медицине // Генетика, 2016. Т. 52. № 12. С. 1371-1387. DOI: 10.7868/S001667581612002X. БД «Y-base». URL: http://www.genofond.ru/genofond.ru/ default26d8e.html?s=0&p=346. (дата обращения - 20.06.2019). Биобанк Северной Евразии. URL: Ь|йр:/Л№№«.биобанк.рф (дата обращения - 20.06.2019).
Генофонд и геногеография народонаселения / Под ред. Ю.Г. Рыч-кова: Т. 1. Генофонд населения России и сопредельных стран. СПб.: Наука. 2000. 611 с.
История и культура дальневосточных эвенков: историко-этногра-фические очерки / Ред. В.А. Тураев. СПб.: Наука, 2010. 334 с. Итоги Всероссийской переписи населения 2010 г. Т. 4. Национальный состав и владение языками, гражданство. Электронный ресурс. URL: http://www.gks.ru/ free_ doc/new_ site/ perepis2010/croc/perepis_ itogi1612.htm. (дата обращения -20.06.2019).
Кошель С.М. Геоинформационные технологии в геногеогра-фии // Современная географическая картография. М: Дата+, 2012. С. 158-166.
Этническая история народов Севера / Отв. ред. И.С. Гурвич. М.: Наука, 1982. 269 с.
Сведения об авторах
Агджоян Анастасия Торосовна, к.б.н.; ОРСЮ Ю 0000-0002-8776-2934; [email protected]; Богунов Юрий Васильевич, к.б.н.; [email protected]; Богунова Анна Александровна, к.б.н.; [email protected]; Каменщикова Елена Николаевна, к.п.н.; [email protected]; Запорожченко Валерий Владимирович, [email protected]; Пылёв Владимир Юрьевич, [email protected];
Короткова Надежда Александровна,
ORCID ID 0000-0003-1760-4700; [email protected];
Утриван Сергей Александрович, ORCID ID 0000-0002-4459-7211;
Схаляхо Роза Арамбиевна, к.б.н.; ORCID ID 0000-0002-4005-1887;
Кошель Сергей Михайлович, к.г.н.;
ORCID ID 0000-0002-4540-2922; [email protected];
Балановский Олег Павлович, д.б.н., профессор РАН;
ORCID ID 0000-0003-4218-6889; [email protected];
Балановская Елена Владимировна, д.б.н., профессор;
ORCID ID 0000-0002-3882-8300; [email protected].
Поступила в редакцию 24.07.2019, принята к публикации 01.08.2019.
Agdzhoyan AT.U), Bogunov Y.V.2a), Bogunova А.А.2), Kamenshikova E.N.3), Zaporozhchenko VV.2^, Pylev V.Yu.2), Korotkova NA34), Utrivan SAy), Skhalyakho RA.1A5),
Koshel S.M.6), Balanovsky O.P.1A4), Balanovska E.V.2,4)
1 Vavilov Institute of General Genetics, RAS, Gubkina st, 3, Moscow, 119991, Russia; 2) FSBI«Research Centre for Medical Genetics», Moskvorechie st, 1, Moscow, 115522, Russia; 3) Department of Pedagogics and Professional Education, Amur State University of Humanities and Pedagogy,
Kirova st, 17/2, Komsomolsk-on-Amur, 681000, Russia; 4) Biobank of Northern Eurasia, Kotlyakovsaya st, 3, Moscow, 115201, Russia; 5) Science Park, St. Petersburg State University, St. Petersburg, 198504; 6 Lomonosov Moscow State University, Faculty of Geography, Department of Cartography and Geoinformatics,
GSP1, Moscow, 119991, Russia
THE MOSAIC OF THE EVENKS GENE POOL: TRANSBAIKALIAN
AND AMUR SEGMENTS
Materials and methods. The gene pool of two eastern Evenk groups, Transbaikalian and Amur ones, has been studied using detailed panel including 60 Y-chromosomal SNP markers. The frequencies of haplogroups identified have been analyzed by multidimensional statistics and genegeography methods. The venous blood was sampled in the course of expeditions in 2010-2015 from unrelated males. No less than three generations of male ancestors identified themselves as belonging to the given Evenk population and were born in the studied region. The sampling was performed with the voluntary written informed consent.
Results. In spite of considerable geographic distance between the ranges of Transbaikalian and Amur Evenks there is a pronounced similarity in their Y-chromosomal gene pools: haplogroups C2-M48(xSK1066), N-M2118 and R1a-M198 account for more than 2/3 of their gene pools. In the Siberian multidimensional genetic space, Transbaikalian and Amur Evenks, together, with the Evenks of Okhotsk coast and Kamchatka, formed the Amur region cluster, considerably removed from Buryats, Mongols, Yakuts and South Siberian groups - Tuvans, Altaians and Khakass. Maps of genetic distances from both the Transbaikalian and Amur Evenks display their close relation to Amur region and Kamchatka gene pools, and a noticeable remoteness from gene pools of Siberia, East and Central Asia. When the data from literature are included into consideration, the differences become clear between West and East Evenk gene pools, which may be attributed to their geographic position in different regions of Siberia: the West Evenks are closer to Yakuts and partly to South Siberian peoples, and the Eastern Evenks - to populations of the Far East.
Conclusion. No genetic traces of assimilation of Evenk tribes by Buryats, nor Buryat tribe-inclusion into Evenks have been found. The results obtained give grounds for the hypothesis about "East Tungusic" genetic component being preserved in Transbaikalian and Amur Evenks. The component occurs mostly between Baikal Lake and Okhotsk coast and is included to the Far East gene pool system either as a part of a wider "Amur region" constituent, or along with the latter.
Keywords: human population genetics; gene pool; Far East; Evenks; Y-chromosome; haplogroups
References
Agdzhoyan A.T., Bogunov Yu.V., Bogunova A.A., Kamenshchikova E.N., Kagazezheva Zh.A. et al. Geneticheskii portret okhotskikh i kamchatskikh evenov [Genetic portrait of the Okhotsk and Kamchatka Evens], Moscow University Anthropology Bulletin [Vestnik Moskovskogo universiteta. Seriya 23. Antropologiya], 2019, 2, pp. 116-125. (In Russ.).
Afanas'eva E.F. Fonologiya i fonetika evenkiiskogo yazyka: ucheb. posobie [Phonology and phonetics of the Evenki language: studies. allowance]. Ulan-Ude, Buryat State University Publ., 2010. 116 p. (In Russ.).
Balanovska E.V., Balanovsky O.P. Russkii genofond na Russkoi ravnine [Russian gene pool on the Russian Plain]. Moscow, Luch Publ., 2007. 416 p. (In Russ.).
Balanovskaya E.V., Zhabagin M.K., Agdzhoyan A.T., CHuhryaeva M.I., Markina N.V. s soavt. Populyacionnye biobanki: principy organizacii i perspektivy primeneniya v genogeografii i personalizirovannoj medicine [Population biobanks: Organizational models and prospects of application in genegeography and personalized medicine]. Genetika [Russian Journal of Genetics], 2016, 52 (12), pp. 1227-1243. DOI: 10.1134/S1022795416120024. (In Russ.). BD "Y-base". [DB Y-base]. Available at: http://www.genofond.ru (Accessed: 20.06.2019). (In Russ.).
Biobank Severnoj Evrazii. [Biobank of Northern Eurasia]. Available at: www.6Mo6aHK.p$. (Accessed: 20.06.2019). (In Russ.). Genofond i genogeografiya narodonaseleniya: Tom 1. Genofond naseleniya Rossii i sopredel'nykh stran. [Gene Pool and Genegeography of population: Gene Pool of population of Russia and contiguous countries] Ed. YU.G. Rychkova. St. Petersburg, Nauka Publ., 2000, 611 p. (In Russ.).
Istoriya i kul'tura dal'nevostochnykh evenkov: istoriko-etno-graficheskie ocherki [The history and culture of the Far Eastern Evenks: historical and ethnographic essays]. St. Petersburg, Nauka Publ., 2010. 334 p. (In Russ.).
Itogi Vserossiiskoi perepisi naseleniya 2010 goda. T. 4. Nacional'nyj sostav i vladenie yazykami, grazhdanstvo. [The Results of population census in Russia in 2010]. Available at: http://www.gks.ru/ free_doc/new_site/perepis2010/croc/perepis_itogi1612.htm. (Accessed:- 20.06.2019). (In Russ.).
Koshel S.M. Geoinformatsionnye tekhnologii v genogeografii [Geoinformation technologies in genogeography]. In: Sovremennaya geograficheskaya kartografiya [Modern geographical cartography]. Moscow, Data+ Publ., 2012, pp. 158-166. (In Russ.). Etnicheskaya istoriya narodovSevera [Ethnic history of the peoples of the North]. Moscow, Nauka Publ., 1982. 269 p. (In Russ.). Balanovsky O., Dibirova K., Dybo A. et al. Parallel evolution of genes and languages in the Caucasus Region. Mol. Biol. Evol., 2011, 28, pp. 2905-2920. DOI: 10.1093/molbev/msr126. Balanovsky O., Rootsi S., Pshenichnov A. et al. Two sources of the Russian patrilineal heritage in their Eurasian context. Am. J. Hum. Genet., 2008, 82, pp. 236-250. DOI: 10.1016/j.ajhg.2007.09.019. Derenko M., Malyarchuk B., Grzybowski T., Denisova G., Dambueva I. et al. Phylogeographic analysis of mitochondrial DNA in northern Asian populations. Am. J. Hum. Genet., 2007, 81 (5), pp. 1025-1041. DOI: 10.1086/522933.
Ilumae A.M., Reidla M., Chukhryaeva M., Jarve M., Post H. et al. Human Y Chromosome Haplogroup N: A Non-trivial Time-Resolved Phylogeography that Cuts across Language Families. Am. J. Hum. Genet., 2016, 99 (1), pp. 163-173. DOI: 10.1016/j.ajhg.2016.05.025. Kong Q.P., Yao Y.G., Liu M., Shen S.P., Chen C. et al. Mitochondrial DNA sequence polymorphisms of five ethnic populations from northern China. Hum. Genet., 2003, 113 (5), pp. 391-405. DOI: 10.1007/s00439-003-1004-7.
Mallick S., Li H., Lipson M., Mathieson I., Gymrek M. et al. The Simons Genome Diversity Project: 300 genomes from 142 diverse populations. Nature, 2016, 538 (7624), pp. 201-206. DOI: 10.1038/ nature18964.
Pagani L., Lawson D.J., Jagoda E., Mörseburg A., Eriksson A. et al. Genomic analyses inform on migration events during the peopling of Eurasia. Nature, 2016, 538 (7624), pp. 238-242. DOI: 10.1038/ nature19792.
Pakendorf B., Novgorodov I.N., Osakovskij V.L., Danilova A.P., Protod'jakonov A.P. et al. Investigating the effects of prehistoric migrations in Siberia: genetic variation and the origins of Yakuts. Hum. Genet., 2006, 120 (3), pp. 334-353. DOI: 10.1007/s00439-006-0213-2.
Pakendorf B., Novgorodov I.N., Osakovskij V.L., Stoneking M. Mating patterns amongst Siberian reindeer herders: inferences from mtDNA and Y-chromosomal analyses. Am. J. Phys. Anthropol., 2007, 133 (3), pp. 1013-1027. DOI: 10.1002/ajpa.20590. Starikovskaya E.B., Sukernik R.I., Derbeneva O.A., Volodko N.V., Ruiz-Pesini E. et al. Mitochondrial DNA diversity in indigenous populations of the southern extent of Siberia, and the origins of Native American haplogroups. Ann. Hum. Genet., 2005, 69 (Pt 1), pp. 67-89. DOI: 10.1046/j.1529-8817.2003.00127.x. Statistica Software, version 7.1. (StatSoft Inc., Tulsa, OK, USA, 2005).
Tamm E., Kivisild T., Reidla M., Metspalu M., Smith D.G. et al. Beringian standstill and spread of Native American founders. PLoS One, 2007, 2 (9), p. e829. DOI: 10.1371/journal.pone.0000829. Torroni A., Sukernik R.I., Schurr T.G., Starikorskaya Y.B., Cabell M.F. et al. mtDNA variation of aboriginal Siberians reveals distinct genetic affinities with Native Americans. Am. J. Hum. Genet., 1993, 53 (3), pp. 591-608.
Underhill P.A., Poznik G.D., Rootsi S., Jarve M., Lin A.A. et al. The phylogenetic and geographic structure of Y-chromosome haplogroup R1a. Eur. J. Hum. Genet., 2015, 23 (1), pp. 124-131. DOI: 10.1038/ejhg.2014.50.
Y-DNA Haplogroup Tree 2019. URL: https://isogg.org/tree/ (date of the application - 20.06.2019).
Information about Authors
Agdzhoyan Anastasiya T, PhD in Biology;
ORCID ID 0000-0002-8776-2934; [email protected];
Bogunov Yurij V., PhD in Biology; [email protected];
Bogunova Anna A., PhD in Biology; [email protected];
Kamenshikova Evgeniya N, PhD in Pedagogy;
Zaporozhchenko Valery V., [email protected]; Pylev Vladimir Yu.; ORCID ID: 0000-0001-9541-8319; [email protected];
Korotkova Nadezhda A, ORCID ID 0000-0003-1760-4700; [email protected];
Utrivan Sergey A, ORCID ID 0000-0002-4459-7211;
Skhalyakho Roza A., PhD in Biology;
ORCID ID 0000-0002-4005-1887; [email protected]; Koshel Sergey M., PhD in Geography; ORCID ID 0000-0002-4540-2922; [email protected]; Balanovsky Oleg P., Doctor in Biology, Professor RAS; ORCID ID 0000-0003-4218-6889; [email protected]; Balanovska Elena V., Doctor in Biology, Professor; ORCID ID 0000-0002-3882-8300; [email protected].
