Влияние островной изоляции и эффекта "бутылочного горлышка" на генетический полиморфизм командорских песцов

Плошница, Анна Ивановна

Бесплатный автореферат и диссертация по биологии на тему
Влияние островной изоляции и эффекта "бутылочного горлышка" на генетический полиморфизм командорских песцов
ВАК РФ 03.02.04, Зоология

Автореферат диссертации по теме "Влияние островной изоляции и эффекта "бутылочного горлышка" на генетический полиморфизм командорских песцов"

На правах рукописи

ПЛОШНИЦА Анна Ивановна

ВЛИЯНИЕ ОСТРОВНОЙ изоляции и ЭФФЕКТА «БУТЫЛОЧНОГО ГОРЛЫШКА» НА ГЕНЕТИЧЕСКИЙ ПОЛИМОРФИЗМ КОМАНДОРСКИХ ПЕСЦОВ (АЮРЕХЫвОРиЗ ВЕШМСЕКШ И А. Е ЯЕМЕАЮУ!)

03.02.04-зоология 03.00.03 - молекулярная биология

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата биологических наук

Москва-2010

004603500

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

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

Гольцман Михаил Ефимович

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

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

доктор биологических наук Холодова Марина Владимировна Институт проблем экологии и эволюции им. А.Н. Северцова РАН

кандидат биологических наук Семенова Серафима Константиновна Институт биологии гена РАН

Зоологический институт РАН

Защита диссертации состоится 24 мая 2010 г. в 15:30 на заседании диссертационного совета Д 501.001.20 при Московском государственном университете им. М.В.Ломоносова по адресу: 119991, Москва, Ленинские горы, д. 1, стр. 12, Биологический факультет, ауд. М-1. Факс: 8 (495) 939 17 46; эл. почта: tadorna@mail.ru

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

Автореферат диссертации разослан 22 апреля 2010 г.

Ученый секретарь диссертационного совета, кандидат биологических наук

А.Б. Поповкина

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

Актуальность темы. Дифференциация популяций и видообразование на океанических островах служит предметом изучения зоологов, начиная с Дарвина и Уоллеса (см., например, MacArthur, Wilson, 1967, van der Geer et al., 2010). Изоляция, простые экосистемы, небольшие размеры популяций и высокая скорость эволюционных изменений делают островные виды привлекательными для исследователей. Однако островные популяции особенно уязвимы к любым средовым изменениям, и гораздо беззащитнее материковых перед антропогенными угрозами (Aleover et al., 1998, Frankham, 1998). Исследование изменений, происходящих при сокращении небольших изолированных популяций, тем более актуально, что все больше фрагментированных материковых популяций превращается в «островные» из-за непреодолимых антропогенных преград. Радикальные изменения и исчезновения популяций происходят не на эволюционных временах, а на протяжении лишь десятков лет, а иногда значительно быстрее. Процессы, происходящие на этих временах, то есть временах поколений, изучены очень слабо (Johnson, 2003). Морфологические индикаторы слишком инерционны, а экологические характеристики, напротив, слишком изменчивы и бывают редко известны для более чем десяти поколений.

Командорские песцы представляют исключительно удобную модель для исследования этих процессов. Популяции изолированы на океанических островах, по крайней мере, в течение нескольких десятков тысяч лет, на краю ареала, в условиях относительно простой экосистемы, но с исключительно стабильными и богатыми кормовыми ресурсами (Гольцман и др., 2003; Goltsman et al., 2005).

Одна из командорских популяций, популяция о-ва Медный в 1970-х годах прошла через демографическое «бутылочное горлышко» из-за внезапной эпизоотии ушной чесотки (Goltsman et al., 1996). После этого популяция восстановилась до 1015% былой численности и стабильно сохраняет этот уровень в течение 30 лет. Сейчас эта популяция занесена в Красную Книгу РФ, в категорию редкости 1, как находящаяся под угрозой исчезновения (Гольцман, Крученкова, 2001). История командорских песцов хорошо документирована, а в музее хранится коллекция черепов, добытых около 100 лет назад. Кроме того, в начале прошлого века были предприняты успешные попытки интродукции командорских песцов на острова Ушишир Курильской гряды (Воронов, 1963; Костенко и др., 2004). Это дает уникальную возможность испытать применимость молекулярно-генетического подхода в исследовании популяционных процессов в масштабе времени десятков поколений, и изучить влияние островной изоляции и прохождения популяции через «бутылочное горлышко».

в формировании адекватного иммунного ответа на разнообразные патогены (Klein, 1986; Parham & Ohta, 1996; Apanius et ai, 1997; Ploegh & Watts 1998; Hedrick, 1999; Sommer, 2005). Предполагается, что балансирующий отбор, действующий в генах ГКГ и направленный на поддержание высокого уровня разнообразия, способен снижать эффект дрейфа генов и замедлять скорость фиксации аллелей (Richman, 2000).

Цель и задачи. Целью работы было исследовать влияние островной изоляции и эффекта «бутылочного горлышка» на нейтральный и адаптивно значимый генетический полиморфизм командорских песцов.

Для достижения этой цели были поставлены следующие задачи:

1.. Описать нейтральный полиморфизм D-петли мтДНК и микросателлитов у песцов предковой популяции о-ва Медный и микросателлитов у современных островных и континентальных песцов.

2. Определить полиморфизм адаптивно значимых последовательностей DRB и DQB генов II класса ГКГ у предковых и современных островных, а также континентальных песцов.

3. Оценить интенсивность балансирующего отбора в антигенсвязывающем регионе генов II класса ГКГ.

4. Оценить влияние островной изоляции и демографического «бутылочного горлышка» на полиморфизм адаптивно значимых генов II класса ГКГ и нейтральных маркеров песцов.

5. На основании анализа мтДНК выяснить происхождение песцов о-ва Янкича (о-ва Ушишир Курильской гряды).

Научная новизна. Впервые описан полиморфизм DRB и DQB генов II класса ГКГ у песца (Alopex lagopus). Получены данные о разнообразии нейтральных (D-петля мтДНК и микросателлиты) и селективных (гены II класса ГКГ) маркеров в исторической популяции песцов о-ва Медный до последнего демографического «бутылочного горлышка». Впервые исследовано влияние прохождения популяции через «бутылочное горлышко» на генетическое разнообразие генов ГКГ на основе прямого сравнения данных по разнообразию в современных образцах с таковым в исторических образцах животных, добытых до спада численности. С помощью молекулярно-генетических методов проведено сравнение командорских популяций песцов между собой и с материковыми популяциями. Исследован уровень полиморфизма в генетических маркерах, находящихся под действием разных селективных сил. Найдены генетические подтверждения того, что современные песцы о-ва Янкича происходят с о-ва Медный.

Теоретическое и практическое значение результатов работы.

так и в селективных маркерах, чем меньшая популяция о-ва Медный. Полученные данные важны для понимания того, как демографические события (эффект «бутылочного горлышка») отражаются на разнообразии нейтральных генетических локусов и локусов, находящихся под действием балансирующего отбора. Понимание микроэволюционных процессов в генах ГКГ, которые определяют способность организма формировать иммунный ответ, - важный этап в исследовании эволюции иммунной системы. Полученные данные являются важными сведениями об адаптивных способностях островных популяций и имеют большое значение в определении дальнейших путей их охраны.

Апробация работы. Основные результаты диссертации доложены на международных симпозиумах «Сохранение разнообразия животных и охотничье хозяйство России» (Москва, 2006), «Молекулярно-генетические основы сохранения биоразнообразия млекопитающих Голарктики» (Черноголовка, 2007), а также на 87ой ежегодной встрече Американской Ассоциации Маммологов (87th annual meeting of American Society of Mammalogists - Albuquerque, USA, 2007), на 6°" и 7™ международных симпозиумах по физиологии, поведению, генетике и охране диких животных (International Zoo and Wildlife Research Conference on Behaviour, Physiology and Genetics - Berlin, Germany, 2007, 2009) и на 1ой междунородной конференции, посвященной генетике сообществ (1st Community Genetics meeting - Manchester, UK, 2010). Результаты работы обсуждались на лабораторных семинарах кафедры зоологии позвоночных биологического факультета МГУ им. М. В. Ломоносова и в исследовательской группе иммуногенетики IZW в Берлине.

Публикации. По материалам диссертации опубликовано 11 научных работ, в том числе 5 статей, из них 1 в журнале из списка ВАК.

Объем и структура работы. Текст диссертации изложен на 114 страницах и включает: введение, обзор литературы, материалы и методы, результаты исследования и их обсуждение, выводы и список литературы. Работа иллюстрирована 10 рисунками и 10 таблицами. Список литературы включает 228 источников, в том числе 200 - на иностранных языках.

Благодарности. Я признательна всем, кто помогал на разных этапах этой работы. Прежде всего, хочу выразить огромную благодарность моему научному руководителю М.Е.Гольцману, а также моим руководителям в лабораториях университета Манчестера (Великобритания) LJ.Kennedy и заведующей исследовательской группой иммуногенетики IZW (Германия) S.Sommer за чуткое руководство, терпение, внимание и помощь в работе. Хочу особо поблагодарить М.В. Холодову, А.А.Банникову, А.А.Колесникова и Н.А.Формозова за ценные рекомендации, полезные замечания и помощь в редакции работы. Я признательна всем тем, кто помогал в организации и участвовал в экспедициях на Командорские острова: С.В.Загребельному, Н.Н.Павлову, Л.О.Дорониной, А.Н.Шиенку. Особая благодарность всем тем, кто осуществлял помощь в проведении лабораторных анализов, а также давал ценные рекомендации и советы, делился опытом: Е.Л.Джикии, П.А.Сорокину, В.А.Орехову, С.В.Найденко, А.С.Карнаухову. Большое

спасибо коллективу зоологического музея МГУ им. М.В.Ломоносова: О.Г.Нановой, С.В.Крускопу и В.С.Лебедеву за доброжелательное отношение, интерес к исследованию и предоставление исторических образцов костей. Хочу поблагодарить

A.Ю.Косорукову и участников проекта "Arctic predators as indicators of tundra ecosystem" за предоставление современных образцов из Сибири. Я особенно признательна D.W .Macdonald и K.Jewgenow за помощь в организации поездок в Великобританию и Германию для проведения там лабораторных работ, а так же

B.К.Боженко за возможность работы в лаборатории ФГУ РНЦРР (Москва). И, конечно, отдельно я хочу поблагодарить моих родных: маму Т.С.Плошницу, брата И.И.Плошницу, тетю Н.С.Захарову, а также моего молодого человека А.de Palma Garrido за помощь и поддержку.

Исследование финансировали: Peoples Trust for Endangered Animals to DWM; DAAD-project Leonard-Euler: A08/01104 и 50077303; Sigma-Xi, Grants-in-Aid: G200631112503926; РФФИ: 07-04-007-45a.

СОДЕРЖАНИЕ РАБОТЫ

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

ГЛАВА 1. ОБЗОР ЛИТЕРАТУРЫ

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

Во второй части обзора рассматриваются эволюционные силы, действующие _на кодирующие и некодирующие последовательности: описывается взаимодействие факторов, снижающих генетическое разнообразие в популяции, и сил балансирующего отбора, направленного на поддержание полиморфизма в генах ГКГ. Рассмотрены исследования естественных популяций, свидетельствующие о влиянии дрейфа генов на полиморфизм. Описаны случаи, когда при малой популяционной численности балансирующий отбор способен поддерживать разнообразие генов ГКГ, несмотря на утрату полиморфизма нейтральных маркеров.

В третьей части обзора описывается специфика биологии популяций песцов на Командорских островах, самых древних изолированных эндемичных тихоокеанских популяций песца. Описывается история популяций после открытия островов, характер основных антропогенных воздействий и проблемы, возникшие в связи с распространением на о-ве Медный эпизоотии ушной чесотки и прохождением популяции через демографическое «бутылочное горлышко».

ГЛАВА 2. МАТЕРИАЛЫ И МЕТОДЫ

2.1. Сбор образцов

В настоящей работе были использованы образцы из популяций песцов Командорских островов (и=137), острова Янкича (и=3), континентальной Аляски (п=24) и Сибири («=33).

Современные образцы тканей песцов о-ва Медный («=41) были получены во время ежегодного мечения (июнь-август) с 1997 по 2006 гг. Выборка проб от современных песцов о-ва Беринга (п=44) включает фрагменты шкур животных, добытых во время зимнего промысла 2001 года, которые любезно предоставлены сотрудниками Командорского природного биосферного заповедника (КПБЗ) C.B. Загребельным и М.В. Жалгаубаевым, а также пробы тканей, собранные автором во время экспедиций 2008 и 2009 гг. Образцы тканей песцов с о-ва Янкича (о-ва Ушишир) собраны в 2007 г. и любезно предоставлены М.Е. Гольцманом и Е.П. Крученковой. Генетические пробы современных континентальных песцов Сибири предоставлены А.Ю. Косоруковой, образцы ДНК песцов Аляски - Л.Дж. Кеннеди (L.J. Kennedy, University of Manchester, UK). Собранные во время отлова и мечения мягкие ткани уха песцов, а также ткани обнаруженных мертвыми или убитых во время охоты животных фиксировали в 96% этиловом спирте.

Исторические образцы (черепа) командорских песцов были любезно предоставлены сотрудниками зоологического музея МГУ им. М. В. Ломоносова (Москва). Для анализа использовали образцы предковых песцов о-ва Медный (я=43), добытых с 1911 по 1946 гг. Кроме этого, были проанализированы 9 исторических образцов с о-ва Беринга, добытых в тот же период времени.

2.2. Выделение ДНК

2.2.1. Выделение современной ДНК из мягких тканей

Современную ДНК из мягких тканей песцов выделяли фенол-хлороформным методом (Sambrook et al., 1989) с модификациями. Перед процедурой экстракции ДНК фрагменты ткани промывали в бидистиллированной воде и измельчали скальпелем. Добавляли 650 мкл буфера NET100, 65 мкл саркозила натрия 35% концентрации и протеиназы К до конечной концентрации 200-400 ед. Смесь инкубировали 4 часа при 65°С и 9-12 часов при 37°С, после инкубации добавляли 1/5 часть объема 5 M хлорида натрия, после чего проводили последовательную экстракцию ДНК сначала равным объемом фенола, потом смесью хлороформа с изоамиловым спиртом (24:1) и в заключении - хлороформом. ДНК осаждали 1/10 объемом 3 M ацетата натрия и 2,2 объемами охлажденного этилового спирта. Осадок ДНК дважды промывали 70% этиловым спиртом и высушивали. ДНК растворяли в стерильной бидистиллированной воде, а полученный раствор хранили при -20°С.

2.2.2. Выделение исторической ДНК из костной ткани

Историческую ДНК выделяли из нижнечелюстных костей (mandibulare) и зубов песцов, хранившихся в Зоологическом музее МГУ им. М. В. Ломоносова. Распил и выделение проводили в специальной лаборатории. Кости тщательно мыли, после чего обрабатывали ультрафиолетовым излучением в течение не менее 60 мин. Распил проводили, используя стерильные стоматологические боры. Костную пудру хранили до выделения ДНК при -20°С.

