Бесплатный автореферат и диссертация по биологии на тему

Филогения и филогеография кожанов рода Eptesicus восточной Палеарктики

ВАК РФ 03.02.04, Зоология

Автореферат диссертации по теме "Филогения и филогеография кожанов рода Eptesicus восточной Палеарктики"

На правах

005059216

АРТЮШИН Илья Витальевич

ФИЛОГЕНИЯ И ФИЛОГЕОГРАФИЯ КОЖАНОВ РОДА ЕРТЕБЮиБ ВОСТОЧНОЙ ПАЛЕ АРКТИКИ

03.02.04 - зоология

АВТОРЕФЕРАТ

диссертации на соискание учёной степени кандидата биологических наук

т

Москва 2013

005059216

Работа выполнена на кафедре зоологии позвоночных биологического факультета МГУ имени М. В. Ломоносова.

Научный руководитель:

кандидат биологических наук Банникова Анна Андреевна

Научный консультант:

кандидат биологических наук Крускоп Сергей Вадимович

Официальные оппоненты:

Ведущая организация:

доктор биологических наук Холодова Марина Владимировна

«Институт проблем экологии и эволюции им. А.Н. Северцова», заведующая кабинетом методов молекулярной диагностики.

доктор биологических наук Ильин Владимир Юрьевич «ФГБОУ Пензенский Государственный Университет» (ПГУ), заведующий кафедрой зоологии и экологии Зоологический институт РАН, Санкт-Перербург

Защита состоится 20 мая 2013г. в 17 часов на заседании диссертационного совета Д 501.001.20 при биологическом факультете МГУ им. М.В. Ломоносова, по адресу: 119991 Москва, Ленинские горы, МГУ, биологический факультет, ауд. М-1

Тел/факс: +7 (495) 939-43-09; E-mail: irbeme@mail.ru

С диссертацией можно ознакомиться в Фундаментальной библиотеке Московского государственного Университета им. М.В. Ломоносова

автореферат разослан апреля 2013г.

Ученый секретарь диссертационного совета доктор биологических наук, профессор И.Р. Беме

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность исследования. Стремительное усиление роли молекулярных исследований во всех отраслях биологической науки обусловливает новые достижения в области филогенетики, филогеографии, популяционной генетики и объясняет появление огромного разнообразия конкурирующих гипотез относительно эволюционной истории различных групп организмов. Одно из заметных последствий процесса экспансии молекулярных методов в зоологию - это дальнейшее развитие и усиление генетической концепции вида (Bateson, 1909; Dobzhansky, 1934; Muller, 1939; Bradley, Baker, 2001; Baker, Bradley, 2006), которой придерживается автор данной работы. Согласно генетической концепции, вид — это совокупность генетически совместимых популяций, которые генетически изолированы от других подобных совокупностей, а видообразование определяется как процесс накопления генетических различий между двумя филогенетическими линиями, приводящий к генетической изоляции двух генных пулов, имеющих каждый свою собственную эволюционную судьбу.

Понимание филогенетических взаимоотношений летучих мышей семейства гладконосых (Mammalia: Chiroptera: Vespertilionidae) всегда вызывало затруднения из-за очевидного эволюционного и экологического успеха группы (с точки зрения видового разнообразия и биогеографии) при ограниченности морфологической диверсификации. Параллельная или конвергентная эволюция многих морфологических признаков приводят к заблуждениям относительно эволюционной истории Vespertilioninae (Arnback-Christie-Linde, 1909; Ellerman, Morrison-Scott, 1951; Hill, Harrison, 1987; Hill, Topal, 1973; Koopman, 1975; Rosevear, 1962; Täte, 1942; Volleth, Heller, 1994; Zima, Horacek, 1985 и др.), что, в свою очередь, затрудняет построение естественной и объективной классификации семейства. Использование молекулярных методов для исследования семейства Vespertilionidae принесло новое понимание филогенетических взаимоотношений между существующими видами. Оно привело к выявлению большого числа критических форм, чей статус, вследствие противоречий между уровнем генетической и морфологической дивергенции, стал предметом дискуссий среди систематиков. (Mayer, Heiversen, 2001; Heiversen, Heller 2001; Ibanez et al., 2006; Matveev et al, 2005; Mayer et al., 2007; Garcia-Mudarra et al., 2008; Juste et al. 2010; Salicini et al., 2011;Salicinietal„ 2012).

Род настоящих кожанов (Eptesicus) - один из центральных родов семейства, особенно важен для понимания систематики Vespertilionidae в целом. Он обладает практически всесветным распространением, включает значительное число видов, некоторые из которых обладают обширными ареалами и выраженной географической изменчивостью (Simmons, 2005). При этом границы и состав рода неоднократно пересматривались (Heller, Volleth, 1984; Kitchener et al., 1987; Hill, Harrison, 1987; Csorba, Lee, 1999; Hoofer, Van Den Bussche, 2003; Roehrs et al., 2010; Roehrs et al., 2011), а родст-

венные взаимоотношения видов в течение длительного времени оставались противоречивыми (Morales et al., 1991).

Для понимания структуры рода важный вопрос представляют собой взаимоотношения широко распространённых в Палеарктике видов: позднего кожана Е. serotinus и северного кожанка Е. nilssonii. Ранее (Mayer, Helversen, 2001) было указано на невозможность разделения по мтДНК этих хорошо различимых морфологически и весьма далёких с точки зрения классической систематики видов в западноевропейской части ареала. Это привело к появлению двух альтернативных гипотез, до сих пор нуждающихся в проверке: (1) интрогрессивная гибридизация или (2) быстрая морфологическая эволюция одного из видов за время, недостаточное для полной сортировки мито-хондриальных линий. Одновременно была показана глубокая митохондри-альная гетерогенность Е. serotinus (Ibanez et al., 2006), однако, популяции Е. serotinus с территории Восточной Европы, Сибири и Средней Азии в генетическом отношении оставались не изученными. Генетические данные об относительно узкоареальных формах Е. bobrinskoi, Е. gobiensis, Е. bottae ognevi, населяющих аридные районы Средней и Центральной Азии (Кузякин, 1950; Крускоп, 2012) также отсутствовали.

