Бесплатный автореферат и диссертация по биологии на тему
Генетическая дифференциация азиатских популяций тихоокеанского лосося-чавычи, Oncorhynchus Tschawytscha (Walbaum)
ВАК РФ 03.02.07, Генетика
Автореферат диссертации по теме "Генетическая дифференциация азиатских популяций тихоокеанского лосося-чавычи, Oncorhynchus Tschawytscha (Walbaum)"
/ '
На правах рукописи
Шпигальская Нина Юрьевна
ГЕНЕТИЧЕСКАЯ ДИФФЕРЕНЦИАЦИЯ АЗИАТСКИХ ПОПУЛЯЦИЙ ТИХООКЕАНСКОГО ЛОСОСЯ - ЧАВЫЧИ, ONCOSHШCHVS Т8СПАШТ8СНА
(\VALBAUM)
03.02.07 - генетика
Автореферат
диссертации на соискание ученой степени кандидата биологических наук
Петропавловск-Камчатский - 2010
004603482
Работа выполнена в лаборатории популяциониой биологии и генетики Камчатского научно-исследовательского института рыбного хозяйства и океанографии (ФГУП «КамчатНИРО»)
Научный руководитель:
доктор биологических наук
Варнавская Наталья Васильевна
Официальные оппоненты:
доктор биологических наук
Балакирев Евгений Станиславович
кандидат биологических наук
Зеленина Дарья Александровна
Ведущая организация:
Институт общей генетики им. Н.И. Вавилова РАН
Защита состоится 4 июня 2010 г. в 11 часов на заседании диссертационного совета Д.005.008.01 при Институте биологии моря им. A.B. Жирмунского ДВО РАН по адресу: 690041, г. Владивосток, ул. Пальчевского, 17. Телефон: (4232) 310-905, факс: (4232) 310-900, e-mail: irmarbio@mail.primorye.ru
С диссертацией можно ознакомиться в библиотеке Института биологии моря им A.B. Жирмунского ДВО РАН (690041, г. Владивосток, ул. Пальчевского, 17). Отзывы просим присылать на e-mail: mvaschenko@mail.ru
Автореферат разослан 29 апреля 2010 г.
Ученый секретарь диссертационного совета, кандидат биологических наук
Ващенко М.А.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы. Тихоокеанские лососи рода Oncorhynchus, являясь ценными объектами промысла, имеют большое экономическое значение для стран Тихоокеанского бассейна. Для всех представителей рода характерна сложная популяционная организация, основными чертами которой являются соответствие принципу иерархии и упорядоченная во времени и пространстве субпопуляционная структура (Altukhov, 1971, 1981, 1990; Алтухов, 1973, 1974; Altukhov, Salmenkova, 1981, 1991; Алтухов и др., 1997; Altukhov et al., 2000). Эмпирическое подтверждение существования в пределах видов генетически отличающихся, репродуктивко изолированных единиц, дает основание для утверждения, что успех краткосрочной стратегии управления ресурсами и долговременных целей сохранения видов зависит от наличия, наряду с демографическими и экологическими данными, информации о генетической внутривидовой структуре (Алтухов и др., 1997; Алтухов и др., 2004). Для решения теоретических проблем популяционной организации в сочетании с долговременным мониторингом популяционных генофондов и прикладных задач генетической идентификации также необходимы данные о генетическом своеобразии внутривидовых структурных компонентов и оценки уровня их дифференциации.
Чавыча, О. tschawytscha (Walbaum), является одним из наиболее ценных видов тихоокеанских лососей, который в азиатской части ареала обитает преимущественно на Камчатке. Популяционно-генетические исследования чавычи Северной Америки, где данный вид распространен практически по всему тихоокеанскому побережью, были начаты во второй половине 20 в. и показали высокий уровень межпопуляционной изменчивости по аллозимным локусам (Kristiansson, Mclntyre, 1976; Gharrett et al., 1987; Reisenbichler, Phelps, 1987; Beacham et al., 1989, 1996; Utter et al., 1989; Bartley, Gall, 1990; Gall et al., 1992; Crane et al., 1996; Teel, 2000, 2004; Templin, Seeb, 2003), структуре митохондриальной ДНК (Wilson et al., 1987; Cronin et al., 1993; Adams et al., 1994), микросателлитам (Moran et al., 2005; Beacham et al., 2006b; Seeb et al., 2007) и единичным нуклеотидным заменам (SNPs) (Smith et al., 2007; Narum et al., 2008). Были созданы объединенные реперные базы данных по частотам аллозимных локусов, на основе которых с успехом осуществлялась идентификация смешанных уловов в прибрежных водах США и Канады (Teel et al., 1999), и базы реперных данных по аллельной изменчивости микросателлитных локусов (Moran et al., 2005; Seeb et al., 2007). Частоты аллозимных генов служили также маркером для оценки успешности деятельности рыбоводных заводов (Waples, 1990a,b; Waples, Do, 1994).
В Азии исследования генетической структуры популяций чавычи на основе биохимической изменчивости и полиморфизма мтДНК находятся на начальном этапе, анализ микросателлитной ДНК и полиморфизма единичных нуклеотидных замен пока не начат.
До недавнего времени для выявления внутри- и межпопуляционной дифференциации и оценки биологического разнообразия азиатских популяций тихоокеанских лососей наиболее широко использовали такие маркеры генетической изменчивости, как полиморфные белковые локусы (Алтухов, 2003; Алтухов и др., 1997; Варнавская, 2001), но в последние десятилетия все большее значение приобретают исследования на уровне ДНК (Avise, 1994; Брыков, 2001; Чуриков, 2001; Алтухов, Салменкова, 2002). Внедрение молекулярно-генетических методов в популяционные исследования позволяет дополнить и расширить уже сформированную картину генетической подразделенности и процессов эволюции видов. Анализ изменчивости нуклеотидной последовательности может дать дополнительное количество высокоинформативных маркеров, что особенно актуально при исследовании видов со сложной генетической структурой и отсутствием четкой корреляции между генетическим и географическим расстояниями.
Исследования популяционной структуры тихоокеанских лососей с использованием молекулярно-генетических методов в силу их больших и принципиально новых возможностей получили в настоящее время широкое распространение. Анализ изменчивости неядерной наследственной структуры - митохондрнальной ДНК является эффективным инструментом оценки генетической дивергенции и выявления филогенетической истории популяций (Полякова и др., 1992; Brykov et al., 1996; Брыков и др., 1999а,б; Брыков и др., 2002; Брыков и др., 2005). Исследования полиморфизма мтДНК чавычи могут дать дополнительную информацию о внутри- и межпопуляционной генетической структуре данного вида, позволят выявить филогенетические связи между внутривидовыми структурными единицами, а также увеличить число информативных генетических маркеров, пригодных для различных аспектов популяционных исследований.
Цель н задачи исследования. Целью настоящей работы является оценка уровня генетического разнообразия азиатских популяций чавычи на основе биохимического полиморфизма и гаплотипической изменчивости митохондрнальной ДНК, а также создание баз генетических данных для последующей идентификации региональной принадлежности чавычи в смешанных морских скоплениях.
Для достижения поставленной цели были определены следующие задачи: . 1. изучить биохимический полиморфизм и изменчивость последовательностей мтДНК в популяциях чавычи азиатской части ареала;
2. выявить максимально возможное число информативных генетических маркеров в популяциях чавычи;
3. выделить набор наиболее дискриминирующих азиатскую чавычу генетических маркеров для идентификации ее в смешанных уловах;
4. оценить внутри- и межпопуляционную изменчивость генетической структуры чавычи Азии;
5. исследовать генетическую дивергенцию азиатских и североамериканских популяционных комплексов.
Научная новизна. Впервые получены результаты исследования полиморфизма митохондриальной ДНК у чавычи Камчатки. Анализ гаплотипической изменчивости позволил получить оценки внутри- и межпопуляционной дифференциации, а также уровня дивергенции популяций азиатской части ареала. Выявлено, что максимальное время независимой эволюции камчатских популяций чавычи соответствует оценке во временном интервале, не превышающем 10-20 тыс. лет. Впервые получены сравнительные оценки генетического разнообразия по аллозимным локусам чавычи различных частей ареала. Характер пространственной изменчивости аллельных частот, уровень внутри- и межпопуляционной дифференциации по исследованным аллозимным локусам позволил выделить ряд наиболее информативных для идентификации азиатских популяций чавычи в смешанных морских скоплениях.
Практическая значимость. Результаты, представленные в работе, могут быть использованы для создания баз реперных данных по двум типам исследованных генетических маркеров в целях идентификации смешанных морских скоплений чавычи, а также генетического мониторинга популяционных процессов и сохранения естественного биоразнообразия природных популяций данного вида в условиях морского и берегового промысла. Помимо этого, полученные результаты на практике могут быть использованы при организации искусственного разведения чавычи для долговременного слежения за состоянием популяционных генофондов, оценки и прогнозирования их динамики во времени и пространстве.
Апробация работы. Основные положения и результаты работы были представлены на Международной конференции «Новые технологии в защите биоразнообразия в водных экосистемах» (Москва, 27-29 мая 2002 г.), Международной рабочей группе Северотихоокеанской комиссии по анадромным видам рыб NPAFC International Workshope on Application of Stock Identification in Defining Marine Distribution and Migration of Salmon (Honolulu, Hawaii, USA, November 1-2, 2003), объединенной отчетной сессии НТК «ТИНРО» (Петропавловск-Камчатский, 2006 г.), Международной конференции «Ихтиологические исследования на внутренних водоемах» (Саранск, 2007), X Всероссийском популяционном семинаре «Современное состояние и пути развития популяционной биологии» (Ижевск, 1722 ноября 2008 г.), а также в виде научных годовых отчетов ФГУП «КамчатНИРО» в 2001 г. и 2003-2006 гг.
Публикации. По материалам диссертации опубликовано 10 работ, в том числе, I статья - в рецензируемом журнале из списка ВАК,
Объём и структура диссертации. Диссертация изложена на 145 страницах машинописного текста, состоит из Введения, Обзора литературы (гл. 1), главы 2 «Материалы и методы», результатов и обсуждения (гл. 3 и 4), Заключения, Выводов, Списка литературы, включающего 212 ссылок, из которых 152 - на английском языке и 4-х Приложений. Работа иллюстрирована 23 рисунками и содержит 19 таблиц.
Благодарности. Автор искренне благодарит своего научного руководителя Наталью Васильевну Варнавскую за конструктивное руководство, создание творческой атмосферы и постоянный научный интерес. Автор выражает глубокую признательность заведующему лабораторией генетики Института биологии моря им. А.В Жирмунского ДВО РАН Владимиру Алексеевичу Брыкову за внимательное отношение, консультации и предоставленную возможность обработки материала, а также всему коллективу сотрудников, работающему под его руководством за помощв и поддержку в работе. Автор искренне признателен P.A. Шапореву за неоценимую помощь в статистической обработке данных и всем коллегам лаборатории популяционной биологии и генетики, а также сотрудникам КамчатНИРО, принимавшим участие в сборе материала в сложных условиях труднодоступных районов Камчатки.
СОДЕРЖАНИЕ РАБОТЫ Глава 1. Обзор литературы
В первом разделе обзора литературы приводится краткий очерк биологии чавычи (Oncorhynchus tschawytscha Walbaum), рассматриваются различные аспекты внутривидовой организации данного вида тихоокеанских лососей.
Во втором разделе приводится описание результатов исследований популяционной структуры тихоокеанских лососей на основе молекулярно-генетических методов, оценивается состояние изученности данного вопроса у чавычи азиатской и североамериканской частей видового ареала.
Глава 2. Характеристика материала и методов исследования
Материал исследования
В работе для анализа изменчивости на основе биохимического полиморфизма было исследовано тринадцать выборок из семи популяций чавычи Камчатки (638 экз.); для анализа популяционной изменчивости на основе полиморфизма мтДНК - тринадцать выборок (579 экз.) чавычи из восьми локальных популяций Западной Камчатки и пяти,
включая бас. Олюторского залива, Восточной Камчатки. Объёмы выборок представлены в табл. 1, а их локализация - на рис. 1-2.
Таблица 1. Локализация сбора и объём выборок чавычи Камчатки
Для
характеристики
генетической
Регион № Локальность Кол-во, экз. Год американской части ареала на
основе аллозимной изменчивости
биохимический полиморфизм использовали опубликованные
Северо-Западная 1 р. Хайрюзова 40 1997
Камчатка 2 р. Воровская 50 1994 ранее данные по частотам
3 р. Воровская 50 1995 белковых локусов в популяциях
Западная 4 р. Утка 42 1997
Камчатка 5 р. Быстрая 50 1998 Юго-Восточной и Западной
б р. Большая 50 1996 Аляски, о-ва Кодьяк, п-ва Аляска
Юго-Восточная 7 р. Большая 50 1998 (Crane et al., 1996), Орегона и
8 р. Авача 40 1997
Камчатка 9 р. Авача 40 1998 Калифорнии (Gall et al., 1992).
Восточная 10 р. Камчатка, поздняя 72 1997 Объединенная база данных
Камчатка 11 р. Камчатка, ранняя 64 1995
12 р. Камчатка, ранняя 50 1998 включает генетические
13 р. Камчатка, поздняя 40 1998 характеристики 7843 экз. чавычи
полиморфизм мтДНК из 100 локальных популяций.
Северо-Западная 1 р. Палана 50 2002
Камчатка 2 р. Тигиль 42 2002 Молекулярно-генетическне
3 р. Воровская 50 2003 методы анализа. Основной метод
Западная 4 р. Коль 50 2003 выявления биохимической
50 2002
Камчатка 5 р. Пымта
б р. Кихчик 33 2003 изменчивости, использованный в
7 р. Утка 50 2002 работе - электрофорез в
8 р. Большая 60 2002
Юго-Восточная 9 р. Авача 44 2003 крахмальном геле с последующим
Камчатка 10 р. Жупанова 30 2003 субстратспецифичным
Восточная Камчатка 11 р. Камчатка 60 2003 окрашиванием (Aebersold et al,
Северо-Восточная 12 Олюторский залив 50 2002
Камчатка 13 р. Пахача 10 2001 1987).
Для анализа гаплотипической изменчивости тотальную ДНК выделяли из ткани сердечной мышцы стандартным способом с использованием метода протеиназного гидролиза в присутствии додецилсульфата натрия с последующим высаливанием белков, удалением их вместе с клеточными обломками центрифугированием и осаждением ДНК из супернатанта изопропанолом (Маниатис и др., 1984; ЗашЬгоок е1 а!., 1989).
Изменчивость мтДНК исследовали, используя рестрикций нный анализ трех фрагментов, которые амплифицировали посредством полимеразной цепной реакции с использованием праймерных последовательностей, разработанных на других видах лососей -микиже, Опсогкупскш туШз, и кижуче, ОпсогИупскиз киШск (2агс1оуа е! а!., 1995; ОЬаггеИ
е! а1., 2001). Обозначение исследуемых в данной работе фрагментов и праймерные последовательности, использованные при амплификации, приведены в табл. 2.
Рис. 1. Локализация сбора проб Рис. 2. Локализация сбора проб
чавычи Камчатки, исследованных при чавычи Камчатки, исследованных при анализе
анализе биохимического полиморфизма (1 -р. Хайрюзова, 2-3 - р. Воровская, 4 - р. Утка, 5 - р. Быстрая, 6-7 - р. Большая, 8-9 -р. Авача, 10-13 - р. Камчатка).
полиморфизма мтДНК (1 - р. Пахача, 2 -Олюторский залив, 3 - р. Палана, 4 - р. Тигиль, 5 - р. Воровская, б-р. Коль, 7 - р. Пымга, 8 - р. Кихчик, 9 - р. Утка, 10 - р. Большая, 11 - р. Авача, 12-р. Жупанова, 13 -р. Камчатка).
Таблица 2. Последовательность праймеров и гены в
Структурные гены Последовательности
A8/A6/COIII/ND3 ND3/ND4L/ND4 COI/COII/A8 5' -CTAGTGACATGCCCCAACTCAACC-3' 5 '-TCATAAGGCGGTCATGGACTTAAACC-3' 5' -TTACGCGTATAAGTGACTTCCAA-3' 5' - TTTTGGTTCCTAAGACCAATGGAT-3' S'-TAATCGTCACAGCCCATGCCTTCGT-З' 5'-GGTCAGTTTCAGGGTTCAGGTTTAGC-3'
Изменчивость нуклео-тидной последовательности анализировали, используя набор
рестриктных ферментов: Я?гМ, ЯгяД МЬо\, Агв1, Ше1, Га?1, Мр1, &ги961.
После реакции эндонуклеазного гидролиза пробы подвергали электрофорезу в 1.8-2% агарозном геле в трис-боратном буфере. Для определения молекулярной массы образующихся фрагментов использовали маркерный набор фрагментов ДНК, кратных 100 парам оснований (BRL Gibco, Grand NY) и ДНК фага X, обработанную рестриктазой Pst\. Фрагменты ДНК в геле окрашивали этидиумбромидом и фотографировали в проходящем ультрафиолетовом свете.
Методы статистического анализа. Для статистической обработки результатов биохимических исследований применяли стандартные методы анализа частот генотипов и аллелей (Ли, 1978). Критерий х2 и G-тест (G-критерий) (Sokal, Rolf, 1981) применяли для
оценки внутри- и межпопуляционной генетической гетерогенности локальных стад и региональных популяционных комплексов. Как мера количественных различий между популяциями вычислялись генетические расстояния с использованием хордового метода (Cavalli-Sforsa, Edwards, 1967) и критерия идентичности М. Нея (Nei, 1987) с последующим построением UPGMA-дендрограмм (Sneath, Sokal, 1973). Другим подходом к оценке генетического сходства популяций был метод многомерного шкалирования (Kruskal, 1964).
Оценку количества замен оснований на нуклеотидный сайт и их стандартные отклонения между гаплотипами рассчитывали, используя пакет программ REAP (McElroy et al., 1992). Нуклеотидную дивергенцию между всеми парами гаплотипов определяли с помощью программы D, а дивергенцию между всеми парами популяций - с помощью программы DA (McElroy et al., 1992). Нуклеотидную и гаплотипическую величины разнообразия рассчитывали по М. Нею (Nei, Tajima, 1983; Nei, 1987) в пакетах программ REAP. Оценку величины дивергенции между выборками и их стандартные ошибки оценивали по М. Нею и Ч. Ли (Nei, Li, 1979). Полученная в REAP матрица количественных значений дивергенции между популяциями использовалась для кластерного анализа (программа UPGMA, пакет программ NTS YS) и построения фенограмм (Rolf, 1990). Для оценки значимости межпопуляционных различий частот гаплотипов использовали методы F-статистики в пакете программ Arlequin2000 (Schneider et al., 2000). Межпопуляционную гетерогенность оценивали с использованием псевдовероятностного теста (Zaykin, Pudovkin, 1993).
Глава 3. Биохимический полиморфизм в популяциях чавычи
Уровень изменчивости аллозимных популяционных маркеров. При анализе чавычи азиатской части ареала был определен уровень аллельной изменчивости 55 аллозимных генов. Доля полиморфных локусов составила меньше половины их общего числа - 20 (36.4%). Количество аллелей в подавляющем большинстве полиморфных локусов не превышало двух и в сумме составило 41. Только при анализе CPIA* в камчатских популяциях было выявлено три аллельных варианта.
При объединении наших данных по аллельной изменчивости чавычи Азии с опубликованными ранее материалами по популяциям чавычи Северной Америки (Gharrett et al., 1987; Bartley, Gall, 1990; Bartley et al., 1992; Gall et al., 1992; Utter et al., 1992; Wilmot et al., 1992; Adams et al., 1994; Crane et al, 1996; Quirin et al, 1996) количество полиморфных локусов оказалось равным 69 (73.4%) и значительно превысило данный показатель в выборках из азиатских популяций. Из общего числа мономорфных локусов (25) большая часть (20) у чавычи Камчатки исследована не была. Только пять локусов (СКА1*, FH*,
G3PDH1*, LDHA2*, LDHB1*) были определены как мономорфные во всех включенных в анализ выборках. Большинство локусов, мономорфных в Азии, оказались изменчивыми в выборках из локальных популяций Северной Америки.
Число аллелей полиморфных локусов в выборках на североамериканской части видового ареала чавычи варьировало от 2 до 7. Общее число аллелей у исследованных 7843 экз. чавычи составило 201, среднее число аллелей на локус оказалось равным 2.91.
Число полиморфных локусов в различных популяциях чавычи Камчатки менялось от 9 до 16, и, соответственно, диапазон варьирования был достаточно широким - от 16.4% до 29.1%. Этот факт согласуется с данными по североамериканским популяциям чавычи, где отмечен еще более широкий диапазон варьирования доли полиморфных локусов в отдельных популяциях (5.8-32.1%%) (Bartley, Gall, 1990).
Следует отметить, что некоторые аллозимные маркеры у азиатской чавычи являются полиморфными лишь в отдельных локальностях и фиксированы в остальных. Так, локус GR* полиморфен только в выборке из р. Хайрюзова, локусы I.DHC* и PEPD2*- в р. Большая, TPI3* - в р. Утка, а локус mAATl* изменчив только в двух популяциях - р. Камчатка и р. Быстрая. Подобные особенности полиморфизма ферментных локусов отмечены и для чавычи американской части ареала (Gharrett et al., 1987; Utter et al,, 1989; Bartley et al., 1992).
Для последующей оценки генетической изменчивости популяций и региональных комплексов азиатской чавычи было отобрано всего 16 полиморфных локусов, но и среди этого количества присутствовали локусы, частота основного аллеля которых в некоторых локальностях была равна 1.0.
Оценки генетического разнообразия, представленные в виде усредненной на локус гетерозиготности по аллозимным генам для разных регионов (HJ и для отдельных популяций региона (Hs), а также оценки относительного межпопуляционного разнообразия GsI приведены в табл. 3. Показатель общего генного разнообразия региональных популяционных систем чавычи, (Я,), усредненный по всем исследованным выборкам и Таблица 3. Оценки генетического разнообразия по аллозимным локусам у чавычи.
Регион Кол-во выборок Число локусов (полиморфные) Я, Hs Gsh %
п-ов Камчатка 13 47 (20) 0.0633 0.0615 2.9
Юго-Восточная Аляска 18 41 (28) 0.0764 0.0725 5.1
Западная Аляска . ' • 32 41 (29) 0.0533 0.0505 5.4
Нортон Саунд, шт. Орегон 6 53 (41) 0.1065 0.1012 5.0
Шт. Южный Орегон и 5 53 (45) 0.0858 0.0840 2.2
Северная Калифорния
Бассейны р. Кламат и 10 53 (48) 0.0613 0.0577 5.9
р. Тринити
Центральная Калифорния 11 53 (47) 0.0701 0.0681 2.8
Бассейн р. Сакраменто 5 53 (29) 0.0884 0.0870 1.6
по всей совокупности аллозимных локусов, оказался равен 0.0784, среднее значение величины - 0.0922. Уровень межпопуляционного разнообразия чавычи Аляски и некоторых южных регионов Северной Америки несколько выше, чем у чавычи Камчатки. Используя собственные и ранее опубликованные данные, были сопоставлены величины генного разнообразия на разных уровнях внутривидовой иерархии (табл. 4).
Как видно из приведенных оценок дифференциация региональных групп популяций (соответствуют табл. 3) заметно выше уровня
В соответствии с распределениями частот аллелей были выявлены достаточно четко выраженные отличия групп популяций по большому числу аллозимных локусов. Так, можно отметить, что по локусам sAATl,2* и mAATI* аллель 100*, фиксирован почти на всем ареале кроме Северной Калифорнии и Орегона, а также Юго-Восточной Аляски. По локусу PEPLT* наиболее изменчивы популяции Юго-Восточной Аляски и Западной Аляски. По локусу sAHl* все регионы изменчивы, кроме Западной Аляски и Центральной Калифорнии, где этот локус почти фиксирован. По локусу GPIB2* изменчивы популяции Калифорнии и Орегона, во всех остальных частях ареала он почти мономорфен. Локус IDDH1 * является наиболее дискриминирующим для популяций Азии: в большинстве популяций Северной Америки он мономорфен, но высокополиморфен как в популяциях Западной, так и Восточной Камчатки. Локусы PEPD2*, sMDHBl,2* sMEP2* PEPA*, РЕРВ1*, sIDHP2* также мономорфны в одних регионах и полиморфны в других. Локусы TPI4* и sAAT3* являются дифференцирующими для популяций Калифорнии, sMEPl* - для популяций Западной и Северной Аляски, а также Азии, PGK2* - для популяций Азии и Аляски, sJDHPl *, FGDH* -для популяций Юго-Восточной Аляски. Таким образом, на основе выраженной аллельной изменчивости отдельных аллозимных локусов, имеет место значительная дифференциация региональных популяционных групп чавычи.
Ранее, в связи с наблюдаемыми закономерностями распределения аллельных частот на ареалах тихоокеанских лососей, Н.В. Варнавской все гены были условно разделены на четыре типа: I - гены, высокополиморфные на всем ареале, II - гены, мономорфные или почти мономорфные в одних регионах и высокополиморфные в других, III - гены,
Таблица 4. Подразделение генетической изменчивости в совокупности выборок из популяций чавычи Азии и Северной Америки (100 выборок, 53 аллозимных локуса; по нашим данным и по Gall et al., 1992; Crane et al., 1996)_
Уровень иерархии Абсолютное разнообразие Относительное разнообразие, %
Общее, Н, 0.0829 100.00
Внутри выборок, 0.0718 86.61
Между всеми выборками 13.39
Между регионами 9.05
Между популяциями в регионах 4.34
межпопуляционной дифференциации в пределах данных регионов.
мономорфные на всем ареале и низкополиморфные в отдельных регионах, IV - гены, почти мономорфные на всем ареале с появляющимися иногда в разных регионах редкими аллелями (Варнавская, 2005, 2006). У чавычи, на основании проанализированных в настоящем исследовании данных, к первому типу можно отнести только один локус - MPI*. Остальные являются мономорфными хотя бы в одном из исследованных регионов.
Ко второму, наиболее дискриминирующему региональные группы популяций типу, были отнесены 17 локусов (табл. 5), причем уровень изменчивости по многим из них очень велик: на протяжении ареала происходит замена доминирующего аллеля на альтернативный (частота изменяется от 1.0 почти до 0.0). Такая картина характерна для локусов sMEPl*, PGK2*, IDDH1*, PEPA*. Локус ADA1*, являясь мономорфным в популяциях чавычи шт. Орегон и Калифорния, имеет высокий уровень полиморфизма на Аляске и Камчатке. Сходная картина наблюдается по локусам sAAT3*, TPI4*, IDDH1 *, mMDH2*. Локусы PEPA *, GPIB2*, sSODl*, sMEPl*, PGK2*, напротив, высокополиморфны в южных районах американской части ареала и мономорфны, или почти мономорфны, в Азии и значительном числе локальностей Аляски.
К третьему типу у чавычи были отнесены еще 11 локусов, для которых характерен четко ограниченный региональными рамками полиморфизм, хотя и с невысокими (0.00-0.15) частотами редких аллелей. Особенным своеобразием отличаются южные локальные стада Североамериканского континента и относительно небольшая часть популяций Аляски.
Степень изменчивости полиморфных аллозимных генов, условно отнесенных к разным типам в зависимости от распределения их аллельных частот на ареале, отражена в величинах полокусной внутри- и межпопуляционной изменчивости (Hs и Fsl , соответственно) (табл. 5). Как видно из представленных данных, уровень дифференциации наиболее значителен по локусам, условно отнесенным ко второй группе генов, для которых характерен наиболее широкий диапазон варьирования частоты основного аллеля и достаточно четко выраженная межрегиональная изменчивость.
Таким образом, в представленной работе у чавычи было выявлено около 30 высокоинформативных, с точки зрения анализа популяционной структуры, аллозимных локусов. Можно подчеркнуть, что наиболее дискриминирующими внутри азиатской части ареала являются 10 локусов - sAHl* sIDHP2*, mMDH2* MPI*, TPI4* sMEPl*, sMEP2*, ALAT* mSODl* и GPIA* (по трем последним локусам нет данных по аллельной изменчивости в локальностях Орегона и Калифорнии); Азию и Северную Америку хорошо дифференцируют 23 локуса-mAATl*, тААТ2*, sAATl,2* sAAT3*, ADAl* sAHI*, GPIB2* IDDH1 *, sIDHPl* siDHP2*, mMDH2* sMDHB¡,2*sMEPl*, sMEP2*, MPI*, PEPLT*. PEPA*, PEPB1 *, PEPD2*, PGK2*, PGM2* sSODl*, TPI4*. Следует отметить сравнительно высокий
Таблица 5. Оценка средней популяционной гетерозиготности (Щ и уровня межпопуляционной
дифференциации (%) по 29
полиморфным локусам чавычи Азии и
уровень полиморфизма по аллозимным локусам у чавычи относительно других исследованных ранее видов тихоокеанских лососей (КоЬауаэЫ е1 а1., 1984; БЬак1ее е! а1.,
Тип Локус Hs 1990Ь,с; Варнавская, 2005).
изменчивости
I MPI* 0.247 11.0 Клинальная изменчивость частот
II sAAT3* 0.131 10.5 аллозимных локусов. Ранее при исследовании
II sAAT4* 0.118 14.9
II ADA1* 0.122 11.0 биохимического полиморфизма у разных видов
II sAHl* 0.114 6.3 рыб было отмечено, что частоты аллелей ряда
II GPIB2* 0.146 33.0
II IDDH1 * 0.141 24.2 аллозимных локусов закономерно меняются по
II sIDHPl * 0.059 12.2 мере продвижения по территории, занимаемой
II sIDHP2* 0.090 7.7 видом (Кирпичников, 1987). Данное явление,
II sMEPJ* 0.239 30.9
II sMEP2* 0.129 15.0 при котором происходит постепенное
II mMDH2* 0.091 35.8
II PEPA* 0.109 12.4 изменение генных частот или других признаков
II PEPB1* 0.084 6.4 в широтном (или в каком-либо другом)
II PGK2* 0.206 27.5 направлении, принято называть клинальной
II sSODl* 0.186 11.2
II TP13* 0.066 9.9 изменчивостью. Обнаружение клинальной
II TP 14* 0.092 9.6 изменчивости, оценка направления и
III sAATl,2* 0.013 5.4
III mAATl* 0.012 3.0 характерных свойств географических клин
III mAAT2* 0.023 4.8 может способствовать выявлению
III ESTD* 0.012 1.2
III FBALD3* 0.012 1.2 функциональных особенностей аллозимов, а
III FGDH* 0.010 6.2 также дать информацию о характере расселения
III LDHC* 0.002 1.2
III sMDHBl,2* 0.003 7.1 исследуемых видов и о механизмах
III PEPLT* 0.021 4.7 формирования их современной популяционной
III PEPD2* 0.007 3.1 структуры.
III PGM2* 0.015 6.4
Частоту основного аллеля аллозимных локусов, исследованных в каждой популяции, сопоставляли с географической широтой устья соответствующего речного бассейна, после чего вычисляли коэффициенты корреляции и оценивали их достоверность.
Результаты анализа показали, что частота основного аллеля большинства высокополиморфных локусов, условно отнесенных в зависимости от типа географической изменчивости аллельных частот к II группе, достоверно коррелирует с географической широтой устья нерестовых рек. Исключением являются локусы лАН1*, $МЕР2*, РЕРВ1*. Для III группы аллозимов оказались характерны относительно более низкие значения коэффициентов корреляции, хотя достоверность связи была также подтверждена в большинстве случаев (кроме локусов тААТI *, тААТ2*, РСИН*, РЕРЭ2*). Из числа локусов,
не исследованных в Азии, было выявлено еще шесть, аллельные частоты которых достоверно коррелировали с широтой на Северо-Американском континенте.
В целом у чавычи обнаружено 28 локусов, географическая изменчивость которых имеет характер клинальной, что, вероятно, отражает особенности миграционных потоков в период колонизации современных речных систем в период отступления ледников после Плейстоценового оледенения. При исследовании географических клин по биохимическим локусам у других видов рыб были описаны различные механизмы их появления (Кирпичников, 1987), но главной причиной возникновения определен отбор (Аронштам и др., 1977), т.е. клина, в широком понимании этого явления, - это ответ вида на различия в среде обитания в разных участках занятого им ареала.
Внутри- и межпопуляционная генетическая гетерогенность чавычи Азии. Сезонная изменчивость. Во многих нерестовых водоемах у чавычи была выявлена временная гетерогенность нереста, т.е. весенняя, летняя и осенняя расы (Вронский, 1983, 1984; Горшков и др., 1985; Burger et al., 1985). Известно, что для анадромных лососей характерна сильная генетическая детерминированность времени возврата на нерест, маркируемая аллозимными различиями между сезонными расами (Алтухов, 1997). Для чавычи Америки установлено, что темпоральная (или межрасовая) генетическая дифференциация невелика - она составляет 7.3% всей межпопуляционной компоненты генетического разнообразия (Utter et al., 1989; Алтухов, 1997).
Внутрипопуляциоиная генетическая гетерогенность в настоящей работе оценивалась на примере чавычи р. Камчатка. Уровень гетерогенности оценивали для ранней и поздней рас чавычи 1998 г. (что соответствует по общепринятой градации весенней и летней расам), а также для объединенных выборок 1995 и 1998 гг. (ранняя раса) и 1997, 1998 гг. (поздняя раса). Выявлено только два локуса (sSODl* и sMEP2*), гетерогенность по которым значима для исследованных выборок (р<0.05). Суммарные значения G-теста и теста х2 оказались недостоверны. В целом, полученные результаты свидетельствуют об отсутствии достоверных генетических различий между темпорально изолированными группировками чавычи данной локальности.
Межгодовая изменчивость. Межгодовая изменчивость изучалась на примере чавычи двух популяций - р. Большая и р. Камчатка. Величины теста %2 и G-теста, оценивающих межвыборочную гетерогенность генотипических распределений, показали, что для р. Большая статистически значимый результат отмечен только по локусу TPI4*. Оценки суммарной генетической гетерогенности между выборками 1996 и 1998 гг. (ранняя раса) показали ее недостоверность.
Анализ межгодовой изменчивости чавычи р. Камчатка показал однородность выборок 1997 и 1998 гг. (поздняя раса) и наличие значимой гетерогенности по локусам slDHP2*, mSODl*, sMEPl* и PEPA* для выборок 1995 и 1998 гг. (ранняя раса). Для второго случая суммарный результат теста у^ был недостоверен, а результат более чувствительного G-теста был значимым - 38.08 (d.f. = 17,р<0.01).
Кроме того, оценка генетической гетерогенности для всех четырех выборок из р. Камчатка позволила выявить значимые различия только по трем локусам - SIDHP2*, mSODl * и sMEPl *. Суммарный результат показал отсутствие достоверных внутрипопуляционных различий при использовании теста - 58.73 (d.f. = 43, р>0.05) и статистически достоверные различия при использовании G-теста - 65.41 (d.f. = 43, р<0.05), причем как в первом, так и во втором случаях результаты находились на границе уровня значимости.
Таким образом, полученные нами оценки генетической гетерогенности популяции чавычи р. Камчатка не позволяют сделать окончательный вывод о наличии межгодовой изменчивости, оставляя этот вопрос открытым и требующим дальнейшего изучения. Можно предположить, что исследования 1995 и 1998 гг. были проведены на выборках из разных субпопуляций чавычи, мигрирующих на нерест в р. Камчатка, и показатели межгодовой изменчивости в данном случае могут отражать пространственную неоднородность данной речной системы.
Межпопуляционная гетерогенность. Межпопуляционная гетерогенность по 16 полиморфным генам оценивалась для чавычи всех исследуемых азиатских популяций и на примере только двух самых крупных стад - pp. Большая и Камчатка. И в первом, и во втором случаях в анализ были включены суммарные выборки из каждой локальной популяции. Оценки генетической гетерогенности выборок, выполненные с помощью критерия %2и G-теста, оказались сходными.
Гетерогенность всех исследованных камчатских выборок была достоверной по частотам 10 генов - ALAT*, sAHl*, slDHP2* mMDH2*, MPI*, mSODl*, TPI4*, sMEPl* sMEP2*, GPIA*. Суммарные оценки гетерогенности по всем полиморфным локусам были высокодостоверными (р<0.001), что говорит о значительной репродуктивной изоляции исследованных популяций.
Сравнение двух самых крупных популяций чавычи Азии - рек Большая и Камчатка показало их достоверную генетическую неоднородность. Межпопуляциокные различия выявлены по 8-ми полиморфным локусам: ALAT* sAHl*, IDDH1*, s1DHP2*, MPI*, sMEPl*, sMEP2*, GPIA*.
Региональная изменчивость азиатской чавычи. Для оценки региональной изменчивости использовали два варианта анализа. В первом из них региональная генетическая гетерогенность оценивалась для четырех групп локальных популяций - выборки из рек Северо-Западной Камчатки (pp. Хайрюзова и Воровская), Западной Камчатки (pp. Утка, Быстрая и Большая), Юго-Восточной (р. Авача) и Восточной Камчатки (р. Камчатка). Во втором варианте исследования проводили, выделив всего два региона - Восточная и Западная Камчатка.
Для четырех региональных групп выявлены достоверные отличия по большинству исследованных локусов. Суммарные результаты тестов показали генетическую гетерогенность этих популяционных систем на уровне значимости /><0.001. Аналогичные результаты получены и для второго варианта анализа. Региональные популяционные. комплексы Западной и Восточной Камчатки достоверно отличались по 8 из 16, включенных в анализ локусов {ALAT*, sAHl*, IDDH1 *, sIDHP2* MPI*, sMEPl*, sMEP2* GPIA*). Суммарные отличия между выборками этих регионов были достоверны также на уровне pO.OOl. Таким образом, можно отметить, что уровень генетической дифференциации близко соответствует географической подразделенности региональных групп популяций. Оценка генетических расстояний между камчатскими популяциями чавычи. В результате проведенных исследований была создана база данных на основе аллозимной изменчивости, включающая частоты аллелей 10 выборок из 7 локальных популяций чавычи Камчатки. Вычисленные попарные стандартные генетические расстояния методом Нея (Nei, 1987) и генетические расстояния в соответствии с хордовым методом (Cavally-Sforsa, Edwards, 1967) далее были исследованы методами кластерного анализа и многофакторного шкалирования.
Анализ стандартных генетических расстояний Нея (рис. 3) показал, что весь массив данных разделился на хорошо обособленные кластеры: Восточная Камчатка, центральная Западная Камчатка (р. Утка и р. Воровская), и Западная Камчатка (р. Большая и р. Хайрюзова). Сходная дендрограмма была получена при вычислении хордовых расстояний (рис. 4), но в данном случае наблюдается более выраженная дифференциация выборок из популяций Западной и Восточной Камчатки.
Вычисления расположения точек в системе трехмерных координат по методу многомерного шкалирования позволили получить ещё более четкую картину генетических отношений между исследуемыми популяциями (рис. 5). Точки, соответствующие выборкам из локальных стад крупных речных бассейнов (pp. Камчатка и Большая), легли близко друг к другу, тогда как точки, соответствующие более северным рекам Западной Камчатки, образовали обособленную, но более разреженную область. Ближе всех к выборкам р.
Большая оказалась выборка из р. Утка, затем р. Воровская и р. Хайрюзова, что соответствует их географическому расположению по отношению друг к другу.
Западная Камчатка
Восточная Камчатка
р. Хяйрюаова |>. Быстрая р. Большая, 1996 р. Большая, 1998 р. Авача
р. Камчатка, 1997 р. Камчатка, 1995 р. Камчатка, 1998 р. Воровская р. Утка
Западная Камчатка
Восточная Камчатка
р. Хайрюзова р. Быстрая р. Большая, 1996 р. Большая, 1998 р. Воровская р. Утка р. Авача
р. Камчатка, 1997 р. Камчатка, 1998 р, Камчатка, 1У95
Рис. 3. UPGMA-дендрограмма, построенная на основе стандартных генетических расстояний (Nei, 1987), вычисленных по 20 полиморфным аллозимным локусам в популяциях чавычи Камчатки (sAATJ* mAATl*, AD Al*, ALAT* s AHI *, mMDH2* GP1A*, GR* IDDH1 *, sIDHP2*, LDHC* sMEPl*, sMEP2*, MPI*, PEPA*, PEPD2*, mSODl*, sSODl* ТРИ* ТРИ*).
Рис. 4. UPGMA-дендрограмма, построенная на основе хордовых генетических расстояний (Cavally-Sforsa, Edwards, 1967), вычисленных по 20 полиморфным аллозимным локусам в популяциях чавычи Камчатки (набор локусов тот же, что указан для рис. 3.).
[р^БодьшаяТшЙ
р. Большая, 1996) |р. Быстрая] '
»jSsf
;|р- ХаУфтоова}--,
Рис. 5. Генетическое сходство популяций чавычи, вычисленное методом главных компонент по 20 полиморфным аллозимным локусам (набор локусов тот же, что указан для рис. 3.).
Q Восточная Камчатка I I бас. р. Большая (Западная Камчатка) PI Западная Камчатка
Генетическая дифференциация чавычи Азии и Северной Америки. Для оценки уровня дивергенции популяций чавычи в различных регионах ареала вида вычисляли генетические расстояния хордовым методом Кавалли-Сфорца (Cavally-Sforsa, Edwards, 1967). Результаты кластеризации показали, что весь массив данных разделился на две дифференцированных ветви - южные и северные популяции. Южная ветвь включает популяции Калифорний и Орегона. Северная ветвь включает три обособленных кластера -Центральная и Западная Аляска, Юго-Восточная Аляска и Камчатка. Кластер камчатских
популяций разделился, в свою очередь, на два кластера - Восточная и Западная Камчатка. В целом, группы популяций, выделенные на основе генетических расстояний, соответствуют географическим региональным группам.
Анализ генетического сходства методом факторного шкалирования позволил выделить семь хорошо обособленных областей, которые полностью соответствуют региональным группам популяций - 1) бас. р. Сакраменто, 2) бас. рек Кламат и Тринити, 3) Калифорния и Орегон, 4) Юго-Восточная Аляска, 5) Центральная и Северо-Западная Аляска, .6) п-ов Камчатка (рис. 6).
I
Рис. 6. Генетическое сходство популяций чавычи Азии и Северной Америки, вычисленное методом главных компонент по аллельным частотам 23 локусов.
На основе проделанного анализа был сформирован набор наиболее информативных дискриминирующих генов и созданы базы реперных данных, которые были использованы для идентификации смешанных уловов чавычи в Тихом океане и сопредельных морях (Varnavskaya et al., 2003; Varnavskaya, Shpigalskaya, 2004).
Глава 4. Изменчивость митохондриальной ДНК в популяциях чавычи Камчатки
Длина амплифицированных в результате полимеразной цепной реакции (ПЦР) фрагментов мтДНК в сумме составляет почти 7000 пар нуклеотидов (пн), или около 40% всего митохондриального генома чавычи (Wilhelm et al., 2003). В результате ПДРФ-анализа
было выявлено 53 рестриктных сайта (=1.5% митохондриалького генома), из которых только 6 оказались полиморфными.
Объединённые данные, полученные в результате анализа трёх фрагментов мтДНК по всем сайтам рестрикции и использованным эндонуклеазам, позволили выявить 11 комбинированных гаплотипов у 579 особей чавычи Камчатки, обозначенные как 01-011 (табл. 6). Как явно доминирующий может быть охарактеризован гаплотип 01, средняя частота других невелика и только в одном случае (02) достигала значения 0.03.
Таблица 6. Выявленные комбинированные гаплотипы митохондриальной ДНК чавычи Азии
Комбинированный гаплотип A8/C0III ND3/ND4 COI/COII
BstN I Hinfl Mbo I Rsa I Dde I Taq I Hinfl Rsa I Sau96 1 Mspl
G1 В В в А А А А А А А
G2 А В в А А А А А В А
G3 В В в А В А А А А А
G4 А В в в А В А А А А
G5 В в в В А А А А А А
G6 А в в А А В А А А А
G7 В в в В А В А А А А
G8 В в в А А А А А В А
G9 А в в А А А А А А А
G10 В в в А А В А А А А
G11 В в А А А А А А А А
Сеть гаплотипов, построенная по принципу минимального числа нуклеотидных
различий между ними, приведена на рис. 7. Из рисунка видно, что распределение гаплотипов чавычи всех исследованных локальностей, в целом формирует картину «звёздчатой» радиации гаплотипов в результате, как правило, единичных нуклеотидных замен. Такое распределение свидетельствует об общности происхождения исследованных популяций и,
Рис. 7. Генеалогическая сеть гаплотипов чавычи Камчатки, построенная по принципу минимального числа нуклеотидных замен и показывающая мутационные различия между гаплотипами. Размер окружности приблизительно отражает частоту встречаемости данного гаплотипа во всех исследованных популяциях. Число поперечных штрихов на ветвях обозначает число нуклеотидных замен. Пунктиром обозначены вероятные альтернативные связи между гаплотипами.
Кроме того, «звездчатая» структура гаплотипической сети свидетельствует о недавнем происхождении популяций на исследованной части ареала (Rogers, Harpending,
соответственно, о едином предковом генофонде.
1992; Rogers, 1995). Наличие возможных альтернативных связей между гаплотипами (петлевые структуры в сетке) может быть результатом гомоплазии, т.е. повторных и обратных мутаций в мтДНК (Churikov et al., 2001) с достаточной вероятностью происходящих на значительной части видового ареала.
Внутри- и межпопуляционная изменчивость частот гаплотипов мтДНК чавычи.
Частота встречаемости комбинированных гаплотипов, выявленных в популяциях чавычи, приведена в табл.7.
Таблица 7. Встречаемость комбинированных гаплотипов митохондриальной ДНК чавычи
Локальность N G1 G2 G3 G4 G5 G6 G7 G8 G9 G10 G11
р. Палана 50 42 8 0 0 0 0 0 0 0 0 0
р. Тигиль 42 38 0 3 1 0 0 0 0 0 0 0
р. Пымта 50 37 0 4 1 4 4 0 0 0 0 0
р. Утка 50 41 2 2 5 0 0 0 0 0 0 0
Олюторский залив 50 42 0 0 2 0 0 3 0 3 0 0
р. Коль 50 49 0 0 0 0 0 0 1 0 0 0
р. Воровская 50 48 0 0 0 0 1 0 0 1 0 0
р. Кихчик 33 32 0 1 0 0 0 0 0 0 0 0
р. Авача 44 43 0 0 0 0 0 0 0 1 0 0
р. Жупанова 30 29 0 0 0 0 0 0 1 0 0 0
р. Камчатка 60 52 5 0 0 0 0 0 0 0 3 0
р. Большая 60 49 2 4 2 0 2 1 0 0 0 0
р. Пахача 10 7 0 0 1 0 0 0 0 1 0 1
Всего 579 509 17 14 12 4 7 4 2 6 3 1
Доминирующий гаплотип G1 представлен во всех азиатских популяциях. В 13
исследованных выборках его частота варьировала от 0.70 (р. Пахача) до 0.98 (р. Коль), средняя частота составила 0.88. В шести камчатских популяциях отмечен гаплотип G4. Шесть гаплотипов (G2, G3, G6, G7, G8, G9) с низкой частотой присутствовали в меньшем количестве локальных популяций. Гаплотипы G5, G10, G11 обнаружены только в одной из 13 выборок, причем последний был выявлен у одной из исследованных особей. Вероятно, об указанных гаплотипах можно говорить как об уникальных, хотя в данном случае необходимо учитывать относительно небольшие объемы исследованных выборок камчатской чавычи. Редкие гаплотипы, являясь дискриминирующими, могут служить молекулярными маркерами для особей соответствующих локальностей, что определяет необходимость продолжения и расширения поисковых исследований в данной области.
Значения гаплотипического и нуклеотидного разнообразия в популяциях камчатской чавычи приведены в табл. 8. На азиатской части ареала наиболее высокие значения данных показателей характерны для чавычи самой северной восточнокамчатской популяции - р. Пахача, несколько ниже эти значения у особей в выборках из рек Пымта, Большая и Утка. Как было отмечено ранее, в выборках из первых двух популяций присутствуют редкие
гаплотипы. Для всех остальных локальностей свойственны заметно более низкие величины гаплотипической и нуклеотидной изменчивости.
Таблица 8. Гаплотипическое и нукяеотидное Анализ количества
нуклеотидных замен во всех гаплотипических вариантах при попарном сравнении популяций чавычи Камчатки позволил выявить максимальное число нуклеотидных замен - девять - в выборке из р. Большая, минимальное - одна нуклеотидная замена между двумя выявленными
гаплотипами - у чавычи рек
Ранее было показано, что в изолированных популяциях или группах локальных популяций конфигурация распределения числа нуклеотидных замен определяется демографической историей (Rogers, Harpending, 1992; Rogers, 1995). Размер выборок, использованных в нашем анализе, недостаточен для подобных заключений, но позволяет сделать первичную оценку эволюционного разнообразия. В соответствии с вышеизложенным, количество нуклеотидных замен в митохондриальных гаплотипах камчатских популяций чавычи варьирует от преимущественно одной до девяти, что позволяет говорить об очень небольшом эволюционном разнообразии азиатских популяций и о небольшой глубине дивергенции между ними.
Оценка географической гетерогенности частот гаплотипов. Частоты и варианты гаплотипов в исследованных выборках заметно различались. Среди особей из р. Большая удалось выявить наибольшее количество гаплотипов - шесть, в остальных популяциях азиатской части ареала их количество варьировало от двух (pp. Авача, Жупанова, Палана, Коль, Кихчик) до пяти (р. Пымта).
Анализ частот гаплотипов с использованием псевдовероятностного теста на гомогенность позволил оценить уровень неоднородности чавычи исследованных популяций. При сравнении всех 13 выборок п-ова Камчатка тест показал значимые различия (х2=314.62; d.f.~ 120; pO.OOl). При раздельном анализе гетерогенности региональных групп популяций Западной и Восточной Камчатки значения ^2-теста также оказались статистически
разнообразие в популяциях чавычи
Популяция Гаплотипическое разнообразие Нуклеотидное разнообразие
р. Пахача 0.5333 ±0.18007 0.002454
Олюторский залив 0.2914 ±0.08160 0.001430
р. Палана 0.2743 ± 0.07093 0.002867
р. Тигиль 0.1800 ±0.07676 0.000712
р. Воровская 0.0792 ±0.05213 О.ОООЗОЗ
р. Коль 0.0400 ±0.03803 0.000310
р. Пымта 0.4416 ±0.08286 0.001796
р. Кихчик 0.0606 ±0.05612 0.000152
р. Утка 0.3208 ±0.08147 0.002400
р. Большая 0.3305 ± 0.07737 0.001958
р. Авача 0.0455 ± 0.04292 0.000117
р. Жупанова 0.0667 ±0.06125 0.000517
р. Камчатка 0.2435 ±0.07001 0.001873
Кихчик и Авача.
значимыми. Для восьми западнокамчатских популяций %2=125.88 (с/./.=56, р<0.001), для пяти локальностей Восточной Камчатки - х2=66.91 (с//=28, р<0.001). Вероятность гомогенности сравниваемых выборок при всех вариантах анализа оказалась равна нулю.
В результате применения методов ^статистики для оценки уровня генетической дифференциации азиатских популяций чавычи (оценка проводилась на основе изменчивости числа нуклеотидных замен) было выявлено, что между большинством сравниваемых пар значимые различия отсутствуют. На азиатской части ареала пять выборок (рр. Палана, Пымта, Коль, Воровская, Авача) при попарном сравнении значимо отличаются более чем от половины включенных в анализ локальностей.
Таким образом, полученные данные свидетельствуют о том, что, несмотря на достоверный уровень гетерогенности чавычи на изучаемой части ареала по структуре мтДНК, межпопуляционная дивергенция в целом невелика даже с учётом того, что в отдельных случаях величина выраженности отличий между популяциями достигает значимого уровня.
Популяционная дивергенция чавычи на основе анализа изменчивости мтДНК.
Расчет матрицы дивергенции по нуклеотидным последовательностям между популяциями чавычи Камчатки позволил выявить, что максимальное значение различий (около 0.02% нуклеотидных замен) имеет место между выборкой из р. Палана и остальными. Фенограмма
Рис. 8. иРОМА-фенограмма, иллюстрирующая генетические различия между локальными популяциями чавычи Камчатки.
Камчатские популяции, в соответствии с уровнем нуклеотидной дивергенции, сформировали три кластера, каждый из которых включает выборки, имеющие относительно близкие значения данного показателя. В наиболее дифференцированный кластер вошли популяции рек Палана и Камчатка. Кластер с несколько меньшими генетическими
на основе матрицы представлена на рис. 8.
- р. Палана
- р. Камчатка _ р. Тигиль
_ р. Кихчнк
_ р. Коль
_ р. Жупанова
_ р. Воровская
_ р. Авача
р. Пымта
Олюторский I залив
р. Пахача
- р. Утка
- р. Большая
0.00000 .0.00001 нуклеотидная дивергенция
расстояниями между образовавшими его выборками, включает популяция рек Тигиль, Кихчик, Коль, Жупанова, Воровская и Авача. В третьем кластере объединены популяции, нуклеотидная дивергенция между которыми оказалась минимальной - pp. Большая, Утка, Пымта, а также выборки из Олюторского залива. Необходимо отметить, что кластеризация популяций на дендрограмме не показывает четкой картины дифференциации популяций чавычи Западной и Восточной Камчатки.
Полученные данные были проанализированы также с использованием программы AMOVA (Analysis of Molecular Variance) в пакете программ Arlequin2000. В первом варианте анализировали всю совокупность выборок. Расчеты показали, что на долю внутрипопуляционной генетической дисперсии приходится 95.98%, на долю межпопуляционной - 4.02%. В другом варианте расчета популяции объединяли в две группы в соответствие с географической принадлежностью - Западная Камчатка и Восточная Камчатка. При таком объединении выявляется близкая к нулю величина изменчивости, приходящаяся на межгрупповую компоненту. Доля межпопуляционной дисперсии внутри групп составляет 3.95%, а доля внутрипопуляционной изменчивости (96.01%) незначительно увеличивается по сравнению с предыдущим вариантом анализа.
Вероятной причиной неглубокой дивергенции камчатских популяций чавычи по структуре мтДНК может являться их исторически недавнее происхождение в сочетании с относительно более низкой, чем у других тихоокеанских лососей, скоростью смены генераций, обратная зависимость между которой и уровнем изменчивости мтДНК других видов лососей была показана ранее (Брыков, 2001). Недавнее происхождение большинства популяций чавычи в реках Камчатки подтверждается также «звездчатой» структурой сети мтДНК. Можно предположить, что относительно небольшая глубина межпопуляционной дивергенции свидетельствует о независимой эволюции камчатских популяций в течение исторически небольшого временного периода, что подтверждается также низкой (менее 5%) долей генетической изменчивости, приходящейся на межпопуляционную компоненту. Если принять во внимание существующие оценки скорости накопления нуклеотидных замен в мтДНК наземных позвоночных - 2% замен нуклеотидов в течение 1 млн. лет (Avis, 1994; Avis et al. 1998), и допустить, что у тихоокеанских лососей приблизительная калибровка «молекулярных часов» соответствует еще меньшей скорости дивергенции пары нуклеотидных последовательностей, равной 1% за 1 млн. лет (Smith, 1992), то с учетом того, что выявленная величина различий между популяциями камчатской чавычи составляет около 0.02%, максимальное время дивергенции соответствует оценке во временном интервале, не превышающем 10-20 тыс. лет. Вероятно, отступление позднеплейстоценового ледника (от 10 до 18 тыс. лет назад) (Брайцева и др., 1968) обусловило начало процессов
расселения, реколонизации и формирования современной картины подразделенное™ вида на структурные компоненты - популяции.
Заключение
В результате проведенных исследований на основе двух генетических методов -анализа изменчивости аллельных частот аллозимных локусов и анализа изменчивости структуры митохондриальной ДНК, описана генетическая дифференциация азиатских популяций чавычи, и проведена оценка популяционной дивергенции на всем ареале данного вида.
Сравнительный анализ распределений генных частот на протяжении всего ареала, включающего его азиатскую и американскую части, показал, что по многим аллозимным локусам наблюдается высокий уровень изменчивости региональных комплексов популяций, при этом часто происходит замена преобладающего аллеля. В исследованных популяциях обнаружено 28 лок)гсов, географическая изменчивость которых имеет характер клинальной, что, вероятно, отражает особенности миграционных потоков в период колонизации современных речных систем в период отступления ледников после Плейстоценового оледенения.
Выявлена достоверная генетическая межпопуляционная гетерогенность азиатской чавычи. Региональная изменчивость генетической структуры популяций четко выражена и позволяет выделить два региональных популяционных комплекса - Западная и Восточная Камчатка.
Кластерный и многофакторный анализы генетических расстояний, определенны* на основе изменчивости аллельных частот полиморфных аллозимных локусов, показали, что генетическая дифференциация соответствует пространственной подразделенности видового ареала. В структуре иРйМА-дендрограмм, также как и при многомерном шкалировании, кластеризация по генетической близости согласуется с географической приуроченностью популяций, что свидетельствует о генетическом своеобразии региональных популяционных комплексов. Создана реперная база данных, включающая мультилокусные характеристики по аллозимным генам большинства крупных популяционных систем ареала вида, на основе которых может осуществляться идентификация региональной принадлежности смешанных морских уловов в период нагульных и преднерестовых миграций чавычи.
Получены результаты исследования полиморфизма мтДНК у чавычи азиатских популяций. ПДРФ-анализ с применением восьми эндонуклеаз рестрикции позволил выявить изменчивость нуклеотидных последовательностей в трех амплифицированных участках мтДНК по 6 полиморфным рестриктным сайтам. Анализ распределения частот комбинированных гаплотипов чавычи Камчатки и изменчивости числа нуклеотидных замен
позволил оценить внутри- и межполуляционные генетические различия. При сравнении выборок из локальных популяций чавычи Азии статистически подтверждена их гетерогенность по структуре митохондриального генома, что свидетельствует о внутривидовой генетической подразделенности. При попарном сравнении выборок, между большинством сравниваемых пар обнаружены значимые различия. Полученные результаты свидетельствуют о возможности использования митохондриального генома в различных аспектах популяционных исследований в качестве информативного генетического маркера.
Выявленные отличия в числе нуклеотидных замен в мтДНК камчатских популяций свидетельствуют о том, что азиатской чавыче свойственно небольшое эволюционное разнообразие, что подтверждается относительно малыми межпопуляционными генетическими расстояниями и, соответственно, невысоким уровнем генетической дивергенции. Вероятной причиной относительно неглубокой дивергенции камчатских популяций по мтДНК может являться их исторически недавнее происхождение из предковой популяции с единым генным пулом. В этом случае небольшие межпопуляционные различия могут быть обусловлены небольшим временем самостоятельной эволюции.
Высказано предположение, что максимальное время дивергенции для камчатских популяций чавычи, величина различий между которыми составляет около 0.02% нуклеотидных замен, соответствует оценке во временном интервале, не превышающем 10-20 тыс. лет. Вероятно, отступление позднеплейстоценового ледника обусловило начало процессов расселения, реколонизации и формирования современной картины подразделенности вида на структурные компоненты. Представленные ранее палеогеографические данные показывают, что описанные особенности межпопуляционной генетической дифференциации согласуются с историей позднеплейстоценовых оледенений, которые были ясно выражены на Камчатке и побережье Охотского моря (Баранова, Бискэ, 1964; Мелекесцев и др., 1974; Глушкова, 1984; Беспалый, Глушкова, 1987). После отступления ледника, сохранившиеся в различных рефугиумах популяции, занимали экологически соответствующие ниши (Алтухов, 1997). Свидетельствами определенных сценариев расселения и процессов формирования современной внутривидовой структуры являются, в первую очередь, генетические особенности популяций.
Выводы
1. Выявлена статистически значимая межпопуляционная генетическая гетерогенность чавычи Азии по аллельной изменчивости исследованных аллозимных локусов и гаплотипам мтДНК.
2. Для региональных групп популяций Восточной Камчатки, бассейна р. Большая, центральной Западной и Северо-Западной Камчатки характерно генетическое
своеобразие, позволяющее по изменчивости частот аллозимных локусов идентифицировать данные регионы в смешанных морских скоплениях.
3. Создана основа для формирования базы данных по частотам аллелей аллозимных локусов чавычи. Выявлены наиболее информативные, с точки зрения анализа популяционной изменчивости, аллозимные локусы. Наиболее дискриминирующими в пределах азиатской части ареала являются 10 локусов - sAHl*, ALAT*, GPIA*, sIDHP2*, mMDH2*, sMEPl*, sMEP2*, MPI*, mSODl*, TPI4*; Азию и Северную Америку хорошо дифференцируют 23 локуса - mAATI* тААТ2*, sAATl,2*, sAATS*, ADA!*, sAHl*, GPIB2*, IDDH1* sIDHPl*, sIDHP2*, mMDH2*, sMDHBl.2* sMEPl*, sMEP2* MPI* PEPLT*, PEPA*, PEPB1 *, PEPD2* PGK2* PGM2*, sSODi* TPI4*.
4. Географическая изменчивость 28 аллозимных локусов чавычи имеет характер клинальной, что, вероятно, отражает особенности миграционных потоков в период колонизации современных речных систем.
5. Наиболее высокие оценки межпопуляционного разнообразия по аллозимньш локусам свойственны североамериканским популяциям чавычи.
6. Невысокий уровень изменчивости митохондриальных гаплотипов, выявленных в камчатских популяциях, свидетельствует о небольшом эволюционном разнообразии азиатской чавычи.
7. Несмотря на достоверный уровень гетерогенности азиатской чавычи по структуре мтДНК, межпопуляционная дивергенция на этой части ареала в целом невелика.
8. Вероятной причиной относительно неглубокой дивергенции камчатских популяций чавычи по структуре мтДНК может являться их исторически недавнее происхождение. Максимальное время самостоятельной эволюции азиатской чавычи соответствует оценке во временном интервале, не превышающем 10-20 тыс. лет.
СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ
1. Варнавская Н.В., Шпнгальская Н.Ю. Генетическое разнообразие азиатских популяционных комплексов тихоокеанского лосося - чавычи, Oncorhynchus tschawytscha (Walbaum) / Тез. докл. Межд. конф. «Новые технологии в защите биоразнообразия в водных экосистемах», Москва, 27-29 мая 2002 г. - М.: МАКС Пресс, 2002. С. 87.
2. Варнавская Н.В., Шпнгальская Н.Ю., Савин В.А. Оценка дивергенции искусственно-природных и естественных популяций чавычи и кеты Камчатки / Тез. докл. Межд. конф. «Новые технологии в защите биоразнообразия в водных экосистемах», Москва, 27-29 мая 2002 г. - М.: МАКС Пресс, 2002. С. 86.
3. Varnavskaya N., Templin W., Shpigalskaya N. Genetic stock identification of chinook salmon, Oncorhynchus tschawytscha (Walbaum) / Abstracts of NPAFC Int. Workshope «Application of Stock Identification in Defining Marine Distribution and Migration of Salmon» (Honolulu, Hawaii, USA, November 1-2,2003). P. 53.
4. Varnavskaya N., Shpigalskaya N. Genetic stock identification of chinook salmon, Oncorhynchus tschawytscha (Walbaum) / Thesis of NPAFC Int. Workshop «Application of Stock Identification in Defining Marine Distribution and Migration of Salmon» (Honolulu, Hawaii, USA, November 1-2,2003). -NPAFC Tech. Rep. 2004. №. 5. P. 76-78.
5. Brykov V.V., Shpigalskaya N.Yu., Kukhlevsky A.D., Varnavskaya N.V. Mitochondrial DNA variation in populations of Asian chinook salmon / Proc. NPAFC Intern. Workshop «BASIS-2004: Salmon and Marine Ecosystems in the Bering Sea and Adjacent Waters». Sapporo, October 30-31,2004. P. 71.
6. Варнавская H.B., Шпнгальская Н.Ю., Шапорев P.A. Изменчивость частот ферментных генов в популяциях чавычи, Oncorhynchus tschawytscha (Walbaum), Азии / Популяционная биология, генетика и систематика гидробионтов. Сб. научных трудов КамчатНИРО. Т. 1. - Петропавловск-Камчатский: КамчатНИРО, 2004. С. 84-96.
7. Шпнгальская Н.Ю. Межпопуляционные различия по структуре митохондриальной ДНК у чавычи // Материалы отчетной сессии КамчатНИРО по итогам НИР в 2004 г. -Петропавловск-Камчатский: КамчатНИРО, 2005. С. 83-92.
8. Шпнгальская Н.Ю., Варнавская Н.В., Брыков В.А. Генетическая изменчивость чавычи, Oncorhynchus tschawytscha (Walbaum), Азии на основе анализа частот аллозимных генов и структуры митохондриальной ДНК II Ихтиологические исследования на внутренних водоемах: Материалы Межд. науч. конф. / А.Б. Ручин, отв. ред. - Саранск: Мордов. гос. ун-т, 2007. С. 192-194.
9. Шпнгальская Н.Ю., Варнавская Н.В. Генетическая дифференциация популяций чавычи Oncorhynchus tschawytscha (Walbaum) // Современное состояние и пути развития популяционной биологии: Материалы X Всероссийского популяционного семинара. Ижевск, 17-22 ноября 2008 г. Ижевск: КиигоГрад, 2008. С. 382-384.
10. Шпигальская Н.Ю., Брыков Вл.А., Гарретт А.Дж., Кухлевский А.Д., Шапорев P.A., Варнавская Н.В. Межпопуляционная изменчивость митохондриальной ДНК чавычи, Oncorhynchus tschawytscha (Walbaum), Камчатки // Генетика. 2008. Т. 44, № 7. С. 972-982.
Нина Юрьевна Шпигальская
Генетическая дифференциация азиатских популяций тихоокеанского лосося - чавычи, ОпсогсЬупсЬм исАошуИсЛа ^а1Ьаит)
03.02.07 - генетика
АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата биологических наук
Подписано к печати 21 апреля 2010 г. Заказ от 22 апреля 2010 г. Издательство ФГУП «КамчатНИРО» 683602, г. Петропавловск-Камчатский, ул. Набережная, 18 Объем 24 стр. А5. Тираж 100 экз.
Содержание диссертации, кандидата биологических наук, Шпигальская, Нина Юрьевна
ВВЕДЕНИЕ.
ГЛАВА 1. ОБЗОР ЛИТЕРАТУРЫ.
1.1. Ареал, краткая биологическая характеристика и внутривидовая организация чавычи.
1.2. Исследования популяционной организации тихоокеанских лососей на основе молекулярно-генетических методов и состояние изученности данного вопроса у чавычи.
1.2.1. Биохимический полиморфизм.
1.2.2. Полиморфизм ДНК.
ГЛАВА 2. ХАРАКТЕРИСТИКА МАТЕРИАЛА И
МЕТОДОВ ИССЛЕДОВАНИЯ.
2.1. Материалы.
2.2. Методы.
2.2.1. Методы выявления биохимической изменчивости.
2.2.2. Методы выявления изменчивости мтДНК.
2.2.3. Методы статистического анализа.
ГЛАВА 3. БИОХИМИЧЕСКИЙ ПОЛИМОРФИЗМ В ПОПУЛЯЦИЯХ ЧАВЫЧИ.
3.1. Уровень изменчивости аллозимных популяционных маркеров.
3.2. Пространственная изменчивость аллельных частот.
3.3. Клинальная изменчивость частот аллозимных локусов.
3.4. Внутри- и межпопуляционная генетическая гетерогенность чавычи Азии.
3.4.1. Сезонная изменчивость.
3.4.2. Межгодовая изменчивость.
3.4.3. Межпопуляционная гетерогенность.
3.4.4. Региональная изменчивость азиатской чавычи.
3.4.5. Оценка генетических расстояний между камчатскими популяциями чавычи.
3.5. Генетическая дифференциация чавычи Азии и Северной Америки.
3.6. Создание базы реперных данных для популяционной идентификации азиатской чавычи.
ГЛАВА 4. ИЗМЕНЧИВОСТЬ МИТОХОНДРИАЛЬНОЙ ДНК В ПОПУЛЯЦИЯХ ЧАВЫЧИ КАМЧАТКИ.
4.1. Результаты ПЦР-ПДРФ-анализа мтДНК.
4.2. Внутри- и межпопуляционная изменчивость частот гаплотипов мтДНК чавычи.
4.3. Оценка географической гетерогенности частот гаплотипов.
4.4. Популяционная дивергенция чавычи на основе анализа изменчивости мтДНК.
Введение Диссертация по биологии, на тему "Генетическая дифференциация азиатских популяций тихоокеанского лосося-чавычи, Oncorhynchus Tschawytscha (Walbaum)"
Актуальность темы. Тихоокеанские лососи рода Oncorhynchus, являясь ценными объектами промысла, имеют большое экономическое значение для стран Тихоокеанского бассейна. Для всех представителей рода характерна сложная популяционная организация, основными чертами которой являются соответствие принципу иерархии и упорядоченная во времени и пространстве субпопуляционная структура (Altukhov, 1971, 1981, 1990; Алтухов, 1973, 1974; Altukhov, Salmenkova, 1981, 1991; Алтухов и др., 1997; Altukhov et al., 2000). Эмпирическое подтверждение существования в пределах видов генетически отличающихся, репродуктивно изолированных единиц, дает основание для утверждения, что успех краткосрочной стратегии управления ресурсами и долговременных целей сохранения видов зависит от наличия, наряду с демографическими и экологическими данными, информации о генетической внутривидовой структуре (Алтухов и др., 1997; Алтухов и др., 2004). Для решения теоретических проблем популяционной организации в сочетании с долговременным мониторингом популяционных генофондов и прикладных задач генетической идентификации также необходимы данные о генетическом своеобразии внутривидовых структурных компонентов и оценки уровня их дифференциации.
Чавыча, О. tschawytscha (Walbaum), является одним из наиболее ценных видов тихоокеанских лососей, который в азиатской части ареала обитает преимущественно на Камчатке. Популяционно-генетические исследования чавычи Северной Америки, где данный вид распространен практически по всему тихоокеанскому побережью, были начаты во второй половине 20 в. и показали высокий уровень межпопуляционной изменчивости по аллозимным локусам (Kristiansson, Mclntyre, 1976; Gharrett et al., 1987; Reisenbichler, Phelps, 1987; Beacham et al., 1989, 1996; Utter et al., 1989; Bartley, Gall, 1990; Gall et al., 1992; Crane et al., 1996; Teel, 2000, 2004; Templin, Seeb, 2003), структуре митохондриальной ДНК (Wilson et al., 1987; Cronin et al., 1993; Adams et al.,
1994), микросателлитам (Moran et al., 2005; Beacham et al., 2006b; Seeb et al., 2007) и единичным нуклеотидным заменам (SNPs) (Smith et al., 2007; Narum et al., 2008). Были созданы объединенные реперные базы данных по частотам аллозимных локусов, на основе которых с успехом осуществлялась идентификация смешанных уловов в прибрежных водах США и Канады (Teel et al., 1999), и базы реперных данных по аллельной изменчивости микросателлитных локусов (Moran et al., 2005; Seeb et al., 2007). Частоты аллозимных генов служили также маркером для оценки успешности деятельности рыбоводных заводов (Waples, 1990a,b; Waples, Do, 1994).
В Азии исследования генетической структуры популяций чавычи на основе биохимической изменчивости и полиморфизма мтДНК находятся на начальном этапе, анализ микросателлитной ДНК и полиморфизма единичных нуклеотидных замен пока не начат.
До недавнего времени для выявления внутри- и межпопуляционной дифференциации и оценки биологического разнообразия азиатских популяций тихоокеанских лососей наиболее широко использовали такие маркеры генетической изменчивости, как полиморфные белковые локусы (Алтухов, 2003; Алтухов и др., 1997; Варнавская, 2001), но в последние десятилетия все большее значение приобретают исследования на уровне ДНК (Avise, 1994; Брыков, 2001; Чуриков, 2001; Алтухов, Салменкова, 2002). Внедрение молекулярно-генетических методов в популяционные исследования позволяет дополнить и расширить уже сформированную картину генетической подразделенности и процессов эволюции видов. Анализ изменчивости нуклеотидной последовательности может дать дополнительное количество высокоинформативных маркеров, что особенно актуально при исследовании видов со сложной генетической структурой и отсутствием четкой корреляции между генетическим и географическим расстояниями.
Исследования популяционной структуры тихоокеанских лососей с использованием молекулярно-генетических методов в силу их больших и принципиально новых возможностей получили в настоящее время широкое распространение. Анализ изменчивости неядерной наследственной структуры — митохондриальной ДНК является эффективным инструментом оценки генетической дивергенции и выявления филогенетической истории популяций (Полякова и др., 1992; Brykov et al., 1996; Брыков и др., 1999а,б; Брыков и др., 2002; Брыков и др., 2005). Исследования полиморфизма мтДНК чавычи могут дать дополнительную информацию о внутри- и межпопуляционной генетической структуре данного вида, позволят выявить филогенетические связи между внутривидовыми структурными единицами, а также увеличить число информативных генетических маркеров, пригодных для различных аспектов популяционных исследований.
Цель и задачи исследования. Целью настоящей работы является оценка уровня генетического разнообразия азиатских популяций чавычи на основе биохимического полиморфизма и гаплотипической изменчивости митохондриальной ДНК, а также создание баз генетических данных для последующей идентификации региональной принадлежности чавычи в смешанных морских скоплениях.
Для достижения поставленной цели были определены следующие задачи:
1. изучить биохимический полиморфизм и изменчивость последовательностей мтДНК в популяциях чавычи азиатской части ареала;
2. выявить максимально возможное число информативных генетических маркеров в популяциях чавычи;
3. выделить набор наиболее дискриминирующих азиатскую чавычу генетических маркеров для идентификации ее в смешанных уловах;
4. оценить внутри- и межпопуляционную изменчивость генетической структуры чавычи Азии;
5. исследовать генетическую дивергенцию азиатских и североамериканских популяционных комплексов.
Научная новизна. Впервые получены результаты исследования полиморфизма митохондриальной ДНК у чавычи Камчатки. Анализ гаплотипической изменчивости позволил получить оценки внутри- и межпопуляционной дифференциации, а также уровня дивергенции популяций азиатской части ареала. Выявлено, что максимальное время независимой эволюции камчатских популяций чавычи соответствует оценке во временном интервале, не превышающем 10-20 тыс. лет. Впервые получены сравнительные оценки генетического разнообразия по аллозимным локусам чавычи различных частей ареала. Характер пространственной изменчивости аллельных частот, уровень внутри- и межпопуляционной дифференциации по исследованным аллозимным локусам позволил выделить ряд наиболее информативных для идентификации азиатских популяций чавычи в смешанных морских скоплениях.
Практическая значимость. Результаты, представленные в работе, могут быть использованы для создания баз реперных данных по двум типам исследованных генетических маркеров в целях идентификации смешанных морских скоплений чавычи, а также генетического мониторинга популяционных процессов и сохранения естественного биоразнообразия природных популяций данного вида в условиях морского и берегового промысла. Помимо этого, полученные результаты на практике могут быть использованы при организации искусственного разведения чавычи для долговременного слежения за состоянием популяционных генофондов, оценки и прогнозирования их динамики во времени и пространстве.
Апробация работы. Основные положения и результаты работы были представлены на Международной конференции «Новые технологии в защите биоразнообразия в водных экосистемах» (Москва, 27-29 мая 2002 г.), Международной рабочей группе Северотихоокеанской комиссии по анадромным видам рыб NPAFC International Workshope on Application of Stock Identification in Defining Marine Distribution and Migration of Salmon (Honolulu, Hawaii, USA, November 1-2, 2003), объединенной отчетной сессии НТК «ТИНРО» (Петропавловск-Камчатский, 2006 г.), Международной конференции «Ихтиологические исследования на внутренних водоемах» (Саранск, 2007), X Всероссийском популяционном семинаре «Современное состояние и пути развития популяционной биологии» (Ижевск, 17-22 ноября 2008 г.), а также в виде научных годовых отчетов ФГУП «КамчатНИРО» в 2001 г. и 2003-2006 гг.
Публикации. По материалам диссертации опубликовано 10 работ, в том числе, 1 статья - в рецензируемом журнале из списка ВАК.
Объём и структура диссертации. Диссертация изложена на 145 страницах машинописного текста, состоит из Введения, Обзора литературы (гл. 1), главы 2 «Материалы и методы», результатов и обсуждения (гл. 3 и 4), Заключения, Выводов, Списка литературы, включающего 212 ссылок, из которых 152 - на английском языке и 4-х Приложений. Работа иллюстрирована 23 рисунками и содержит 19 таблиц.
Заключение Диссертация по теме "Генетика", Шпигальская, Нина Юрьевна
ВЫВОДЫ
1. Выявлена статистически значимая межпопуляционная генетическая гетерогенность чавычи Азии по аллельной изменчивости исследованных аллозимных локусов и гаплотипам мтДНК.
2. Для региональных групп популяций Восточной Камчатки, бассейна р. Большая, центральной Западной и Северо-Западной Камчатки характерно генетическое своеобразие, позволяющее по изменчивости частот аллозимных локусов идентифицировать данные регионы в смешанных морских скоплениях.
3. Создана основа для формирования базы данных по частотам аллелей аллозимных локусов чавычи. Выявлены наиболее информативные, с точки зрения анализа популяционной изменчивости, аллозимные локусы. Наиболее дискриминирующими в пределах азиатской части ареала являются 10 локусов - sAHl* ALAT* GPIA* sIDHP2* mMDH2* sMEPl* sMEP2* MPI* mSODl * TPI4*\ Азию и Северную Америку хорошо дифференцируют 23 локуса - mAATl* тААТ2* sAATl,2* sAATS*, ADA1* sAHl* GPIB2*, IDDH1 * sIDHPl* sIDHP2*, mMDH2* sMDHBl,2* sMEPl* sMEP2*, MPI*, PEPLT* PEP A* PEPB1 * PEPD2* PGK2* PGM2* sSODl*, TPI4*
4. Географическая изменчивость 28 аллозимных локусов чавычи имеет характер клинальной, что, вероятно, отражает особенности миграционных потоков в период колонизации современных речных систем.
5. Наиболее высокие оценки межпопуляционного разнообразия по аллозимным локусам свойственны североамериканским популяциям чавычи.
6. Невысокий уровень изменчивости митохондриальных гаплотипов, выявленных в камчатских популяциях, свидетельствует о небольшом эволюционном разнообразии азиатской чавычи.
7. Несмотря на достоверный уровень гетерогенности азиатской чавычи по структуре мтДНК, межпопуляционная дивергенция на этой части ареала в целом невелика.
8. Вероятной причиной относительно неглубокой дивергенции камчатских популяций чавычи по структуре мтДНК может являться их исторически недавнее происхождение. Максимальное время самостоятельной эволюции азиатской чавычи соответствует оценке во временном интервале, не превышающем 10-20 тыс. лет.
ЗАКЛЮЧЕНИЕ
В результате проведенных исследований на основе двух генетических методов - анализа изменчивости аллельных частот аллозимных локусов и анализа изменчивости структуры митохондриальной ДНК, описана генетическая дифференциация азиатских популяций чавычи, и проведена оценка популяционной дивергенции на всем ареале данного вида.
Сравнительный анализ распределений генных частот на протяжении всего ареала, включающего его азиатскую и американскую части, показал, что по многим аллозимным локусам наблюдается высокий уровень изменчивости региональных комплексов популяций, при этом часто происходит замена преобладающего аллеля. В исследованных популяциях обнаружено 28 локусов, географическая изменчивость которых имеет характер клинальной, что, вероятно, отражает особенности миграционных потоков в период колонизации современных речных систем в период отступления ледников после Плейстоценового оледенения.
Выявлена достоверная генетическая межпопуляционная гетерогенность азиатской чавычи. Региональная изменчивость генетической структуры популяций четко выражена и позволяет выделить два региональных популяционных комплекса - Западная и Восточная Камчатка.
Кластерный и многофакторный анализы генетических расстояний, определенных на основе изменчивости аллельных частот полиморфных аллозимных локусов, показали, что генетическая дифференциация соответствует пространственной подразделенности видового ареала. В структуре UPGMA-дендрограмм, также как и при многомерном шкалировании, кластеризация по генетической близости согласуется с географической приуроченностью популяций, что свидетельствует о генетическом своеобразии региональных популяционных комплексов. Создана реперная база данных, включающая мультилокусные характеристики по аллозимным генам большинства крупных популяционных систем ареала вида, на основе которых может осуществляться идентификация региональной принадлежности смешанных морских уловов в период нагульных и преднерестовых миграций чавычи.
Получены результаты исследования полиморфизма мтДНК у чавычи азиатских популяций. ПДРФ-анализ с применением восьми эндонуклеаз рестрикции позволил выявить изменчивость нуклеотидных последовательностей в трех амплифицированных участках мтДНК по 6 полиморфным рестриктным сайтам. Анализ распределения частот комбинированных гаплотипов чавычи Камчатки и изменчивости числа нуклеотидных замен позволил оценить внутри- и межпопуляционные генетические различия. При сравнении выборок из локальных популяций чавычи Азии статистически подтверждена их гетерогенность по структуре митохондриального генома, что свидетельствует о внутривидовой генетической подразделенности. При попарном сравнении выборок, между большинством сравниваемых пар обнаружены значимые различия. Полученные результаты свидетельствуют о возможности использования митохондриального генома в различных аспектах популяционных исследований в качестве информативного генетического маркера.
Выявленные отличия в числе нуклеотидных замен в мтДНК камчатских популяций свидетельствуют о том, что азиатской чавыче свойственно небольшое эволюционное разнообразие, что подтверждается относительно малыми межпопуляционными генетическими расстояниями и, соответственно, невысоким уровнем генетической дивергенции. Вероятной причиной относительно неглубокой дивергенции камчатских популяций по мтДНК может являться их исторически недавнее происхождение из предковой популяции с единым генным пулом. В этом случае небольшие межпопуляционные различия могут быть обусловлены небольшим временем самостоятельной эволюции.
Высказано предположение, что максимальное время дивергенции для камчатских популяций чавычи, величина различий между которыми составляет около 0.02% нуклеотидных замен, соответствует оценке во временном интервале, не превышающем 10-20 тыс. лет. Вероятно, отступление позднеплейстоценового ледника обусловило начало процессов расселения, реколонизации и формирования современной картины подразделенности вида на структурные компоненты. Представленные ранее палеогеографические данные показывают, что описанные особенности межпопуляционной генетической дифференциации согласуются с историей позднеплейстоценовых оледенений, которые были ясно выражены на Камчатке и побережье Охотского моря (Баранова, Бискэ, 1964; Мелекесцев и др., 1974; Глушкова, 1984; Беспалый, Глушкова, 1987). После отступления ледника, сохранившиеся в различных рефугиумах популяции, занимали экологически соответствующие ниши (Алтухов, 1997). Свидетельствами определенных сценариев расселения и процессов формирования современной внутривидовой структуры являются, в первую очередь, генетические особенности популяций.
102
Библиография Диссертация по биологии, кандидата биологических наук, Шпигальская, Нина Юрьевна, Петропавловск-Камчатский
1. Алтухов Ю.П. Генетические процессы в популяциях. М.: ИКЦ «Академкнига». 2003.431 с.
2. Алтухов Ю.П. Локальные стада рыб как генетически стабильные популяционные системы // Биохимическая генетика рыб. JL: Ин-т цитологии АН СССР. 1973. С. 43-53.
3. Алтухов Ю.П. Популяционная генетика рыб. М.: Пищ. пром-сть. 1974. 274 с.
4. Алтухов Ю.П., Салменкова Е.А. Полиморфизм ДНК в популяционной генетике //Генетика. 2002. Т. 38. № 9. С. 1173-1195.
5. Алтухов Ю.П., Салменкова Е.А., Омельченко В.Т. Популяционная генетика лососевых рыб. М.: Наука. 1997. 287 с.
6. Алтухов Ю.П., Салменкова Е.А., Сачко Г.Д. Дупликация и полиморфизм генов лактатдегидрогеназы у тихоокеанских лососей // ДАН СССР. 1970. Т. 195. №3. С. 711-714.
7. Аронштам Л.А., Боркин Л.Я., Пудовкин А.И. Изоферменты в популяционной и эволюционной генетике / Генетика изоферментов. М. 1977. С. 199— 249.
8. Атлас распространения в море различных стад тихоокеанских лососей в период весенне-летнего нагула и преднерестовых миграций. Под ред. О.Ф. Гриценко. 2000. М.: ВНИРО. 190 с.
9. Баранова Ю.П., Бискэ С.Ф. Северо-Восток СССР (История развития рельефа Сибири и Дальнего Востока). М.: Наука. 1964. 289 с.
10. Беспалый В.Г., Глушкова О.Ю. Северо-Восток // Четвертичные оледенения на территории СССР. М.: Наука. 1987. С. 62-69.
11. Бирман И.Б. Морской период жизни и вопросы динамики стада тихоокеанских лососей. М.: Изд-во ФГУП «Национальные рыбные ресурсы». 2004. 172 с.
12. Брайцева О.А., Мелекесцев И.В., Евтеева И.С., Лупикина Е.Г.
13. Стратиграфия четвертичных отложений и оледенения Камчатки. М.: Наука. 1968. 266 с.
14. Брыков В.А. Эволюция генома, изменчивость и дивергенция ДНК у морских животных. Автореф. дис. . д-ра биол. наук / Ин-т общей генетики им. Н.И. Вавилова РАН. Москва. 2001. 48 с.
15. Брыков Вл.А., Полякова Н.Е., Подлесных А.В. Дивергенция митохондриальной ДНК в популяциях нерки {Oncorhynchus пегка Walbaum) озера Азабачьего // Генетика. 2002. Т. 38. № 11. С. 1524-1529.
16. Брыков Вл.А., Полякова Н.Е., Подлесных А.В., Голубь Е.В., Голубь А.П., Жданова О.Л. Влияние биотопов размножения на генетическую дифференциацию популяций нерки Oncorhynchus пегка II Генетика. 2005. Т. 41. №5. С. 1-10.
17. Бугаев В.Ф., Вронский Б.Б., Заварина Л.О., Зорбиди Ж.Х., Остроумов А.Г., Тиллер И.В. Рыбы реки Камчатка. Петропавловск-Камчатский: Изд-во КамчатНИРО, 2007. 459 с.
18. Варнавская Н.В. Генетическая дифференциация популяций тихоокеанских лососей. Петропавловск-Камчатский: Изд-во КамчатНИРО. 2006. 488 с.
19. Варнавская Н.В. Геногеография по аллозимным локусам у тихоокеанских лососей / В кн. «Популяционная биология, генетика и систематика гидробионтов». Сб. науч. трудов КамчатНИРО. Т. 1. Петропавловск-Камчатский: Изд-во КамчатНИРО. 2005. С. 28-62.
20. Варнавская Н.В. Принципы генетической идентификации популяций тихоокеанских лососей рода Oncorhynchus spp. в связи с задачами рационального промысла. Автореф. дис. . д-ра биол. наук / Ин-т общей генетики им. Н.И. Вавилова РАН. Москва. 2001. 48 с.
21. Варнавская Н.В., Варнавский B.C., Мидяная В.В., Ростомова С.А., Збоева Е.Н. Анализ генетической изменчивости в популяциях нерки Oncorhynchus nerka (Walbaum) Азии и Северной Америки методами многомерной статистики//Генетика. 1996. № 7. С. 962-977.
22. Варнавская Н.В., Шпигальская Н.Ю., Шапорев Р.А. Изменчивость частот ферментных генов в популяциях чавычи, Oncorhynchus tschawytscha
23. Вронский Б.Б. Сезонные расы чавычи в бассейне р. Камчатки. Петропавловск
24. Камчатский: Архив КамчатНИРО. 1979. 42 с. Вронский Б.Б. Сезонные расы чавычи Oncorhynchus tschawytscha (Walbaum) в бассейне р. Камчатка // Тез. докл. X Всесоюз. Симп. Магадан. 1983. С. 159-160.
25. М.: Наука. 1977. 278 с. Ли Ч. Введение в популяционную генетику. М.: Мир. 1978. 526 с. Маниатис Т., Фрич Э., Сэмбрук Дж. Молекулярное клонирование. М.: Мир. 1984. 479 с.
26. Чуриков Д.Ю. Генеалогия гаплотипов митохондриальной ДНК у нескольких видов тихоокеанских лососей. Автореф. дис. . канд. биол. наук // Санкт-Петербургский Государственный университет. СПб. 2001. 17 с.
27. Шпигальская Н.Ю., Брыков Вл.А., Гарретт А.Дж., Кухлевский А.Д., Шапорев Р.А., Варнавская Н.В. Межпопуляционная изменчивость митохондриальной ДНК чавычи, Oncorhynchus tschawytscha (Walbaum), Камчатки // Генетика. 2008. Т. 44. № 7. С. 972-982.
28. Aebersold P.B., Winans G.A., Teel D.G., Milner G.B., Utter F.M. Manual for starch gel electrophoresis: a method for the detection of genetic variation // NOAA Tech. Rep. NMFS. 1987. V. 61. 19 p.
29. Allendorf F.W., Ryman N. Utter F.M. Genetics and fisheries management: Past, present and future // In: «Population Genetics and Fishery Management» Univ. Wash. Press. Seattle. (Eds. N. Ryman & F. Utter). 1987. P. 1-19.
30. Altukhov Yu.P. Population genetics: Diversity and stability. L.: Harwood. 1990. 3671. P
31. Altukhov Yu.P. Recent physiological, biochemical and immunological studies on the problem of intraspecific differentiation of marine fish // Rapp. et proc. verb. Cons, intern, explor. mer. 1971. V. 161. P. 103-108.
32. Altukhov Yu.P. The stock concept from the viewpoint of population genetics // Canad. J. Fish. Aquat. Sci. 1981. V. 38. № 12. P. 1523-1538.
33. Altukhov Yu.P., Salmenkova E.A. Applications of the stock concept to fish populations in the USSR // Canad. J. Fish. Aquat. Sci. 1981. V. 38. № 12. P. 1591-1600.
34. Altukhov Yu.P., Salmenkova E.A. The genetic structure of salmon populations // Aquaculture. 1991. V. 98. P. 11-40.
35. Altukhov Yu.P., Salmenkova E.A., Omelchenko V.T. Salmonid fishes: Population biology, genetics and managment. Oxford: Blackwell. 2000. 354 p.
36. Anderson S., Bankier A.T., Barrel B.G. et al. Sequence and organization of the human mitochondrial genome //Nature. 1981. V. 190. P. 457-465.
37. Avise J.C. Molecular Markers, Natural History and Evolution. N.Y.; London: Chapman and Hall. 1994. 511 p.
38. Avise J.C. Phylogeography. The history and formation of species. Harvard Univ. Press. Cambridge. MA, London, England. 2000. 447 p.
39. Avise J.C., Lansman R.A. Polymorfism of mitochondrial DNA in population of higher animals // In «Evolution of Genus & Proteins». Sinauer Ass. Inc. Publ. Sounderland (Mass.). 1983. P. 147-184.
40. Avise J.C., Neigel J.E., Arnold J. Demographic influences on mitochondrial DNA lineage surviorship in animalpopulations // J. Mol. Evol. 1984. V. 20. P. 90105.
41. Avise J.C., Walker D., Johns G.C. Speciation durations and Pleistocene effects on vertebrate phylogeography // Pr. R. Soc. Lond. B. 1998. V. 265. P. 17071712.
42. Bailey G.S., Tsuyuki H., Wilson A.C. The number of genes for lactate dehydrogenase in salmonid fishes // J. Fish. Res. Bd. Can. 1976. V. 33. P. 760-767.
43. Beacham T.D., Jonsen K.L., Supernault J., Wetklo M., Deng L., Varnavskaya N.
44. Pacific rim population structure of chinook salmon as determined from microsatellite analysis // Trans. Am. Fish. Soc. 2006b. V. 135. P. 1604-1621.
45. Beacham T.D., Murrey C.B., Withler R.E. Age, morphology and biochemical genetic variation of Yukon river chinook salmon // Trans. Amer. Fish. Soc. 1989. V. 118. P. 46-63.
46. Beacham T.D., Withler R.E., Murrey C.B., Barner L.W. Variation in body size, morphology, egg size, and biochemical genetics of pink salmon in British Columbia // Trans. Am. Fish. Soc. 1988. V. 117. № 2. P. 109-126.
47. Beacham T.D., Withler R.E., Stevens T.A. Stock identification of chinook salmon {Oncorhynchus tschawytscha) using minisatellite DNA variation // Can. J. Fish. Aquat. Sci. 1996. V. 53. P. 380-394.
48. Billington N., Hebert P.D.N. Mitochondrial DNA diversity in fish and its implications for introductions // Can. J. Fish. Aquat. Sci. 1991. V. 48. P. 8094.
49. Brodziak J., Bentley В., Bartley D., Gall G.A., Gomulkiewicz R., Mangel M. Test of genetic stock identification using coded wire tagged fish // Can. J. Fish. Aquat. Sci. 1992. V. 49. P. 1507-1517.
50. Brookes A.J. The essence of SNPs // Gene. 1999. V. 234. № 2. P. 177-186.
51. Brown W.M., George M., Wilson A.C. Rapid evolution of animal mitochondrial DNA // Proc. Nat. Acad. Sci. US. 1979. V. 76. P. 1967-1971.
52. Brykov V.A., Polyakova N.E., Skurikhina L.A., Kukhlevsky A.D. Geographical and temporal mitochondrial DNA variability in populations of pink salmon // J. Fish. Biol. 1996. V. 48. P. 899-909.
53. Burger C.V., Wilmot R.L., Wangaard D.B. Comparison of spawning areas and times for two runs of chinook salmon (Oncorhynchus tschawytscha) in the Kenai river, Alaska // Can. J. Fish. Aquat. Sci. 1985. V. 42. № 4. P. 693-700.
54. Cavalli-Sforsa L.L., Edwards A.W.E. Phylogenetic analysis: models and estimation procedures // Am. J. Hum. Genet. 1967. V. 19. P. 233-257.
55. Churikov D.U., Gharrett A.J. Comparative phylogeography of the two pink salmon brodlines: An analysis based on a mitochondrial DNA genealogy // Mol. Ecol. 2002. V. 11. P. 1077-1101.
56. Crane P.A., Templin W.D., Seeb L.W. Genetic stock identification of Alaska chinook salmon. Final Rep., Fish & Game, NOAA, Anchorage, No. NA46FD0356. 1996. 71 p.
57. Crow J.F. Spontaneus mutation as a risk factor // Exp. Clin. Immunogenet. 1995. V. 12. №3. P. 121-128.
58. Excoffier L., Smouse P.E., Quattro J.M. Analyses of molecular variance inferred from metric distance among DNA gaplotypes: Application to human mitochondrial DNA restriction data // Genetics (US). 1992. V. 131. P. 479491.
59. Ferris S.D., Berg W.J. The utility of mitochondrial DNA in fish genetics and management // Population genetics and fishery management / Ed. N. Ryman, F. Utter. Seattle; L.: Univ. Wash, press. 1987. P. 277-301.
60. Ferris S.D., Sage R.D., Prager E.M. et al. Mitochondrial DNA evolution in mice // Genetics. 1983. V. 105. P. 681-721.
61. Ford E.B. Polymorphism and taxonomy // The new systematics. Oxford: Clarendon press. 1940. P. 493-513.
62. Gall G.A.E., Bartley D.M., Bentley В., Brodziak J., Gomulkiewicz R., Mangel M. Geographic variation in population genetic structure of chinook salmon from California and Oregon // Fish. Bull. U.S. 1992. V. 90. P. 77-100.
63. Gharrett A.J., Gray A.K., Brykov V.A. Phylogeographical analysis of mitochondrial DNA variation in Alaskan coho salmon, Oncorhynchus kisutch //Fish. Bull. 2001. V. 99. P. 528-544.
64. Gharrett A. J., Shirley S.M., Tromble G.R. Genetic relationship among populations of Alaskan chinook salmon {Oncorhynchus tschawytscha) // Can. J. Fish. Aquat. Sci. 1987. V. 44. P. 165-114.
65. Gustafson R.G., Winans G.A. Distribution and population genetic structure of river-and sea-type sockeye salmon in western North America // Ecology of Freshwater Fish. 1999. V. 8. P. 181-193.
66. Hawkins S.L., Wilmot R.L. Development of Pacific Rim even-and odd-year pink salmon baselines / Proc. NPAFC Int. Workshope on Appl. of Stock Identification in Defining Marine Distribution and Migration of Salmon. Honolulu, Hawaii, USA. Nov. 2003. P. 67.
67. Healey M.C. Coastwide distribution and ocean migration patterns of stream- and ocean-type chinook salmon, Oncorhynchus tschawytscha II Can. Field-Naturalist. 1983. V. 97. P. 427-433.
68. Healey M.C. The life history of chinook salmon (Oncorhynchus tschawytscha) // Pacific salmon life histories / Ed. C. Groot, L. Margolis. Vancouver: UBC Press, 1991. P. 311-393.
69. Hedrick P.W., Rashbrook V.K., Hedgecock D. Effective population size of winter-run chinook salmon based on microsatellite analysis of returning spawners // Can. J. Fish. Aquat. Sci. 2000. V. 57. P. 2368-2373.
70. Hodgins H., Ames W.E., Utter F.M. Variants of lactate dehydrogenase isosymes in sera of sockeye salmon (Oncorhynchus nerka) II J. Fish. Res. Board Canada. 1969. V. 26. № i.p. 15-19.
71. Hunter R.L., Markert C.L. Histochemical demonstration of enzymes separated by zone electrophoresis in starch gels // Science. 1957. V. 125. P. 1294-1295.
72. Jeffreys A.J., Wilson V., Thein S.L. Hypervariable "minisatellite" regions in human DNA //Nature. 1985a. V. 316. P. 67-73.
73. Jeffreys A.J., Wilson V., Thein S.L. Individual specific "fingerprints" of human DNA//Nature. 1985b. V. 316. P. 76-79.
74. Jin L., Macaubas C., Hallmayer J., Kimura A., Mignot E. Mutation rate varies among alleles at a microsatellite locus: Phylogenetic evidence // Proc. Natl. Acad. Sci. USA. 1996. V. 93. P. 15285-15288.
75. Kobayashi Т., Milner G.B., Teel D.J., Utter F.M. Genetic basis for electrophoretic variation of adenosine deaminase in chinook salmon // Trans. Am. Fish. Soc. 1984. V. 113. P. 86-89.
76. Kondzela C.M. Guthrie C.M., Hawkins S.L., Russell C.D., Helle J.H. Gharrett A.J. Genetic relationship among chum salmon population in southeast Alaska and northern British Columbia // Canad. J. Fish. Aquat. Sci. 1994. V. 51. suppl. 1. P. 50-64.
77. Kristiansson A.C., Mclntyre D. Genetic variation in chinook salmon (Oncorhynchus tshawytscha) from the Columbia River and three Oregon coastal rivers //Trans. Amer. Fish. Soc. 1976. V. 105. P. 620-623.
78. Kruskal J.B. Nonmetric multidimensional scaling: a numerical method // Psychometrika. 1964. V. 29. P. 28^12.
79. W.H., Ellsworth D.L., Krushkal J., Chang B.H.J., Hewett-Emmett D. Rates of nucleotide substitution in primates and rodents and the generation-time effect hypothesis //Mol. Phylogenet. Evol. 1996. V. 5. № 1. P. 182-187.
80. Major R.L., Ito J., Ito S., Godfrey H. Distribution and origin of chinook salmon (Oncorhynchus tschawytscha) in offshore waters of the North Pasific Ocean // Bull. INPFC. 1978. № 38. P. 1-54.
81. Manchenko G.H. Handbook of detection of enzymes on electrophoretic gels. Boca Raton: CRC press. 1994. 341 p.
82. McElroy D.M., Moran P., Bermingham E., Kornfield I. REAP: An integrated environment for the manipulation and phylogenetic analysis of restriction data // J. Heredity. 1992. V. 83. P. 153-158.
83. McPhail J.D., Lindsey C.C. Freshwater fishes of Northwestern Canada and Alaska //Bull. Fish. Res. Board Canada. 1970. № 173. P. 1-375.
84. McPhail J.D., Lindsey C.C. Zoogeography of the freshwater fishes of Cascadia (the Columbia system and rivers North to the Stikine) // Zoogeography of North American freshwater fishes. Copyright (c). N.Y.: Wiley. 1986. P. 615-637.
85. Milner G.B., Teel D.J., Utter F.M., Winans G.A. A genetic method of stock identification in mixed populations of Pacific salmon, Oncorhynchus spp. // US Nat. Mar. Fish. Serv. Mar. Fish. Rev. 1985. V. 47. № 1. P. 1-8.
86. Molecular Systematics (ed. Hills D.M., Moritz C., Mable B.K.). Sinauer Ass. Inc. Publ. Sounderland (Mass.). 1996. 655 p.
87. Moritz C., Dowling Т.Е., Brown W.M. Evolution of animal mitochondrial DNA: Relevance for populatuon biology and sistematics // Annu. Rev. Ecol. Syst. 1987. V. 18. P. 269-292.
88. Mullis K., Faloona F., Scharf S., Saiki R.K., Horn G., Elrich H. Specific enzymatic amplification of DNA in vitro: the polymerase chain reaction // Cold Spring Harbor Symp. Quant. Biol. 1986. V. 51. P. 263-273.
89. Nei M. Molecular evolutionary genetics. N.Y.: Columbia Univ. Press. 1987. 512 p.
90. Nei M., Li W.-H. Mathematical model for shifting genetic variation in terms of restriction endonucleases // Proc. Natl Acad. Sci. USA. 1979. V. 76. P. 52695273.
91. Nei M., Tajima F. DNA polymorphism detectable by restriction endonucleases // Genetics. 1981. V. 105. P. 207-217.
92. Nei M., Tajima F. Maximum likelihood of the number of nucleotide substitution from restriction sites data // Genetics. 1983. V. 105. P. 207-217.
93. Nelson R.J., Small M.P., Beacham T.D., Supernault K.J. Population structure of Fraser River chinook salmon (iOncorhynchus tschawytscha): an analysis using microsatellite DNA markers // Fish. Bull. 2001. V. 99. P. 94-107.
94. O'Malley K.G., Camara M.D., Banks M.A. Candidate loci reveal genetic differentiation between temporally divergent migratory runs of chinook salmon {Oncorhynchus tschawytscha) // Mol. Ecol. 2007. V. 23. № 16. P. 4930-4941.
95. Park L.K., Moran P. Developments in molecular genetic techniques in fisheries // Rev. Fish Biol. 1994. V. 4. P. 272-299.
96. PCR Technology. Principles and Applications for DNA Amplification (ed. H.A. Erlich). Stockton Press. 1989. 246 p.
97. Pella J.J., Milner G.B. Use of genetic marks in stock compositions analysis // Population genetics and fishery management / Eds. N. Ryman, F. Utter. Seattle; L.: Univ. Wash, press. 1987. P. 247-276.
98. Phelps S.R., Le Clair L.L. Young S., Blankenship H.L. Genetic diversity of chum salmon in the Pasific northwest // Canad. J. Fish. Aquat. Sci. 1994. V. 51. suppl. l.P. 65-83.
99. Quinn T.P. A review of homing and straying of wild and hatchery-produced salmon //Fish. Res. 1993. V. 18. P. 29^14.
100. Quinn T.P., Fresh К. Homing and straying in chinook salmon (Oncorhynchus tschawytscha) from Cowlitz River Hatchery, Washington // Canad. J. Fish. Aquat. Sci. 1984. V. 41. P. 1078-1082.
101. Reisenbichler R.R., Phelps S.R. Genetic variation in chinook, Oncorhynchus tschawytscha, and coho O. kisutch, salmon from the north coast of Washington // Fish. Bull. 1987. V. 85. № 4. P. 681-701.
102. Rogers A.R. Genetic evidence for Pleistocene population explosion // Evolution. 1995. V. 49. P. 608-615.
103. Rogers A.R., Harpending H. Population growth makes waves in the distribution of pairwise genetic distances // Mol. Biol. Evol. 1992. V. 5. P. 552-569.
104. Rolf F.J. NTSYS-ps: Numerical taxonomy and multivariate analysis system. Version. 1.60. N.Y.:Exter publishing, Ltd. 1990.
105. Saiki R.K., Gelfand D.H., Stoffel S. et al. Premier-directed enzymatic amplification of DNA with termostable DNA polymerase // Science. 1988. V. 239. P. 487491.
106. Sambrook J., Fritsch E.F. Maniatis T. Molecular Cloning: A Laboratory Manual. N.Y.: Cold Spring Harbor Lab. Press. 1989. 1626 p.
107. Schneider S., Roessli D., Excoffier L. Arlequin ver. 2.000: A software for population genetics data analysis. Genetics and Biometry Laboratory, Univ. Geneva, Switzerland. 2000.
108. Shaklee J.B., Allendorf F.W., Morizot D.C., Whitt G.S. Gene nomenclature for protein-coding loci in fish // Trans. Amer. Fish. Soc. 1990a. V. 119. P. 2-15.
109. Shaklee J.B., Phelps S.R. Chinook salmon NADP+-dependent cytosolic isocitrate dehydrogenase: Electrophoretic and genetic dissection of a complex isozyme system and geographic patterns of variation // Biochem. Genet. 1992. V. 30. №9/10. P. 455-489.
110. Smith G.R. Introgression in fishes: significance for paleontology, cladistics, and evolutionary rates // Systematic Biology. 1992. V. 41. P. 41-57.
111. Smith S.T., Seeb J.E., Schwenke P., Seeb L.W. Use of the 5' nuclease reaction for single nucleotide polymorphism genotyping in chinook salmon // Trans. Amer. Fish. Soc. 2005. V. 134. P. 207-215.
112. Smith S.T., Seeb L.W., Seeb J.E. Single nucleotide polymorphism (SNP) genotyping using the 5' nuclease reaction provides high-throughput/high-resolution stock identification of chinook, sockeye, and chum salmon // NPAFC Doc. 733. 2003. 11 p.
113. Sneath P.H.A., Sokal R.R. Numerical Taxonomy. Freeman W.H. ed.. San-Francisco, Ca. 1973. 573 p.
114. Sokal R.R., Rohlf F.G. Biometry // 2nd ed. W.H. Freemen & Co., San Francisco, CA. 1981. 859 p.
115. Springer M. SNPs a great catch for salmon genotyping // American Laboratory: 2006. V. 38. № 10. P. 34-38.
116. Tautz D. Hypervariability of symple sequences as a general sourse for polymorphic DNA markers //Nucl. Acids Res. 1989. V. 17. P. 6463-6471.
117. Tautz D., Trick M., Dover G.A. Cryptic simplicity in DNA is major source of genetic variation // Nature. 1986. V. 322. P. 652-656.
118. Teel D.J. Genetic mixed stock analysis of juvenile chinook salmon in coastal areas of western North America // NPAFC Tech. Rep. 2004. № 5.
119. Teel D.J., Crane P.A., Guthrie Ch., Marshall A.R., Van Doornik D.M., Templin W.D., Varnavskaya N.V., Seeb L.W. Comprehensive allozyme database discriminates chinook salmon around the Pasific Rim / NPAFC Doc. 440. 1999. 25 p.
120. Teel D.J., Milner G.B., Winans G.A., Grant W.S. Genetic population structure and origin of life history types in chinook salmon in British Columbia, Canada // Trans. Amer. Fish. Soc. 2000. V. 129. P. 194-209.
121. Templin W.D., Smith S.T., Seeb L.W., Seeb J.E. SNP provide high-throughput resolution for migratory studies of Chinook salmon // NPAFC Doc. 908. 2005. 10 p.
122. Thomas W.K., Beckenbach A.T. Variation in salmonid mitochondrial DNA: evolutionary constraints and mechanisms of substitution // J. Mol. Evol. 1989. V. 29. P. 233-245.
123. Utter F.M., Aebersold P.B., Winans G.A. Interpreting genetic variation detected by electrophoresis / In: «Ryman N. & Utter F. eds. Population Genetics And Fisheries Management». Seattle: University of Washington Press. 1987a. P. 21-46.
124. Utter F.M., Allendorf F.W., May B. Genetic basis of creatine kinase isozymes in skeletal muscle of salmonid fishes // Biochem. Genet. 1979. V. 17. № 11/12. P. 1079-1091.
125. Utter F.M., Chapman D.W., Marshall A.R. Genetic population ctructure and history of chinook salmon of the Upper Columbia river // Amer. Fish. Soc. Symp. 1995. V. 17. P. 149-165.
126. Utter F.M., Milner G.B., Stahl G., Teel D.J. Genetic population structure of chinook salmon, Oncorhynchus tschawytscha, in the Pasific Norhwest U.S. // Fish. Bull. 1989. V. 87. № 2. P. 239-264.
127. Utter F.M., Teel D.J., Milner G.B., Mclsaac D. Genetic estimates of stock composition of 1983 chinook salmon, Oncorhynchus tschawytscha, harvests off the Washington coast and Columbia River /U.S. Nat. Mar. Fish. Serv. Fish. Bull. 1987b. V. 85. P. 13-23.
128. Utter F.M., Waples R.S., Teel D.J. Genetic isolation of previously indistinguishable chinook salmon populations of the Snake and Klamath Rivers: Limitations of negative data//Fish. Bull. U.S. 1992. V. 90. P. 770-777.
129. Van Doornik D.M., Winans G.A. Inheritance of allozymes in coho salmon // Trans. Amer. Fish. Soc. 1998. V. 127. № 5. P. 833-837.
130. Varnavskaya N.V., Beacham T.D. Biochemical genetic variation in pink salmon {Oncorhynchus gorbuscha) from Kamchatka // Can. J. Zool. 1992. V. 70. P. 1788-1792.
131. Waples R.S., Do C. Genetic risk associated with supplementation of Pacific salmonids: captive broodstock programs // Can. J. Fish. Aquat. Sci. 1994. V. 51. P. 310-329.
132. Waples R.S., Teel D.J., Myers J.M., Marshall A.R. Life-history divergence in chinook salmon: Historic contingency and parallel evolution // Evolution. 2004. V. 58. P. 386-403.
133. Waples R.S., Winans G.A., Utter F.M., Mahnken C. Genetic approaches to the management ofPasific salmon//Fisheries. 1990a. V. 15. № 5. P. 19-25.
134. Waples R.S., Winans G.A., Utter F.M., Mahnken C. Genetic monitoring of Pacific salmon hatcheries // NOAA Tech. Rep. 1990b. № 92. P. 33-37.
135. Wilhelm V., Villegas J., Miquel A., Engel E., Bernales S., Valenzuela P.D.T., Burzio L.O. The complete sequence of the mitochondrial genome of the chinook salmon, Oncorhynchus tschawytscha II Biological Research. 2003. V.36. P. 223-231.
136. Wilmot R.L., Everett R., Spearman W.J., Baccus R. Genetic stock identification of Yukon River chum and chinook salmon 1987-1990 / U.S. Fish and Wildlife Service. Prog. Rep. Anchorage. 1992. 132 p.
137. Wilmot R.L., Everett R.J. Spearman W.J., Baccus R., Varnavskaya N.V., Putivkin S.V. Genetic stock structure of western Alaska chum salmon and comparison with Russian Far East stock // Can. J. Fish. Aquat. Sci. 1994. V. 51. suppl. l.P. 84-94.
138. Wilson G.M., Thomas W.K., Beckenbach A.T. Mitochondrial DNA analysis of Pasific Nortthwest population of Oncorhynchus tschawytscha I I Can. J. Fish. Aquat. Sci. 1987. V. 44. P. 1301-1305.
139. Winans G.A., Aebersold P.B., Waples R.S. Allozyme variability of Oncorhynchus nerka in the Pasific Northwest with special consideration to populations of Redfish Lake, Idaho // Trans. Amer. Fish. Soc. 1996. V. 125. № 5. P. 645663.
140. Winans G.A., Urawa S. Allozyme variability of Oncorhynchus nerka in Japan // Ichthyol. Res. 2000. V. 47. № 4. P. 343-352.
141. Wright J.M., Bentzen P. Microsatellites: Genetic markers for the future // Rev. Fish. Biol. Fish. 1994. V. 4. P. 384-388.
142. Zardoya R., Garrido-Pertierra A., Bautista J.M. The complete nucleotide sequence of the mitochondrial DNA genome of the rainbow trout, Oncorhynchus mykiss//J. Mol. Evol. 1995. V. 41. P. 942-951.
143. Zaykin D.V., Pudovkin A.I. Two programs to estimate significance of Chi-square values using pseudo-probability test // J. Heredity. 1993. V. 84. P. 152.
- Шпигальская, Нина Юрьевна
- кандидата биологических наук
- Петропавловск-Камчатский, 2010
- ВАК 03.02.07
- Экология чавычи и кижуча азиатских стад в морской период жизни
- Особенности воспроизводства кеты Oncorhynchus кета в связи с ее размерно-возрастной структурой
- Экология и структура популяции симы Oncorhynchus masou (Brevoort) на севере ареала
- Биология и динамика численности горбуши зал. Анива
- Биология, динамика численности и рациональное использование рыб Камчатки и прилегающих морских акваторий