Балановский О.П. , Фролова С.А. , Пшеничнов А.С. , Кузнецова М.А.
advertisement
Основные черты митохондриального генофонда восточных славян Балановский О.П.1, Пшеничнов А.С. 1, Фролова С.А.1, Кузнецова М.А.1, Дибирова Х.Д. 1, Васинская О.А.1, Кошель С.М., Запорожченко В, Чурносов М.И., Атраментова Л.А., Утевская О., Евсеева И.В., Тегако О., Почешхова Э.А., Микулич А.И., Виллемс Р., Балановская Е.В.1 1 Медико-генетический научный центр РАМН, 115478, г. Москва, ул. Москворечье, д.1, факс 007(095)3240702, balanovsky@inbox.ru АННОТАЦИЯ Впервые по маркерам митохондриальной ДНК проведено подробное изучение обширных выборок из генофонда восточнославянских народов: образцы ДНК 2396 тщательно отобранных представителей коренного населения из 23 популяций русских, украинцев и белорусов охарактеризованы по изменчивости ГВС1 и набору филогенетически информативных SNP маркеров. Результаты статистического анализа показали большую степень сходства всех восточнославянских популяций друг с другом. К ним оказались очень близки и западные славяне, а также южнославянские и неславянские народы Балканского полуострова. Что же касается неславянских народов Восточной Европы, то популяции Волго-Уральского региона (волго-финские и тюркские народы), а также западно-финские народы демонстрируют отличия от генофонда восточных славян. Мы трактуем это как результат недавнего (в генетическом масштабе времени) появления восточных славян на своем нынешнем ареале. Показана гомогенность генофонда белорусов, при выявлении небольших различий между северными и южными белорусами. Обнаружено сходство генофонда украинцев с населением более северных территорий (Польши, Белоруссии и юга России). Северные русские популяции за счет интенсивной ассимиляции автохтонного населения отличаются от основной массы восточных славян. Ключевые слова: восточнославянские народы, митохондриальная ДНК, геногеография, русские, украинцы, белорусы, генофонд 1 ВВЕДЕНИЕ Комплексное изучение популяций восточных славян представляет важную фундаментальную проблему, новые возможности для решения которой появились благодаря бурному развитию ДНК-технологий. Изучение этногенеза и этнической истории народов традиционно проводится рядом взаимодействующих дисциплин – археологией, историей, лингвистикой, этнографией, физической антропологией. Популяционно-генетическое изучение изменчивости ДНК маркеров неотделимо от антропологических исследований, поскольку также опирается на биологическую изменчивость популяций, только имеет дело не с фенотипическими (внешними) признаками, а с ДНК маркерами. Молекулярногенетическое изучение популяционных генофондов с каждым десятилетием занимает все более важное место в системе наук о народонаселении и решает, на наш взгляд, две задачи. При решении первой – наиболее популярной - задачи сведения о генофонде становятся еще одним историческим источником. Например, сведения о генетическом сходстве двух соседних популяций свидетельствуют об их общем происхождении, а бесспорно доказанное генетическое сходство двух удаленных популяций является основанием для выдвижения гипотезы об имевшей место миграции. Вторая задача является уже собственной задачей популяционной генетики: изучить структуру генофонда и закономерности ее изменения. Логически эта задача является первостепенной, поскольку только после получения и осмысления собственных данных популяционная генетика может предоставлять свои результаты и выводы в пользование смежным наукам. В данном исследовании мы попытались решить обе задачи применительно к восточнославянским народам, подробное изучение которых является традиционным для российской науки. Но генетические данные о них на момент начала нашего исследования были представлены в мировой литературе лишь фрагментарно, и в результате огромный славянский мир, занимающий почти половину Европы, исчезал из современных знаний о генофонде населения. Из всего спектра генетических маркеров мы выбрали митохондриальную ДНК как один из наиболее эффективных инструментов современной генетики. Кроме того, мы провели сравнительный анализ и с имеющимися в литературе данными по второй информативной генетической системе – Y хромосоме, а также и генетическими системами третьего типа – аутосомными ДНК-маркерами. В ходе предварительной работы, потребовавшей многих лет и неоднократных экспедиций, нами собраны уникальные по объему (2396 человек) и географическому охвату коллекции биологических образцов от коренного населения - русских, украинцев и 2 белорусов. Экспериментальный анализ митохондриальной ДНК в созданных ДНК коллекциях выполнен во многом при поддержке Российского гуманитарного научного фонда и благодаря сотрудничеству с Эстонским биоцентром. Это позволило сделать митохондриальные генофонды восточнославянских народов одними из наиболее подробно изученных этнических генофондов мира. Анализ полученных данных проведен как статистическими, так и картографическими методами (созданный картографический Атлас опубликован на сайте нашего коллектива www.genofond.ru). Он позволил не только выявить строение генофонда, но и воочию увидеть его закономерности на географической карте. Полученные данные находят широкий спектр приложений - от эколого-генетического мониторинга и судебно-медицинской экспертизы до преподавательской деятельности популяризации научных исследований. и В данной же статье анализ полученных данных проводится в основном аспекте - исторической геногеографии. Его цель -приблизиться к пониманию того, какие исторические процессы в народонаселении (в первую очередь, соотношение исторических миграций и частичной изоляции популяций) сформировали наблюдаемую картину генофонда. МАТЕРИАЛЫ И МЕТОДЫ Изученные выборки В ходе исследования впервые проведено изучение генофонда восточных славян по большим выборкам из основных региональных групп русских, украинцев, белорусов. Нами исследованы тринадцать русских, пять украинских и пять белорусских популяций. В таблице 1 приведена характеристика изученных популяций. В ней (и в последующем анализе) две популяции из Западной Украины рассматриваются суммарно: увеличение объема выборки обеспечивает более высокую статистическую достоверность результатов. Суммарно нами проанализированы почти две с половиной тысячи образцов. 3 Таблица 1. Изученные популяции русских, украинцев и белорусов. условное название народ географическая локализация Бобруйск белорусы Могилевская обл, Бобруйский р-н Столин белорусы белорусы Брестская обл., Столинский р-н Гомельская обл., Светлогорский район Светлогорск объем выборки долгота широта источник 29,20 53,12 [11] 26,77 51,97 наши данные 29,72 52,62 наши данные 29,00 55,00 наши данные 24,00 52,00 наши данные 30,00 52,00 наши данные 40,17 44,50 наши данные 124 43,50 43,50 наши данные 76 38,83 61,83 [10] 46,53 63,43 наши данные 42,00 66,00 [17] 92 89 71 100 104 121 132 Витебск белорусы Витебская обл. Брест белорусы Брестcкая обл. Гомель белорусы Гомельская обл. Кубанские казаки русские Терские казаки русские Ошевенск русские Пинега русские казаки кубанские (Адыгея) казаки терские (КабардиноБалкария) Архангельская обл., Каргопольский рн, Ошевенск Архангельская обл., Пинежский район Поморы русские Архангельская обл. (поморы) 144 81 Русские Белгород русские Белгородская обл., Яковлевский, Красненский, Прохоровский районы 148 36,48 50,78 наши данные Ростов русские 111 39,70 47,20 [12, 14] Смоленск русские 147 32,88 53,93 наши данные Кострома (Унжа) русские Ростовская обл. Смоленская обл., Рославльский и Ершичский районы Костромская обл., Мантуровский и Межевской районы (р. Унжа) 44,77 58,33 наши данные Боровский русские Калужская обл., Боровский район 36,50 55,00 наши данные Барятинский русские Калужская обл., Барятинский район 34,53 54,00 наши данные Болховский русские Орловская обл., Болховский район 36,00 53,43 наши данные 37,25 51,87 наши данные 39,00 54,22 наши данные 40,53 54,40 наши данные 40,25 52,63 наши данные 35,75 50,40 наши данные 32,07 49,43 наши данные 79 70 75 76 62 82 Черемисиновский русские Курская обл., Черемисиновский р-н Михайловский русские Спасск-Рязанский русские Рязанская обл., Михайловский район Рязанская обл., Спасск-Рязанский район Петровский русские Украинцы Белгород украинцы Тамбовская обл., Петровский район Белгородская обл., Красногвардейский и Грайворонский районы Черкассы украинцы Черкасская обл. 95 179 Хмельницкая Украинцы западные украинцы Хмельницкая обл., Староконстантиновский район 179 27,30 49,70 наши данные украинцы Львовская и Ивано-Франковская обл. 157 24,00 49,25 наши данные ИТОГО 3 народа 21 регион в трех странах 86 76 2756 Эта огромная по меркам современной популяционной генетики выборка была дополнительно расширена за счет привлечения опубликованных данных других авторов. Сравнительно с объемами нашего исследования их немного - суммарно 360 образцов. Это: данные по двум северным русским популяциям (Belyaeva et al., 2003; Tonks et al., 2006), по населению Ростовской области (Kornienko et al., 2004; Richards et al., 2000) и белорусам 4 Бобруйска (Behar et al., 2006). С учетом этих данных общая анализируемая выборка по восточным славянам составляет 2756 образцов. В работах Б.А. Малярчука также проанализировано около 500 образцов из различных областей Европейской части России, но, к сожалению, недостаточно тщательная стратегия формирования этих выборок (анкетирование лишь по родителям или даже по месту рождения самого пробанда) делает их несопоставимыми с нашими данными. По сходным причинам мы не рассматриваем и 101 образец, опубликованный в работе Орехова (Orekhov et al., 1999). Две другие имеющиеся в литературе небольшие выборки (русские Магадана 50 чел, и русские Башкирии – 83 чел) не включены в наше исследование, поскольку представляют русские популяции, проживающие за пределами исторического русского ареала. Таким образом, нами проведен тщательный скрининг литературы и отобраны популяции, способные повысить информативность проводимого исследования, но в целом наша работа основана преимущественно на собственных результатах (90% от суммарной анализируемой выборки). Подчеркем, что принципиальной особенностью анализируемого материала является не столько его объем, сколько качество сформированных выборок и их география. Для выявления географической структурированности восточнославянского генофонда было необходимо охватить всю обширную территорию расселения восточных славян, для чего требуется анализ многочисленных популяций. Без детальной характеристики различных региональных групп русского, украинского и белорусского народов изучение их генофонда невозможно. Наше исследование проведено в пределах исторического ареала восточнославянских народов, и намеренно не включило русские популяции Урала, Сибири, Средней Азии, сформировавшихся относительно поздно: их генофонд отражает не столько историю формирования русского народа, сколько историю его позднейшей экспансии. Та же цель – выявить не сиюминутный портрет генофонда, а его исторически устойчивые черты – определила и подход к формированию выборок. В обследование включались только те люди, оба дедушки и обе бабушки которых родились на территории данного изучаемого региона и относили себя к данному народу. Этот принцип «анкетирования до третьего поколения» в последние годы становится общепринятым в мировой науке, а нашим коллективом применялся, начиная с наших первых же экспедиций на Кавказе в 90х годах. Также в обследование не включались родственники (как минимум до третьего поколения) – их заведомое генетическое сходство могло бы сместить и исказить получаемые генетические характеристики популяции. Во всех популяциях проводился забор венозной крови, который всегда сопровождался письменным информированным согласием. Столь строгие требования к сбору первоначального биологического материала 5 потребовали значительных и постоянных усилий для организации и проведения множества экспедиций, которые для восточных славян начались в 2000 году и проводятся нашим коллективом ежегодно по настоящее время. Лабораторный анализ Генотипирование митохондриальной ДНК выполнено на уровне общепринятых мировых стандартов: проанализирован первый гипервариабельный регион митохондриальной ДНК (ГВС1 мтДНК) и проведен анализ отдельных SNP маркеров кодирующей части мтДНК. Секвенирование ГВС1 выполнено в интервале от 16024 до 16400 позиции скорректированной последовательности мтДНК [8]. Определение основных гаплогрупп проведено по данным о SNP маркерах согласно общепринятой классификации [18], тогда как отнесение образца к филогенетической ветви внутри гаплогруппы выполнено по данным о характерных мутационных заменах в ГВС1. Генотипирование пяти русских, двух белорусских и пяти украинских популяций выполнено в Эстонском биоцентре О.П. Балановским и А.С. Пшеничным; генотипирование остальных популяций выполнено в Медико-генетическом научном центре РАМН С.А. Фроловой (секвенирование), О.А. Васинской, М.А. Кузнецовой и Х.Д. Дибировой (анализ SNP-маркеров). Секвенирование в обеих лабораториях проводилось одинаковым методом на секвенаторах Applied Biosystems 377/3130/3130xl с использованием наборов Applied Biosystems BigDye 1.1./3.1, следуя инструкциям производителя. Анализ SNP маркеров в Эстонском биоцентре выполнен методом ПЦР-ПДРФ, а в МГНЦ РАМН – частично методом ПЦР-ПДРФ, но в основном с использованием технологии TaqMan (на приборе для ПЦР в реальном времени AB 7900HT). Методы анализа данных Методы, примененные для анализа полученного обширного массива данных, можно подразделить на две группы: статистические и картографические. Статистические методы включили: расчет матрицы попарных генетических расстояний (по Нею) между изученными популяциями, выполненный с помощью программы DJ, созданной нашим коллективом, и визуализацию созданной матрицы методом многомерного шкалирования, выполненного с помощью пакета Statistica 6.0. Картографический анализ включил создание карт распространения основных гаплогрупп в ареале восточных славян и соседних популяций. Определение географических координат изученных популяций проведено с помощью электронного атласа Microsoft Encarta, а построение геногеографических карт – с помощью разрабатываемого нашим коллективом геногеографического программного пакета GeneGeo. Часть карт созданного Атласа размещена на сайте www.genofond.ru. 6 РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ Спектр гаплогрупп у восточных славян Нами изучен состав гаплогрупп мтДНК в 23 популяциях восточнославянских народов. Средние частоты гаплогрупп мтДНК у восточных славян представлены в таблице 2. Большое число обнаруженных гаплогрупп принято в первом приближении объединять в два класса:: восточно-евразийские гаплогруппы (наиболее распространенные у народов Восточной Азии и Сибири) и западно-евразийские (свойственные народам Европы, Передней и Южной Азии). Таблица 2. Частоты гаплогрупп мтДНК у восточных славян. Гаплогруппа A B C D F G Z ВосточноЕвразийские суммарно L M1 Восточные Русские славяне 0,001 0,001 0,001 0,002 0,005 0,007 0,004 0,009 0,001 0,003 0,002 0,003 0,003 0,006 0,017 0,003 0,002 0,030 0,01 0 preHV H HV J 0,000 0,284 0,015 0,052 0 0,40 0,02 0,08 T* T1 T2 I N1a N1b R K V W X 0,004 0,023 0,044 0,013 0,002 0,004 0,005 0,031 0,030 0,017 0,007 0 0,03 0,06 0,02 0 0,01 0 0,03 0,05 0,03 0,01 7 U1a U1b U2 U3 U4 U5a U5b U6 U7 U8 0,001 0,002 0,013 0,008 0,029 0,047 0,025 0,000 0,002 0,004 0 0 0,02 0,02 0,04 0,08 0,03 0 0 0,01 other ЗападноЕвразийские суммарно 0,010 0,01 Как следует из таблицы 2, восточные славяне целиком относятся к западному кругу популяций, поскольку западно-евразийские гаплогруппы составляют подавляющее большинство в их генофонде, а суммарная доля восточных гаплогрупп не достигает даже двух процентов (табл. 2). В пределах западно-евразийских популяций по спектру гаплогрупп можно выделить две основных региональных группы - переднеазиатские и европейские народы. При таком делении восточные славяне полностью примыкают к генофонду европейских популяций (что, например, очевидно по сравнительно невысокой частоте типично «переднеазиатской» гаплогруппы J). Далее в пределах Европы простое рассмотрение частот гаплогрупп неинформативно, поскольку все европейские народы практически идентичны по спектру и очень похожи по частотам гаплогрупп мтДНК: такая высокая гомогенность митохондриального генофонда Европы была даже предметом специальных публикаций [15, 16]). Поэтому для выявления взаимного сходства и различий отдельных европейских популяций необходимо проводить специальный анализ, что сделано нами с использованием как картографических, так и статистических методов. География митохондриальной ДНК в восточнославянских популяциях Одной из задач данного исследования было представить его результаты не только в виде данных, понятных узкому кругу специалистов, но и в виде универсального электронного картографического Атласа. Такой Атлас реализован на страницах нашего сайта www.genofond.ru, и содержит карты распространения основных гаплогрупп мтДНК в масштабе Евразии. Это делает полученные результаты доступными не только генетикам, но и представителям смежных наук о человеке (медицины, археологии, лингвистики, 8 физической антропологии), которые, в силу трудности чтения узкоспециальной генетической литературы для неспециалиста, обычно вынуждены черпать сведения о генофондах из устаревших и не всегда достоверных источников. Рассмотрим для примера только три карты из 35 карт, опубликованных на сайте, причем применительно к географическому масштабу данной работы приведем карты не для Евразии в целом, а для региона Восточной Европы (рис. 1-3). На всех картах изученные популяции восточных славян отмечены звездочками. При построении карт учитывались и данные по неславянским народам Восточной Европы (эти популяции показаны треугольниками). Такое рассмотрение восточных славян на фоне остальных генофондов Восточной Европы делает картографический анализ наиболее информативным. Карта распространения гаплогруппы H (рис. 1) показывает ее высокую частоту у восточных славян (более 35%) и резкое снижение к востоку (20% и ниже в районе Урала). Внутри восточнославянского ареала более южные популяции (украинцы, часть белорусов, самые южные русские популяции) характеризуются несколько более низкими частотами гаплогруппы H, чем северные белорусы и большинство русских популяций. Отметим, что без картографирования, опираясь лишь на табличные данные, обнаружить эту географическую закономерность было невозможно. Другой пример – распространение гаплогруппы С (рис. 2). Эта типичная «сибирская» гаплогруппа обладает противоположным генетическим ландшафтом. Частоты гаплогруппы С высоки за Уралом, снижаются в Волго-Уральском регионе и низки у восточных славян. Еще одна гаплогруппа V представлена на рис. 3. Она имеет более сложный характер распространения. Ее частота повышена в южных районах и снижена в северных. Исключением – и очень ярким - являются саамы, у которых частота этой гаплогруппы резко возрастает, достигая своего мирового максимума [2]. В последующем изложении мы будет пользоваться этими результатами, но важно отметить, что созданный Атлас имеет самостоятельную научную ценность: освещая структуру восточноевропейского генофонда, он позволяет изучать целый ряд вопросов, не затрагиваемых нашим исследованием. 9 Рис. 1. Карта распространения гаплогруппы H мтДНК в Восточной Европе 10 Рис. 2. Карта распространения гаплогруппы С мтДНК в Восточной Европе 11 Рис. 3. Карта распространения гаплогруппы V мтДНК в Восточной Европе Генетические взаимоотношения восточных славян и соседних групп Применение разработанной нашим коллективом базы данных, содержащей практически всю опубликованную информацию по генофондам народов Европы, позволило получить четкую картину их генетических взаимоотношений, представленную на рис. 4. При построении этого графика мы избегали его чрезмерной перегрузки точками-популяциями и потому намеренно оставили лишь те, которые имеют непосредственное отношение к этногенезу и этнической истории восточной ветви славян. Эти группы перечислены в легенде к рисунку и занимают обширный ареал от Балканского полуострова на юго-западе до Уральского хребта на востоке. Общая структура генофонда, выявляемая этим анализом, сводится к трем основным кластерам (группам генетически сходного населения): Приуральскому, Восточно- европейскому и Северно-европейскому. В Приуральский кластер вошли восточные финноугорские (марийцы, коми) и тюркские народы Приуралья (татары, чуваши). В Восточноевропейский кластер вошли западные и восточные славяне, а также народы Балканского региона. В Северно-европейский кластер вошли западные финно-угорские народы и 12 северные русские. Из народов Кавказа два (осетины и армяне) в целом близки к восточноевропейскому кластеру, но другие народы Кавказа оказались ближе к народам Передней Азии, генетически весьма удалены от генофонда Восточной Европы и потому не могут быть отображены в «генетическом масштабе» данного графика. Условные овалы на рис. 4 показывают примерное генетическое пространство, популяции которого относятся к каждому их трех перечисленных кластеров. Легко заметить, что три популяции не относятся ни к одному из них. Первая из них (белорусы Бобруйска) была изучена другими авторами и в частотах гаплогрупп «прочие» гаплогруппы достигают высокой величины - 17%, что видимо, связано с методическими особенностями генотипирования или со сложностями конвертации информации в базу данных. Это означает, что особость этой популяции – кажущаяся, и вызвана она, скорее всего, техническими причинами. Что касается татар-кряшен, то они генетически удалены от всех других изученных нами популяций татар, но приближены к русским, что, возможно, является результатом смешанных браков с единоверным русским населением в течение последних двух-трех столетий. Наконец, особое положение мордовской популяции не находит однозначного объяснения, хотя, по нашему мнению, вызвано особенностями стратегии формирования этой выборки авторами публикации [4]. 13 Рисунок 4. График многомерного шкалирования популяций Восточной Европы - русские Волго-Уральский регион (финно-угры, татары, чуваши) - белорусы западные финно-угры - украинцы западные славяне и балты - северные русские Балканский регион Scatterplot 2D Final Configuration, dimension 1 vs. dimension 2 2,5 OsN 2,0 1,5 Dimension 2 1,0 0,5 0,0 ChuT -0,5 Rbol Bbre BSvet Rspas Bgom RomR Arm TtrK Ubel Rbar Rche Rter BVit RBel Uche TtrO Rmih RSm RKub RomS TtrB Hung Cze Est Slov Mari BStol TtrKr LatvRunj Pol Bulg LitZ UW Uhme Chuv TtrMiRrostRborLitAFinS Mord TtrE Rosh Rpom KarRpetr Rpin -1,0 BBobr Komi -1,5 -2,0 -3,0 FinN -2,5 -2,0 -1,5 -1,0 -0,5 0,0 Dimension 1 14 0,5 1,0 1,5 2,0 2,5 Итак, по результатам анализа митохондриальной ДНК выявляется единая, генетически гомогенная балто-славянская общность, к которой тяготеют также и популяции Балканского региона. От этой общности отличаются финно-угорские и тюркские народы Восточной Европы, генофонд которых, в свою очередь, подразделяется на западную и восточную части. Сравнение расположения популяций в генетическом пространстве и на географической карте показывает сходство в основных чертах (чтобы увидеть это сходство, достаточно допустить, что «восток» в генетическом пространстве находится слева, а «север» - внизу). Лишь восточные славяне, которые географически должны бы тяготеть к Приуральскому кластеру, показывают сходство с расположенными намного западнее группами балтов и западных славян. На наш взгляд, такая картина могла сформироваться в ходе славянской колонизации Восточно-Европейской равнины: двигающееся с запада славянское население несло с собой «восточно-европейский» митохондриальный генофонд, который замещал ранее существовавший на этих территориях автохтонный «приуральский» генофонд. Вопрос об ассимиляции славянами автохтонного финно-угорского населения, многократно обсуждавшийся в антропологической литературе, также находит отражение в данных по митохондриальной ДНК. Речь идет о северных русских популяциях, которые более сходны не с основной массой русских, украинских и белорусских популяций, а с неславянскими народами севера Восточной Европы (финнами и карелами). Можно предполагать, что славянская колонизация на своих поздних этапах (Русский Север) в наибольшей степени включала в себя ассимиляцию местного населения, и в генофонде северных русских до сих прослеживается дославянский компонент. Это и приводит к тому, что их размещение на графике - промежуточное между восточнославянскими и неславянскими популяциями. Вопрос о том, что представлял этот дославянский субстрат, который впитали северные русские, не находит однозначного ответа. С одной стороны, разумно предполагать финно-угорский компонент, однако специальный подробный анализ, проведенный для одной из северно-русских популяций (популяция с верхней Пинеги, находящаяся на границе с республикой Коми [1]) показал отсутствие ярко выраженного сходства с соседними коми, но при этом дал указания на сходство русских Пинеги с более западными популяциями Фенноскандии, причем говорящими на языках не только финноугорской, но и германской группы [1]. Что же касается украинских популяций, то три из четырех изученных популяций (за исключением наиболее восточной) располагаются практически в той же области графика, 15 что и западно-славянские и балтоязычные популяции. Это согласуется с ранее установленным по картам генетических расстояний тяготением украинского митохондриального генофонда к популяциям, расположенным севернее (поляки, белорусы, южные русские) [5]. Генетические взаимоотношения русских, украинских и белорусских популяций Определив место восточнославянских групп в генофонде Европы, перейдем к рассмотрению внутренней структуры восточнославянского генофонда. Для этого построен аналогичный график многомерного шкалирования, но на него выведены только восточнославянские популяции (рис. 5). График наглядно восточнославянских показывает популяций сходство друг с митохондриального другом: все они генофонда образуют всех единый малодифференцированный кластер. Такая генетическая гомогенность восточных славян находит параллели в их антропологическом и лингвистическом сходстве друг с другом. Однако при более пристальном рассмотрении начинает вырисовываться и внутренняя структура генофонда: все-таки существующие различия между популяциями по мтДНК. Так, большинство белорусских популяций (кроме уже обсуждавшейся кажущейся особости белорусов Бобруйска и белорусов Столинского района) образуют скопление в левой части графика. Все три северные русские популяции (Ошевенск, Пинега, поморы), напротив, смещены вправо (особость северных популяций была выявлена нами и ранее на несколько отличающемся наборе исходных данных [3]). Шесть русских популяций образуют единое облако со всеми украинскими популяциями, а также белорусами и кубанскими казаками. Таким образом, именно эта группа оказывается ядром восточнославянского генофонда (к ней относятся 12 из 26 изученных популяций). Остальные шесть русских популяций занимают различное положение на графике, которое в большинстве случаев, видимо, объясняется статистическими случайностями, неизбежными при многомерном шкалировании большого числа популяций (поскольку двумерные геометрические расстояния никогда не могут в точности соответствовать реальным многомерным генетическим расстояниям). Так, парадоксально расположение тамбовских русских (Rpetr) рядом с северными; кажется странным положение Курской и Орловской групп (Rche, Rbol), а также одной из двух рязанских популяций (Rspas) поскольку вторая рязанская популяция вошла в основную восточнославянскую генетическую общность. Лишь особость русских Ростовской области (Rros) и терских казаков (Rter) легко может быть объяснима их периферийным положением на самом юге исторического русского ареала. Любопытно, что группа кубанских казаков (Rkub) неотличима от основного восточнославянского ядра. Аналогичное поведение тех же 16 групп было выявлено нами по маркерам Y хромосомы [9]. Итак, рассмотрение генетических взаимоотношений двадцати шести восточнославянских популяций по маркерам митохондриальной ДНК показывает: - выходящую на первый план гомогенность восточнославянского генофонда; - умеренное своеобразие большинства белорусских популяций, занимающих особое положение в генетическом пространстве (что объясняется, вероятнее всего, наличием балтского компонента в составе их генофонда); - выраженное своеобразие северных русских популяций; - возможное своеобразие терского казачества. Scatterplot 2D Final Configuration, dimension 1 vs. dimension 2 1,4 1,2 Rspas Rter 1,0 Bbre 0,8 Rrost BBobr 0,6 Dimension 2 0,4 Runj BSvet Rbar 0,2 0,0 Bgom RBel Uche -0,2 -0,4 Rbol UW Uhme Rbor RSm Rmih Ubel BVit Rpom BStol Rosh RKub -0,6 -0,8 Rpetr Rche -1,0 Rpin -1,2 -1,4 -1,5 -1,0 -0,5 0,0 0,5 1,0 Dimension 1 17 1,5 2,0 2,5 3,0 Рисунок 5. График многомерного шкалирования русских, украинских и белорусских популяций. - русские - украинцы - белорусы - северные русские Сравнение результатов по мтДНК с другими генетическими системами Сравнение результатов анализа восточнославянских популяций по мтДНК с ранее проведенным нами изучением изменчивости Y хромосомы показывает хорошее совпадение выявляемых картин: по обоим маркерам все восточнославянские популяции оказываются чрезвычайно сходными западнославянскими между собой популяциями. При (малодифференцированный этом по обеим системам кластер) и с прослеживается своеобразие северных русских популяций [9]. Различия в результатах, полученных по мтДНК и Y хромосоме, невелики и связаны главным образом с тем, что мтДНК в целом более гомогенна в Европе и поэтому в меньшей степени, чем Y хромосома, позволяет выявлять различия между генофондами. Например, по мтДНК южные славяне в целом присоединяются к остальной славянской общности, хотя в генетическом пространстве и располагаются на ее окраине. Однако по Y хромосоме южные славяне четко отделены от сходных между собою западных и восточных славян. Данные по маркерам третьего типа (панели 500 000 аутосомных ДНК маркеров [13]) также показывают отличия южных славян и большое сходство восточных и западных славян. Таким образом, незначительное расхождение в результатах между мтДНК и Y хромосомой вызвано лишь большей гомогенностью распределения гаплогрупп мтДНК. Сравнение с результатами проведенного нами анализа пяти аутосомных маркеров (ACE, APOA1, B65, PV92, TPA25) показывает сходство и с этой системой признаков [6, 7]. По аутосомным маркерам также вырисовывается единый гомогенный кластер восточных славян, несколько удаленный от коренного населения Кавказа, но значительно отличающийся от народов Волго-Уральского региона (Соловьева и др, в печати). ВЫВОДЫ Результаты статистического анализа показали большую степень сходства всех восточнославянских популяций друг с другом. К ним оказались очень близки и западные славяне, а также южнославянские и неславянские народы Балканского полуострова. Что же 18 касается неславянских народов Восточной Европы, то популяции Волго-Уральского региона (волго-финские и тюркские народы), а также западно-финские народы демонстрируют отличия от генофонда восточных славян. Мы трактуем это как результат недавнего (в генетическом масштабе времени) появления восточных славян на своем нынешнем ареале. Именно поэтому генофонд, например, русских в целом ближе к далеким западным и южным славянам, чем к соседним финно-угорским народам. При этом такая близость является относительной, и анализируемые данные показывают не только миграцию славян, но и ассимиляцию ими местного дославянского населения, в особенности на севере русского ареала. Другие, более частные результаты, касаются генофондов отдельных народов. Для белорусов показана гомогенность их генофонда, сквозь которую все же пробиваются различия между северными и южными белорусами. Для украинцев обнаружено сходство их генофонда с населением более северных территорий (Польши, Белоруссии и юга России) и некоторое своеобразие западных популяций украинцев. Для русских одним из важных выводов стало то, что северные русские популяции за счет ассимилированного населения отличаются от основной массы восточных славян и приближены, с одной стороны, к генофонду финно-угорских народов, а с другой, проявляют генетическое сходство и с популяциями Скандинавии. БЛАГОДАРНОСТИ: Авторы благодарят 2396 представителей коренного населения различных регионов России, Украины и Белоруссии, предоставивших образцы своей ДНК в ходе экспедиционных обследований, а также медицинский персонал и администрации этих регионов, активно помогавших в проведении этих экспедиций. Программа DJ для расчета генетических расстояний была написана Ю.А. Серегиным. Работа частично поддержана грантами РГНФ 06-06–00640а, РФФИ 07–06-00086а и РФФИ 07–04-00340а. 19 СПИСОК ЛИТЕРАТУРЫ 1. Балановская Е.В., Балановский О.П. Русский генофонд на Русской равнине. М.: Луч, 2007. 416 с. 2. Балановский О.П. Генофонд высоких широт Евразии или откуда пришли cаамы? // Путь на север: окружающая среда и самые ранние обитатели Арктики и Субарктики (материалы международной конференции) / Под ред. А.А. Величко, С.А. Васильева. Москва. Институт географии РАН. 2008. С. 277-282. 3. Балановский О.П., Пшеничнов А.С., Фролова С.А., Чурносов М.И., Атраментова Л.А., Тегако О.В., Микулич А.И., Кузнецова М.А., Васинская О.А. Изоляция северных русских популяций на фоне восточно-славянского единства: рассказывает митохондриальная ДНК // VIII Конгресс этнографов и антропологов России, г. Оренбург, 1 - 5 июля 2009 г. 4. Бермишева M. А., Тамбетс К., Виллемс P., Хуснутдинова Э. К. Разнообразие гаплогрупп митохондриальной ДНК у народов Волго-Уральского региона // Молекулярная биология. 2002. N 6. С. 990-1001. 5. Пшеничнов А.С. Структура генофонда украинцев по данным о полиморфизме митохондриальной ДНК и Y хромосомы. Автореф. дис. канд. биол. наук / ГУ Медикогенетический научный центр РАМН. М., 2007. 29 с. 6. Соловьева Д.С., Балановская Е.В., Кузнецова М.А., Васинская О.А., Фролова С.А., Почешхова Э.А., Евсеева И.В., Болдырева М.Н., Балановский О.П. Русский генофонд: геногеография ALU-инсерций (АСЕ, APOA1, B65, PV92, TPA25) // Молекулярная биология (на рецензировании). 7. Соловьева Д.С., Балановский О.П., Кузнецова М.А., Дибирова Х.Д., Дружинина Е.Г., Пузина Т.А., Сперанская А.С., Фролова С.В. Разнообразие генофонда восточных славян по пяти аутосомным ДНК маркерам (ACE, APOA1, B65, PV92, TPA25) // VIII Конгресс этнографов и антропологов России, г. Оренбург, 1 - 5 июля 2009 г. 8. Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA // Nat Genet. 1999 20 Oct;23(2):147. 9. Balanovsky O, S. Rootsi, A. Pshenichnov, T. Kivisild, M. Churnosov, I. Evseeva, E. Pocheshkhova, M. Boldyreva, N. Yankovsky, E. Balanovska and R. Villems R (2008) ‘Two sources of the Russian patrilineal heritage in their Eurasian context’ American Journal of Human Genetics Vol. 82(1):23650. 10. Belyaeva O., Bermisheva M., Khrunin A., Slominsky P., Bebyakova N., Khusnutdinova E., Mikulich A., Limborska S. Mitochondrial DNA variations in Russian and Belorussian populations // Hum. Biol. 2003. 75(5). Р. 647-660. 11. Behar DM, Metspalu E, Kivisild T, Achilli A, Hadid Y, Tzur S, Pereira L, Amorim A, QuintanaMurci L, Majamaa K, Herrnstadt C, Howell N, Balanovsky O, Kutuev I, Pshenichnov A, Gurwitz D, Bonne-Tamir B, Torroni A, Villems R, Skorecki K. The matrilineal ancestry of Ashkenazi Jewry: portrait of a recent founder event // American Journal of Human Genetics. 2006 Mar. Vol 78(3). P. 487-497. 12. Kornienko I.V., Ivanov P.L. Hypervariable region I-II data from the Russian Federation // URL: www.ncbi.nlm.nih.gov [на 11 ноября 2009]. 13. Novembre J, T. Johnson, K. Bryc, Z. Kutalik, A.R. Boyko, A. Auton, A. Indap, K.S. King, S. Bergmann, M.R. Nelson, M Stephens, and C.D. Bustamante (2008) ‘Genes mirror geography within Europe’, Nature Vol 456(7218):98-101. 14. Richards M., Macaulay V., Hickey E., Vega E., Sykes B., Guida V., Rengo C., Sellitto D., Cruciani F., Kivisild T., Villems R., Thomas M., Rychkov S., Rychkov O., Rychkov Y., Gölge M., Dimitrov D., Hill E., Bradley D., Romano V., Calì F., Vona G., Demaine A., Papiha S., Triantaphyllidis C., Stefanescu G., Hatina J., Belledi M., Di Rienzo A., Novelletto A., Oppenheim A., Nørby S., AlZaheri N., Santachiara-Benerecetti S., Scozari R., Torroni A, Bandelt H. J. Tracing European founder lineages in the Near Eastern mtDNA pool // Am. J. Hum. Genet. 2000. V. 67. P. 1251-1276. 15. Richards M, Macaulay V, A. Torroni, and H.J. Bandelt. (2002) ‘In search of geographical patterns in European mitochondrial DNA’, American Journal of Human Genetics Vol 71(5):1168-74. 16. Simoni L, F. Calafell, D. Pettener, J. Bertranpetit, and G. Barbujani (2000) ‘Geographic patterns of 21 mtDNA diversity in Europe’, American Journal of Human Genetics Vol 66(1):262-278. 17. Tonks S., Winney B.J., Evseeva, I. Comparison of sex-linked and autosomal markers in Orkney and other North European populations // URL: www.ncbi.nlm.nih.gov [на 11 ноября 2009]. 18. Torroni A, Achilli A, Macaulay V, Richards M, Bandelt HJ.Harvesting the fruit of the human mtDNA tree // Trends Genet. 2006 Jun;22(6):339-45. Epub 2006 May 4. Review. Principal features of the mitochondrial gene pool of Eastern Slavonic populations Balanovsky O1, Pshenichnov A1, Frolova S.1, Vasinskaya O.1, Dibirova Kh.1, Kuznetchova M.1, Koshel S.1, Churnosov M.1, Atramentova L.1, Utevskaya O. 1, Tegako O.1, Pocheshkhova E.1, Mikulich A.1, Villems R.2, Balanovska E.2 1 Research Centre for Medical Genetics, Russian Academy of Medical Sciences, 115478, Moscow, Moskvoreche str., 1, a fax +7(495) 3240702, balanovsky@inbox.ru THE SUMMARY We have performed the detailed study of gene pool of the Eastern Slavonic populations. 2396 samples of mitochondrial DNA (mtDNA) from finely selected indigenous individuals representing 23 indigenous Russian, Ukrainian and Belorussian populations have been analyzed by sequencing HVS1 and typing informative coding region SNPs. Statistical analysis demonstrated the great similarity of all Eastern Slavonic populations. This cluster virtually coincides with clusters of Western Slavs and non-Slavonic populations from Balkan region. The different clusters are formed by non-Slavonic population of East Europe (Volgo-Finnic, West-Finnic and Turkic speakers). This finding illustrates the historically well known Slavonic colonization of the East European plain. The Belorussian mitochondrial pool has been shown to be rather homogenous, though slight difference between northern and southern Belorussian populations does exist. The Ukranian mitochondrial pool is close to those of the areas situated north from Ukraine ( Poland, Belorussia, southern Russia). North Russian populations clearly differ from main corpus of East Slavonic populations, which might be caused by more significant assimilation processes in the Russian North. 22 Балановский Олег Павлович 115583, г. Москва, ул. Ясеневая, д.21/2, кв. 139 8-915-321-9399, balanovsky@inbox.ru Пшеничнов Андрей Сергеевич Фролова Светлана Александровна Кузнецова Марина Александровна Дибирова Хадижат Дибировна Васинская Ольга Алексеевна Кошель Сергей Михайлович Запорожченко Валерий Чурносов Михаил Иванович Атраментова Любовь Утевская Ольга Тегако Ольга Почешхова Эльвира Аслановна Микулич Алексей Игнатьевич Виллемс Рихард Балановская Елена Владимировна 23