Цель исследования состоит в оценке генетической изменчивости видов рода Eptesicus на основании нескольких независимых генетических локу-сов уделив отдельное внимание генетической структуре Е. serotinus ssp. и Е. nilssonii в восточноевропейской и азиатской частях ареала и установлении филогенетических отношений этих видов с близкородственными видами кожанов.

Для решения этих вопросов поставлены следующие задачи:

1. На основе данных анализа мтДНК оценить видовое и внутривидовое разнообразие рода Eptesicus, проанализировать географическую изменчивость широко распространенных видов.

2. Выявит!) несоответствия между принятой таксономической системой и наблюдаемой структурой генетического разнообразия рода Eptesicus.

3. На основании изучения митохондриальных и ядерных маркеров проанализировать уровень дивергенции и характер взаимоотношений между признаваемыми географическими формами Е. serotinus, ревизовать их таксономический статус.

4. На основании митохондриальных и ядерных генов проанализировать возможные причины сходства мтДНК позднего кожана Е. serotinus и северного кожанка Е. nilssonii в европейской части ареала.

5. На основе всех имеющихся генетических данных предложить гипотезу о филогенетических взаимоотношениях между видами рода Eptesicus.

Научная новизна. Предложена новая молекулярно-филогенетическая гипотеза взаимоотношений внутри рода Eptesicus на основе анализа митохон-дриальной и ядерной ДНК. Результаты работы указали на ненадежность филогенетических реконструкций в данной группе, основанных исключительно на мтДНК. На основе анализа как митохондриальных, так и ядерных генов

впервые получены данные в пользу конспецифичности Е. bobrinskoi и Е. gobiensis, выдвинуто предложение о необходимости ревизиии Е. bottae ognevi и Е. serotinus andersoni, показано, что выделение подрода Amblyotus не поддерживается генетическими данными. Новые данные о генетической изменчивости подвидов Е. serotinus не подтверждают видового статуса E.s. turcomanus. Впервые на основе митохондриальной и ядерной ДНК получены убедительные подтверждения древней межвидовой гибридизации в пределах изучаемого рода, что дает материал для интерпретации ряда других таксоно-мически неоднозначных ситуаций.

Теоретическая и практическая ценность работы. Проведенные исследования финансировались Российским Фондом Фундаментальных Исследований (№08-04-00029, №10-04-00683, №11-04-000200) и соответствуют задачам этого фонда. Работа имеет теоретическое значение для познания процессов формирования биоразнообразия и закономерностей микроэволюции. Результаты данного исследования вносят вклад в понимание многообразных и сложных факторов видообразования и явления интрогрессии у млекопитающих. Внутриродовая и внутривидовая систематика летучих мышей рода Eptesicus, базирующаяся на немногих морфологических признаках, теперь может быть проверена на основе молекулярно-генетических данных. Полученные новые данные по генетическому разнообразию и молекулярной систематике кожанов обусловливают необходимость таксономической ревизии рода Eptesicus и должны быть использованы в этих целях.

Результаты работы могут быть применены в практических мероприятиях, использующих данные таксономических исследований и сведения по инвентаризации фауны (при составлении кадастров, разработке стратегии и тактики охраны генофондов, проведении экологического мониторинга и т.п.)

Апробация работы. Основные результаты работы обсуждались на лабораторных семинарах и заседании кафедры зоологии позвоночных Московского Государственного Университета им. М.В. Ломоносова и доложены на конференциях: «Современные проблемы зоо- и филогеографии млекопитающих» (Пенза, 2009), «Целостность вида у млекопитающих (изолирующие барьеры и гибридизация)» (Петергоф, 2010), «IX Съезд Териологического общества при РАН» (Москва, 2011), 15thIRBC (Прага, 2010), ЕСМ '12 (Париж 2012).

Структура и объем диссертации. Диссертация изложена на 130 страницах машинописного текста и состоит из введения, обзора литературы, материала и методов исследования, результатов исследования, их обсуждения, выводов, заключения и библиографического списка. Работа иллюстрирована 7 таблицами, 35 рисунками. Список литературы включает 160 наименований, в том числе 145 иностранных источников.

СОДЕРЖАНИЕ РАБОТЫ Введение

Во введении обоснована актуальность темы, поставлены цели и задачи исследования.

Глава 1. Изменчивость и систематика и рода Eptesicus и проблема исторической гибридизации (обзор литературы).

В разделе I главы приведен обзор изменения представлений о составе и разнообразии рода Eptesicus и обсуждается современное состояние проблемы на основе морфологических и молекулярных данных. Раздел II посвящен проблемам интрогрессивной гибридизации в природе преимущественно на примере млекопитающих. В разделе III приводится обзор молекулярных маркеров, использовавшихся ранее для решения филогенетических проблем у рукокрылых и краткий обзор методик, применяющихся для изучения интрогрессивной гибридизации.

Глава 2. Материал и методы.

Материал исследования.

Оригинальным материалом для данного исследования послужили 114 особей летучих мышей, из них 66 особей Е. serotinus, 27 особей Е. nilssonii, 2 Е. bottae ognevi, 7 Е. gobiensis, 6 Е. bobrinskoi, 2 Scotimanes ornatus, 2 Arielulus circumdatus, а также по одной особи E. fuscus и E. furinalis. Для исследования использовали мышцы и органы зверьков, фиксированные в 96% этаноле. Заспиртованные тушки исследованных летучих мышей хранятся в коллекциях Зоологического музея МГУ им. М. В. Ломоносова (Москва) и Зоологического института РАН (Санкт - Петербург). Часть образцов в виде спиртовых препаратов тканей была любезно предоставлена Д.Г. Смирновым (1 ЦПУ, Пенза), В.А. Матвеевым, C.B. Крускопом (ЗМ МГУ, Москва), В. Росиной (Палеонтологический Институт РАН, Москва) и A.B. Борисенко (Королевский Музей естественной истории Онтарио), A.C. Голенко (Центр реабилитации рукокрылых, Киев). Образцы из окрестностей посёлка Мирное Туруханского района Красноярского края, посёлка Оленье, хутора Ерецкий, посёлка Клетский Волгоградской области и села Солодники Астраханской области, посёлка Нижнебаканский Краснодарского края, деревни Ляличи, Брянской области, из северо-западного Казахстана и Узбекистана, собраны мною лично в полевые сезоны 2007-2011 г.

