CC BY
131
СЕЛЬСКОХОЗЯЙСТВЕННАЯ БИОЛОГИЯ, 2015, том 50, № 6, с. 776-784
УДК 638.123.54:575.174:[577.2.08+591.4.08 doi: 10.15389/agrobiology.2015.6.776rus
МОРФОМЕТРИЧЕСКАЯ И МОЛЕКУЛЯРНО-ГЕНЕТИЧЕСКАЯ ДИФФЕРЕНЦИАЦИЯ ЛИНИЙ И СЕМЕЙ МЕДОНОСНОЙ ПЧЕЛЫ Apis mellifera caucasica L., РАЗВОДИМЫХ В РАЙОНЕ БОЛЬШОГО СОЧИ*
М.С. ФОРНАРА1, А.С. КРАМАРЕНКО1, 2, С.В. СВИСТУНОВ3, Е.М. ЛЮБИМОВ3,
С.С. СОКОЛЬСКИЙ3, Н.А. ЗИНОВЬЕВА1
Создание специализированных линий — один из приемов генетического совершенствования и сохранения биоразнообразия пород и популяций медоносной пчелы. Серая горная кавказская порода (Apis mellifera caucasica L.) выделяется среди прочих по комплексу хозяйственно полезных признаков (самый длинный хоботок, высочайшая работоспособность и продуктивность, использование широкого видового состава медоносов, чрезвычайно слабая ройливость, миролюбие и др.). Целью настоящего исследования была сравнительная оценка разнообразия и степени дифференциации линий A m. caucasica на основании данных морфометрического анализа и микросателлитов (МС). Материалом для исследований служили рабочие пчелы пяти линий (I — линия № 19, с. Медовеевка; II — линия № 13, ПГТ Красная Поляна; III — линия № 12, ПГТ Красная Поляна; IV — линия № 49, с. Калиновое Озеро; V — линия № 34, с. Аибга; n = 728) с пяти пасек в районе Большого Сочи (Краснодарский край). Для морфометрического анализа измеряли длину хоботка (ДХ, мм), ширину 3-го тергита (ШТ, мм) и рассчитывали кубитальный индекс (КИ). Молекулярно-генетические исследования проводили по семи локусам МС (A024, A88, A113, AP043, HB-C16-05, HB-THE-03, HB-C16-01). Степень межсемейной изменчивости по морфометрическим показателям определяли методом двухфакторного иерархического дисперсионного анализа. Генетические различия между семьями по МС оценивали по показателям Fst при парном сравнении. Полученную матрицу использовали для PCA-анализа (principal component analysis). Для количественной оценки межсемейной изменчивости внутри линий рассчитывали индексы Fst, Rst (AMOVA). Степень дифференциации линий по морфометрическим признакам определяли посредством расчета евклидовых дистанций. Полученные значения использовали для построения дендрограммы подобия методом одиночной связи (single linkage) в алгоритме иерархической кластеризации. Степень дифференциации линий по МС оценивали на основании расчета значений генетических дистанций по М. Nei. Дендрограмму подобия строили с помощью метода UPGMA. Для расчетов использовали программное обеспечение STATISTICA, GenAlEx (v. 6.5.1), PAST (v. 3.03). Результаты морфометрического анализа показали наличие достоверных различий между линиями по ДХ и Ш3Т, тогда как в отношении КИ различий не обнаружили. Установлена большая гетерогенность в отношении изучаемых признаков в линиях II и V и, напротив, большая консолидация в линиях III и IV. Пчелы линии I существенно отличались от остальных по обоим показателям, однако характеризовались значительными межсемейными различиями по ДХ. Анализ МС профилей выявил сходные тенденции при оценке степени внутрисемейной и межсемейной изменчивости. В линии I наблюдается избыток гетерозигот (Fjj = -0,048), что можно рассматривать как указание на высокую гетерогенность. Пчелы этой линии характеризовались минимальной индивидуальной (F,t = 0,052) и максимальной межсемейной (Fst = 0,124) изменчивостью. В линиях II-V наблюдали дефицит гетерозигот (Fjj = 0,062-0,128), относительно высокую индивидуальную изменчивость (Fit = 0,143-0,189) и более низкие по отношению к линии I значения межсемейной изменчивости (Fst = 0,095-0,104). Наименьшие межсемейные различия были отмечены в линиях III и IV (Fst 0,096 и 0,095). Анализ дифференциации изучаемых линий по морфометрическим признакам и МС выявил различия в структуре генеалогического дерева. Дендрограмма, построенная по МС, была отражением географической удаленности изучаемых линий. Структура генеалогического дерева, построенного по морфометрическим признакам, не отражала ни географическую общность (различия) происхождения, ни сходство (различия) в хозяйственно полезных признаках у изучаемых линий. Таким образом, при исследовании морфометрических показателей и МС обнаружены сходные тенденции в оценке внутрисемейной и внутрилинейной изменчивости, однако наблюдались различия при дифференциации линий с использованием двух методов. В дальнейшем комплексный подход позволит с высокой точностью проводить определение не только породы пчел, но и более мелких таксономических групп. Полученные результаты могут быть использованы в селекционной работе по восстановлению чистоты пород медоносных пчел.
* Исследования выполнены при поддержке Российского научного фонда, проект 14-36-00039. При проведении исследований было использовано оборудование ЦКП «Биоресурсы и биоинженерия сельскохозяйственных животных» ВИЖ им. академика Л.К. Эрнста.
776
Ключевые слова: медоносная пчела, Apis mellifera caucasica L., ДНК-маркеры, морфометрия, микросателлиты, биоразнообразие.
Сохранение многообразия пород медоносной пчелы (Apis mellifera L., 1758), а также более мелких таксономических единиц (экотипов и линий) требует глубоких фундаментальных популяционно-генетических исследований состояния их генофонда (1). Внутри вида A. mellifera выделяют около 30 подвидов (2), для обозначения которых используют также термин раса, а в русскоязычной литературе — порода. На территории Российской Федерации разводится несколько пород: среднерусская, краинская (или карника), карпатская (некоторые авторы классифицируют ее как подвид карники 2-го порядка), серая горная кавказская. Последняя по комплексу хозяйственно полезных признаков (самый длинный хоботок, высочайшая работоспособность и продуктивность, использование широкого видового состава медоносов, чрезвычайно слабая ройливость, миролюбие и др.) считается одной из лучших (3). Структура серой горной кавказской породы неоднородна и состоит из популяций, достоверно различающихся по биологическим и хозяйственным признакам (4). Селекционная работа с породой, которая в настоящее время ведется на базе Краснополянской опытной станции (район Большого Сочи), направлена на сохранение генофонда и создание новых специализированных линий с улучшенными хозяйственно полезными качествами (5, 6).
Определение породной принадлежности пчел и их отнесение к той или иной линии — непростая процедура, поскольку различия между линиями в пределах одной породы, как правило, невелики. Кроме того, имеет место спонтанная гибридизация. Преимущественно используется оценка по экстерьерным, физиологическим, этологическим признакам и показателям хозяйственной ценности (2, 7). Наиболее часто измеряют длину хоботка, ширину 3-го тергита и вычисляют кубитальный индекс (8). Была показана возможность дифференциации подвидов медоносной пчелы A. m. mellifera и A. m. carnica с использованием морфометрических методов (9). Однако попытка внутривидовой дифференциации популяций пчел A. m. meda, обитающих в различных географических регионах (на расстоянии 500 км), на основании морфометрических данных оказалась неуспешной (10). Поскольку морфометрические признаки зависят от географического размещения пчел, их возраста, питания и широко варьируют даже в пределах одной породы или популяции, возникает необходимость в разработке подходов, основанных на анализе генотипа.
В последние годы все большую популярность для дифференциации пород, экотипов и линий медоносной пчелы находят молекулярно-генетические методы, включающие исследование аллозимов (11), митохондриальной ДНК (12-14), полиморфизма единичных нуклеотидов (15, 16) и микросателлитов (17-20). Оценка по ДНК-маркерам рассматривается в качестве дополнения к традиционному анализу по морфометрическим критериям (21, 22). Ранние работы, в которых ДНК-анализ использовался для идентификации внутривидовых таксономических единиц медоносной пчелы, базировались на определении уникальных (приватных) для соответствующей линии аллелей (23). Затем применение байесовских методов сделало возможным генетическую кластеризацию индивидуумов (подвидов и более мелких таксонов) на основе генотипа по множеству несцепленных локусов (24). Метод высокоэффективен не только при определении подвидов, относящихся к различным эволюционным ветвям A. mellifera (25, 26), но и для установления происхождения более близких таксонов (27).
В представляемой работе мы показали, что структура генеалогиче-
777
ского дерева, построенного по морфометрическим признакам, не отражала ни географическую общность происхождения, ни сходство в хозяйственно полезных признаках у изученных линий. Следовательно, для их адекватной характеристики и дифференциации необходимо сочетание морфометрических и молекулярно-генетических методов. Использованный нами комплексный подход позволит с высокой точностью определять не только породы пчел, но и более мелкие таксономические группы. Полученные результаты могут быть использованы в работе по восстановлению чистоты пород медоносных пчел.
Целью настоящей работы стала сравнительная оценка разнообразия и степени дифференциации линий Apis mellifera caucasica на основании данных морфометрического анализа и микросателлитов.
Методика. Исследования проводили в 2014 году. В качестве материала использовали рабочих пчел серой горной кавказской породы различных линий, разводимых на пяти пасеках в районе Большого Сочи (Краснодарский край): I — линия № 19 (с. Медовеевка), II — линия № 13 (ПГТ Красная Поляна), III — линия № 12 (ПГТ Красная Поляна), IV — линия № 49 (с. Калиновое Озеро) и V — линия № 34 (с. Аибга). Для анализа отбирали 5 семей каждой линии по 22-30 пчел в семье (25 семей, n = 728), в том числе I — n = 149, II — n = 147, III — n = 147, IV — n = 137, V — n = 148. Пчел консервировали в 96 % спирте и хранили при +4 °С.
Измеряли длину хоботка (ДХ, мм), ширину 3-го тергита (ШТ, мм) и рассчитывали кубитальный индекс (КИ) как отношение между большей и меньшей сторонами третьей кубитальной ячейки переднего крыла (8).
ДНК выделяли с помощью экстракции перхлоратом (28), материалом служила голова рабочей пчелы. Семь локусов микросателлитов (МС) — A024, A88, A113, AP043, HB-C16-05, HB-THE-03, HB-C16-01 амплифици-ровали в одной мультиплексной реакции с использованием флуоресцентно меченных праймеров. Реакции проводили в конечном объеме 15 мкл. В пробирки вносили по 14 мкл реакционной смеси (1 мкл 10* ПЦР-буфера, 0,5 мкл 2 мМ раствора dNTPs, 0,7 мкл 100 мМ раствора DMSO, 4,6 мкл 10 мМ смеси праймеров, 0,1 мкл Таq-полимеразы, 7,1 мкл бидистиллированной воды) и добавляли 1 мкл (50-100 нг) исследуемой геномной ДНК. Состав ПЦР-буфера: 16,6 мМ (NH4)2SO4, 67,7 мМ Трис-НС1 (рН = 8,8), 0,1 объема Tween 20. ПЦР выполняли в следующем режиме: начальная денатурация 8 мин при 94 °С; отжиг 2 мин при 55 °С; полимеризация 2 мин при 72 °С; 30 циклов — денатурация 30 с 94 °С, отжиг 35 с при 55 °С, элонгация 40 с при 72 °С; элонгация 20 мин при 72 °С. Разделение продуктов ПЦР выполняли на генетическом анализаторе ABI3130xl («Applied Biosystems», США) с использованием внутреннего стандарта ROX. Элек-трофореграммы обрабатывали с помощью программного обеспечения Ge-neScan v. 3.7 и Genotyper v. 3.7 («Applied Biosystems», США).
Межсемейную изменчивость по морфометрическим показателям оценивали методом двухфакторного иерархического дисперсионного анализа с помощью программного обеспечения Statistica v. 6 (29). В качестве фиксированных эффектов использовали принадлежность к линии и семье (иерархический подчиненный фактор). Рассчитывали Mst:
Mst = F - 1
F - 1 + n*
где F — дисперсионное отношение, n* — среднее геометрическое объема выборок. Величина Mst может принимать значения от 0 до 1 и оценивать меру морфологической дифференциации между отдельными популяциями и, по сути, представляет собой аналог Fst (30).
778
Генетические различия между семьями медоносной пчелы в отношении МС оценивали по показателю Fst (31) при парном сравнении и на основании генетических дистанций по M. Nei (32). Полученные попарные матрицы использовали для построения распределения центроидов семей в пространстве первых двух главных координат (PCA, principal component analysis). Для количественной оценки степени межсемейной изменчивости внутри линий рассчитывали идексы Fst, Rst (AMOVA) (33), используя программное обеспечение GenAlEx v. 6.5 (34).
Степень дифференциации пяти изучаемых линий A m. caucasica по морфометрическим признакам оценивали посредством расчета евклидовых дистанций. Полученные значения использовали для построения дендрограммы подобия методом одиночной связи (single linkage) в алгоритме иерархической кластеризации (29). Оценку степени дифференциации линий по МС проводили на основании расчета значений генетических дистанций по M. Nei (32). Дендрограмму подобия строили методом UPGMA при помощи программного обеспечения PAST (v. 3.03). Дендрограмму проверяли на устойчивость топологии кластеров посредством проведения 999 перму-таций bootstrap-методом.
Результаты. При формировании каждая из изученных линий пчел была строго ориентирована на получение материала с заданными хозяйственно полезными признаками. Матки линии III обладают высокой яйценоскостью в период подготовки к основному медосбору, а семьи — повышенной продуктивностью по меду. Линии IV и V характеризуются ранним развитием по сравнению с линиями I, II и III. Пчелы линии IV менее зимостойки, в то время как линия II выделяется повышенной, а линия I — высокой зимостойкостью.
Были выявлены достоверные различия между линиями по двум из трех оцененных признаков — длине хоботка и ширине 3-го тергита, тогда как по значениям кубитального индекса различий между линиями мы не обнаружили (табл. 1). У пчел линии I по сравнению с другими изучаемыми линиями величина ДХ оказалась на 0,05-0,07 мм меньше (p < 0,001), а Ш3Т — на 0,04-0,11 мм больше (p < 0,001). Максимальную длину хоботка отмечали у пчел линии III. У пчел линий III и IV Ш3Т была на 0,08-0,11 мм меньше, чем у других линий (p < 0,001). Обнаруженные различия могут быть следствием разной стратегии при селекционно-племенной работе с этими линиями.
1. Морфометрические показатели у различных линий медоносных пчел серой горной кавказской породы (Apis mellifera caucasica L.) (X+x, Краснодарский край, 2014)
Линия Морфометрический признак
ДХ, мм I Ш3Т, мм | КИ
I 6,88±0,006a, ь c d 4^0,0^ ь c d 0,502±0,008
II 6,94±0,005ь 4,67±0,010ь, e f 0,509±0,007
III 6^0,004^ i 4^0,008^ f g 0,513±0,007
IV 6,94±0,004^ i 4,61±0,009<1, e, h 0,515±0,008
V 6,93±0,005a 4^0,0^ g, h 0,524±0,008
Примечание. ДХ — длина хоботка, Ш3Т — ширина 3-го тергита, КИ — кубитальный индекс. Описание линий см. в разделе «Методика».
а, ь, c, d, e, f, g, h p < 0,001, i p < 0,01 (достоверность различий между значениями, обозначенным одной и той же буквой).
Результаты двухфакторного иерархического дисперсионного анализа подтвердили наличие достоверных различий по признакам ДХ и Ш3Т как между линиями, так и между семьями в пределах каждой из линий.
Формирование кластеров по линейному принципу наблюдалось для трех из пяти исследованных линий (рис. 1). Только семьи линии I распре-
779
делились внутри одного обособленного 1-го кластера, однако они характеризовались наибольшей изменчивостью по показателю ДХ (Сдх = 1,11 %) по сравнению с другими линиями (Сдх = 0,69-0,86 %) при средней гетерогенности по Ш3Т (СкщэТ = 2,49 %). Семьи линий III и IV формировали два консолидированных кластера (соответственно Сдх 0,74 и 0,69 % и Сшзт 2,15 и 2,33 %), которые частично перекрывались между собой. У пчел линий II и V сильно варьировали оба признака (Сдх — 0,79 и 0,86 %, Сшзт — 2,66 и 2,52 %). Достоверные различия (р < 0,05) между семьями внутри линий наблюдались по признаку ДХ в линии II (Mst = 0,1359) и признаку Ш3Т в линиях II (Mst = 0,2869) и V (Mst = 0,1525). Межлинейные различия были достоверны (р < 0,05) для показателей ДХ (Mst = 0,1453) и Ш3Т (Mst = 0,1097), в то время как по КИ достоверных различий мы не выявили.
В целом следует отметить большую гетерогенность в отношении изучаемых признаков у пчел линий II и V и, напротив, большую консолидацию у пчел линий III и IV. Пчелы линии I существенно отличались от остальных по обоим показателям, однако характеризовались значительными межсемейными различиями по ДХ.
Анализ МС профилей выявил сходные тенденции в определении степени внутрисемейной и межсемейной
изменчивости (табл. 2).
В линии I (единственной из исследованных) наблюдался избыток гетерозигот (отрицательное значение Fis), что можно рассматривать как
Рис. 1. Распределение 25 семей медоносной пчелы серой горной кавказской породы (Apis meUifsra caucasica L.) в пространстве двух измерений по длине хоботка и ширине 3-го тергита: I-V — изучаемые линии медоносной пчелы; 1-3 — формирующиеся кластеры (Краснодарский край, 2014). Описание линий см. в разделе «Методика».
2. Степень внутрисемейной и межсемейной (внутри линий) генетической изменчивости у медоносной пчелы серой горной кавказской породы (Apis mellifera caucasica L.) по семи локусам микросателлитов (Краснодарский край, 2014)
Линия Показатели внутрисемейной и межсемейной изменчивости
Fis Fit Fst (freq) Rst (AMOVA)
I -0,048 0,052 0,124 0,405
II 0,128 0,189 0,104 0,300
III 0,090 0,162 0,096 0,266
IV 0,062 0,143 0,095 0,141
V 0,077 0,169 0,102 0,462
I-V 0,132 0,162 0,040 0,094
Примечание. Расчеты проводили внутри каждой из линий с использованием семьи в качестве субпопуляции, а также для выборки в целом (I-V) с использованием линии в качестве субпопуляции. Описание линий и локусов микросателлитов см. в разделе «Методика».
указание на высокую гетерогенность. Пчелы этой линии характеризовались минимальной индивидуальной изменчивостью по отношению к выборке в целом (значение Fit было в 2,75-3,63 раза ниже по сравнению с линиями II-V) и максимальной межсемейной изменчивостью. У пчел линии II-V имелся дефицит гетерозигот (положительные значения индекса Fis), отмечалась относительно высокая индивидуальная изменчивость по
780
сравнению с выборкой в целом и более низкая по отношению к линии I межсемейная изменчивость. Наименьшие межсемейные различия, оцененные по показателям Fst и R^ (AMOVA), выявили в линиях III и IV.
Результаты PCA-анализа показали наибольшее сходство семей внутри линии IV, в то время как линии I и II характеризуются наибольшими межсемейными различиями (рис. 2).
Анализ дифференциации изучаемых линий выявил различия в структуре генеалогического дерева, построенного для семи локусов МС и трех морфометрических признаков (рис. 3). Дендрограмма, построенная по МС, отражала географическую удаленность изучаемых линий (см. рис. 3, А). Формирование двух главных кластеров могло быть следствием географической изолированности линий I, II и III, расположенных в окрестностях Красной Поляны (первая пасека находится на южном склоне хребта Ачишхо высотой 2391 м, вторая и третья расположены в ущелье реки Мзымта), от линии IV (пасека расположена на расстоянии 15-25 км от трех перечисленных и отделена от них хребтом Алек высотой 1000 м) и линии V (находится в верховьях реки Псоу, удалена от первой, второй и третьей пасек на 10-16 км и отделена от них хребтом Аибга высотой 2460 м). Существенные генетически различия между линиями IV и V могли быть обусловлены их географической изолированностью. Пасеки, на которых разводят эти линии, расположены на расстоянии около 30 км друг от друга и отделены хребтом Алек, который протянулся на 20 км между реками Мзымта и Сочи на расстоянии 15 км от берега Черного моря, и хребтом Аибга, протянувшимся с северо-запада на юго-восток на 23 км и разделяющим среднюю часть долины Мзымты и верховья Псоу.
Рис. 3. Дендрограмма подобия пяти линий медоносной пчелы серой горной кавказской породы (Apis mellifera caucasica L.), рассчитанная для семи локусов микросателлитов (А) и трех морфометрических признаков (Б) (Краснодарский край, 2014). Попарные генетические дистанции рассчитаны по M. Nei (32). В обоих случаях применен метод кластеризации UPGMA. Описание линий см. в разделе «Методика».
Структура генеалогического дерева, построенного по морфометри-
781
□ □ о> >>• о
□ о + о □ + • + • + +
Координата 1
Рис. 2. Распределение центроидов семей медоносной пчелы серой горной кавказской породы (Apis mellifeia caucasica L.) в пространстве первых двух главных координат, рассчитанных на основе матрицы Fst при парном сравнении для семи локусов микросателлитов: о — линия I, □ — ЛИНИЯ II, + — ЛИНИЯ III, А — линия IV, • — линия V (Краснодарский край, 2014). Описание линий см. в разделе «Методика». Расчеты проводили для выборки в целом с использованием семьи в качестве субпопуляции (число субпопуляций = 25).
ческим признакам, не отражала ни географическую общность (различия) происхождения, ни сходство (различия) в хозяйственно полезных признаках изучаемых линий.
Таким образом, нами обнаружены сходные тенденции при оценке внутрисемейной и внутрилинейной (между семьями) изменчивости на основе морфометрических показателей и микросателлитов. Однако данные о степени дифференциации и характере генеалогических связей между линиями, полученные с использованием двух систем анализа, существенно различались. Возможно, подобные расхождения обусловлены тем, что большинство различий по традиционным экстерьерным признакам (морфологические промеры) укладываются в пределы сезонных колебаний и малопригодны для идентификации локальных популяции. Кроме того, они в большей степени подвержены давлению отбора факторов внешней среды, не связанных с воздействием направленного искусственного отбора. Микросателлиты в большинстве случаев селекционно нейтральны и отражают более отдаленные эволюционные события. Следовательно, для получения наиболее полной информации о состоянии и динамике популяций необходимо учитывать результаты, полученные с использованием обоих методов.
1ФГБНУ Всероссийский НИИ животноводства Поступила в редакцию
им. академика Л. К. Эрнста, 14 августа 2015 года
142132 Россия, Московская обл., Подольский р-н, пос. Дубровицы, e-mail: [email protected];
Николаевский национальный аграрный университет,
54020 Украина, г. Николаев, ул. Парижской коммуны, 9, e-mail: [email protected];
3ФГБНУ Краснополянская опытная станция пчеловодства,
354393 Россия, г. Сочи-А, пос. Молдовка, ул. Пчеловодов, 4, e-mail: [email protected]
Sel’skokhozyaistvennaya biologiya [Agricultural Biology], 2015, V. 50, № 6, pp. 776-784
MORPHOMETRIC AND MOLECULAR GENETIC DIFFERENTIATION OF Apis mellifern caucasica L. HONEY BEE LINES REARED IN SOCHI REGION
M.S. Fomarn* 1, A.S. Kramarenko12, S.V. Svistunov3, E.M. Lyubimov3,
S.S. Sokol’skii3, N.A. Zinovieva1
1L.K. Ernst All-Russian Research Institute of Animal Husbandry, Federal Agency of Scientific Organizations, pos. Dubrovitsy, Podolsk Region, Moscow Province, 142132 Russia, e-mail [email protected];
2Mykolayiv State Agrarian University, vul. Paris Commune, 9, Mykolayiv, 54020 Ukraine, e-mail [email protected]; 3Krasnopolyanskii Experimental Station of Beekeeping, Federal Agency of Scientific Organizations, 4, ul. Pchelo-vodov, pos. Moldovka, Sochi-A, 354393 Russia, e-mail [email protected] Acknowledgements:
The equipment of Bioresources and Bioengineering Center of L.K. Ernst All-Russian Research Institute of Animal Husbandry was used.
Supported by the Russian Science Foundation, project 14-36-00039
Received August 14, 2015 doi: 10.15389/agrobiology.2015.6.776eng
Abstract
Creating specialized lines is one of the techniques of genetic improvement and conservation of breeds and populations of the honeybee. The aim of this study was a comparative assessment of the diversity and differentiation degree of A. m. caucasica lines based on morphometric analysis and microsatellites (MS). Material for this study was the worker bees of Gray Mountain Caucasian breed (lines I-V, n = 728) which were selected in five apiaries in the Greater Sochi of Krasnodar krai. Morphometric analysis included such measurements as the length of the proboscis (LP, mm), width of third tergite (W3T, mm) and cubital index (CI). Molecular genetic studies were based on seven MS loci (A024, A88, A113, AP043, HB-C16-05, HB-THE-03, HB-C16-01). The level of variation between families for morphometric parameters was determined by two-way hierarchical analysis of variance. Genetic differences between families for MS were estimated by paired comparison of Fst values. Fst matrix was used for PCA-analysis. To determine the quantitative estimation of
782
variation between families within lines we calculated Fst, Rst (AMOVA). The degree of line’s differentiation for morphometric characters was evaluated by calculating the Euclidean distances. The obtained values were used to construct a dendrogram of similarity by a single bond (single linkage) of the hierarchical clustering algorithm. The differentiation of lines for MS was based on calculating the values of Nei genetic distances. Similarity dendrogram was constructed using the method of UP-GMA. We performed summary statistic using the software STATISTICA, GenAlEx (v. 6.5.1), PAST (v. 3.03). Morphometric analysis showed the presence of significant differences between the lines for LP and W3T whereas there was no difference between the lines for CI. The greater heterogeneity concerning studied traits in the lines II and V was revealed, and on the contrary, there was more consolidation in the lines III and IV. Bees of line I differed significantly from the rest of the lines on both traits, but they were characterized by significant differences between families in LP. Analysis of MS profiles showed similar trends in assessing the level of intra- and interfamily variability. We observed an excess of heterozygotes in the line I (Fjs = -0.048), which can be considered as an indication of the high heterogeneity. Bees of this line were characterized by a minimal individual (Fjt = 0.052) and the maximal interfamily variability (Fst = 0,124). Lines II-V were characterized by a deficiency of heterozygotes (Fjs = 0.062-0.128), a relatively higher individual variability (Fjt = 0.143-
0.189) and lower values of interfamily variability comparing to line I (Fst = 0.095-0.104). The lowest interfamily differences were observed in the lines III and IV (Fst = 0.096 and 0.095, respectively). Analysis of the differentiation of the studied lines for morphometric characteristics and MS revealed differences in the structure of the family tree. The dendrogram based on MS data is a reflection of the geographical origin of these lines. The structure of the family tree, based on morphometric characters, does not reflect the geographic closeness (differences) of origin or similarities (differences) in the economically useful traits of studied lines. Thus, the results of our studies of the morphometric parameters and MS show similar trends in assessing intra- and interline variation, but there are differences in assessing differentiation of lines using two methods. In the future complex approach will allow to identify not only breeds of bees with high accuracy, but also smaller taxonomic units. It is hoped that the research results in general can be used in breeding work to restore the purity of the honeybee breeds.
Keywords: honeybee, Apis mellifera caucasica L., DNA marker, morphometry, microsatellites, biodiversity.
REFERENCES
1. Brandorf A.Z. Pchelovodstvo, 2012, 7: 14-16.
2. Ruttner F. Biogeography and taxonomy of honey bees. Germany, Berlin, Springer, 1988.
3. Bilash G.D., Krivtsov N.I. Pchelovodstvo, 1984, 2: 12-13.
4. Alpatov V.V. Porody medonosnoi pchely [Breeds of honey bees]. Moscow, 1948.
5. Sokol'skii S.S. Materialy Mezhdunarodnoi nauchnoi konferentsii «Sovremennye napravlenii nauchno-tekhnicheskogo progressa v pchelovodstve» [Proc. Int. Conf. «Recent progress in honeybee rearing»]. Rybnoe, 2007: 80-82.
6. Sokol'skii S.S., Savushkina L.H. Pchelovodstvo, 2006, 2: 17-18.
7. Ruttner F., Tassencourt L., Louveaux J. Biometrical-statistical analysis of the geographic variability of Apis mellifera L. Apidologie, 1978, 9: 363-381 (doi: 10.1051/apido:19780408).
8. Ruttner F. Tekhnika razvedeniya i selektsiorrnyi otbor pchel [Breeding and selection of honeybees]. Moscow, 2006.
9. Oleksa A., Tofilski A. Wing geometric morphometrics and microsatellite analysis provide similar discrimination of honey bee subspecies. Apidologie, 2015, 46: 49-60 (doi: 10.1007/s13592-014-0300-7).
10. Rahimi A., Mirmoayedi A., Kahrizi D., Abdolshahi R., Kazemi E., Ya-ri K. Microsatellite genetic diversity of Apis mellifera meda skorikov. Mol. Biol. Rep., 2014, 41(12): 7755-7761 (doi: 10.1007/s11033-014-3667-7).
11. Ivanova E., Bienkowska M., Panasiuk B., Wilde J., Staykova T., Stoy-ano v I. Allozyme variability in populations of A. mellifera mellifera (Linnaeus 1758.), A. m. carnica (Pollman, 1879) and A. m. caucasica (Gorbachev, 1916) from Poland. Acta Zool. Bulg. Suppl, 2012, 4: 81-88.
12. Kalashnikov A.E., Maslennikov I.V., Kolbina L.M., Udina I.G.
Sel’skokhozyaistvennaya biologiya [Agricultural Biology], 2013, 4: 88-92 (doi:
10.15389/agrobiology.2013.4.88rus, doi: 10.15389/agrobiology.2013.4.88eng).
13. Ilyasov R.A., Kutuev I.A., Petukhov A.V., Poskryakov A.V., Nikolenko A.G. Phylogenetic relationships of dark European honeybees Apis mellifera mellifera L. from the Russian Ural and West European populations. J. Apic. Sci., 2011, 55: 67-76.
14. Gruber K., Schyning C., Otte M., Ki nut hi a W., Hasselmann M. Distinct subspecies or phenotypic plasticity? Genetic and morphological differentiation of mountain honey bees in East Africa. Ecol. Evol, 2013, 3(10): 3204-3218 (doi: 10.1002/ece3.711).
783
15. Whitfield C.W., Behura S.K., Berlocher S.H., Clark A.G., Johnston J.S., Sheppard W.S., Smith D.R., Suarez A.V., Weaver D., Tsutsui N.D. Thrice out of Africa: Ancient and recent expansions of the honey bee Apis mellifera. Science, 2006, 314: 642-645 (doi: 10.m6/science.1132772).
16. Harpur B.A., Kenta C.F., Molodtsova D., Lebonb J.M.D., Alqarnic A.S., O w a y s s c A.A., Zayed A. Population genomics of the honey bee reveals strong signatures of positive selection on worker traits. PNAS USA, 2013, 111(7): 2614-2619 (doi: 10.1073/pnas.1315506111).
17. Jensen A.B., Palmer K.A., Boomsma J.J., Pedersen B.V. Varying degrees of Apis mellifera ligustica introgression in protected populations of the black honeybee, Apis mellifera mellifera, in northwest Europe. Mol. Ecol., 2005, 14: 93-106 (doi: 10.1111/j. 1365-294X.2004.02399.x).
18. Bodur C., Kence M., Kence A. Genetic structure of honeybee, Apis mellifera L. (Hy-menoptera: Apidae) populations of Turkey inferred from microsatellite analysis. J. Apic. Res., 2007, 46: 50-56.
19. Krivtsov N.I., Goryacheva I.I., Udina I.G., Borodachev A.V., Monakhova M.A. Sel'skokhozyaistveimaya biologiya [Agricultural Biology], 2010, 6: 26-29.
20. Zinov'eva N.A., Krivtsov N.I., Fornara M.S., Gladyr' E.A., Borodachev A.V., Lebedev V.I. Sel'skokhozyaistveimaya biologiya [Agricultural Biology], 2011, 6: 75-79.
21. Krivtsov N.I., Gladyr' E.A., Volkova V.V., Fornara M.S., Lebedev V.I., Zinov'eva N.A. Problemy biologii produktivrykh zhivotrykh, 2011, 1: 41-45.
22. O s t r o v e rk h o v a N.V., Konusova O.L., Kucher A.N., Pogorelov Yu.L., Belykh E.A., Vorotov A.A. Biologiya, 2013, 1(21): 161-172.
23. Garnery L., Franck P., Baudry E., Vautrin D., Cornuet J.M., Solig-nac M. Genetic diversity of the west European honey bee (Apis mellifera mellifera and A. m. iberica). II. Microsatellite loci. Genet. Sel. Evol., 1998, 30: 49-74 (doi: 10.1186/1297-9686-30-S1-S49).
24. Pritchard J.K., Stephens M., Donnelly P. Inference of population structure using multilocus genotype data. Genetics, 2000, 155: 945-959.
25. Oleksa A., Chybicki I.J., Tofilski A., Burczyk J. Nuclear and mitochondrial patterns of introgression into native dark bees (Apis mellifera mellifera) in Poland. J. Apic. Res., 2011, 50: 116-129 (doi: 10.3896/IBRA.1.50.2.03).
26. Alburaki M., Bertrand B., Legout H., Moulin S., Alburaki A., Sheppard W.S., Garnery L. A fifth major genetic group among honeybees revealed in Syria. BMC Genetics, 2013, 14: 117 (doi: 10.1186/1471-2156-14-117).
27. Mu n oz I., Dall’Olio R., Lodesani M., De La R й a P. Population genetic structure of coastal Croatian honeybees (Apis mellifera carnica). Apidologie, 2009, 40: 617-626 (doi: 10.1051/apido/2009041).
28. Zinov'eva N.A., Popov A.N., Ernst L.K., Marzanov N.S., Bochkarev V.V., Strekozov N.I., Brem G. Metodicheskie rekomendatsii po ispol’zovaniyu me-toda polimeraznoi tsepnoi reaktsii v zhivotnovodstve [Manual of PCR analysis technique in animal husbandry]. Dubrovitsy, 1998.
29. K h a l a f y a n A.A. STATISTICA 6. Statisticheskii analiz dannykh [STATISTICA 6. Statistical analysis]. Moscow, 2007.
30. Kramarenko S.S. Vestnik zoologii, 2009, 43(1): 51-58.
31. Weir B.S., Cocke rh am C.C. Estimating F-statistics for the analysis of population structure. Evolution, 1984, 38: 1358-1370 (doi: 10.2307/2408641).
32. Nei M., Koehn R. Genetic polymorphism and the role of mutation in evolution. In: Evolution Genes and Proteins. Sinauer, Sunderland, MA, 1983: 165-190.
33. Slatkin M. A measure of population subdivision based on microsatellite allele frequencies. Genetics, 1995, 139: 457-462.
34. Peakall R., Smouse P.E. GenAlEx 6.5: genetic analysis in Excel: population genetic software for teaching and research — an update. Bioinformatics, 2012, 28: 2537-2539 (doi: 10.1093/bioinformatics/bts460).
784