Историческую ДНК выделяли по методу, предложенному Yang et al. (1998) и модифицированному Холодовой и др. (2005) с использованием набора mini elute PCR purification kit (Qiagen). К каждому образцу костной пудры объемом примерно 300400 мкл добавляли 700 мкл Е буфера (ЕДТА - 0,5 М, Tris-HCl - 10 мМ, саркозил Na -3,5%, pH 8,5) и 200-300 ед. протеиназы К, инкубировали 24 часа при температуре 55°С и постоянном перемешивании. После инкубации к надосадочной жидкости добавляли 5 объемов (примерно 3,5 мл) РВ буфера (поставляется в наборе), тщательно перемешивали и пропускали через колонку (поставляется в наборе), центрифугируя при 14 000 об/мин. в течение 1 мин. После этого в колонку добавляли 40 мкл нагретого до 50°С ЕВ буфера (поставляется в наборе). ДНК с колонки снимали центрифугированием. Процедуру проводили дважды. Полученный раствор ДНК хранили при -20°С.

2.3. Амплификация фрагмента D-петли мтДНК

Для амплификации полного I гипервариабельного региона D-петли мтДНК длиной примерно 300 нуклеотидных пар (н.п.) использовались праймеры PexlF и H3R (Dalén et al., 2002). Если из-за деградации исторической ДНК не удавалось получить ПЦР-продукт всего фрагмента целиком, использовали дополнительные праймеры PexlR, Pex2F, pex2R и pex3R (Nyström et al., 2006), а также CytF (Джикия и др., 2007), которые позволяют получить ПЦР-продукты трех перекрывающихся фрагментов исследуемого участка ДНК меньшей длины (примерно 85-155 н.п.).

Полимеразную цепную реакцию (ПЦР) фрагмента D-петли проводили в объеме 20 мкл, содержащем 5 мкл исторической ДНК (2 мкл ДНК в случае современных образцов), 0,3 мМ смеси дезоксирибонуклеотидтрифосфатов (dNTP's, Fermentas), однократный ПЦР-буфер без детергентов (Qiagen), 2,25 ед. Hotstar Taq полимеразы (Qiagen), прямой и обратный праймеры по 0,5 мкМ, 1,8 мМ MgCl2 (Qiagen) и деионизованную воду. Реакцию проводили в следующих условиях: предварительная денатурация 95°С - 15 мин. Далее следовало 43 цикла (30 для современной ДНК) со стандартной денатурацией при 95°С - 30 с; отжигом праймеров на протяжении 30 с при 50°С и элонгацией 40 с при 72°С. После окончания программы проводили заключительный цикл достройки при 72°С в течение 5 мин. Для оценки концентрации, специфичности полученного ПЦР-продукта и отсутствия контаминации проводили электрофорез в 2% агарозном геле в присутствии бромистого этидия.

2.4. Определение последовательности ДНК

Перед анализом нуклеотидной последовательности проводили очистку 1-5 мкл полученного ПЦР-продукта с помощью 10 ед. щелочной фосфатазы теленка (CIAP, Fermentas) и 2 ед. экзонуклеазы I (Exol, Fermentas), инкубируя при 37°С в течение 15 мин., потом при 85°С - 15 мин. Секвенирование каждого фрагмента проводили как минимум дважды (с прямого и обратного праймеров) используя набор для секвенирования BigDye® Terminator vl.l или v3.1 Cycle Sequencing Kit (Applied Biosystems). Для секвенирования применяли те же праймеры, что и для амплификации фрагментов. Реакцию секвенирования проводили в объеме 20 мкл, в состав которого входили: 2 мкл очищенного ПЦР-продукта, 0,1 мкМ одного из праймеров, 0,5 мкл набора для секвенирования BigDye®, 3,75 мкл 5-кратного буфера для секвенирования (Applied Biosystems) и деионизованная вода. Режим реакции секвенирования: 96°С в течение 1 мин., после чего следовали 25 циклов с денатурацией 96°С - 15 с, отжигом праймеров 50°С - 15 с и стадией элонгации 60°С в течение 1 мин. 40 с. Очистку продукта после реакции секвенирования проводили, используя 3 мкл ацетата натрия (3 М, рН 4.6), 14,5 мкл воды и 62,5 мкл этилового спирта. Смесь тщательно перемешивали и 15 минут инкубировали при комнатной температуре. После центрифугирования аккуратно удаляли надосадочную жидкость, а осадок дважды промывали 70% этанолом. После высушивания осадок растворяли в 20 мкл Hi-DÍ формамида (Applied Biosystems). Раствор денатурировали 2 минуты при 96°С и сразу охлаждали во льду. Электрофорез проводили на автоматическом генетическом анализаторе ABI PRISM 3100-Avant (Applied Biosystems) в соответствии с протоколом производителя.

2.5. Анализ микросателлитов

Для генотипирования песцов были протестированы и отобраны две группы микросателлитов: локусы, расположенные во II классе главного комплекса гистосовместимости (ГКГ) - CFA12-4, CFA12-9, CFA12-17, CFA12-21 и CFA12-22 и локусы в интронах других генов (СХХ140, СХХ173, СХХ250, 377 и 771) (Ostrander et ai, 1993; Ostrander et al., 1995; Mellersh et al., 1997; Debenham et al., нсопубл. данные). Использованные локусы изначально были выделены и описаны для собаки (Canis familiaris).

Для исследования микросателлитов, связанных с ГКГ, использовались прямые праймеры, меченые флуоресцентными красителями. Для анализа микросателлитов в составе других генов применяли систему из трех праймеров, содержащую удлиненный прямой праймер, на 5' конце которого добавлялась последовательность праймера М13, обратный праймер и флуоресцентно-меченый праймер М13.

Полимеразную цепную реакцию каждого локуса проводили отдельно, в объеме 10 мкл, содержащем 1 мкл ДНК (3 мкл ДНК в случае исторических образцов), 0,3 мМ смеси dNTP's (Fermentas или Amersham), однократный ПЦР-буфер без детергентов (Qiagen), 1,25 ед. Hotstar Taq полимеразы (Qiagen) и деионизованную воду. Для микросателлитных локусов ГКГ в смесь для ПЦР добавляли по 0,5 мкМ прямого и обратного праймеров и 0,2 мМ MgCl2 (Qiagen). В случае маркеров в составе других

генов смесь содержала 0,02 мкМ удлиненного прямого праймера и по 0,5 мкМ обратного и М13 праймеров, а также 1,5 мМ MgCI2 (Qiagen). Реакцию проводили используя метод "touch down 58" в следующих условиях: предварительная денатурация 95°С - 15 мин. Далее следовало 35 циклов (45 циклов для исторических образцов ДНК) со стандартной денатурацией при 95°С - 30 с; отжигом праймеров на протяжении 30 с при начальной температуре 65°С, которая в течение первых 15 циклов каждый последующий цикл снижалась на 0,5°С до температуры 58°С; элонгация 1 мин. при 72°С. После окончания программы проводили заключительный цикл элонгации при 72°С в течение 5 мин. Для предварительной оценки концентрации полученного ПЦР-продукта проводили электрофорез в 2% агарозном геле.

Анализ микросателлитных локусов делали на генетическом анализаторе ABI PRISM 3100-Avant (Applied Biosystems) в соответствии с протоколом производителя, используя стандарты длины ROX-400HD или ROX-500 (Applied Biosystems).

2.6. Главный комплекс гистосовместимости

2.6.1. Амплификация фрагмента fi-цепи II класса ГКГ

Для амплификации антигенсвязывающего региона |3-цепи (DRB и DQB гены) II класса ГКГ применяли следующие праймеры: DM-l/DM-2 (Hedrick et al., 2000) и DQB1/DQB2 (Fernandez-Vina etal, 1997).

Амплификацию фрагмента проводили В объеме 20 мкл, содержащем 2-5 мкл ДНК, 0,2 мМ смеси dNTP's (QBIOgene), однократный ПЦР буфер без детергентов (Qiagen), 4 ед. Hotstar Tag полимеразы (Qiagen), прямой и обратный праймеры по 0,5 мкМ, 2 мМ MgCb (Q BlOgene) и деионизованную воду. ПЦР проводили в следующих условиях: предварительная денатурация и активация фермента 95°С - 15 мин., после чего следовало 35-45 циклов со стандартной денатурацией при 94°С - 60 с; отжигом праймеров на протяжении 60 с при соответствующей праймерам температуре и элонгацией 60 с при 72°С. После окончания программы проводили заключительный цикл достройки при 72°С в течении 10 мин. Для оценки концентрации, специфичности полученного ПЦР-продукта и отсутствия контаминации проводили электрофорез в агарозном геле. Далее определяли генетический полиморфизм исследуемых фрагментов, применяя два независимых подхода: (а) прямое секвенирование полученных ПЦР-продуктов и (б) анализ полиморфизма методом дифференциального дисплея (SSCP) с последующим секвенированием отдельных вырезанных из геля аллелей. Результаты прямого секвенирования сопоставляли с профилями SSCP-гелей и последовательностями отдельных аллелей.

2.6.2. Определение полиморфизма ГКГ методом дифференциального дисплея (SSCP)

Для получения профилей дифференциального дисплея ПЦР-продукт анализируемого фрагмента (2-5 мкл) смешивали с 8-10 мкл денатурирующего буфера для нанесения (подготовленный по протоколу ETC-manual, ETC Elektrophoresetechnik). Смесь денатурировали при 95°С, а после охлаждали на льду. После этого б мкл смеси наносили на 15% акриламидный гель (CleanGel DNA-HP, ETC. Elektrophoresetechnik), который готовили, следуя инструкциям производителя. Электрофорез проводили при температуре 12°С на горизонтальной охлаждаемой системе (Amersham Pharmacia Biotech) при условиях: 200 V, 20 mA, 10 W в течение 20 мин., после чего 450 V, 30 mA, 20 W в течение 3 час. 30 мин. Затем гель фиксировали и окрашивали серебрением (DNA Plus One Silver Staining Kit, Amersham Pharmacia Biotech). Каждый образец анализировали как минимум два раза. На каждый SSCP-гель для сравнения рядом с новыми образцами всегда наносили все известные аллели. Как минимум три образца каждой аллели были вырезаны из геля, повторно амплифицированы и секвенированы.

2.7. Историческая ДНК: аутентичность и контроль контаминации

Все процедуры подготовки образцов, распила, выделения ДНК и амплификации проводились в строгом соответствии с правилами работы с историческим материалом (Cooper, Poinar, 2001). Меры предосторожности по предотвращению контаминации были предприняты на всех этапах работы с историческим материалом. Распил костей проводили в специализированном боксе, оборудованном УФ лампой. Перед распилом бокс и оборудование вымывали хлорсодержащим дезинфицирующим раствором и облучали ультрафиолетом не менее 60 минут.

Выделение ДНК из музейных образцов и последующий анализ проводили в трех разных лабораториях (Россия, Великобритания, Германия). Для работы использовали специальный набор пипеток. Комнаты для выделения исторической ДНК и постановки ПЦР были пространственно разделены с комнатами, где проводились выделение современной ДНК и работа с ПЦР-продуктами. Музейные пробы, современные образцы и реагенты для их анализа хранили раздельно. Контроль контаминации проводили во всех случаях, где это было возможно. Для многих образцов делали повторное выделение ДНК. И для таких образцов результаты анализа совпадали. В очень редких случаях удавалось получить ПЦР-фрагмент длиной более 300 н.п.

Чтобы исключить вероятность учета ложных гомозигот при анализе ядерных маркеров, проводили многократные повторные исследования (фрагментный анализ микросателлитов, секвенирование и дифференциальный дисплей генов ГКГ) со стадии ПЦР или экстракции ДНК. Для исключения возможных артефактов секвенирования для тех образцов, которые обладали новыми последовательностями, проводили два дополнительных повторных секвенирования.

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

Редактирование и выравнивание нуклеотидных последовательностей проводили с помощью программы MEGA 4.0 (Tamura et al., 2007). Определение длин фрагментов микросателлитов и генотипирование проводили, используя программы Genotyper v.3.6. и Genemapper v.3. (Applied Biosystems).

Число вариабельных позиций, среднее число попарных нуклеотидных и аминокислотных замен, а также среднее по всем сайтам нуклеотидное разнообразие аллелей ГКГ определяли в программе MEGA 4.0. Действие балансирующего отбора на последовательности второго экзона исследованных генов ГКГ оценивали по количеству синонимичных (d¡) и несинонимичных (d¡i) замен в кодонах, участвующих (ABS) и не участвующих в антигенсвязывании (non-ABS). Для оценки достоверности различий применяли Z-тест (в программе MEGA 4.0), исходя из предположения нейтральности исследуемых кодонов (ds = ds).

Построение Neighbor Joining филогенетического дерева выполняли также в программе MEGA 4.0 (Saitou, 1987), используя модель Tamura - Nei. Для статистической поддержки топологии применяли бутстреп анализ (1000 реплик). Построение медианных сетей проводили с помощью программы NETWORK 4.5.1.0 (Bandelt et al., 1999).

Программа ARLEQUIN v3.1 (Schneider et al., 2000) использовалась для оценки нуклеотидного и гаплотипного разнообразия D-петли мтДНК согласно Nci (1987), проведения FSt теста по методу Tamura-Nei (Tamura, Nei, 1993), а также для определения среднего числа аллелей в микросателлитах - Аа и средней гетерозиготности - НА (Tajima, 1983). Исходя из уравнения Харди-Вайнберга, наблюдаемую (Н„) и ожидаемую (Не) гетерозиготности микросателлитов оценивали вместе с отклонениями от самого уравнения, используя алгоритм цепи Маркова (100 000 шагов) и dememorisation в 100 000 шагов (Levene, 1949). Оценку сцепленности микросателлитных локусов (Slatkin, 1994) проводили, используя ЕМ алгоритм с 1000 пермутаций.

ГЛАВА 3. РЕЗУЛЬТАТЫ 3.1. D-петля митохондриалъной ДНК

Консенсусные последовательности I гипервариабельного региона D-петли мтДНК длиной 297 н.п. были получены для 31 образца исторической популяции о-ва Медный, 1 образца исторической популяции о-ва Беринга, 3 образцов песцов о-ва Янкича (о-ва Ушишир) и 10 образцов песцов современной популяции континентальной Аляски.

В исторической популяции о-ва Медный у 26 проанализированных животных был обнаружен гаплотип, последовательность которого идентична гаплотипу MEDN (610 н.п., Genbank: DQ630748), ранее описанному для современных песцов о-ва Медный (Джикия и др., 2007). У остальных 5 образцов была зафиксирована гетероплазмия в 37-ой нуклеотидной позиции от начала последовательности D-петли. В указанной позиции наравне с характерным для гаплотипа MEDN аденином наблюдали также присутствие цитозина (рис. 1).

* * i i 1 п 2 2 2 4 4

2 2 2 3 3 3 3 3 3 3 3 3 9 2 4 5 9 2 3 8 з 4

п 2 4 6 1 2 3 4 5 6 7 8 9 8 8 4 5 5 2 0 8 9 8

MEDN Í 610) 30 с С С - - с с С с А А А Т Т С А Т Т Т Т Т G