Также были использованы последовательности из базы «GenBank» и данные проекта «BOLD», любезно предоставленные C.B. Крускопом.

Выделение, амплификация и секвенирование ДНК.

Для молекулярно-генетического анализа использовали митохондри-альный ген цитохрома b (cytö), три ядерных экзона: экзон 29 аполипопротеи-на В (АРО), dentin matrix acidic phosphoprotein gene (DMP), фрагмент гена активации рекомбинации (RAG2), и два интрона: интрон 2 с фрагментом экзона 3 ß-цепи гена предшественника тиротропина (Thy), 12-й интрон с фрагментом 13-го экзона неэритроидного бета-спектрина 1 (sptbn).

Геномную ДНК выделяли из фиксированных этанолом мышц и печени стандартным методом фенол-хлороформной депротеинизации (Sambrook и др, 1989). Для амплификации и секвенирования были использованы прайме-ры: для фрагментов гена cytè L14734 - (Ohdachi и др., 2001), L486, Н602, H15395_pip (Artyushin et al., 2009); для ядерных генов Thy и sptbn соответствующие праймеры из статьи Eick et al. (2005). Для RAG2 — RAG2 F1, RAG2_R2 (Baker, 2000). Для генов АРО и DMP1 (и дополнительно для sptbn и RAG2) были разработаны оригинальные праймеры. ПЦР амплификацию проводили на приборах My Cycler BioRAD и ТЕРЦИК. Автоматическое секвени-рование проведено на секвенаторе ABI 3100-Avant с использованием ABI PRISM® BigDyeTM Terminator v. 3.1 в лаборатории ЦКП «ГЕНОМ».

Для ограниченной выборки проводили inter-SINE PCR анализ с олиго-нуклеотидами mirl7, mill7, omirl7 и omill7 (Jurka 1995).

Обработка данных молекулярно-генетического анализа.

Полученные последовательности выравнивали вручную. Реконструкция филогенетических деревьев была выполнена по методу максимального правдоподобия (ML) и максимальной экономии (MP) в программах PAUP* version 4.0b 10 (Swofford, 1998) и TREEFINDER (Jobb, G. TREEFINDER version of March 2011. 2011. Munich, Germany). Также проводили байесовский анализ с помощью программы MrBayes 3.1.2 (Huelsenbeck and Ronquist, 2001; Ronquist and Huelsenbeck, 2003; Altekar et al., 2004). Схемы разбиения на партиции и модели для ML и В А были определены в PartitionFinder vi.0.1 (Lan-fear et al., 2012). Для построения супердерева использовался метод MRP (Baum, 1992; Ragan, 1992; Baum, Ragan, 1993) на основе результатов ML анализа по отдельным генам.

Реконструкцию аллельных фаз для генов Thy и sptbn проводили с использованием Байесова метода в программе PHASE 2.1 (Stephens et al., 2001; Stephens and Donnelly, 2003). Для полученных аллелей строили парсимони-альную сеть в программе ТС S vl.21 (Clement et al., 2000). Был проведён NJ дистанционный анализ матрицы реконструированных аллелей в программе MEGA 5 (Tamura, 2011). В среде R (R Core Team, 2012) проводили многомерное шкалирование матрицы дистанций для аллелей и усреднённых дистанций между особями.

Анализ главных координат по матрице NET-дистанций между популяциями был выполнен в NTSYSpc 2.20 (Rohlf, 2008).

Глава 3. Результаты

Анализ митохондриального гена cytb.

В итоговое выравнивание вошло 146 последовательностей. Его длина составила 707 п.н. после приведения к длине последовательностей, доступных в БД «GeneBank». Также для анализа использовали 66 полных последовательностей длиной 1140 п.н.

По результатам анализа митохондриальной ДНК выделяется не менее семи группировок (рис. 1). Последовательности Е. s. serotinus вместе с Е. s. turcomamis распадаются на две группы: (1) Е. s. serotinus с территории Запад-

ной Европы объединяются с образцами Е. nilssonii, причем Е. nilssoni пара-филетичен относительно западно-европейских Е. serotinus: в одну из его подгрупп попадают гаплотипы из Западной Европы и с запада Европейской части России, в другую - гаплотипы с Урала и из Сибири. В этот западноевропейский кластер nilssonii!'serotinus также попадают три из шести гапло-типов из Брянской области, один из Калининградской и три из четырёх из окрестностей Киева (2) Остальные последовательности Е. serotinus из Восточной Европы и все, полученные с территории Средней Азии образуют отдельную ("южную") ветвь.

'00 J id

1 E.bottae ognevi

E. serotinus

L*.

100 L

E.s.mirza E.b.hingstoni E.b.taftanimontis

E.gobiensis + E.bobrinskoi E.nilssonii Зап. Европа E.nilssonii Сибирь E.serotinus Зап. Европа E.s.andersoni

Amblyotus

100

I кожаны Нового Света

-4 E.isi

О '

■ E.hottentottus

I Внешняя группа

Рис.1. Филогенетическое древо, построенное в результате Байесова анализа по данным гена cyt Ь (707 п.н.). В узлах указаны апостериорные вероятности для соответствующих ветвей.

Е. s. mirza является сестринской группой к "южной" ветви Е. serotinus.

Е. s. andersoni является сестринской группой Е. nilssoni и отделён значительными дистанциями от остальных гаплотипов Е. serotinus ssp.

Комплекс Е. bottae не монофилетичен: последовательности Е. anatolicus (ранее Е. Ъ. anatolicus) оказываются базальными относительно всех палеарктических кожанов, за исключением Е. isabellinus, подвиды Е. Ъ. ognevi, Е. b. hingstoni+E. b. taftanimontis образуют сестринскую группу с Е. s.

яегоНп^+Е. ж. ппгга. Дистанции между Е. Ь. ognevi, и ветвью Е. Ь. Ит^от+Е. Ь. 1сф.аттопи$ высоки и достигают видового уровня (7,7%). Положение монофилетической группировки Е. %оЫет1$+Е. ЪоЪптког внутри клады европейских кожанов не определено, причем оба вида в составе этой группировки не монофилетичны, их гаплотипы разделены крайне низкими дистанциями.

