doi: 10.15389/agrobiology.2022.2.283rus
УДК 638.123:577.2
Данное научное исследование выполнено при поддержке Программы повышения конкурентоспособности ТГУ.
НЕКОТОРЫЕ ПРОБЛЕМЫ ИДЕНТИФИКАЦИИ ПОДВИДОВ МЕДОНОСНОЙ ПЧЕЛЫ И ИХ РЕШЕНИЕ НА ПРИМЕРЕ ИЗУЧЕНИЯ Apis mellifera В СИБИРИ
Н.В. ОСТРОВЕРХОВА1, 2 ✉, О.Л. КОНУСОВА1
При изучении медоносной пчелы Apis mellifera L. применяются как классический морфометрический, так и молекулярно-генетические методы, включая полногеномное секвенирование. С использованием морфометрического анализа тридцать подвидов A. mellifera были распределены по четырем эволюционным линиям, которые соответствовали их географическому происхождению. Исследование полиморфизма мтДНК, например локуса COI-COII, позволило выделить три основные эволюционные линии — А, М и С. Вместе с тем метод имеет ограничение, связанное с материнским наследованием митохондриального генома. Поскольку у медоносной пчелы нет половых хромосом, информацию о наследовании по линии трутней (а также по линии матки) можно получить только на основании анализа аутосомных локусов, таких как SNP (single nucleotide poly-morphism) и SSR (simple sequence repeats) маркеры. В настоящей работе впервые проведена оценка информативности морфометрического и молекулярно-генетического методов для идентификации подвидов A. mellifera, обитающих на пасеках Сибири, и показано, что анализ изменчивости основных параметров крыла (кубитальный индекс, гантельный индекс, дискоидальное смещение) и локуса COI-COII мтДНК позволяет точно определить происхождение пчелиной семьи. Впервые изучено генетическое разнообразие среднерусской породы сибирских популяций с использованием микросателлитных локусов. Выявлены породо- и экоспецифичные диагностические аллели, позволяющие дифференцировать среднерусскую породу и ее экотипы, а также подвиды южного происхождения (карпатскую породу, карнику). Цель работы — оценить информативность морфометрического и молекулярно-генетических методов для идентификации подвидов A. mellifera при изучении медоносной пчелы в Сибири. Объектом исследования были медоносные пчелы, полученные с 92 пасек 69 населенных пунктов, расположенных в пяти регионах Сибири: в Томской и Кемеровской областях, в Красноярском крае, Алтайском крае и Республике Алтай. Чтобы изучить породный состав медоносных пчел на пасеках Сибири, первоначально рабочие особи из 414 пчелиных семей были исследованы с использованием морфометрического метода и анализа вариабельности локуса COI-COII мтДНК. На втором этапе выявленные пчелиные семьи среднерусской породы A. m. mellifera были изучены с использованием комплекса SSR-маркеров. Исследована ва-риабельность 31 микросателлитного локуса. С целью поиска уникальных или специфичных для разных подвидов медоносной пчелы SSR-маркеров также изучено генетическое разнообразие пчел двух подвидов A. m. carpathicaи A. m. carnica, имеющих южное происхождение (группа сравнения). Расчет популяционно-генетических параметров (частота аллелей, показатели наблюдаемой и ожидаемой гетерозиготности) проводили с использованием программы GenAlEx 6.5 (https://biology-assets.anu.edu.au/GenAlEx/). Оценка степени интрогрессии генов эволюционной линии С в линию М была выполнена на основе данных полиморфизма микросателлитных локусов с использованием программы STRUCTURE 2.3.4 (https://web.stanford.edu/group/pritchardlab/home.html). Выявлено, что три показателя крыла, а именно кубитальный индекс, гантельный индекс и дискоидальное смещение, совместно с данными по полиморфизму мтДНК необходимы и достаточны для дифференциации подвидов A. mellifera, причем показатель «дискоидальное смещение» одним из первых отклоняется от стандартных значений породы в процессе гибридизации пчел. В результате микросателлитного анализа выявлены локусы, которые дифференцируют как подвиды разных эволюционных линий (линии М и С), так и различные экотипы среднерусской пчелы. В качестве диагностических ДНК-маркеров могут быть рассмотрены локусы А043, Аp081, Аp049, AT139, A113, mrjp3, состав и частота преобладающих аллелей которых (породоспецифичные аллели) различались у подвидов A. m. mellifera (линия М) и двух подвидов южного происхождения (карпатская порода A. m. carpathica и карника A. m. carnica, линия С). У пчел среднерусской породы A. m. mellifera аллели размером 128 п.н. локуса А043, 124 п.н. локуса Ap081, 127 п.н. локуса Ap049, 190 п.н. локуса AT139, 218 п.н. локуса A113 и 529 п.н. — mrjp3 регистрировали с высокой частотой (0,54-0,99), тогда как у пчел южного происхождения A. m. carpathica и A. m. carnicaэти аллелиоказались более редкими (0,01-0,27). Микросателлитный локус А008 рассматривается как наиболее перспективный молекулярный маркер для дифференциации различных экотипов темной лесной пчелы A. m. mellifera, обитающих на территории Сибири, Урала и Европы (выявлены экоспецифичные аллели). На основании изучения генетического разнообразия медоносных пчел, обитающих на территории Сибири, с использованием микросателлитных локусов разработана диагностическая панель молекулярных маркеров для дифференциации подвидов и экотипов медоносной пчелы, принадлежащих к эволюционным линиям М и С (A. m. mellifera, A. m. carpathica и A. m. carnica).
Ключевые слова: Apismellifera, медоносная пчела, морфометрические признаки, мтДНК, микросателлитные локусы, COI-COII, ДНК-маркеры, Сибирь.
N.V. Ostroverkhova1, 2 ✉, O.L. Konusova1
Studies of the honeybee Apis mellifera L. are carried out by both classical morphometric and molecular methods, including whole genome sequencing. Morphometric analysis has revealed thirty subspecies of A. mellifera for four evolutionary lineages that corresponded to their geographical origin. mtDNA analysis, e.g., for the variability of COI-COII locus (the sequence between the cytochrome oxidase I and cytochrome oxidase II genes), has identified three major evolutionary lineages, A, M, and C. However, this method has a limitation associated with maternal inheritance of the mitochondrial genome. The honeybee does not have sex chromosomes, so information about inheritance in the paternal line (as well as in the maternal line) can only be obtained from the analysis of autosomal loci, such as SNP (single nucleotide polymorphism) and SSR (simple sequence repeats) markers. The present work, for the first time, evaluates the informativity of morphometric and molecular methods for the identification of A. mellifera subspecies inhabiting Siberian apiaries. We have shown that the analysis of the variability of the main parameters of the wing (cubital index, hantel index, and discoidal shift) and the mtDNA COI-COII locus accurately detect the origin of the bee colony. We also studied for the first time the genetic diversity of the A. mellifera mellifera Siberian populations for microsatellite loci. Diagnostic alleles specific for subspecies and ecotypes of the honeybees have been identified to differentiate the A. m. mellifera subspecies and its ecotypes, as well as subspecies of southern origin (Carpathian bee, Carnica). This work aimed to evaluate prospects for morphometric and molecular analysis methods in differentiation of A. mellifera subspecies reared in Siberia. Honeybees from 92 apiaries in 69 settlements located in five regions of Siberia (Tomsk and Kemerovo regions, Krasnoyarsk Territory, Altai Territory, and the Altai Republic) were studied. The first stage was the investigation of worker bees from 414 bee colonies using the morphometric method and analysis of the mtDNA COI-COII locus variability. The second stage was the identification of the A. m. mellifera colonies based on a complex of SSR markers. We examined the variability of 31 microsatellite loci. To search for unique or specific SSR markers for different honeybee subspecies, we also examined the genetic diversity of two southern subspecies, A. m. carpathica and A. m. carnica (a comparison group). Population genetic parameters (allele frequency, observed and expected heterozygosity Ho and He) were calculated using GenAlEx 6.5 software (https://biology-assets.anu.edu.au/GenAlEx/). Introgression of the genes of the evolutionary lineage C into the lineage M was assessed based on the microsatellite loci polymorphism data using the STRUCTURE 2.3.4 program (https://web.stanford.edu/group/pritchardlab/home.html). It is shown that three wing parameters, i.e., cubital index, hantel index, and discoidal shift, together with mtDNA polymorphism analysis data, are necessary and sufficient for differentiation of A. mellifera subspecies. The discoidal shift parameter is one of the first morphometric trait to deviate from the breed standard values in the honeybee hybridization. Microsatellite analysis revealed loci that differentiate both subspecies of different evolutionary lineages (M and C) and different A. m. mellifera ecotypes. Loci A043, Ap081, Ap049, AT139, A113, mrjp3, etc. can be considered as diagnostic (subspecies-specific) loci the composition and frequency of the prevailing alleles of which differ in A. m. mellifera subspecies (lineage M) and two subspecies of southern origin (A. m. carpathica and A. m. carnica, lineage C). In A. m. mellifera honeybees, alleles 128 bp at the A043 locus, 124 bp at the Ap081 locus, 127 bp at the Ap049 locus, 190 bp at the AT139 locus, 218 bp at the A113 locus, and 529 bp at the mrjp3 occur in high frequency (0.54-0.99). In honeybees of southern origin (A. m. carpathica and A. m. carnica), these alleles are rarer (0.01-0.27). The microsatellite locus A008 is the most promising molecular marker to differentiate A. m. mellifera ecotypes from Siberia, the Urals and Europe (eco-specific locus). Based on the genetic diversity of Siberian honeybees for microsatellite loci, a diagnostic panel of molecular markers has been developed to differentiate subspecies and ecotypes of honeybees belonging to the evolutionary M and C lineages (A. m. mellifera, A. m. carpathica, and A. m. carnica).
Keywords: honeybee, Apis mellifera, morphometric signs, molecular genetic methods, mtDNA, microsatellite loci, COI-COII, DNA markers, Siberia.
1ФГАОУ ВО Национальный исследовательский |
Поступила в редакцию |