MEDN (297) 29 - - п п

МЕШ- -Н 5 - - М п п

ВА1 1 - - С с п п

В1 2 - - с

32 9 - - А с А G

ВЗ В4 2 4 - - д А с с G

Б5 4 - - с G

86 7 - - А г т

Е7 2 - - А А С А G

S2 4 т А А с т т А А А с п п

Н2 2 т А А с т т А А А с С п п

С1 1 т G А с т т А А А с п п

N1 2 т А А - т т А А А с С п п

Н5 1 т А А - т т А А А с с п п

Рисунок i. Вариабельные позиции фрагмента D-петли мтДНК песцов Командорских островов и Аляски. Гаплотипы MEDN (610) (современная популяция о-ва Медный) и В1-В7 (современная популяция о-ва Беринга) опубликованы Джикия и др. (2007); п - количество животных с данным гаплотипом. Мономорфные позиции обозначены звездочкой.

У изученных песцов о-ва Янкича был выявлен такой же, как и у песцов с о-ва Медный, гаплотип MEDN.

В историческом образце с о-ва Беринга был зарегистрирован новый, ранее не описанный, гаплотип В Al, отличающийся от гаплотипа В1, описанного для современных песцов о-ва Беринга (Джикия и др., 2007), наличием цитозина в 37-й позиции. Примечательно, что именно эта позиция является полиморфной в случае гетероплазмы, выявленной у исторических образцов популяции о-ва Медный.

У материковых песцов Аляски было выявлено 5 разных гаплотипов (S2, Н2, С1, N1 и Н5), В обнаруженных последовательностях было зарегистрировано 7 полиморфных позиций (6 замен и 1 вставка). Наиболее распространенным гаплотипом в изученной выборке популяции Аляски оказался гаплотип S2, который был выявлен у 4 животных. Интересно, что этот гаплотип ранее был обнаружен только у песцов Сибири (Dalén et al., 2005). Гаплотип С1 был найден у песцов Канады, гаплотип N1 - распространен по всей Неоарктике, а Н2 и Н5 - по всей Голарктике (Dalén et ai, 2005).

3.2. Микросателлиты

Тридцать шесть исторических и 27 современных образцов ДНК песцов с о-ва Медный. 12 современных образцов с о-ва Беринга и 24 континентальных образца с Аляски были успешно генотипированы по 10 микросателлитным локусам. Исключение составил связанный с ГКГ микросателлит CFA12-21, который не удалось проанализировать в исторической популяции песцов о-ва Медный.

Общее количество аллелей в 10 микросателлитах для песцов о-ва Беринга составило 43, для песцов с Аляски - 80, а у современных животных с о-ва Медный обнаружено всего 14 аллелей. Даже, несмотря на то, что локус СРА12-21 для исторических образцов о-ва Медный не удалось проанализировать, общее количество аллей по оставшимся 9 локусам составило 17. Таким образом, для 9 микросателлитов, изученных в обеих выборках песцов о-ва Медный, все аллели, встреченные в современной выборке, фиксировали у предковых образцов. В то же время, 5 аллелей, выявленных у предковых песцов, не были обнаружены у современных животных. В исторической популяции о-ва Медный 3 локуса уже были мономорфны, а эффект «бутылочного горлышка» привел к фиксации мономорфизма еще в 3 микросателлитах.

Командорские песцы характеризуются значительно более низким полиморфизмом микросателлитов, чем континентальные животные Аляски. Генетическое разнообразие в предковой популяции песцов о-ва Медный оказалось немного выше, чем в современной популяции, но все же обе выборки с о-ва Медный менее вариабельны по сравнению с современной популяцией песцов о-ва Беринга (табл. 1).

Таблица 1. Количество выявленных аллелей (Ал) в 10 микросателлитных локусах, а также средние показатели количества аллелей (Ал) и гетерозиготности (НА) для песцов о-ва Медный, о-ва Беринга и Аляски. * - без данных по локусу СРА12-21.

Популяция Количество аллелей ^ Среднее число аллелей на локус Аа±80 Средняя гетерозиготность На±80 Количество образцов N

Истор. о. Медный 17* 1.89 ±0.64 0.26 ±0.17 36

Совр. о. Медный 14 1.40 ±0.49 0.19 ± 0.12 27

12* 1.33 ±0.47* 0,16±0.11*

Совр. о. Беринга 43 4.30 ± 1.00 0.60 ±0.33 12

Аляска 80 8.00 ± 2.00 0.79 ±0.41 24