Кожаны Нового Света (Е. /тст, Е. /,игтаШ), кожаны из Южной Испании (Е. [ьаЪеШпш), и южноафриканский вид Е. ЪоПепШт образуют три самостоятельные ветви, на митохондриальном дереве отделившиеся от общего ствола раньше палеарктических Ер1е$и:и$ эрр.

Для уточнения положения ветви кожанов Нового Света был проведён дополнительный анализ объединённых данных для митохондриальных генов су\Ъ и N01 (последний из БД «СепеВапк»), который показал наиболее раннее отделение этой ветви в радиации группы. Отделение Е. ЬоМепШш — наиболее ранняя дивергенция среди кожанов Старого Света.

с

о

□=L

п

0.05 0.10

о - дистанция

Рис. 2. Гистограмма распределения попарных невзвешенных р-дистанций (ген cyt b)

Распределение попарных р-дистанций (рис. 2) имеет выраженный мультимодальный характер. Первая группа (0-4%) соответствует внутривидовым дистанциям, также в эту группу попадают дистанции между западноевропейской группой Е. serotinus и обеими ветвями Е. nilssonii и дистанция между Е. gobiensis и Е. bobrinskoi. Во вторую (4-11 %) попадают дистанции между большинством видов, кроме Е. isabellinus, Е. hottentottus и видов Нового Света. Третью группу составляют дистанции между Е. isabellinus, Е. hottentottus, видами Нового Света и остальными группами. Дистанции между Е. s. mirza и "южной" гаплогруппой Е. serotinus не укладываются в диапазон внутривидовых и находятся на нижней границе межвидовых.

Анализ ядерных данных THY

Для анализируемой выборки (24 экз. Е. nilssonii, 6 - Е. gobiensis, 4 - Е. bobrinskoi, 54 — Е. serotinus, 5 - Е. s. andersoni), программа PHASE реконструировала наибольшее количество аллелей для Е. nilssonii —11, и для Е. serotinus - 10. Из трёх аллелей Е. gobiensis - два общие с Е. bobrinskoi. Три аллеля из пяти у Е. bobrinskoi - уникальные. У Е, s. andersoni найдено 3 алле-ля. Топология парсимониальной сети, построенной по реконструированным аллелям (рис. 3), лишена циклов и большинство связей соответствует единичным заменам. Аллели Е. serotinus составляют два подграфа (гаплогруппы ТА и ТВ), которые связаны между собою и с группой гаплотипов Е. nilssonii через гаплотип Е. s. anderssoni.

I ЕодпВебЗТ |

Рис. 3. Парсимониальная сеть, отображающая разнообразие аллелей THY. Частоты аллелей не указаны, каждое ребро соответствует одной замене. Числа над ребрами графа показывают положение замен в выравнивании вариабельных позиций.

Для анализа географической изменчивости гена THY выборка Е. serotinus была подразделена на 8 географических подвыборок. Распределение частот встречаемости аллелей демонстрирует переход от доминирования группы ТА на востоке исследованной территории к преобладанию группы ТВ на западе (рис. 4).

Метод главных координат (РСА), применённый к матрице NET-дистанций между географическими подвыборками Е. serotinus располагает их

вдоль оси (ТНУ_е1), соответствующей долготному распределению локалите-тов (рис. 5). Расстояние в пространстве главных координат между восточными подвыборками ("Астрахань", "Казахстан") и подвыборкой "Волгоград" оказывается непропорционально большим относительно географических дистанций.

УКР БРЯ КРС САМ ВОЛ ACT КАЗ УЗБ

Рис. 4. Географическое распространение гаплогрупп THY. Ось категорий - географические подвыборки (слева направо: «Узбекистан», «Казахстан», «Астрахань», «Волгоград», «Самара», «Краснодар», «Брянск», «Украина»), Ось значений - относительная частота гаплоплогрупп. 0.008 см 0.004 0.000 I- -0.004 -0.008 -0.012

0.03 0.02 0.01 0.00 -0.01 -0.02 -0.03 -0.0'

THY е1

Рис. 5. Размещение выборок в главных координатах по данным NET-дистанций между последовательностями гена THY. Обозначения подгрупп соответствуют обозначениям на рис. 4.

SPTB

Для анализируемой выборки (Е. nilssonii - 18 экз., Е. serotinus - 47 экз., Е. bobrinskoi и Е. gobiensis - по 6 экз.) в программе PHASE было реконструировано 23 аллеля Е. nilssonii, 30 аллелей Е. serotinus - 4 и 5 аллелей для Е. bobrinskoi и Е. gobiensis соответственно, причём один аллель Е. gobiensis совпадает с аллелем Е. bobrinskoi.

Парсимониальная сеть гаплотипов SPTB имеет более сложную топологию, чем для гена THY, и отличается от неё наличием циклических связей (рис. 6). Граф распадается на несколько несвязанных сегментов: сегмент, включающий большинство аллелей Е. nilssonii, три отдельных сегмента, состоящих из аллелей Е. serotinus (SA, SB, SC), два сегмента для Е. bobrinskoi+E. gobiensis. Эти сегменты соответствуют слабо- и умеренно-поддержанным кластерам на NJ дереве (не приводится).

-11-

in- О

ьт

КРС

о

УКР о

САМ

Рис. 6. Парсимониальная сеть, отображающая разнообразие аллелей SPTB. Частоты аллелей не указаны, каждое ребро соответствует одной замене. Числа над ребрами графа показывают положение замен в выравнивании вариабельных позиций.

Географическое распределение (рис. 7.) аллелей из гаплогрупп Е. serotinus так же демонстрирует переход от преобладания одной из групп (SC) в восточных популяциях к преобладанию других (SA, SB) в западных.

УКР БРЯ КРС САМ ВОЛГ ACT КАЗ УЗБ

Рис. 7. Географическое распространение гаплогрупп SPTB. Ось категорий - географические подвыборки, ось значений - относительная частота гаплогрупп.

Расположение выборок при анализе последовательностей интрона SPTB в пространстве двух осей хуже соответствует географическим дистанциям между подвыборками, чем таковое для интрона THY. Однако остаётся верным утверждение об относительно удалённом расположении друг относительно друга восточной ("Казахстан", "Узбекистан", "Астрахань") и западной ("Украина", "Краснодар", "Брянск") группы выборок. Выборка "Волгоград"

оказывается внутри западной группы, что противоречит её географической близости к подгруппе "Астрахань" (рис. 8.).

0.020

Я о.ою

0)

m' 0.000 н

§5 -0.010 -0.020

С

БРЯ УК р

О хз кр<

зол у; )Ь

КАЗ

САМ

О

0.10

0.08

0.06

0.04

0.02 0.00 SPTB е1

-0.02

-0.04

-0.06

-0.08

Рис. 8. Размещение выборок в главных координатах по данным ЫЕТ-дистанций между последовательностями гена ЭРТВ. Обозначения подгрупп соответствуют обозначениям на рис. 4.

№ анализ с бутстрэп-процедурой матрицы ЫЕТ-дистанций (рис. 9) с высокой степенью статистической поддержки восстанавливает восточную ("Казахстан", "Узбекистан", "Астрахань") и западную ("Украина", "Краснодар", "Брянск") группу подвыборок, объединяя "Волгоград" с западной.

971

99

981

88

98

■ КАЗ

■ УЗБ ' ACT

■ САМ

■ БРЯ

■ ВОЛ

■ УКР

■ КРС

Рис. 9. NJ дерево по матрице NET-дистанций между географическими подвыборками Е. serotinus. В узлах указаны бутстрэп поддержки (1000 реплик). Обозначения терминальных ветвей соответствуют кодам географических подвыборок.

Анализ данных по пяти ядерным генам.

В объединённое выравнивание по ядерным генам (4058 п.н.) вошли 56 последовательностей. (АРОВ — 705 п.н., DMP1 — 962 п.н., RAG2 — 1285 п.н., SPTBN — 544 п.н., Thy — 562 п.н.)

Все деревья, построенные методами MP, ML, MB на объединённых данных по ядерным генам, содержат 5 основных клад:

1. Кожаны Нового Света, включая род Histiotus (клада "fuscus") ;

2. Е. serotinus ssp. вместе с Е. s. andersoni и Е. isabellinus, ранее входивший в качестве подвида в Е. serotinus (клада "serotinus");

3. Е. bottae ognevi и Е. hottentotus - (клада "hottentotus/bottae")

4. Группа Е. gobiensis/E. bobrinskoi - клада "gobiensis";

5. Е. nilssonii - клада "nilssonii";

Две последние клады, не объединяются в общий кластер, однако разделены наименьшими генетическими дистанциями.

Общий анализ оригинальных и доступных в GeneBank (Juste et al 2010) последовательностей гена RAG2 поддерживает объединение Е. Ъ. bottae и Е. anatolicus в общую группу, а также, хотя и с низкими поддержками, но объединяет Е. serotinus ssp.

В целом филогенетические отношения между найденными на ядерном дереве пятью группами не являются разрешёнными.

84/79/1/72

74/81/0.99/92

62/54/0.95/69

82/95/0.99/58

97/95/1/93

—E.s.andersoni

98/0.97/1/100

Е.serotinus

97/95^1/93

E.gobiensis 4 E.nilssonii

98/99/1/100

E.isabellinus

96/95/1/77

98/97/1/88

E.bottae ognevi

E.hottentotus

кожаны Нового Света

100/100/-/10jfr

Внешняя группа

Рис. 10. Филогенетическое древо, построенное методом ML по конкатенированным ядерным данным (3.9 тыс. п.н.). В узлах указаны поддержки ML/MP/BPP/бутстрэп поддержки для супердерева (MRP).

Глава 4. Обсуждение

Митохондриальная филогения рода.

Полученная нами топология митохондриального филогенетического дерева (рис. 1) не полностью соответствует принятому таксономическому делению рода на уровне видов.

Общая структура дерева и распределение попарных дистанций заставляют предположить два этапа радиации в пределах рода. Во время первого из них происходит обособление южно-африканской (Е. hottentotus) и американской (Е. fuscus, Е. furinalis) ветвей, а также Е. isabellinus, от клады, включающей остальные виды рода; во время второго этапа последняя разделяется на 7-8 линий, соответствующих формам, относимым к Е. serotinus, Е. bottae, Е. nilssonii, Е. gobiensis.

Тем самым, наши митохондриальные данные поддерживают монофи-лию палеарктических кожанов, при исключении Е. isabellinus. По данным дополнительного анализа двух митохондриальных генов (ND1 и cytb) ветвь, объединяющая кожанов Нового Света, отделяется несколько раньше других, что было ранее продемонстрировано при использовании других молекулярных маркеров (Stadelmann 2004, Roehrs et al 2010, Roehrs et al. 2011).

Расположение вида E. nilssonii внутри основной радиации рода на митохондриальных филогенетических деревьях не согласуется с традиционным пониманием общей морфо-эволюционной тенденции кожанов, выражающейся в большей редукции I2 и М3 у более морфологически продвинутых видов. Исходя из большей морфологической дифференциации I2 и М3 у Е. nilssonii, Е. bobrinskoi/E. gobiensis, иногда выделяемых в подрод (или даже род) Amblyotus, эти таксоны должны были бы занимать базальное положение на филогенетическом дереве. Более того, для северного кожанка описано факультативное присутствие в зубном ряду малого предкоренного зуба (Кру-скоп, Фадеева, 2008), что, несомненно, является крайне примитивной чертой в контексте рода Eptesicus. Несомненная архаичность морфологического облика кожанков даже заставляла исследователей, не признававших Amblyotus в качестве отдельного таксона, выделять их в особую группу видов (напр., Koopman, 1994). Тем не менее, по данным анализа митохондриальных генов положение северного и гобийского кожанков далеко от базального.

Этому можно предложить два объяснения. Первое из них связано с вероятностью древней интрогрессии мтДНК предка Е. serotinus в геном предка современных Amblyotus. Другое объяснение предполагает ошибочность современных представлений о преобладающих тенденциях в эволюции морфологических признаков кожанов.

Филогения рода по ядерным данным и события интрогрессии в эволюции рода.

Ядерные данные противоречат митохондриальной филогении рода и лучше последней соответствуют традиционным взглядам на систематику кожанов, поддерживая существование видового комплекса Е. bottae s.l. (включая anatolicus) и Е. serotinus s.l. (включая формы isabellinus и andersoni). Оче-