В исторической и современной популяциях песцов о-ва Медный ни для одного из локусов не было выявлено отклонений от равновесия Харди-Вайнберга. Однако для песцов о-ва Беринга достоверное отклонение было зафиксировано в 377 локусе (Яо=0.417, #£=0.493; Р=0.025), а для микросателлита СРА12-4 в составе ГКГ отклонение было достоверным как для популяции о-ва Беринга (Яо=0.333, Яе=0.547; Р=0.029), так и для популяции Аляски (Яо=0.458, Я^0.854; /М).000).

Результаты фрагментного анализа были также протестированы на наличие сцепленности между исследуемыми локусами. Случаи сцепления аллелей были зарегистрированы среди микросателлитов, связанных с ГКГ, что и ожидалось, так как эти микросателлиты находятся в близко расположенных генах II класса ГКГ в составе одной (12-й) хромосомы.

3.3. Главный комплекс гистосовместимости

3.3.1. Филогенетическая связь и полиморфизм DRB и DQB генов II класса ГКГ

песцов

Генетический полиморфизм генов (З-цепи II класса ГКГ был проанализирован для 184 образцов ДНК песцов, 47 из которых были взяты у континентальных животных. Среди 137 образцов с Командорских островов были проанализированы исторические образцы о-ва Медный (л=43), и о-ва Беринга (п=9), а также современные образцы о-ва Медный (п=41) и о-ва Беринга («=44). Всего было выявлено 25 уникальных аллелей фрагмента (З-цепи II класса ГКГ. Несмотря на то, что праймеры DM-1/DM-2 и DQB1/DQB2 сконструированы для DRB и DQB локусов, соответственно, многие аллели были выявлены обоими парами одновременно. Основываясь на результатах сопоставления аллелей ГКГ песцов с генетической базой данных (GenBank), а также на основании филогенетического анализа (см. ниже) выявленные у песцов аллели были названы как Alia-ТЖВ*01 - *13 и yl//a-DQB*01 -*12 согласно предложенной Klein et al. (1990) схеме номенклатуры аллелей ГКГ. Большинство аллелей было описано на основании анализа 207 н.п., но последовательности /i//o-DQB*01, *08 и * 12 описаны только по 172 н.п.

Филогенетическая связь аллелей ГКГ песцов. Для того чтобы оценить принадлежность выявленных у песцов аллелей к DRB или DQB локусу, было проведено построение медианных сетей, а также филогенетического дерева с использованием neighbor-joining алгоритма (рис. 2). Оба построен™ отражают взаимосвязь между аллелями сходным образом. При анализе топологии полученных построений четко прослеживалось разделение всех исследуемых аллелей на две ветви (кластера), одна из которых включала DRB последовательности, другая - DQB. Последовательность /l//a-DRB*05 расположена относительно близко к области разделения DRB и DQB ветвей медианных сетей. Однако в случае построения филогенетического дерева указанная аллель располагалось в составе DRB кластера.

BLASTn. Сравнение с генетической базой данных подтвердило, что все описанные аллели являются фрагментом Р-цепи II класса ГКГ и показало высокий процент сходства с другими DRB или DQB аллелями представителей семейства псовых. При анализе первых десяти наиболее сходных последовательностей из базы данных /4//a-DRB*01 - 13 обладали 100% сходством с DRB геном, анализ Alla-DQB*04, 05, 07 и 09-12 выявил 100% сходства с DQB геном. Аллели Alia-DQB*02, 03 и 06 также обладали высоким уровнем сходства с DQB локусом, но 10-40% последовательностей среди первой десятки были аллелями DRB гена. Среди первых десяти наиболее сходных с yi//a-DQB*01 и 08 аллелями 80% принадлежало DRB гену, тогда как только 20% - DQB локусу.

Alla-DRB*03 AI!a-DRB*08 Alla-DRB*10 Alla-DRB*09 Alla-DRB*13

Alla-DRB*01

AI!a-DRB*02 Alla-DRB*12 .

• Alla-DRB*04

■ Alla-DRB*11

■ Alla-DRB*05

■ Alla-DQB*09 - Alla-DQB*10

j— Alla-DQB*03 »1— Alla-DQB*06

- Alla-DQB*12

7i Alla-DQB*01 Alla-DQB*08 Alla-DQB*02 A!la-DQB*05 Alla-DQB*07

DQB*04 AIIa-DQB*ll

1— aib-UI, S^r— Alla-I m- Alla-DC

Рисунок 2. Филогенетическое дерево выявленных у песцов 25 аллелей р-цепи II класса ГКГ (¿//e-DRB*01-*13, zl//a-DQB*01-*12), построенное по алгоритму neighbor-joining. В узлах - результаты бутстреп анализа.

Нуклеотидная и аминокислотная вариабельность. Для последовательностей DRB локуса число вариабельных нуклеотидных позиций составило 32 (16 аминокислот). Для DQB аллелей выявлено 34 нуклеотидных 18 аминокислотных замен (табл. 2). Среднее число попарных нуклеотидных замен составило 13,6 ± 2,3 (от 2 до 25 н.п.) для DRB аллелей и 14,6 ± 2,3 (1-27 н.п.) для DQB аллелей.

Среди описанных аллелей не было выявлено инсерций/делеций или стоп-кодонов. Все описанные для песцов аллели обладали уникальной нуклеотидной и аминокислотной последовательностями. Исключение составили yi//«-DQB+01 и Alla-DQB*08, которые различались одной синонимичной заменой и обладали идентичной аминокислотной последовательностью.

Число различных аллелей ¡J-цепи II класса ГКГ, регистрированных одновременно в одном образце, варьировало от 2 до 5. В подавляющем большинстве случаев (91,3%) для каждого из DRB или DQB генов фиксировали 1 или 2 аллели. Однако в 8,7% случаев фиксировали до 3 аллелей каждого гена. Наблюдаемое распределение аллелей указывает на то, что исследуемые локусы у некоторых животных могут быть дуплицированы. Анализ гаплотипов исследованных образцов подтвердил, что дупликация DRB или DQB генов встречается не у всех образцов и связана с определенными сочетаниями аллелей.

Таблица 2. Основные характеристики генетического разнообразия 1ЖВ и 0<ЗВ аллелей II класса ГКГ песцов. Показатели приведены со стандартной ошибкой (вЕ).

DRB локус

DQB локус

Вариабельные позиции

32 16

Нуклеотиды Аминокислоты

Среднее число попарных замен Нуклеотиды Аминокислоты Нуклеотидное разнообразие, среднее по всем сайтам

I 0.07 ±0.01

13.6 ±2.3 (2-25) 8.2 ± 1.9 (2-13)

34 18

14.6 ±2.3 (1-27) 8.1 ± 1.8 (0-17)

0.08 ±0.02

3.3.2. Балансирующий отбор ß-цепи II класса ГКГ

Исследуемый участок ß-цепи II класса ГКГ участвует в распознавании антигенов и формировании иммунного ответа, однако не все сайты, кодирующие этот пептид, а только антигенсвязывающие сайты принимают в этом непосредственное участие (Brown et al., 1988; 1993). Поэтому в составе исследованного региона разные кодоны испытывают разное селективное давление. В результате анализа аминокислотных последовательностей было выявлено, что большинство вариабельных позиций аллелей ГКГ песцов расположено именно в антигенсвязывающих сайтах (ABS, рис. 3).

Соотношение синонимичных (ds) и несинонимичных (4v) замен было вычислено для входящих и не входящих в состав антигенсвязывающего участка позиций (ABS и non-ABS, соответственно, таблица 3). Для антигенсвязывающих сайтов (ABS), число несинонимичных мутаций оказалось значительно выше, чем синонимичных, а их соотношение составило 3,3 (/"=0.009) и 3,4 (Р=0.007), для DRB и DQB последовательностей, соответственно. Напротив, для сайтов, не участвующих в узнавании антигенов (non-ABS), соотношение между несинонимичными заменами и синонимичными оказалось близким к единице. Такой характер распределения нуклеотидных замен указывает на действие селективных сил в антигенсвязывающем регионе, направленных на поддержание высокого уровня аминокислотной вариабельности.

Alla Alla Alla' Alla' Alla Alla-Alla-Alla' Alla' Alla' Alla Alla Alla Alla Alla Alla Alla Alla Alla Alla Alla Alla Alla' Alla Alla

-DRB'Ol -DRB*02 -ORB* 0 3 -DR3*04 •DRB*05 ■DRB*06 -DRB*07 -DRB*C8 -DRB*09 -DRS*10 -DRB'll -DRBU2 -DRB*13 -DQB*Ol -DQ8*02 -DQB*03 -DQB*04 -DQB*05 -DQB*06 -DQB*07 -DQB*08 -DCB'09 -DQB*10 -DQ8*11 -DQB*12

NGTERVRiLlRHItoREEb^RFDSD ■¿Q..V......ГА

, .B.:.-. VQ..X-.

p:i B.jf.l

fco,.s.t VQ.'.S.t. i.p.f.l' a.ftK...

J.kid.,

ft.jSK^.J."

. A-.-' ■ -r.

Ä.ak,.; Ä.ftsi.'..

,.SK,. .

• Fft

GE^RAVTELGRtDAEYiSNGQKEILEDARAAVDTSCRilNYGVCE

.Z. .

. ZI .

■ S,1

6 f,

•ili

■ S .

■ S..

■ р..

. .... Л..

• R. .

.R. .

R.(i P. ;

-S.

.R.'. • F..

•.VI;.

.v!

r?:

R. R. R. R.'. R. . R. .

,.№. Л.....

. '..... RV

9.EK..4..RV.

.....c.. ..

. ГЗ .....

: RH .E.....

. ......

Л ..< 1.

y' EK. .£. .P?77?^? . ES..V..P7...

V'.to«. .р.. .я.ка. ; .рк. .v. ,R,V. , F.<. .V. .FV. .'. f. .iR..g...|?.K. |.

.EK. .V. .RV . . L .E. .Ri??'???

DEVDRV!..s..,y..1;. DEV3CT. fe . V. .'. . .

SK..». .ву. •s- V-DEVDi». .E. . .'???./■?;

Рисунок 3. Аминокислотный полиморфизм выявленных у песцов DRB и DQB аллелей ГКГ с указанием сайтов, участвующих в антигенсвязывании, по Brown et al., 1988; 1993 (выделено серым цветом). Точками отмечены позиции, идентичные с верхней последовательностью.

Таблица 3. Соотношение синонимичных (ds) и несинонимичных (dN) замен, вычисленных для сайтов, участвующих в связывании антигенов (ABS) и сайтов, не участвующих в антигенсвязывании (non-ABS). п - количество кодоиов; Р указывает на достоверность различий и определено при помощи Z-теста.

Сайты_и_dN±SE_dS±SE dN/dS_Р

DRB аллели

ABS 21 0.22±0.06 0.07±0.05 3.3 0.009

Non-ABS 48 0.03±0.02 0.03±0.02 0.9 0.880

ВСЕ 69 0.08±0.02 0.04±0.02 2.1 0.029

DQB аллели

ABS 21 0.19±0.06 0.06±0.04 3.4 0.007

Non-ABS 48 0.05±0.02 0.06±0.03 0.9 0.809

ВСЕ 69 0.09±0.02 0.06±0.02 1.6 0.216

3.3.3. Сопоставление генетического разнообразия ГКГ в командорских и материковых популяциях песцов

В настоящем исследовании ГКГ было выявлено отсутствие дифференциации материковых популяций (^7=0.051, Р=0.071). Из 25 аллелей, выявленных у песцов, 20 различных последовательностей регистрировали у континентальных особей: 9 БИВ аллелей Ы//л-ОКВ+02 и *06-*13) 11 ВрВ аллелей (АПа-0(}В*01, *02 и *04-*12; Таблица 4). И только 8 аллелей было выявлено у песцов Командорских островов (Л//а-01Ш*01-*05 и Л//а-БС>В*01-*03). Три из всех описанных аллелей (А11а-ВШ*02,

Л//я-0(ЗВ*01 и *02) оказались общими для популяций командорских и континентальных песцов.

Таблица 4. Относительная частота ВГШ и аллелей II класса ГКГ у песцов. В скобках указано число особей, у которых регистрировали данную аллель; п - размер выборки.

Континентальные песцы Сибирь Аляска ВСЕГО »=33 /1=14 «=47

Песиы Командорских островов Истор. Истор. Совр. Совр.

о.Беринга о.Медный о.Беринга о.Медный /1=09 л=43 п=44 //=41

Л11а-оав*01 - - - 0.44 (4) 0.81 (35) 0.50 (22) 1.00(43)

АИа-ОЮЗ'Ог 0.73 (24) 0.93 (13) 0.79 (37) 0.56 (5) - 0.52 (23) -

А11а-ОИ.В*03 - - - 0.11(1) - 0.07(3) -

А11а-С11В*04 - - - 0.33 (3) - 0.23 (10) -

А11а-0ав*05 - - - 0.44 (4) 0.42(18) 0.45 (20) -

А11а-0!1В*06 - 0.21 (3) 0.06 (3) - - - -

АИа-01Ш*07 0.39(13) 0.36(5) 0.38(18) - - - -

А!!а-ГЖВ*08 0.03(1) - 0.02 (1) - - - -

А11а-0ав*09 0.27(9) 0.07(1) 0.21(10) - - - -

АИа-ОЯВЧО 0.03(1) - 0.02(1) - - - -

АИа-ОЯВ*! 1 0.03 (1) 0.07(1) 0.04 (2) - - - -

А1ЫЭЯВ*12 0.12(4) - 0.09 (4) - - - -

А11а-ОЯВ*!3 0.09 (3) - 0.06 (3) - - - -

А11а-0(}В*01 0.61(20) 0.29 (4) 0.51 (24) 0.11(1) 0.95(41) 0.34(15) 1.00 (43)

А11а-одВ*02 0.64(21) 0.93(13) 0.72 (34) 0.78 (7) 0.33(14) 0.84 (37) -

АНа-ООВ'ОЗ - - - 0.67 (6) 0.33(14) 0.52 (23) -

А11а-0дв*04 0.18(6) - 0.13(6) - - - -

А1ЫХ5В*05 - 0.07(1) 0.02(1) - - - -

А11а-В(}В»06 0.06 (2) 0.14(2) 0.09 (4) - - - -

А11а-ОС>В»07 0.03 (1) 0.14(2) 0.06(3) - - - -

АНз-ООВ'ОЕ - 007(1) 0.02(1) - - - -

А!!а-0<ЗВ»09 0.12(4) - 0.09 (4) - - - -

АНа-РОВ* 10 0.06 (2) - 0.04 (2) - - - -

А11а-ОС)В*11 0.06 (2) - 0.04 (2) - - - -

А11а-ОС>В»12 0.12(4) - 0.09 (4) - - - -

У всех проанализированных образцов песцов о-ва Медный в каждом из исследованных генов было выявлено по одной аллели (Л//а-011В*01 и Л//а-ОС?В*01), которые находились в гомозиготном состоянии. Однако среди проанализированных исторических образцов с о-ва Медный, добытых до начала эпизоотии, было выявлено 5 аллелей ГКГ. Кроме упомянутых атлеяей для современных песцов, в исторической популяции были обнаружены Л&ЯЖВ*05, Л//д-0<ЗВ*02 и *03. Наиболее распространенными последовательностями были Л/и-ОКВ*01 и /Ша-Б<ЗВ*01, которые встречались в 81% и 95% образцов исторической популяции о-ва Медный, соответственно.

В исторической популяции песцов о-ва Беринга, как и у современных песцов, были зарегистрированы 8 аллелей (^//д-ВКВ*01-*05 и Л//я-ВС?В*01-*03). Изменения частот аллелей за последние 80-100 лет в популяции песцов о-ва Беринга не обнаружено (Рзт = -0.009, /М).544). Последовательности Л//а-В11В*02-*04 были уникальными для популяции о-ва Беринга, тогда как остальные пять аллелей встречались в обеих командорских популяциях. Таким образом, популяция предковых медновских песцов характеризовалась подвыборкой аллелей ГКГ, выявленных в популяции о-ва Беринга.

ГЛАВА 4. ОБСУЖДЕНИЕ

4.1. Гетероплазмия D-петли мтДНК песцов о-ва Медный

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

Ранее считалось, что гетероплазмия мтДНК - явление крайне редкое, однако в последнее время публикуется все большее количество результатов исследований, свидетельствующих об обратном (Calloway et al., 2000; Grzybowski, 2000; Kmiec et ai, 2006). Гетероплазмия, рассматриваемая как следствие мутационного процесса, представляет собой промежуточное состояние от возникновения полиморфизма до того момента, как новый вариант мтДНК будет фиксирован в популяции как отдельный гаплотип или же утерян (Kmiec et al., 2006; Millar et al., 2008). Темп сегрегации разных вариантов гетероплазмы может варьировать в широких пределах в зависимости от вида организма и типа гетероплазмы (Bendall, Sykes, 1995; Gocke et al, 1998).

4.2. Разнообразие микросателлитов: локусы, связанные с ГКГ и другими

генами

В настоящей работе в анализе нейтрального полиморфизма использовали микросателлиты, расположенные в некодирующей части генов, находящихся под разными селективными силами. Гены ГКГ находятся под влиянием балансирующего отбора, остальные микросателлиты расположены в генах, не испытывающих такого селективного воздействия. При сравнении исследованных локусов не было выявлено более высокой вариабельности в локусах, связанных с генами ГКГ, по сравнению с микросателлитами других генов (табл. 5).

Наиболее вероятно, что дефицит наблюдаемой гетерозиготности, выявленный в локусе 377 и в локусе, связанном с ГКГ (CFA12-4), является следствием наличия нулевых аллелей. Праймеры, используемые в настоящей работе, были сконструированы на основании ДНК собаки (С. familiaris). Возможная мутация в области отжига праймеров, которая могла возникнуть у песцов, не позволяет фиксировать мутантный аллель. Отклонение от равновесия Харди-Вайнберга были зарегистрировано только в популяции о-ва Беринга и Аляски. Присутствие не выявленных аллелей микросателлитов в этих популяциях только подтверждает более высокое их аллельное разнообразие по сравнению с песцами о-ва Медный.

Таблица 5. Средняя гетерозиготность (IIА) микросателлитов в составе ГКГ и других генов (всего 10 локусов) в исследованных популяциях песцов (с указанием стандартного отклонения, Б Б). * - без данных по локусу СРА12-21; п -число проанализированных проб.

Микросателлиты Истор. о. Медный Совр. о. Медный Совр. о. Беринга Аляска

Другие гены 0.32±0.21 0.29±0.20 0.65±0.39 0.84±0.47

Связанные с ГКГ 0.18±0.15* 0.09±0.09 0.00* 0.54±0.33 0.73±0.42

ВСЕ - 0.19±0.12 0.60±0.33 0.79±0.41

0.26±0.16* 0.16±0.11*

4.3. Сходство DRB и DQB генов II класса ГКГ песцов

Исследованные гены р-цепи II класса ГКГ выполняют общую функцию распознавания, связывания антигенов и последующей репрезентации Т-лимфоцитам и характеризуются сходной структурой (Klein, 1986). О высокой степени сходства между аллелями указанных генов ранее сообщалось в других исследованиях представителей семейства псовых (Seddon, Ellegren, 2002). Объяснением сходства DRB и DQB аллелей может быть действие отбора, конвергенция, но более вероятно причиной является межгенная рекомбинация (Andersson, Mikko, 1995; Seddon, Ellegren, 2002).

В случае филогенетического построения разделение аллелей на два кластера прослеживается весьма четко. Каждый из кластеров содержал /í//a-DRB*01 или Alla-DQB*01 аллель. Все исследованные образцы о-ва Медный характеризовались одинаковым генотипом, представленным одновременным присутствием указанных последовательностей, что неоспоримо свидетельствует о гомозиготном состоянии каждого из локусов. Принадлежность аллелей к DRB или DQB гену в большинстве случаев также подтверждалась результатами сравнения с генетической базой данных (BLASTn).

Для 8,7% образцов была выявлена дупликация DRB или DQB гена. События дупликации генов ГКГ неоднократно фиксировали для некоторых видов и отдельных популяций (Slierendregt et al, 1997; Baker et al., 2006; Van Oosterhout et al., 2006; Rosal-Sánchez et al., 2008). Считается, что дупликация играет важную роль в эволюции ГКГ, одним из следствий чего является повышение полиморфизма, так как увеличивается число возможных аллелей дуплицированного локуса в отдельном организме.

4.4. Влияние островной изоляции на генетический полиморфизм

4.4.1. Нейтральные маркеры

Проведенный анализ выявил значительную утрату нейтрального генетического разнообразия в островных формах песцов по сравнению с континентальными. Кроме того, уровень полиморфизма коррелирует с размером острова. Исследование исторических образцов подтвердило крайне низкий нейтральный полиморфизм песцов о-ва Медный еще до эпизоотии. За исключением выявленной гетероплазмы О-петли (МБЮМ-Н), не было зафиксировано гаплотипного разнообразия мтДНК. Также был зарегистрирован низкий уровень разнообразия в микросателлитах. Предковые песцы о-ва Медный характеризовались значительно более низким разнообразием по сравнению с более крупной популяцией о-ва Беринга. Таким образом, результаты проведенного исследования подтверждают негативное влияние продолжительной островной изоляции на генетический полиморфизм командорских песцов.

Уточнение данных по нейтральному генетическому разнообразию современных песцов о-ва Медный. Полученные в настоящем исследовании данные по вариабельности мтДНК в исторической популяции о-ва Медный свидетельствуют об отсутствии гаплотипного полиморфизма. Единственный гаплотип, зарегистрированный в исторической популяции, является идентичным последовательности, выявленной в современной популяции гаплотипа МЕОМ Отсутствие вариабельности в исторических образцах подтверждает мономорфизм современных песцов, обнаруженный Джикией и др. (2007). Кроме того, в трех исследованных образцах песцов с о-ва Янкича из популяции, интродуцированной с о-ва Медный еще в начале прошлого века, был выявлен тот же единственный гаплотип МЕБ>), что и в исторических образцах о-ва Медный. Опубликованное Ое1Геп е< а!., (2007) митохондриальное разнообразие современных медновских песцов (в 7 пробах было выявлено 3 гаплотипа) не согласуется с мономорфизмом мтДНК в исторической популяции.

Происхождение песцов о-ва Янкича.

Результаты нуклеотидного анализа фрагмента О-петли мтДНК подтверждают происхождение песцов о-ва Янкича из популяции о-ва Медный. В 3 образцах из популяции о-ва Янкича выявлен тот же уникальный гаплотип МЕОКГ, что и в популяции о-ва Медный. Исследованная выборка слишком мала для каких-либо популяционных выводов, однако она составляет примерно 10-15% от всей популяции (Гольцман, персональное сообщение). Стоит отметить, что мтДНК наследуется только по материнской линии, поэтому нельзя с уверенностью утверждать об отсутствии самцов-предков из других популяций.

4.4.2. Селективно значимый полиморфизм ГКГ

Полученные в данной работе результаты свидетельствуют о том, что командорские песцы характеризуются значительно меньшим генетическим разнообразием генов II класса ГКГ, чем континентальные животные. В случае сравнения вариабельности ГКГ между командорскими популяциями те песцы, которые населяли меньший по площади остров (о-в Медный), характеризовались

более низким разнообразием, чем большая популяция о-ва Беринга. Понижение уровня селективно значимого генетического разнообразия в островных популяциях песцов по сравнению с материковыми сопоставимо с таковым в нейтральных локусах.

Выявленная общность аллелей ГКГ командорских песцов подтверждает их общее происхождение и указывает на влияние эффекта основателя в момент отделения популяции о-ва Медный от песцов о-ва Беринга. Сохранение сходного аллельного разнообразия генов ГКГ также может быть результатом эволюционирования популяций в сходных экологических условиях.

Фиксирование трех общих аллелей ГКГ у командорских и континентальных песцов указывает на длительную их устойчивость в природе, по крайней мере, до момента, когда командорские песцы отделились от предковой популяции. Одним из объяснений продолжительной устойчивости указанных последовательностей может быть следствие балансирующего отбора (1Мщеш1а-Мога1ез е( а/., 2007). Можно также предположить, что указанные аллели обладают более высокой адаптивной ценностью по сравнению с другими аллелями.

4.5. Влияние демографического «бутылочного горлышка» на генофонд полированной популяции песцов

4.5.1. Нейтральное генетическое разнообразие

Столь низкий полиморфизм мтДНК, выявленный у песцов о-ва Медный еще до эпизоотии, не позволил использовать этот маркер в описании влияния прохождения популяции через «бутылочное горлышко» на её генетическую изменчивость.

Использование более вариабельных микросателлитов позволило провести сравнение нейтрального разнообразия до и после эпизоотии. Следствием эпизоотии была утрата 5 из 17 аллелей в 9 проанализированных локусах и мономорфизм в трех локусах. Средняя гетерозиготность (НЛ) уменьшилась в популяции с 0,26±0,17 до 0,16±0,11 (рассчитано для 9 локусов). Таким образом, были зарегистрированы четкие признаки утраты генетического разнообразия популяцией о-ва Медный в связи с прохождением ей демографического «бутылочного горлышка».

4.5.2. Полиморфизм генов ГКГ

После эпизоотии вариабельность была утеряна в обоих исследованных локусах ГКГ песцов о-ва Медный, в то время как до эпизоотии оба гена сохраняли полиморфизм. Таким образом, интенсивность потери генетического разнообразия генов II класса ГКГ сопоставима или даже превышает потерю полиморфизма в нейтральных локусах (микросателлиты). Важно отметить, что две аллели, которые были наиболее распространены в исторической популяции (Л//<я-01Ш*01 и АИа-0(}В*01), сохранились у современных особей, в то время как три более редкие последовательности были утеряны. Песцы о-ва Беринга не испытали столь значительного сокращения численности в прошлом веке. В исторических образцах популяции о-ва Беринга, как и в современных, фиксировали те же 8 аллелей ГКГ, а их частоты за последние 80-100 лет не изменились.

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

4.6. Генетический полиморфизм и приспособленность

Популяция о-ва Медный с сильно редуцированным нейтральным полиморфизмом, но сохранявшая вариабельность генов р-цепи II класса ГКГ, была способна успешно существовать продолжительное время в условиях изоляции и значительного антропогенного воздействия. Однако эпизоотия, приведшая к незначительному уменьшению и без того предельно низкого нейтрального полиморфизма, кардинально повлияла на гены ГКГ. В результате прохождения популяцией «бутылочного горлышка» оба исследованных локуса ГКГ утратили полиморфизм полностью. Современная популяция о-ва Медный остается на критически низком уровне, не восстанавливаясь до прежних размеров уже более 30 лет, несмотря на благоприятные условия среды (отсутствие конкурентов, антропогенного влияния, стабильность кормовых ресурсов, мягкий климат).

Возможно, что полиморфизм генов ГКГ в популяции песцов о-ва Медный до эпизоотии был недостаточным для формирования адекватного иммунного ответа на эктопаразита - ушного клеща, скорее всего, интродуцированного на Командорские острова с собаками, и это оказалась решающим фактором в высокой чувствительности к заболеванию.

Современные песцы о-ва Беринга обладают большим числом аллелей, чем предковые песцы о-ва Медный, но значительно меньшим, чем континентальные песцы. Стоит отметить, что частые случаи инфицирования ушной чесоткой, ведущие к смерти щенков, также регистрировали у современных песцов о-ва Беринга (наблюдения автора). Учитывая историю популяции о-ва Медный и современное состояние песцов о-ва Беринга, численность которых находится в пределах 600-800 особей, мы рекомендуем обратить особое внимание на охрану этой популяции. Популяция о-ва Беринга может оказаться уязвимой в случае значительного снижения численности, ведущего к потере имеющегося полиморфизма, или проникновения новых патогенов.

выводы

1. Впервые проведенное описание антигёнсвязывающего региона DRB и DQB генов II класса ГКГ у песцов выявило 25 уникатьных последовательностей, в которых обнаружены свидетельства действия балансирующего отбора, направленного на поддержание высокого аминокислотного разнообразия.

2. Результаты анализа ГКГ и нейтральных маркеров свидетельствуют о потоке генов между популяциями песцов Северной Америки и Азии и о значительно более высоком генетическом полиморфизме у континентальных, чем у островных песцов.

3. Продолжительная изоляция песцов ведет к снижению генетического полиморфизма. На о-ве Медный еще до последнего демографического «бутылочного горлышка» было полностью утрачено разнообразие гаплотипов D-петли мтДНК, и в значительной степени разнообразие микросателлигов и генов ß-цепи II класса ГКГ (Alla-DRB*01,*Q5 и AIla-DQB*01-*03).

4. Демографическое «бутылочное горлышко» на о-ве Медный привело к еще большему снижению полиморфизма микросателлитов (исчезновению 5 из 17 аллелей в 9 локусах) и к полной утрате разнообразия антигенсвязывающего региона II класса ГКГ.

5. Большая по размерам популяция о-ва Беринга сохранила более высокий уровень полиморфизма, чем популяция о-ва Медный,: большее число аллелей в микросателлитах и генах ß-цепи II класса ГКГ (Alla-DRB*01-*05 и Alla-DQB*01-*03). При этом в генах ГКГ с начала прошлого века изменений в частотах и аллельном разнообразии не произошло.

6. Анализ D-петли мтДНК подтверждает происхождение песцов о-ва Янкича из популяции о-ва Медный.

Список работ, опубликованных по теме диссертации

Джикия Е.Л., Плошница А.И., Колесников A.A., Гольцман М.Е. 2008. Генетический мономорфизм второго экзона гена DRB1 главного комплекса гистосовместимости у песцов острова Медный (Alopex lagopus semenovi Ognev, 1931). ДАН. T. 418. №2. С. 277-278.

Джикия E.Jl., Чудакова ДА., Плошница А.И., Загребельный C.B., Колесников A.A., Гольцман М.Е. 2006. Генетический полиморфизм Командорских популяций песца // Сохранение разнообразия животных и охотничье хозяйство России: сборник материалов международной конференции, Москва. С. 47-48.

Джикия Е.Л., Плошница А.И., Колесников A.A., Гольцман М.Е. 2007. Генетический мономорфизм второго экзона гена DRB1 главного комплекса гистосовместимости у песцов (Alopex lagopus semenovi Ognev, 1931) острова Медный (Командорские острова). // Сохранение биоразнообразия Камчатки и прилегающих морей: Материалы VIII международной научной конференции, посвященной 275-летию с начала Второй Камчатской экспедиции (1732-1733 гг.). С.312.

Джикия Е.Л., Плошница А.И., Колесников А.А., Голъцман М.Е. 2007. Генетический мономорфизм в популяции песцов острова Медный // Молекулярно-генетические основы сохранения биоразнообразия млекопитающих Голарктики: сборник материалов международной конференции (Черноголовка). М.: Товарищество научных изданий КМК. С. 57-60.

Плошница А.И., Джикия E.JI., Голъцман М.Е. 2007. Анализ контрольного региона исторической мтДНК в изолированной популяции песцов до прохождения через «бутылочное горлышко» // Молекулярно-генетические основы сохранения биоразнообразия млекопитающих Голарктики: сборник материалов международной конференции (Черноголовка). М.: Товарищество научных изданий КМК. С. 191-195.

Ploslmitsa A., Dzhikiya Е., Goltsman М. 2007. The MHC-DRB1 monomoiphism in the island subspecies of Arctic foxes (Alopex lagopus semenovi, Ognev 1931) // 87th Annual Meeting of American Society of Mammalogists: conference contributions, Albuquerque, New Mexico, USA. P. 68.

Ploshnitsa A. 2007. Genetic monomorphism (control region mtDNA and MHC-DRB1) in the Mednyi Arctic fox // 6th International Zoo and Wildlife Research Conference on Behavior, Physiology and Genetics: conference contributions, Berlin, Germany. P. 172.

Bocharova N., Ploshnitsa A., Doronina L., Solomatnikova A., Goltsman M. 2007. Preliminary study of the blood in the endemic insular Arctic fox population // 6lh International Zoo and Wildlife Research Conference on Behavior, Physiology and Genetics: conference contributions, Berlin, Germany. P.44.

Bocharova N.. Stefanski V., Ploshnitsa A., Doronina L„ Goltsman M. 2009. Preliminary immunologic study in endangered Arctic fox subspecies (Alopex lagopus semenovi, Ognev. 1931) // 7th International Zoo and Wildlife Research Conference on Behavior, Physiology and Genetics: conference contributions, Berlin, Germany. P. 67.

Ploshnitsa A., Goltsman M, Kennedy L., MacdonaldD., SommerS. 2009. Effects of historical founder effects and a recent bottleneck on MHC variability in Commander arctic foxes // 7th International Zoo and Wildlife Research Conference on Behavior, Physiology and Genetics: conference contributions, Berlin, Germany. P. 194.

Ploshnitsa A., Goltsman M, Kennedy L„ Macdonald D. W., Sommer S. 2010. The P-chain MHC class II variability in Arctic foxes // 1st Community Genetics meeting: conference contributions, Manchester, UK. P. 12.

Отпечатано на ризографе вОНТИГЕОХИРАН Тираж 150 экз.

Содержание диссертации, кандидата биологических наук, Плошница, Анна Ивановна

СПИСОК СОКРАЩЕНИЙ.

ВВЕДЕНИЕ.

ГЛАВА 1. ОБЗОР ЛИТЕРАТУРЫ.

1.1. Маркеры генетического разнообразия популяции

1.1.1. Нейтральные маркеры.

1.1.2. Иммунные гены главного комплекса гистосовместгшости.

1.1.2.1. Структура и функция

1.1.2.2. Селективные механизмы

1.1.2.3. Эмпирические данные, свидетельствующие о значимости разнообразия ГКГ в устойчивости к инфекциям.

1.1.2.3.1. Исследования значимости вариабельности ГКГ и селективных механизмов у человека и лабораторных животных.

1.1.2.3.2. Исследования значимости вариабельности ГКГ и селективных механизмов в естественных популяциях.

1.2. Эволюционные силы, действующие на генетическое разнообразие кодирующих и не кодирующих последовательностей.

1.2.1. Процессы, ведущие к снижению генетического разнообразия.

1.2.2. Свидетельства действия балансирующего отбора в популяциях с отсутствием или крайне низким нейтральным полиморфизмом.

1.3. Специфика биологии островных тихоокеанских популяций песца.

1.3.1. Распространение и среда'обитания.

1.3.2. Генетическая дифференциация популяций.

1.3.3. Особенности биологии островных тихоокеанских популяций.

1.3.4. История командорских популяций и влияние деятельности человека.2%

1.3.5. Эпизоотия чесотки в популяции о-ва Медный и состояние демографического «бутылочного горлышка».

ГЛАВА 2. МАТЕРИАЛЫ И МЕТОДЫ.

2.1. Сбор образцов.

2.2. Выделение ДНК.

2.2.1. Выделение современной ДНК из мягких тканей.

2.2.2. Выделение исторической ДНК из костной ткани.

2.3. Амплификация фрагмента D-петли мтДНК.

2.4. Определение последовательности ДНК.

2.5. Анализ микросателлитов.

2.6. Главный комплекс гистосовместимости.

2.6.1. Амплификация фрагмента р-цепи II класса ГКГ.

2.6.2. Определение полиморфизма ГКГ методом дифференциального дисплея (SSCP).

2.7. Историческая ДНК: аутентичность и контроль контаминации.

2.7.1. Меры предосторожности и контроль контаминации.

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

ГЛАВА 3. РЕЗУЛЬТАТЫ.

3.1. D-петля митохондриальной ДНК.

3.2. Микросателлиты.

3.3. Главный комплекс гистосовместимости.

3.3.1. Филогенетические связи и полиморфизм DRB и DQB генов II класса ГКГ песцов.

3.3.2. Балансирующий отбор /3-цепи II класса ГКГ.

3.3.3. Сопоставление генетического разнообразия DRB и DQB ГКГ в командорских и материковых популяциях песцов.

ГЛАВА 4. ОБСУЖДЕНИЕ.

4.1. Гетероплазмия D-петли мтДНК песцов о-ва Медный.

4.2. Разнообразие микросателлитов: микросателлиты, связанные с ГКГ и другими генами.

4.3. Сходство DRB и DQB генов II класса ГКГ песцов.

4.4. Влияние островной изоляции на генетический полиморфизм.

4.4.1. Нейтральные маркеры.

4.4.2. Селективно значимый полиморфизм ГКГ.

4.5. Влияние демографического «бутылочного горлышка» на генофонд изолированной популяции песцов.

4.5.1. Нейтральное генетическое разнообразие.

4.5.2. Полиморфизм генов ГКГ.

4.6. Генетический полиморфизм и приспособленность.

ВЫВОДЫ.Г.

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

Введение Диссертация по биологии, на тему "Влияние островной изоляции и эффекта "бутылочного горлышка" на генетический полиморфизм командорских песцов"

Дифференциация популяций и видообразование на океанических островах служит предметом изучения зоологов, начиная с Дарвина и Уоллеса (см., например, MacArthur, Wilson, 1967, Van der Geer et al., 2010). Изоляция, простые экосистемы, небольшие размеры популяций и высокая скорость эволюционных изменений делают островные виды привлекательными для исследователей. Однако островные популяции особенно уязвимы к любым средовым изменениям, и гораздо беззащитнее материковых перед антропогенными угрозами (Alcover et al., 1998, Frankham, 1998). Исследование изменений, происходящих при сокращении небольших изолированных популяций тем более актуально, что все больше фрагментированных материковых популяций превращается в «островные» из-за непреодолимых антропогенных преград. Радикальные изменения и исчезновения популяций происходит не на эволюционных временах, а на протяжении лишь десятков лет, а иногда значительно быстрее. Процессы, происходящие на этих временах, то есть временах поколений, изучены очень слабо (Johnson, 2003). Морфологические индикаторы слишком инерционны, а экологические характеристики, напротив слишком изменчивы и обычно неизвестны для более, чем десяти поколений.

Одна из командорских популяций, популяция о-ва Медный в 1970-х годах прошла через демографическое «бутылочное горлышко», из-за внезапной эпизоотии ушной чесотки (Goltsman et al., 1996). После этого популяция восстановилась до 10-15% былой численности и стабильно сохраняет этот уровень в течение 30 лет. Сейчас эта популяция занесена в Красную Книгу РФ, в категорию редкости 1, как находящаяся под угрозой исчезновения (Гольцман, Крученкова, 2001). История командорских песцов хорошо документирована, а в музее храниться коллекция черепов, добытых около 100 лет назад. Кроме того, в начале прошлого века не безуспешно были предприняты попытки интродукции командорских песцов на острова Ушишир Курильской гряды (Воронов, Воронов, 1963; Костенко и др., 2004). Это дает уникальную возможность испытать применимость молекулярно-генетического подхода в исследовании популяционных процессов в масштабе времени десятков поколений, и изучить влияние островной изоляции и прохождения популяции через «бутылочное горлышко».

Нейтральные маркеры, такие как D-петля мтДНК или микросателлиты, информативны в исследовании филогенетических связей, популяционной истории, молекулярных часов, потока генов и родства (Blouin et al., 1996; Nabata et al., 2004; Bruford, Wayne, 1993; Selkoe, Toonen, 2006; Balloux, Lugon-Moulin, 2002; Kasapidis et al., 2005; Lemay, Bouldin, 2009). Но они не могут дать полноценную информацию о селективных процессах, протекающих, в популяции и ее адаптивных возможностях. В последнее время особый интерес привлекают к себе гены ГКГ, полиморфизм которых важен в формировании адекватного иммунного ответа на разнообразные патогены (Klein, 1986; Parham, Ohta, 1996; Apanius et al., 1997; Ploegh, Watts 1998; Hedrick, 1999; Sommer, 2005). Предполагается, что балансирующий отбор, действующий в генах ГКГ и направленный на поддержание высокого уровня разнообразя, способен снижать эффект дрейфа генов и замедлять скорость фиксации аллелей (Richman, 2000).

Цель работы

Исследовать влияние островной изоляции и эффекта «бутылочного горлышка» на нейтральный и адаптивно значимый генетический полиморфизм командорских песцов.

Задачи

1. Описать нейтральный полиморфизм D-петли мтДНК и микросателлитов у песцов предковой популяции о-ва Медный и микросателлитов у современных островных и континентальных песцов.

3. Оценить интенсивность балансирующего отбора в антигенсвязывающем регионе генов II класса ГКГ.

5. На основании анализа мтДНК выяснить происхождение песцов о-ва Янкича (о-ва Ушишир, Курильской гряды).

Научная новизна. Впервые описан полиморфизм DRB и DQB генов II Класса ГКГ у песца (.Alopex lagopus). Получены данные о разнообразии нейтральных (D-петля мтДНК и микросателлиты) и селективных (гены II класса ГКГ) маркеров в исторической популяции песцов о-ва Медный до последнего демографического «бутылочного горлышка». Впервые исследовано влияние прохождения популяции через «бутылочное горлышко» на генетическое разнообразие генов ГКГ на основе прямого сравнения данных по разнообразию в современных образцах, с таковым в исторических образцах животных, добытых до спада численности. С помощью молекулярно-генетических методов проведено сравнение командорских популяций песцов между собой и с материковыми популяциями. Исследован уровень полиморфизма в генетических маркерах, находящихся под действием разных селективных сил. Найдены генетические подтверждения того, что современные песцы о-ва Янкича происходят с о-ва Медный.

Теоретическое и практическое значение результатов работы.

Сравнительное исследование островных и континентальных песцов позволило оценить влияние островной изоляции на формирование современной генетической изменчивости популяций песцов Командорских островов. Показано, что генетический полиморфизм в популяциях островных песцов значительно ниже, чем у континентальных, при этом более крупная по размерам популяция о-ва Беринга способна сохранять более высокий генетический полиморфизм как в нейтральных, так и в селективных маркерах, чем популяция о-ва Медный. Полученные данные важны для понимания того, как демографические события (эффект «бутылочного горлышка») отражаются на разнообразии нейтральных генетических маркеров и локусов, которые находиться под действием балансирующего отбора. Понимание микроэволюционных процессов в генах ГКГ, которые определяют способность организма формировать иммунный ответ, важный этап в исследовании эволюции иммунной системы. Полученные данные является важными сведениями об адаптивных способностях островных популяций, и имеют большое значение в определении дальнейших путей их охраны.

Заключение Диссертация по теме "Зоология", Плошница, Анна Ивановна

выводы

1. Впервые проведенное описание антигенсвязывающего региона DRB и DQB генов II класса ГКГ у песцов выявило 25 уникальных последовательностей, в которых обнаружены свидетельства действия балансирующего отбора, направленного на поддержание высокого аминокислотного разнообразия.

3. Продолжительная изоляция песцов ведет к снижению генетического полиморфизма. На о-ве Медный еще до последнего демографического «бутылочного горлышка» было полностью утрачено разнообразие гаплотипов D-петли мтДНК, и в значительной степени разнообразие микросателлитов и генов Р-цепи II класса ГКГ (Alla-DRB*01,*05 и АПа-DQB*01-*03).

5. Большая по размерам популяция о-ва Беринга сохранила более высокий уровень полиморфизма, чем популяция о-ва Медный,: большее число аллелей в микросателлитах и генах р-цепи II класса ГКГ (Alla-DRB*01-*05 и Alla-DQB*01-*03). При этом в генах ГКГ с начала прошлого века изменений в частотах и аллельном разнообразии не произошло.

6. Анализ D-петли мтДНК подтверждает происхождение песцов о-ва Янкича из популяции о-ва Медный.

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

Я признательна всем, кто помогал на разных этапах этой работы. Прежде всего, хочу выразить огромную благодарность моему научному руководителю: М.Е.Гольцману, а так же моим руководителям в лабораториях университета Манчестера в Великобритании L.J.Kennedy и заведующей исследовательской группой иммуногенетики Leibniz-IZW в Берлине S.Sommer за чуткое руководство, терпение, внимание и помощь в работе. Хочу особо поблагодарить М.В. Холодову, А.А.Банникову, А.А.Колесникова и Н.А. Формозова за ценные рекомендации, полезные замечания и помощь в редакции настоящей работы. Я признательна всем тем, кто помогал в организации и учувствовал в экспедициях на Командорские острова: С.В.Загребельному, Н.Н.Павлову, Л.О.Дорониной, А.Н.Шиенку. Особая благодарность всем тем, кто осуществлял помощь в проведении лабораторных анализов, а так же давал ценные рекомендации и советы, делился опытом: Е.Л.Джикии, П.А.Сорокину, В.А.Орехову, С.В.Найденко и А.С.Карнаухову. Большое спасибо коллективу зоологического музея МГУ: О.Г.Нановой, С.В.Крускопу и В.С.Лебедеву за доброжелательное отношение, интерес к исследованию и предоставление исторических образцов костей командорских песцов. Кроме того, хочу поблагодарить А.Ю.Косорукову и участников проекта "Arctic predators as indicators of tundra ecosystem" за предоставление современных континентальных образцов. Я особенно признательна D.W.Macdonald и K.Jewgenow за помощь в организации поездок в Великобританию и Германию, с целью проведения там лабораторных работ, а так же В.К.Боженко за поддержку и возможность работы в лаборатории РНЦРР г. Москвы. И конечно, отдельно я хочу поблагодарить моих родных маму Т.С.Плошницу, брата И.И.Плошницу, тетю Н.С.Захарову, а так же моего молодого человека A.de Palma Garrido за помощь и поддержку.

Библиография Диссертация по биологии, кандидата биологических наук, Плошница, Анна Ивановна, Москва

1. Архивы Командорского зверосовхоза. 1928 - 1960. с. Никольское.

2. Батурин А.Д. 1928. Островное хозяйство Дальневосточного Края // Производительные силы Дальнего Востока. Вып.4. Хабаровск-Владивосток: Книжное дело. С.337-383.

3. Бочарова Н.А. Гольцман М.Е. 2008. Микроскопическое исследование кожного покрова песца о-ва Медный (Alopex lagopus semenovi Ognev 1931) // Материалы IX международной научной конференции 25-26 ноября 2008 г. Петропавловск-Камчатский. С. 250-253.

4. Васильев В.А. 1931. Голубой песец на острове Кильдине.// Союзпушнина. Вып.5.

5. Воронов В.Г., Воронов Г.А. 1963. Голубой песец на островах Ушишир (Курильская гряда) // Тр.Сах. НИИ АН СССР, Сибирское отд. Т. 14. С. 79-81.

6. Гептнер В.Г., Наумов Н.П., Юргенсон П.Б., Слутский А.А., Чиркова А.Ф., Банников А.Г. Млекопитающие Советского Союза. М.: Высшая школа. 1967. Т.2. ч.2. С.200-265.

7. Гольцман М.Е., Крученкова Е.П. 1995. Родительское поведение песца {Alopex lagopus semenovi) на острове Медном. Факторы, определяющие пространственно-временной режим активности // Зоологический журнал РАН. С. 88-103.

8. Гольцман М.Е., Крученкова Е.П., Сергеев С.Н., Володин И.А. 2003. Песец острова Медного {Alopex lagopus semenovi). Особенности экологии островной популяции // Зоологический журнал. Т. 82. № 4. С. 514-524.

9. Гончаров И. 1929. Опыт акклиматизации голубого песца на Шантарских острвоах.// Пушное дело. Вып.5.

10. Гончаров И. 1931. Об акклиматизации голубого песца на шантарских островах. // Союзпушнина. Вып. 6.

11. Данилина Н.Р. 1987. Состояние популяции беринговского песца и проблемы охраны природы Командор // Рациональное природопользование на Командорских островах. М.: МГУ. С.84-87.

12. Загребельный С.В. 2000. Экология питания песцов с о-ва Беринга (Alopex Iagopus beringensis) и о-ва Медного (A. I. semenovi) (Carnivora, Canidae), Командорские о-ва // Зоологический журнал. Т. 79, 595-607.

13. Зоотехнический отчет Медновского зверосовхозаза 1934 г. С. 1-21.

14. Ильина Е.Д. Островное звероводство. М.: Международная книга. 1950. 302 с.

15. Костенко В.А., Нестеренко В.А., Трухин A.M. 2004. Млекопитающие Курильского архипелага. Владивосток: Дальнаука. 186 с.

16. Лавров П.Н. 1932. Песец. М.: Внешторгиздат. 56 с.

17. Мараков С.В. Природа и животный мир Командор. М.: Наука, 1972. 185 с.

18. Нанова О.Г. 2009. Структура морфологического разнообразия признаков черепа и зубов трех видов хищных млекопитающих (Mamalia: Carnivora). Автореферат канд. диссертации. М: МГУ им. М.В. Ломоносова.

19. Овсяников Н.Г. 1981. Медновский песец может исчезнуть // Природа. С. 59-61.

20. Огнев С.И. Звери Восточной Европы и Северной Азии. Л.: Главнаука. 1931.776с.

21. Суворов Е.К. Командорские острова и пушной промысел на них. СПб. 1912. 324 с.

22. Холодова М.В., Сорокин П.А., Лущекина А.А. 2005. Изменение генетического разнообразия европейской популяции сайгаков (Saiga tatarica) в период резкого снижения численности // Доклады Академии наук. т. 404. №2. С. 1-2.

23. Цалкин В.И. 1944. Географическая изменчивость в строении черепа песцов Евразии // Зоологический журнал. Т. 23. № 4. С. 156-169.

24. Черский А.И. 1920. Командорский песец. Материалы по изучениюрыболовства и пушного промысла на Дальнем Востоке. Токио // Изд. Управления Рыб. и Мор. Звер. промыслами. Вып.1. С.60-107.

25. Эрлих Э.Н., Мелекесцев И.В. 1974. Командорские острова. История развития рельефа Сибири и Дальнего Востока. Камчатка, Курильские и Комадорские острова. М.: Наука. С. 337-344.

26. Якушкин Г.Д. 1966. Песец Таймыра и рациональное использование его запасов. Автореферат канд. диссертации. Норильск: НИИ Сельского Хозяйства Крайнего Севера.

27. Aguilar A., Roemer G., Debenham S., Binns M., Garcelon D., Wayne R.K. 2004. High MHC diversity maintained by balancing selection in an otherwise genetically monomorphic mammal // Proceedings of the National Academy of Sciences. V. 101. P. 3490-3494.

28. Alberts S. C., Ober C. 1993. Genetic variability in the major histocompatibility complex: a review of nonpathogen-meditated selective mechanisms. // Yearbook of Physical Anthropology. V. 36. P. 71-89.

29. Alcover LA., Sans A., Palmer M. 1998. The extent of extinctions of mammals on islands //Journal of Biogeography. V. 25. P. 913-918.

30. Altizer S., Harvell D., Friedle E. 2003. Rapid evolutionary dynamics and disease threats to biodiversity. // Trends Ecol Evol. V. 18. P. 589-596.

31. Amills M., Jimenez N., Jordana J., Riccardi A. et al. 2004. Low diversity in the major histocompatibility complex class II DRB1 gene of the Spanish ibex

32. Capra pyrenaica). // Hered. V. 93. P. 266-272.

33. Anderson L., Mikko S. 1995. Generation of Mhc class II diversity by Intra-and lntergenetic recombination. // Immunol Rev. V. 14. N. 3. P. 5-12.

34. Apanius V., Penn D., Slev P.R., Ruff L.R. 1997. The nature of selection on the major histocompatibility complex. // Crit Rev Immunol. V. 17. P. 179224.

35. Baker C. S., Vant M.D., Dalebout M.L., Lento G.M., O'Brien S.J., Yuhki N. 2006. Diversity and duplication of DQB and DRB -like genes of the MHC in baleen whales (suborder: Mysticeti). II Immunogenetics. V. 58. N. 4. P. 283296.

36. Balloux F., Goudet J., Perrin N. 1998. Breeding system and genetic variance in the monogamous, semi-social shrew, Crocidura russula. II Evol. V. 52. P. 1230-1235.

37. Balloux F., Lugon-Moulin N. 2002. The estimation of population differentiation with microsatellite markers. // Molecular Ecology. V. 11. P. 155-165.

38. Bandelt H.J., Forster P., Rohl A. 1999. Median-joining networks for inferring intraspecific phylogenies. // Mol Biol Evol. V. 16. N. 1. P. 37-48.

39. Behnke J.M., Wahid F.N. Immunological relationships during primary infection with Heligmosomoides polygyrus (Nematospiroides dubius): H-2 linked genes determine worm survival. // Parasitol. V. 103. P. 157-164.

40. Bendall K.E., Sykes B.C. 1995. Length Heteroplasmy in the First Hypervariable Segment of the Human mtDNA Control Region Am. J. Hum. // Genet. V. 57. P. 248-256.

41. Berger-Wolf T.Y., Sheikh S.I., DasGupta В., Ashley M.V., Caballero I.C., Chaovalitwongse W., Putrevu S.L. 2007. Reconstructing sibling relationships in wild populations. //Bioinformatics. V. 23. N.13. P. 149-156.

42. Bergstrom Т., Gyllensten U. 1995. Evolution of Mhc class II polymorphism: the rise and fall of class II gene function in primates. // Immunol Rev. V. 143. P. 13-31.

43. Bernatchez L., Landry C. 2003. MHC studies in nonmodel vertebrates: what have we learned about natural selection in 15 years? // Journal of evolutionary biology. V. 16. N.3. P. 363-77.

44. Berube M., Urban J., Dizon A.E., Brownell R.L., Palsboll P.J. 2002. Genetic identification of a small and highly isolated population of fin whales (Balaenoptera physalus) in the Sea of Cortez, Mexico. // Conserv. Genet. V. 3.N. 2. P. 183-190.

45. Blouin M.S., Parsons M., Lacaille V., Lotz S. 1996. Use of microsatellite loci to classify individuals by relatedness // Molecular Ecology. V. 5. P. 393-401.

46. Bodmer W.F. 1972. Evolutionary significance of the HLA-system. // Nature. V. 237. P. 139-183.

47. Briles W.E., Stone H.A., Cole R.K. 1977. Mareks-disease effects of B-histocompatiblity alloalleles in resistant and susceptible chicken. // Science. V. 195. P. 193-195.

48. Brown J.H., Jardetzky T.S., Gorga J.C., Stern L.J., Urban R.G., Strominger J.L., Wiley D.C. 1993. 3-Dimensional structure of the human class-II histocompatibility antigen HLADR1. //Nature. V. 364. P. 33-39.

49. Brown J.H., Jardetzky Т., Saper M.A., Samraoui В., Bjorkman P.J., Wiley D.C. 1988. A hypothetical model of the foreign antigen-binding site of class-II histocompatibility molecules. // Nature. V. 332. P. 845-850.

50. Brown W.M., George M., Wilson A.C. 1979. Rapid evolution of animal mitochondrial DNA. // Proc Natl Acad Sci USA. V. 76. P. 1967-1971.

51. Brown W.M., Prager M., Wang A., Wilson A.C. 1982. Mitochondrial DNA sequences of primates: tempo and mode of evolution. // J Mol Evol. V. 18. P. 225-239.

52. Bruford M.W., Wayne R.K. 1993. Microsatellites and their application to population genetic studies. // Current Opinion in Genetics & Development. V. 3. N. 6. P. 939-943.

53. Calloway C.D., Reynolds R.L., Herrin G.L.J, Anderson W.W. 2000. The frequency of heteroplasmy in the HVII region of mtDNA differs across tissue types and increases with age. // Am J Hum Genet. V. 66. N.4. P. 1384-1397.

54. Carrington M., Nelson G.W., Martin M.P., Kissner Т., Vlahov D., Goedert

55. J.J., Kaslow R., Buchbinder S., Hoots K., O'Brien S.J. 1999. HLA and HIV-1: heterozygote advantage and B*35-Cw*04 disadvantage. // Science. V. 283. P. 1748-1752.

56. Chen C.Y., Cohen S.A., Zaleski M.B., Albini B. 1992. Genetic control of streptococcus-induced hepatic granulomatous lesions in mice. // Immunogenetics. V. 36. P. 28-32.

57. Clarke В., Kirby D.R.S. 1966. Maintenance of histocompatibility polymorphisms. //Nature. V. 211. P. 999-1000.

58. Cooper A., Poinar H.N. 2001. Ancient DNA: do it right or not at all. // Science. V. 289. P. 1139.

59. Cooper, A. and Wayne, R. 1998. New uses for old DNA. // Current Opinion in Biotechnology. V. 9. P. 49-53.

60. Dalen L., Gotherstrom A., Tannerfeldt M. and Angerbjorn A. 2002. Is the endangered Fennoscandian Arctic fox {Alopex lagopus) population genetically isolated?//Biological Conservation. V. 105. P. 171-178.

61. Decamposlima P.O., Gavioli R., Zhang Q.J., Wallace L.E., Dolcetti R., Rowe M., Rickinson A.B., Masucci M.G. 2003. HLA-A11 epitope loss isolates of Epstein-Barr-virus from a highly A11+ population. // Science. V. 260 P. 98100.

62. Dengjel J., Schoor O., Fischer R., Reich M., Kraus M. et al. 2005. Autophagy promotes MHC class II presentation of peptides from intracellular source proteins. // Proc Nat Acad Sci USA. V. 102. P. 7922-7927.

63. Doherty P.C., Zinkernagel R.M. 1975. Enhanced immunological surveillance in mice heterozygous at the H-2 locus // Nature. V. 256. P. 50-52.

64. Eberhardt L.E., Garrott R.A., Hanson W.C. 1983. Winter movements of Arctic foxes, Alopex Iagopus, in a petroleum development area // Canadian Field-Naturalist. V. 97. P. 66-70.

65. Eizirik E., Kim J., Menotti-Raymond M., Crawshaw P.G. Jr, O'Brien S.J., Johnson W.E. 2001. Phylogeography, population history and conservation genetics of jaguars (Panthera onca, Mammalia, Felidae). // Molecular Ecology. V. 10. P. 65-79.

66. Ejsmond M.J., Radwan J. 2009. MHC diversity in bottlenecked populations: a simulation model. // Conservation genetics. V. 20. P. 1-09.

67. Ellegren H., Miklco S., Wallin K., Andersson L. 1996. Limited polymorphism at major histocompatibility complex (MHC) loci in the Swedish moose A. alces.ll MolEcol. V. 5. P. 3-9.

68. Else K.J., Wakelin D., Wassom D.L., Hauda K.M. 1990. The influence of genes mapping within the major histocompatibility complex on resistance to Trichuris muris infections in mice. // Parasitol. V. 101. P. 61-67.

69. Enriquez F.J., Brooks B.O., Cypess R.H., David C.S., Wassom D.L. 1988. Nematospiroides dubius: two Н-2-linked genes influence levels of resistance to infection in mice. // Exp Parasitol. V. 67. P. 221-226.

70. Fernando P., Pfrender M., Encalada S., Lande R. 2000. Mitochondrial DNA variation, phylogeography and population structure of the Asian elephant. // Heredity. V. 84. P. 362-372.

71. Fernandez-Vina M., Lazaro A., Marcos C., Nulf C., Raimondi E., Haas E., Stastny P. 1997. Dissimilar evolution of B-locus versus A-locus and class II loci of the HLA region in South American Indian tribes. // Tissue Antigens. V. 50. P. 233-250.

72. Frafjord К. 1993. Circumpolar size variation in the skull of the arctic fox Alopex lagopus II Polar Biology. V. 13. P. 235-238.

73. Frankham R. 1997. Do island populations have less genetic variation than mainland populations? // Heredity. V. 78. P. 311-327.

74. Frankham R. 1998. Inbreeding and extinction island populations // Conservation Biology. V 12. N.3. P. 665-675.

75. Freeman A.R., Machugh D.E., McKeown S., Walzer C., McConnell D.J., Bradley D.G. 2001. Sequence variation in the mitochondrial DNA control region of wild African cheetahs (Acynonyx jubatus). // Heredity. V. 86. P. 355-362.

76. Froeschke G., Sommer S. 2005. MHC Class II DRB constitution and parasite load in the striped mouse, Rhabdomys pumilio, in the Southern Kalahari. // Mol Biol Evol. V. 22. P. 1254-1259.

77. Gardenfors U. 2005. The 2005 red list of Swedish species. // Art Databanken, SLU, Uppsala.

78. Garrott R.A., Eberhardt L.E. 1982. Mortality of Arctic fox pups in northern Alaska // Journal of Mammalogy. V. 63. P. 173-174.

79. Geffen E., Waidyaratne S., Dalen L., Angerbjom A., Vila C., Hersteinsson P., Fuglei E., White P.A., Goltsman M.E., Kapel C.M.O., Wayne R.K. 2007. Sea ice occurrence predicts genetic isolation in the Arctic fox // Molecular Ecology. V. 16. P. 4241-4255.

80. Gocke C.D., Benko F.A., Rogan P.K. 1998. Transmission of mitochondrial DNA heteroplasmy in normal pedigrees. // Hum Genet. V. 102. P. 182-186.

81. Goda N., Mano Т., Masuda R. 2009. Genetic Diversity of the MHC Class-II DQA Gene in Brown Bears (Ursus arctos) on Hokkaido, Northern Japan // Zoological science. V. 26. P. 530-535.

82. Goldstein D.B., Roemer G.W., Smith D.A., Reich D.E., Bergman A., Wayne R.K. 1999. The Use of Microsatellite Variation to Infer Population Structure and Demographic History in a Natural Model. // System Genetics. V. 151. P. 797-801.

83. Goltsman-M.E., Kruchenkova E.P., Macdonald D.W. 1996. The Mednyi Arctic Foxes: treating a population imperiled by disease // Oryx. V. 30. N 4. P. 251-258.

84. Goltsman M., Kruchenkova E.P., Sergeev S., Volodin I.A., Macdonald D.W. 2005. "Island syndrome" in a population of Arctic foxes (Alopex lagopus) from Mednyi Island // Journal of Zoology. V. 267 N.4. P. 405-418.

85. Grzybowski T. 2000. Extremely high levels of human mitochondrial DNA heteroplasmy in single hair roots. // Electrophoresis. V. 21. N. 3. P. 548-53.

86. Harf R., Sommer S. 2005. Association between MHC Class II DRB alleles and parasite load in the hairy-footed gerbil, Gerbillurus paeba, in the Southern Kalahari. // Mol Ecol. V. 14. P. 85-91.

87. Hedrick P.W. 2004a. Comment on "parasite selection for immunogenetic optimality". // Science. V. 303. P. 957.

88. Hedrick P.W. 2004b. Evolutionary genomics: Foxy MHC selection story. // Heredity. V. 93. P. 237-238.

89. Hedrick P.W. 1994. Evolutionary genetics of the major histocompatibility complex // Am. Nat. V. 143. P. 945-964.

90. Hedrick P.W. 1999. Balancing selection and MHC // Genetica. V. 104. P. 207-214.

91. Hedrick P.W., Gutierrez-Espeleta G.A., Lee R.N. 2001. Founder effects in an island population of bighorn sheep. // Mol Ecol. V. 10. P. 851-857.

92. Hedrick P.W., Lee R.N., Parker K.M. 2000. Major histocompatibility complex (MHC) variation in the endangered Mexican wolf and related canids. // Heredity. V. 85. P. 617-624.

93. Hedrick P.W., Parker K.M., Gutierrez-Espeleta G.A., Rattink A., Lievers K. 2000. Major histocompatibility complex variation in the Arabian oryx. // Evol. V. 54. P. 2145-2151.

94. Hedrick P.W., Thomson G. 1988. Maternal-fetal interactions and the maintenance of HLA polymorphism. // Genetics. V. 119. P. 205-212.

95. Hess C.M., Edwards S.V. 2002. The evolution of the major histocompatibility complex in birds. // Bioscience. V. 52. P. 423-431.

96. Hingston M., Goodman S.M., Ganzhorn J.U., Sommer S. 2005. Reconstruction of the colonization of southern Madagascar by introduced Rattus rattus. //J Biogeography. V. 32. P. 1549-1559.

97. Hoelzel A.R., Stephens J.C., O'Brien S.J. 1999. Molecular genetic diversity and evolution at the MHC DQB locus in four species of pinnipeds. // Mol Biol Evol. V. 16. P. 611-618.

98. Horai S., Hayasaka K., Kondo R., Tsugane K., Takahata N. 1995. Recent African origin of modern humans revealed by complete sequences of hominoid mitochondrial DNAs. // Proc Natl Acad Sci USA. V. 92. P. 532536.

99. Hoss M., Dilling A., Currant A., Paabo S. 1996. Molecular phylogeny of the extinct ground sloth Mylodon darwinii. I I Proc. Natl. Acad. Sci. V. 93. P. 181185.

100. Huang S.W., Yu H.T. 2003. Genetic Variation of Microsatellite Loci in the Major Histocompatibility Complex (MHC) Region in the Southeast Asian House Mouse (Mus musculus castaneus). II Genetica. V. 119. N. 2. P. 201218.

101. Hughes A.L., Hughes M.K. 1995. Natural selection on the peptide-binding regions of major histocompatibility complex molecules. // Immunogenet. V. 42. P. 233-243.

102. Hughes A.L., Nei M. 1988. Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection. // Nature. V. 335. P. 167-170.

103. Hughes A.L., Nei M. 1989. Nucleotide substitution at major histocompatibility complex class II loci: Evidence for overdominant selection. // Proc Natl Acad Sci USA. V. 86. P. 948-962.

104. Hughes A.L., Yeager M. 1998. Natural selection at major histocompability complex loci of vertebrates // Annu. Rev. Genet. V. 32. P. 415-435.

105. Jarvi S.I., Tarr C.L., Mcintosh C.E., Atkinson C.T., Fleischer R.C. 2004. Natural selection of the major histocompatibility complex (MHC) in Hawaiian honeycreepers (Drepanidinae). // Mol Ecol. V. 13. P. 2157-2168.

106. Jeffery K.J., Bangham C.R. 2000. Do infectious diseases drive MHC diversity?//Microb Infect. V. 2. P. 1335-41.

107. Johnson. 2003. Island biogeography and evolution genetic divergence and speciation of island taxa // Comments on Theoretical Biology. V. 8. N. 1. P. 339-356. •

108. Klein J. 1986. Natural History of the Major Histocompatibility Complex. // Wiley & Son, New York.

109. Klein J., Horejsi V. 1997. // Immunology. Oxford, Blackwell Science.

110. Kmiec В., Woloszynska M., Janska H. 2006. Heteroplasmy as a common state of mitochondrial genetic information in plants and animals. // Curr Genet. V. 50. P. 149-159.

111. Koskinen M.T., Nilsson J., Veselov A.J., Potutkin A.G., Rantal E., Primmer C.R. 2002. Microsatellite data resolve phylogeographic patterns in European grayling, Thymallus thymallus, Salmonidae. //Heredity. V. 88. P. 391-401.

112. Kruchenkova E.P., Goltsman M.E., Sergeev S.N., Macdonald D.W. 2009. Is alloparenting helpful for Mednyi Island arctic foxes, Alopex lagopus semenovi? //Naturwissenschaften. V. 96. N. 4. P. 457-466.

113. Lebas N.R. 2001. Microsatellite determination of male reproductive success in a natural population of the territorial ornate dragon lizard, Ctenophorus ornatus. //Molecular Ecology. V. 10. N. 1. P. 193-203.

114. Lemay M.A., Bouldin E.G. 2009. Microsatellite pedigree analysis reveals high variance in reproductive success and reduced genetic diversity in hatchery-spawned northern abalone. // Aquaculture, Blackwell Science, Ltd. V. 295. N. 1-2. P. 22-29.

115. Levene H. 1949. On a matching problem arising in genetics. // Ann. Math. Stat. V. 20. P. 91-94.

116. Levinson G., Gutman G.A. 1987. Slipped-Strand Mispairing: A Major Mechanism for DNA Sequence Evolution. // Mol. Biol. Evol. V. 4. N. 3. P. 203-221.

117. MacArthur R.H., Wilson E.O. 1967. The theory of island biogeography. Princeton, NJ: Princeton University Press.

118. Matson C.W., Baker R.J. 2001. DNA sequence variation in the mitochondrial control region of red-backed voles (Clethrionomys). // Molecular Biology and Evolution! V. 18. P. 1494-1501.

119. Matsuda M., Yonekawa H., Hamaguchi S., Sakaizumi M. 1997. Geographic Variation and Diversity in the Mitochondrial DNA of the Medaka, Oryzias latipes, as Determined by Restriction Endonuclease Analysis. // Zoological Science. V. 14. N. 3. P. 517-526.

120. McClelland E.E., Granger D.L., Potts W.K. 2003. Major histocompatiblity complex-dependent susceptibility to Cryptococcus neoformans in mice. // Infect Immun. V. 71. P. 4815-4817.

121. Meagher S., Penn D.J., Potts W.K. 2000. Male male competition magnifies inbreeding depression in wild house mice. // Proc Natl Acad Sci USA. V. 97. P. 3324-3329.

122. Mellersh C.S., Langston A.A., Acland G.M. Fleming M.A., Ray K., Wiegand N.A., Francisco L.V., Gibbs M., Aguirre G.D., Ostrander E.A.I997. A linkage map of the Canine genome. // Genomics. V. 46. P. 326-336.

123. Meyer-Lucht Y., Sommer S. 2005. MHC diversity and the association to nematode parasitism in the yellow necked mouse (Apodemus flavicollis). И Mol Ecol. V. 14. P. 2233-2243.

124. Mikko S., Andersson L. 1995. Low major histocompatiblity complex class II diversity in European and North American moose. // Proc Natl Acad Sci USA. V. 92. P. 4259-4263.

125. Mikko S., Roed K., Schmutz S., Andersson L. 1999. Monomorphism andpolymorphism at MHC DRB loci in domestic and wild ruminants. // Immunol Rev. V. 167. P. 169-178.

126. Millar C.D., Dodd A., Anderson J., Gibb G.C., Ritchie P.A., Baroni C., Woodhams M., Hendy M.D., Lambert D.M. 2008. Mutation and evolutionary rates in Adelie penguins from the Antarctic. // PLoS Genet. V.4. P. 20-29.

127. Miller H.C., Lambert D.M. 2004. Genetic drift outweighs balancing selection in shaping post-bottleneck major histocompatiblity complex variation in New Zealand robins (Petroicidae). // Mol Ecol. V. 13. P. 3709-3721.

128. Moum Т., Вакке I. 2001. Mitochondrial control region structure and single , site heteroplasmy in the razorbill (Alea torda\ Aves). // Current Genetics. V.39. P. 198-203.

129. Nabata D., Masuda R., Takahashi O. 2004. Bottleneck Effects on the Sika Deer Cervus nippon Population in Hokkaido, Revealed by Ancient DNA Analysis. //Zoological science. V. 21. P. 473-481.

130. Nei M. 1987. Molecular evolutionary genetics. // Columbia University Press, New York.

131. Nevo E., Kirzher V., Beiles A., Karol A. 1997. Selection versus random drift: long-term polymorphism persistence in small populations (evidence and modelling). // Phil Trans R Soc Lond В. V. 352. P. 381-389.

132. Nigenda-Morales S. Flores-Ramirez S., Urban-R. J., Vazquez-Juarez R. 2008. MHC DQB-1 Polymorphism in the Gulf of California Fin Whale (Balaenoptera physalus) Population. // Journal of Heredity. V. 99. N. 1. P. 1421.

133. Norris A.T., Bradley D.G., Cunningham E.P. 2000. Parentage and relatedness determination in farmed atlantic salmon (Salmo salar) using microsatellite markers. // Aquaculture. V. 182. N. 1-2. P. 73-83.

134. Nystrom V., Angerbjorn A., Dalen L. 2006. Genetic consequences of a demographic bottleneck in the Scandinavian Arctic fox. // OIKOS. 00. P. 0000.

135. O'Brien S.J., Wildt D.E., Goldman D., Merril C.R., Bush M. 1985. The cheetah is depauperate in genetic variation. // Science. V. 221. P. 459-462.

136. Ostrander E.A., Mapa F.A., Yee M., Rine J. 1995. One hundred and one simple sequence repeate-based markers for the canine genome. // Mamm. Gen. V. 6. P. 92-195.

137. Ostrander E.A., Sprague J.R., Rine G.F.J. 1993. Identification and characterization of dinucleotide repeat (CA)n markers for genetic mapping in dog. // Genomics. V. 16. P. 207-213.

138. Paabo S. 1989. Ancient DNA: extraction, characterisation, molecular cloning and enzymatic amplification. // Proc. Natl. Acad. Sci. USA. V. 86. P. 19391943.

139. Paetkau D., Calvert W., Stirling I., Strobeck C. 1995. Microsatellite analysis of population structure in Canadian polar bears. // Mol Ecol. V. 4. N. 3. P. 347-54.

140. Paterson S., Wilson K., Pemberton J.M. 1998. Major histocompatibility complex variation associated with juvenile survival and parasite resistance in a large unmanaged ungulate population (Ovis aries L.). // Evol. V. 95. P. 3714-3719.

141. Parham P., Ohta T. 1996. Population biology of antigen presentation by MHC class I molecules // Science. Y. 272. P. 67-74.

142. Penn, D.J. 2002. The scent of genetic compatibility: sexual selection and the major histocompatibility complex. //Ethology. V. 108. P. 1-21.

143. Penn D.J., Damjanovich K., Potts W.K. 2002. MHC heterozygosity confers a selective advantage against multi-strain infections. // Proc Natl Acad Sci. V.99. P. 11260-11264.

144. Penn D.J., Potts W.K. 1998. Chemical signals and parasite-mediated sexual selection. // Trends in Ecology & Evolution. V. 13. P. 391-396.

145. Ploegh H., Watts C. 1998. Antigen recognition. // Curr. Opin. Immunol. V. 10. P. 57-58.

146. Potts W.K., Wakeland E.K. 1990. Evolution of diversity at the major histocompatibility complex // Trends Ecol. Evol. V. 5. P. 181-187.

147. Queller D.C. 1993. Genetic relatedness and its components in polygynous colonies of social insects. In: Keller L (ed) Queen number and sociality. // Oxford University Press, Oxford. P. 132-152.

148. Radwan J., Biedrzycka A., Babik W. 2010. Does reduced MHC diversity decrease viability of vertebrate populations? // Biological Conservation. V. 143. P. 537-544.

149. Radwan J., Kawalko A., Wyjcik J.M., Babik W. 2007. MHC-DRB3 variation in a free-living population of the European bison, Bison bonasus I I Molecular Ecology.V. 16. P. 531-540.

150. Richman A. 2000. Evolution of balanced genetic polymorphism. // Mol Ecol. V. 9. P. 1953-1963.

151. Robertson A. 1962. Selection for heterozygotes in small populations. // Genetics. V. 47. P. 1291-1300.

152. Saitou N., Nei M. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees // Mol. Biol. Evol. V. 4. P. 406-425.

153. Sambrook J., Fritsch E.F., Maniatis T. 1989. Molecular cloning: A laboratory manual. 2 ed. // Cold Spring Harbor laboratory Press. Prainview. NY. 320 p.

154. Santucci F., Ibrahim K.M., Bruzzone A., Hewit G.M. 2007. Selection on МНС-linked microsatellite loci in sheep populations. // Heredity. V. 99. P. 340-348.

155. Sarmiento U.M., Storb R.F. 1988. Restriction fragment length polymorphism of the major histocompatibility complex of the dog. // Immunogenetics. V. 28. P. 117-124.

156. Schad J., Ganzhorn J.U., Sommer S. 2005. MHC constitution and parasite burden in the Malagasy mouse lemur, Microcebus murinus. II Evol. V. 59. P. 439-450.

157. Schad J., Sommer S., Ganzhorn J.U. 2004. MHC variability of a small lemur in the littoral forest fragments of southeastern Madagascar. // Conserv Genet. V. 5. P. 299-309.

158. Schmitt L.H. 1978. Genetic variation in isolated populations of the Australian bush rat Rattus fuscipes. // Evolution. V. 32. P. 1-14.

159. Schneider S., Roessli D., Excoffier L. 2000. ARLEQUIN. Version 2.000 // University of Geneva. Switzerland.

160. Seddon J.M., Baverstock P.R. 1999. Variation on islands: major histocompatibility complex (Mhc) polymorphism of the Australian bush rat. // Molecular Ecology. V. 8. P. 2071-2079.

161. Seddon J.M., Ellegren H. 2002. MHC class II genes in European wolves: a comparison with dogs. // Immunogenetics. V. 54. P. 590-600.

162. Selkoe K.A., Toonen R.J. 2006. Microsatellites for ecologists: a practical guide to using and evaluating microsatellite markers. // Ecology Letters. V. 9. P. 615-629.

163. Senseney H.L., Briles W.E., Abplanalp H., Taylor R.L. 2000. Allelic complementation between MHC haplotypes B-Q and B-17 increases regression of Rous sarcomas. // Poultry Sci. V. 79. P. 1736-1740.

164. Sher A., Hieny S., James S. 1984. Mechanisms of protective immunity against S. mansoni infection in mice vaccinated with irradiated cercariae: VI. Influence of the major histocompatibility complex. // Parasit Immunol. V. 6. P. 319-328.

165. Shriver M.D., Jin L., Ferrell R.E., Deka R. 1997. Microsatellite Data Support an Early Population Expansion in Africa. // Genome Res. V. 7. P. 586-591.

166. Simonsen B.T., Siegismund H.R., Arctander P. 1998. Population structure of African buffalo inferred from mtDNA sequences and microsatellite loci: high variation but low differentiation. // Molecular Ecology. V. 7. N. 2. P. 225-237.

167. Slamovits C.H., Rossi M.S. 2002. Satellite DNA: agent of chromosomal evolution in mammals. A review. // Mastozoologia Neotropical. V. 9. P. 297308.

168. Slatkin M. 1994. Linkage disequilibrium in growing and stable populations. Genetics. V. 137. P. 331-336.

169. Slierendregt B.L., Otting N., Van Besouw N., Jonker M., Bontrop R.E. 1997. Expansion and contraction of rhesus macaque DRB regions by duplication and deletion. // The Journal of Immunology. V. 152. N. 5. P. 2298-2307.

170. Smulders M.J.M., Snoek L.B., Booy G., Vosman B. 2003. Complete loss of

171. MHC genetic diversity in the Common hamster (Cricetus cricetus) population in the Netherlands. Consequences of conservation strategies. // Conserv Genet. V. 4. P. 441-451.

172. Sommer S. 2005a. МНС-dependent mate choice in a monogamous rodent. // Behav Ecol Sociobiol. V. 58. P. 181-189.

173. Sommer S. 2005b. The importance of immune gene variability (MHC) in evolutionary ecology and conservation. // Frontiers in Zoology. V. 2. P. 16-29.

174. Sommer S., Tichy H. 1999. Major histocompatibility complex (MHC) class II polymorphism and paternity in the monogamous Hypogeomys antinema, the endangered, largest endemic Malagasy rodent. // Mol Ecol. V. 8. P. 12591272.

175. Sommer S., Schwab D., Ganzhorn J.U. 2002. MHC diversity of endemic Malagasy rodents in relation to range contraction and social system. // Behav Ecol Sociobiol. V. 51. P. 214-221.

176. Spencer C.C., Neigel J.E., Leberg P.L. 2000. Experimental evaluation of the usefulness of microsatellite DNA for detecting demographic bottlenecks. // Molecular Ecology. V. 9. P. 1517-1528.

177. Sugg D.W., Chesser R.K., Dobson F.S., Hoogland J.L. 1996. Population genetics meets behavioural ecology. // Trends Ecol Evol. V. 11. P. 338-342.

178. Sunnucks P. 2000. Efficient genetic markers for population biology. // TREE. V. 15. P. 199-203.

179. Tajima F. 1983. Evolutionary relationship of DNA sequences in finite populations // Genetics. V. 105. P. 437-460.

180. Takahata N. 1990. A simple genealogical structure of strongly balanced allelic lines and trans-species evolution of polymorphism. // Proc. Natl. Acad. Sci. USA. V. 87. P. 2419-2423.

181. Takahata .N., Nei M. 1990. Allelic genealogy under overdominant and frequency-dependent selection and polymorphism of major histocompatibility complex loci. // Genetics. V. 124. P. 967-978.

182. Tannerfeldt M., Angerbjorn A. 1998. Fluctuating resources and the evolution of litter size in the Arctic fox // Oikos. V. 83. P. 545-559.

183. Tamaki K., Jeffreys A.J. 2005. Human Tandem Repeat Sequences in Forensic DNA Typing. // Journal of Legal Medicine. V. 7. P. 244-250.

184. Tamura K., Dudley J., Nei M., Kumar S. 2007. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. // Molecular Biology and Evolution. V. 24. P. 1596-1599.

185. Tamura K., Nei M. 1993. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. // Molecular Biology and Evolution. V. 10. P. 512-526.

186. Tannerfeldt M. 1997. Population fluctuations and life history strategies in the arctic fox. // PhD Thesis, Stockholm University, Sweden.

187. Tautz D., Schlotterer C. 1994. Simple sequences. // Curr. Opin. Genet. Dev. V. 4. P. 832-837.

188. Thomas M.I., Harger J.H., Wagener D.K., Rabin B.S., Gill T.J. 1985. HLA sharing and spontaneous abortion in humans. // Am. J.Obslet. Gynecol. V. 151. P. 1053-1058.

189. Thursz M.R., Kwiatkowski D., Allsopp C.E. Greenwood B.M., Thomas H.C., Hill АЛЛ 1995. Association between an MHC class II allele and clearance of hepatitis В virus in the Gambia. // New Engl J Med. V. 332. P. 1065-1069.

190. Thursz M.R., Thomas H.C., Greenwood B.M., Hill A.V. 1997. Heterozygote advantage for HLA class-II type in hepatitis В virus infection. // Nat Genet. V. 17. P. 11-12.

191. Van der Geer A., Lyras G., de Vos J., Dermitzakis M. 2010. Evolution of Island Mammals: Adaptation and Extinction of Placental Mammals on Islands. Wiley-Blackwell, 384 p.

192. Van Eden W., Persijn G.G., Bijkerk H., de Vries R.R.P., Schuurman R.K.B., van Rood J.J. 1983. Differential resistance to paralytic poliomyelitis controlled by histocompatibility leukocyte antigens. // J. Infect. Dis. V. 147. P. 422-426.

193. Van Oosterhout C.5 Joyce D.A., Cummings S.M. 2006. Evolution of MHC class IIB in the genome of wild and ornamental guppies, Poecilia reticulate. II Heredity. V. 97. P.l 11-118.

194. Wagner J.L. 2003. Molecular Organization of the Canine Major Histocompatibility Complex. // Journal of Heredity. V. 94. P. 23-26.

195. Wagner J.L., Burnett R.C., Storb R. 1996. Molecular analysis of the DLA-DR region. // Tissue Antigens. V. 48. P. 549-553.

196. Wagner J.L., Hayes-Lattin В., Works J.D., Storb R. 1998. Molecular analysis and polymorphism of the DLA-DQB genes. // Tissue Antigens. V. 52. P. 242250.

197. Waits L.P., Talbot S.L., Ward R.H., Shields G.F. 1998. Mitochondrial DNA Phylogeography of the North American Brown Bear and Implications for Conservation // Conservation Biology. V. 12. N. 2. P. 408-417.

198. Wakeland E.K., Boehme S., She J.X. 1990. The Generation and Maintenance of MHC Class II Gene Polymorphism in Rodents. // Immunological Reviews. V. 113. P.207-226

199. Wassom D.L., Brooks B.O., Babisch J.G., David C.S. 1983. A gene mapping between the S and D regions of the H-2 complex influences resistance to Trichinella spiralis infections in mice. // J Immunogenetics. V. 10. P. 371378.

200. Wassom D.L., Krco C.J., David C.S. 1987. 1-Е expression and susceptibility to parasite infection. // Immunol Today. V. 8. P. 39-43.

201. Weber D.S., Stewart B.S., Schienman J., Lehman N. 2004. Major histocompatibility complex variation at three class II loci in the northern elephant seal. // Mol Ecol. V. 13. P. 711-718.

202. Willerslev E., Cooper A. 2005. Ancient DNA. // Proc. R. Soc. Lond. V. 272. P. 3-16.

203. White P.A. 1992. "Social organization and activity patterns of Arctic foxes on St. Paul Island, Alaska". M.Sc. dissertation, University of California. Berkeley, California, USA.

204. Yang D., Eng В., Waye J., Dudar J., Saunders S. 1998. Technical note: improved DNA extraction from ancient bones using silica based spin columns. // American Journal of Physical Anthropology. V. 105. P. 539-543.

205. Yin В., Yang S., Wei W., Zhang Y. 2009. Male reproductive success in plateau pikas (Ochotona curzoniae): A microsatellite analysis. // Mammalian Biology Zeitschrift fur Saugetierkunde. V. 74. N. 5. P. 344-350.

206. Yoder A.D., Delefosse T. 2002. "The Rise and Fall and Rise of Ancient DNA Studies" in Ancient DNA, McGraw-Hill, // Yearbook of Science and Technology. P. 9-14.

Информация о работе

Плошница, Анна Ивановна
кандидата биологических наук
Москва, 2010
ВАК 03.02.04

Диссертация

Влияние островной изоляции и эффекта "бутылочного горлышка" на генетический полиморфизм командорских песцов - тема диссертации по биологии, скачайте бесплатно

Автореферат

Похожие работы