видные несоответствия митохондриальной и ядерной филогении касаются, во-первых, комплекса видов Е. bottae (включая Е. anatolicus) и, во-вторых, положения Е. s. andersoni.

Положение на митохондриальном дереве видов комплекса Е. bottae отличается от положения на дереве по яДНК. На дереве по объединенным ядерным генам Е. Ъ. ognevi объединяется с Е. hottentottus, в то время как на митохондриальном дереве он ближе к Е. serotinus. Кроме того, дополнительный анализ данных по гену RAG2 показал, что и другие представители видового комплекса Е. bottae, включая Е. anatolicus, объединяются при высоких поддержках, Можно предположить, что в эволюции этого видового комплекса митохондриальный геном был независимо заимствован как минимум трижды.

Е. s. andersoni, будучи близким по ядерным данным к Е. serotinus, по митохондриальным занимает обособленное положение и отличается от Е. serotinus примерно в той же степени, что от Е. nillssoni и Е. gobiensis (-7%), Это так же свидетельствует о возможных интрогрессивных событиях в недавней истории таксона.

Присутствие на дереве по объединенным ядерным генам высоко поддерживаемой клады «hottentotus/bottae» и клады, объединяющей Е. isabellinus, Е. s. serotinus, Е. s. andersoni и Е. s. mirza (последнее обнаруживается в анализе расширенного набора таксонов по гену RAG2), противоречит традиционной трактовке группы видов «serotinus» (Simmons 2005; Koopman, 1993): виды комплекса Е. bottae s.l. следует исключить из него. Вопрос о правомочности выделения подрода Cnephaeus Каир, 1829, объединяющего палеарктических кожанов и южно-африканского Е. hottentotus (Hoofer & Van Den Bussche, 2003), остаётся открытым до установления филогенетических отношений между кладами «hottentotus/bottae» и «fuscus».

Еще одно вероятное противоречие митохондриального и ядерного деревьев касается положения Е. isabellinus. На митохондриальном дереве положение этого вида не определяется наверняка, однако он оказывается внешней ветвью относительно других палеарктических Eptesicus, в то время как по результатам анализа ядерных генов Е. isabellinus образует сестринскую ветвь к Е. serotinus, хотя и с низкими поддержками. Для уточнения его филогенетического положения требуется дальнейшее накопление данных по яДНК.

Взаимоотношения Е. isabellinus и Е. serotinus

Выявленные зарубежными авторами различия в мтДНК между поздними кожанами Западной Европы и Пиренейского полуострова позволили выделить Е. isabellinus в качестве самостоятельного вида. Предположение, что и эта форма, и поздние кожаны Европейской части России - популяции одного вида, избежавшие, в отличие от западно-европейских, недавней ин-трогрессии мтДНК Е. nilssonii, не подтвердилось (Ibanez et al., 2006). Результаты проведенного нами анализа демонстрируют наличие значимых генетических различий между Е. isabellinus и Е. serotinus как из Западной Европы, так и из Европейской части России и Казахстана. Это поддерживает выделе-

ние Е. isabellinus в качестве отдельного вида и указывает на ошибочность представления о синонимичности Е. s. turcomanus и Е. s. isabellinus (Ellerman, Morrison-Scott, 1966).

Положение E. bobrinskoi и E. gobiensis

Как ядерные, так и митохондриальные данные продемонстрировали отсутствие различий видового уровня между формами Е. gobiensis и Е. bobrinskoi. Таким образом, видовая самостоятельность Е. bobrinskoi нашими данными не подтверждается. Вероятно, кожанка Бобринского следует рассматривать как подвид гобийского кожанка, от которого он отличается более мелкими размерами (Artyushin et al., 2012).

Е. bottae

Внутри этого видового комплекса традиционно признается наличие выраженной географической изменчивости (Ellerman, Morrison-Scott, 1966), однако таксономический статус выделяемых форм неоднозначен, в частности, статус формы anatolicus недавно был повышен до видового (Mayer et al., 2007). По данным нашего анализа митохондриальной ДНК не только Е. anatolicus образует далеко отстоящую от других палеарктических кожанов ветвь, но и остальные формы видового комплекса Е. bottae не составляют монофилетической группы. Выделяется два отчетливых кластера, один из которых образован Е. b. ognevi, второй — Е. b. hingstoni, Е. b. innesi, Е. Ь. taftanimontis. Различия между этими двумя группами подвидов достигают видового уровня (р-дистанция 7.7%). По-видимому, статус Е. b. ognevi должен быть пересмотрен на более обширном материале.

Статус Е. serotinus andersoni

Восточные формы Е. serotinus, населяющие Китай (кроме Синдзяня), Корею, север Индокитая, хорошо обособлены от номинативного подвида по данным анализа мтДНК, и, возможно, их таксономический статус при изучении более обширного материала, может быть повышен до видового. В этом случае вид должен получить название Е. andersoni Dobson, 1871.

Имевшиеся в нашем распоряжении пробы поздних кожанов из Кореи, Тайваня и Лаоса группируются в компактный кластер с низкими внутригруп-повыми дистанциями. Однако надёжность группирования Е. s. andersoni вместе с Е. serotinus по данным яДНК и относительно низкие внутригрупповые дистанции свидетельствуют об их достаточно недавней дивергенции.

Интрогрессия мтДНК Е. nilssonii в геном Е. serotinus

На территории Западной Европы обитают популяции весьма несхожих морфологически видов Е. nilssonii и Е. serotinus, при этом крайне слабо различающиеся по мтДНК (1% р-дистанции (Mayer, Helversen, 2001)). Для объяснения этого предлагались две гипотезы: недавнее видообразование с неполной сортировкой линий и быстрой (относительно накопления генетических различий) морфологической эволюцией, либо интрогрессия мтДНК между Е. nilssonii и Е. serotinus (Mayer et al., 2007).

В настоящей работе обнаружено, что территорию Юго-Восточной Европы, Средней и Центральной Азии населяют поздние кожаны с митохондри-альными гаплотипами, принадлежащими к особой группе («южная» гапло-группа), которые не родственны гаплотипам Е. nilssonii и Е. serotinus из западно-европейских популяций. Гаплотипы «южной» группы отделены от гаплотипов Е. nilssonii и западно-европейских гаплотипов Е. serotinus дистанциями, характерными для видового уровня в данной группе (-7.5% ). На территории юго-восточной Европы, в Поволжье, на Северном Кавказе и нижнем Дону, рапространены гаплотипы «южной» группы; в Брянской области и окрестностях Киева присутствуют гаплотипы обеих групп, западнее известны только гаплотипы сходные с Е. nilssonii. Анализ ядерной ДНК показывает, что Е. nilssonii и Е. serotinus обладают уникальными видоспецифичными пулами аллелей. При этом все популяции Восточной Европы (подвид Е. s. serotinus) по яДНК крайне сходны и образуют единую группу вне зависмости от преобладающего варианта мтДНК. Таким образом, ядерные гены не поддерживают гипотезу о недавнем разделении Е. nilssonii и Е. serotinus. Более вероятной представляется гипотеза интрогрессии мтДНК Е. nilssonii в популяции Е. serotinus в результате гибридизации в зоне исторической симпатрии.

По данным анализа мтДНК, западно-европейские формы Е. nilssonii и Е. serotinus образуют весьма близкие, но обособленные монофилетичные группы. Можно предположить, что интрогрессия имела место в прошлом, но не продолжается в настоящее время. В пользу этого говорит отсутствие в настоящее время существенного перекрывания тех частей ареалов этих видов, в которых происходит размножение. При этом Е. nilssonii свойственна выраженная географическая изменчивость — популяции с территории Сибири и Европейские популяции образуют отдельные кластеры с высокой статистической поддержкой. Большее сходство европейского гаплотипа Е. nilssonii с заимствованным в ходе интрогрессии гаплотипом Е. serotinus позволяет предположить более раннее разделение форм Е. nilssonii, чем событие интрогрессии митохондриального генома Е. nilssonii в популяции Е. serotinus.

Соответствие распределения маркеров ядерной ДНК границам видов так же не подтверждает гипотезу современной гибридизации между Е. serotinus и Е. nilssonii.

Проблема взаимоотношений Е. s. turcomanus и номинативного подвида

Традиционно считается, что на юге Европейской части России обитают два хорошо дифференцированных подвида позднего кожана (Огнёв 1927): западный Е. s. serotinus и восточный Е. s. turcomanus («кожан туркменский» или «туранский»). Между этими формами на основании морфологических и генетических данных предполагается гибридизация (Стрелков, Ильин, 1990; Матвеев, 2002).

Результаты анализа мтДНК показывают, что экземпляры из Астраханской области, Казахстана и Узбекистана, морфологически определенные как Е. s. turcomanus, не объединяются в общий кластер.

Анализ ядерных маркеров в нашей работе не позволил выявить однозначных различий между Е. serotinus s.str. и Е. s. turcomanus, однако, разделение по частотам аллелей и NET-дистанциям популяций Е. serotinus на «западные» («Украина», «Брянск») и «восточные» («Узбекистан», «Казахстан», «Астрахань»), а так же существование на деревьях по ядерным генам групп, соответствующих тому же делению, указывает на дифференциацию внутри Е. serotinus. «Восточные» подвыборки Е. serotinus преимущественно включают в себя зверьков с типичным габитусом Е. s. turcomanus, «западные» же наоборот — Е. s. serotinus. В промежуточных по частотам аллелей подвыборках «Краснодар», «Самара», «Волгоград» также доминируют звери схожие по морфологическому описанию с Е. s. serotinus.

Такая ситуация может быть объяснена двумя различными гипотезами. Согласно первому сценарию диверсификация подвидов Е. serotinus - это эво-люционно недавнее событие и тогда ядерные данные отражают результат неполной сортировки генетических линий, а различия в морфологических признаках — результат быстрой морфологической дифференциации.

Второй сценарий предполагает вторичный контакт двух изолированных в прошлом форм с замещающей интрогрессией мтДНК и образованием клины частот ядерных аллелей. Такой сценарий представляется более вероятным, поскольку хорошо согласуется с низким уровнем различий по мтДНК (0.3% р-дистанция по eyt6) по сравнению с яДНК (0.3% и 1.6% р-дистанции между географическими подвыборками по THY и SPTB соответственно).

Филогеография Е. nilssonii

Разделение с высокими индексами бутстрэпа и умеренными генетическими дистанциями, не достигающими видового уровня, европейской и сибирской форм Е. nilssonii, позволяет констатировать относительно давнюю дивергенцию и требует пересмотра традиционного объединения этих форм в единый подвид Е. п. nilssonii.

Ни в одной из проанализированных популяций обе гаплогруппы Е. nilssonii не встречаются одновременно. Популяции Е. nilssonii, в которых присутствуют гаплотипы сибирской группы, распространены не только на территории Сибири и на Урале, но и на Европейской части России (Самарская обл.). Такое положение филогеографической границы достаточно необычно для бореальных видов млекопитающих, для которых филогеографи-ческие границы сосредоточены обычно к востоку от Урала.

Данные по яДНК на приведённом в этой работе материале не выявляют группировок внутри Е. nilssonii, что можно трактовать как более низкий уровень дифференциации этого вида, по сравнению с Е. serotinus.

ЗАКЛЮЧЕНИЕ

Митохондриальные данные успешно используются как инструмент для предварительной оценки видового разнообразия для многих групп рукокрылых и позволяют выдвигать гипотезы относительно филогенетической истории и систематики этих групп. Однако целый ряд процессов и в частности возможность интрогрессивной гибридизации приводят к несоответствию ге-ногенеалогии мтДНК истинной генеалогии для некоторых групп. Наше исследование показывает, что в эволюции рода Eptesicus интрогрессивные процессы происходили многократно и привели к несоответствию ядерной и ми-тохондриальной генеалогии. Вместе с тем митохондриальные данные оказались безусловно полезными для инвентаризации разнообразия (в случае с Е. bobrinskoi, Е. anatolicus и Е. isabellinus) и поиска направлений для пересмотра систематики рода. Филогенетические отношения между найденными на ядерном дереве пятью группами не являются разрешёнными, отчасти вследствие достаточно быстрого разделения некоторых из этих линий за короткий промежуток времени в период базальной радиации рода, отчасти по причине недостаточного количества данных по яДНК. Увеличение объёма ядерных данных позволит уточнить эволюционную историю рода и определить более точно число и порядок событий интрогрессии.

ВЫВОДЫ

1. Полученные молекулярные данные поддерживают выделение следующих Палеарктических видов рода Eptesicus: Е. nilssonii, Е. serotinus, Е. bottae, Е. gobiensis, Е. isabellinus, Е. anatolicus. Таксономический статус хорошо дифференцированных по мтДНК форм Е. bottae ognevi и Е. s. anderssoni подлежит ревизии на более обширном материале. Видовая самостоятельность Е. bobrinskoi не подтверждается, эту форму следует рассматривать в качестве подвида Е. gobiensis.

2. Относительная генетическая близость Е. s. serotinus и Е. s. turcomanus объясняется потоком генов при вторичном контакте либо недавней дифференциацией, которая ещё не привела к накоплению различий в мтДНК. Независимо от принимаемой гипотезы, возведение Е. s. turcomanus в ранг вида не оправдано.

3. Существуют достоверные различия по мтДНК между популяциями Е. nilssonii с территории Сибири и Европы, не достигающие, тем не менее, видового уровня. С помощью использованных ядерных маркеров эта дифференциация не выявлена.

4. Сходство мтДНК популяций Е. serotinus из Западной и Центральной Европы с мтДНК Е.nilssonii - следствие исторической интрогрессии митохондриальных гаплотипов Е. nilssonii в популяции западноевропейских Е. serotinus. Дивергенция европейской и азиатской гаплогрупп Е.

-20-

nilssonii, по-видимому, предшествовала событию интрогрессии мтДНК Е. nilssonii в популяции Е. serotinus.

5. В роде Eptesicus существуют 5 филогенетических линий: линия включающая кожанов Нового Света, а так же представителей рода Histiotus\ линия «hottentotus/bottae» (Е. hottentotus, комплекс Е. bottae, включая Е. anatolicus); линия «serotinus» (Е. serotinus ssp., включая Е. s. andersoni, Е. isabellinus)-, линия «gobiensis» (Е. gobiensis, включая Е. bobrinskoi)', линия «nilssonii» (Е. nilssonii). Эти данные не указывают на обособленное положение мелких архаичных форм Е. nilssoni и Е. gobiensis (подрод Amblyo-tus).

6. Филогенетическая схема, построенная по данным анализа мтДНК некорректно отражает генеологические связи внутри рода Eptesicus, вследствие неоднократных событий интрогрессии митохондриального генома, сопровождавших эволюцию этой группы.

БЛАГОДАРНОСТИ

Автор глубоко признателен всем, кто поддерживал меня и помогал при выполнении этой работы: моим научным руководителям и учителям — Анне Андреевне Банниковой, Сергею Вадимовичу Крускопу; Лебедеву Владимиру Святославовичу за консультации и неоценимую методическую помощь; Бак-лушинской Ирине Юрьевне за рецензирование предварительной версии работы; Борисенко Алексею Владимировичу, Росиной Валентине Владимировне, Матвееву Виталию Александровичу, Голенко Анне Сергеевне, Смирнову Дмитрию Григорьевичу за любезно предоставленные образцы тканей кожанов; коллегам по кабинету молекулярных методов кафедры зоологии позвоночных биологического факультета МГУ; коллективу стационара «Мирное» ИПЭЭ им. А.Н. Северцова; родственникам и друзьям, которые оказывали моральную, материальную и деятельную поддержку при сборе материала и написании работы.

СПИСОК РАБОТ, ОПУБЛИКОВАННЫХ ПО ТЕМЕ ДИССЕРТАЦИИ

Статьи в рецензируемых журналах из перечня ВАК и других изданиях

Artyushin I. V., Bannikova A. A., Lebedev V. S. and Kruskop S. V. Mitochondrial DNA relationships among North Palearctic Eptesicus (Vespertilioni-dae, Chiroptera) and past hybridization between Common Serotine and Northern Bat.// Zootaxa 2262: 40-52, 2009.

Artyushin I. V., Lebedev V. S., Smirnov D. G., and Kruskop, S. V. 2012. Taxo-nomic position of the Bobrinski's serotine (Eptesicus bobrinskoi, Vesper-tilionidae, Chiroptera). Acta Chiropterologica, 14(2), 291-303.

Artyushin I., Lebedev V., Bannikova A. and Kruskop S. 2012. A molecular reappraisal of the taxonomic status of Eptesicus serotinus turcomanus// Vespertilio 16: 31-42.

Тезисы конференций

И.В. Артюшин, А.А. Банникова, B.C. Лебедев, С.В. Крускоп, 2009. Следы гибридизации между Е. serotinus и Е. nilssonii и родственные взаимоотношения палеарктических видов кожанов по данным митохондри-альной ДНК // Современные проблемы зоо- и филогеографии млекопитающих: мат. Всерос. научной конфенции Пенза: ПГПУ. С 8.

С.В. Крускоп, И.В. Артюшин, 2010. Молекулы и морфология в систематике рукокрылых: расхожде-ние во взглядах // Целостность вида у млекопитающих (изолирующие барьеры и гибридизация) (Петергоф: 2010) Материалы конференции. М.: Товарищество научных изданий КМК. С 46.

И.В. Артюшин, А.А. Банникова, С.В. Крускоп, B.C. Лебедев, 2011. Кожаны рода Eptesicus (Vespertilionidae) Северной Палеарктики: молекулярная систематика и филогеография // териофауна россии и сопредельных территорий. Международное совещание (IX Съезд Териологического общества при РАН). М.: Товарищество научных изданий КМК. С 29.

Artyushin I., Bannikova A., Lebedev V., Kruskop S., Mitochondrial DNA Relationships among North Palearctic Eptesicus species // 15th IRBC - the Conference Manual. Prague, Czech. 2010. P 90.

Artyushin I., Bannikova A., Lebedev V., Kruskop S., North palearctic serotines, molecular study // Vlth European Congress of Mammalogy (Abstract volume). Paris, France. 2011. P 77.

Подписано в печать 18.04 Формат А5 Бумага офсетная. Печать цифровая. Тираж ЮОэкз. Заказ № 1688 Типография ООО "Ай-клуб" (Печатный салон МДМ) 119146, г. Москва, Комсомольский пр-т, д.28 Тел. 8(495)782-88-39