越冬西方蜜蜂4个亚种的转录组学差异分析

【目的】为进一步探讨西方蜜蜂Apis mellifera在越冬期基因表达模式的变化规律及抗寒机理,进而更好地指导蜂群越冬管理。【方法】利用转录组学鉴定和比较了西方蜜蜂4个亚种高加索蜂A.m.caucasica、欧洲黑蜂A.m.mellifera、意大利蜜蜂A.m.ligustica和卡尼鄂拉蜂A.m.carnica在越冬期的差异表达基因(differentially expressed genes, DEGs),并进行了GO和KEGG富集分析。利用qRT-PCR验证转录组测序结果的可靠性。【结果】西方蜜蜂4个亚种在越冬期DEGs的表达量变化各具特点,其中意大利蜜蜂在越冬初期到越冬中期DEGs表达量变化活跃,卡尼鄂拉蜂和欧洲黑蜂在越冬中期到越冬末期DEGs表达量变化活跃,高加索蜂在整个越冬期DEGs表达量相对稳定。越冬初期西方蜜蜂4个亚种的DEGs表达模式存在差异,高加索蜂和欧洲黑蜂中表皮蛋白22、卵黄原蛋白、气味结合蛋白14、热休克蛋白70和细胞色素P4509e2基因的表达量比意大利蜜蜂和卡尼鄂拉蜂的高,意大利蜜蜂中几丁质酶、防御素1、视黄醇脱氢酶和细胞色素b5基因的表达量较高,而卡尼...     

p38 MAPK基因在西方蜜蜂越冬期表达模式的研究

[目的]p38 MAPK基因在昆虫的低温响应机制中发挥着重要作用,本研究旨在探究p38 MAPK基因在西方蜜蜂Apis mellifera越冬期内表达规律.[方法]本研究对西方蜜蜂p38 MAPK蛋白序列进行生物信息学分析,并利用荧光定量PCR技术检测该基因在意大利蜜蜂Apis mellifera ligustica、欧洲黑蜂Apis mellifera mellifera、高加索蜂Apis mellifera caucasica和卡尼鄂拉蜂Apis mellifera carnica于不同越冬时期和不同越冬方式下体内的表达量.[结果]西方蜜蜂38 MAPK包含1083 bp的开放阅读框区域,编码360个氨基酸,包含TGY双磷酸化三肽模体序列和磷酸化激活环序列.p38 MAPK在蜜蜂所有组织中均有表达,分别在胸部和腹部处于最高和最低表达丰度.38 MAPK mRNA在不同蜂种不同越冬时期的表达量存在显著差异,4个蜂种于室外越冬时p38 MAPK的表达量均显著高于室内越冬(P＜0.05);随着蜜蜂越冬持续时间的延长,意大利蜜蜂和欧洲黑蜂不同越冬方式下体内p38 MAPK的表达均呈现先上升后下降的表达趋势,高加索蜂在室外与室内越冬时体内该基因的表达量均表现出持续上升的表达趋势,该基因在卡尼鄂拉蜂中的表达趋势与其他3个蜂种呈现一定差异,卡尼鄂拉蜂p38 MAPK在室内越冬方式下表达量保持恒定,但在室外越冬方式下表现出先下降后上升的表达趋势;比较室外越冬情况下4个蜂种于不同越冬时期体内38 MAPK表达量后发现,除在次年的1月外,其他越冬时期4个蜂种体内p38 MAPK的表达量均存在显著差异(P＜0.05).[结论]p38 MAPK基因在西方蜜蜂越冬期内发挥着重要的生理功能,p38 MAPK信号通路可作为蜜蜂抗寒机制研究的候选信号通路.     

西方蜜蜂六个亚种苹果酸脱氢酶Ⅱ基因的遗传差异

研究了西方蜜蜂Apis mellifera 6亚种-浙农大1号意蜂（ZND A. m. ligustica）、东北黑蜂（A. m. ssp.）、卡尼鄂拉蜂（A. m. carnica）、喀尔巴阡蜂（A. m. carpatica）、高加索蜂（A. m. caucasica）和乌克兰蜂（A. m. acervorum）苹果酸脱氢酶Ⅱ的基因型频率、基因频率和杂合纯合度。浙农大1号意蜂、喀尔巴阡蜂和高加索蜂的纯合度较高,但浙农大1号意蜂等位基因c频率最高,喀尔巴阡蜂等位基因b频率最高, 高加索蜂等位基因a频率最高;东北黑蜂、卡尼鄂拉蜂和乌克兰蜂是高度杂合的亚种,但东北黑蜂等位基因a、b、c的频率差异较小,卡尼鄂拉蜂和乌克兰蜂主要存在a、c两个等位基因,b出现频率很小;6亚种的基因型频率、基因频率和杂合纯合度都有极显著差异。这些差异将从遗传和生化角度为西方蜜蜂6个亚种的鉴别提供依据。     

采用微卫星DNA标记分析东北黑蜂小群保种效果

本研究利用16对荧光标记微卫星引物分析国家蜜蜂基因库(吉林)保种群体东北黑蜂的小群保种效果,为进一步保护和利用东北黑蜂遗传资源提供理论依据。分别选取国家蜜蜂基因库(吉林)保种群体和国家级东北黑蜂保护区(原产地)群体东北黑蜂各30头,利用测定优势等位基因频率(Pi)、期望杂合度(He)、多态信息含量(PIC)和群体内近交系数(F_(is))比较基因库保种群体与原产地群体的遗传差异。两群体间东北黑蜂的Pi、He、PIC和F_(is)均没有显著差异(P0.05)。国家蜜蜂基因库(吉林)对东北黑蜂群体保种效果良好,异地小群保种的方法对东北黑蜂具有潜在可行性。     

引入西方蜜蜂对中蜂的危害及生态影响

作者阐述自1896年中国引进西方蜜蜂Apis melliferaL.的优良品种如意大利蜂Apis mellifera ligusticaSpinola和喀尼阿兰蜂Apis mellifera Carnica Pollmann以来,使当地的东方蜜蜂Apis cerana F.受到严重危害,其分布区域缩小75%以上,种群数量减少80%以上。使山林植物授粉总量减少,导致植物多样性减少。文中提出: 建立自然保护区保存本地蜜蜂遗传特性,和采用基因转移等技术,培育具有西方蜜蜂优良生产性能的中蜂新品种,逐步恢复中蜂的种群数量。     

中华蜜蜂mrjp1 cDNA的克隆及其序列分析

构建中华蜜蜂(Apis cerana cerana)8日龄工蜂头部cDNA文库,利用中蜂基因组的mrjp3部分基因片段作为杂交探针,采用DIG标记筛选cDNA文库,获得mrjps阳性克隆120个;对阳性克隆进行PCR扩增和测序,通过NCBI的BLAST序列比对,获得12个与印度蜂(Apis cerana india)、西方蜜蜂(Apis mellifera L.)mrjp1基因同源的中蜂mtjp1 cDNA片段,并进一步对中华蜜蜂mrjp1的cDNA全序列进行测定和分析。序列比对分析表明,东方蜜蜂(Apis cerana)与西方蜜蜂mrjp1的cDNA序列相似性为93．78％,中华蜜蜂与印度蜂的相似性高达99．36％,这一结果从分子水平证实中华蜜蜂与印度蜂有较近的共同祖先,而东方蜜蜂与西方蜜蜂的亲缘关系较远。     

中国意大利蜜蜂微卫星遗传多态性

为了分析我国不同生态条件下意大利蜜蜂Apis mellifera ligustica（简称意蜂）的基因特征,对品种间基因交流或基因渗入现象进行监测,完成对我国主要地区意蜂多样性水平的评价。本研究利用10对微卫星（SSR）引物分析了来自我国黑龙江等14个省市自治区31个采集点的意蜂样品的遗传多态性,并与来自德国的喀尼阿兰蜂A. m. carnica（简称喀蜂）和来自美国的意蜂A. m. ligustica进行比对。通过NTSYS-pc2.1和popGene32软件对SSR数据进行多态性分析。结果表明：10对SSR引物中8对表现出较高的遗传稳定性,可作为意蜂遗传多态性标记。遗传一致度在0.61以上时,浙江、四川、甘肃和云南地区的样品出现较严重分化,其他地区的意蜂样品则保持了美国原种意蜂的典型基因特征。东北和新疆地区的样品与德国喀蜂血统亲缘关系较近（遗传距离在0.18149以内）,分析主要原因可能是该地区的意蜂样品杂有黑蜂A. m. mellifera血统。据此提出应在保护各地意蜂原种特性及种群多态性的基础上,科学有效地进行高产品种的繁育。     

基于线粒体和核基因序列的蜜蜂属系统发育分析

文章测定了中国分布的蜜蜂属(Apis)5种蜜蜂22个样本的线粒体基因ND2、CO2、16S rRNA以及核基因ITPR的序列,对序列的碱基组成和蜜蜂种间的遗传距离进行了分析。结合下载的蜜蜂属其他4个种的相关序列,采用最大简约法、邻接法和最大似然法重建了蜜蜂属系统发育关系。系统发育分析结果支持蜜蜂属划分为3个类群,即小蜜蜂类群(包括小蜜蜂和黑小蜜蜂)、大蜜蜂类群(包括大蜜蜂和黑大蜜蜂)和穴居蜜蜂类群(西方蜜蜂、东方蜜蜂、沙巴蜂、苏拉威西蜂、绿努蜂),且小蜜蜂类群较早分化。结果还显示,我国海南岛的大蜜蜂和大陆的大蜜蜂之间可能存在较大的遗传分歧。     

卡尼鄂拉蜂锌转运蛋白-7相似蛋白（ZnT-7-like）的基因克隆、 抗体制备及差异表达

【目的】本研究旨在克隆卡尼鄂拉蜂Apis mellifera carnica锌转运蛋白 7相似蛋白（zinc transporter 7-like, ZnT-7-like）基因, 制备ZnT-7-like多克隆抗体, 了解该基因在卡尼鄂拉蜂王浆腺中不同时期的差异表达情况。【方法】运用RT PCR法从卡尼鄂拉蜂王浆腺总RNA中扩增ZnT-7-like基因, 进行生物信息学分析, 为避开跨膜结构的干扰, 选择克隆其部分序列（273 bp）作为多肽免疫序列。将其亚克隆入原核表达载体pGEX-4T-1, 转入大肠杆菌Escherichia coli BL21 （DE3）中诱导表达获得融合蛋白, 然后将融合蛋白纯化后免疫新西兰大白兔制备多克隆抗体, 分别用间接ELISA和Western blot检测抗体的效价和特异性, 最后采用实时荧光定量RT-PCR以及Western blot技术检测该基因在不同日龄成虫王浆腺中的相对表达量。【结果】克隆得到卡尼鄂拉蜂ZnT-7-like基因, 大小为1 065 bp。SDS-PAGE电泳结果显示融合蛋白成功表达; 制备的多克隆抗体效价高达1∶64 000, 且具有很高的特异性。ZnT-7-like在卡尼鄂拉蜂5个日龄成虫中的转录情况存在较大差异, 表现为3日龄成虫中的表达量极显著高于其他日龄（P<0.01）, 12日龄表达量最低, 其他日龄间表达量两两差异显著(P<0.05)或极显著（P<0.01）; ZnT-7-like蛋白表达与转录水平基本一致。【结论】成功克隆了卡尼鄂拉蜂ZnT-7-like基因, 制备了兔抗蜂ZnT-7-like多克隆抗体, 并在转录和翻译两个水平上测定了ZnT-7-like在卡尼鄂拉蜂王浆腺中不同时期的相对表达量。这些结果为深入研究卡尼鄂拉蜂ZnT-7-like基因的功能奠定了基础。     

三江白猪与其他5个品种猪的群体遗传关系的分析

应用100条RAPD引物对三江白猪及其杂交亲本长白猪和东北民猪,以及大约克夏、圣特西和法系大白猪进行亲缘关系的DNA多态分析。其中15条引物的扩增结果具有明显的多态,将其扩增结果进行聚类分析,发现三江白猪与长白猪(遗传距离d=0.090)、与东北民猪(d=0.175)的亲缘关系较近,与大约克夏(d=0.263)、法系大白(d=0.223)、圣特西猪(d=0.580)的亲缘关系较远。     

蜜蜂残翅病毒在中华蜜蜂及其蜂箱奇露尾甲的发生与进化分析

蜜蜂残翅病毒(Deformed wing virus,DWV)广泛存在于我国规模化养殖的蜂群中,具有明显的季节流行特征.近年来,DWV不断向野生授粉昆虫及其它昆虫传播,包括蜜蜂新的寄生性病害蜂箱奇露尾甲Aethina tumida.前期本实验室已鉴定来源于中华蜜蜂Apis cerana与意大利蜜蜂Apis mellifera的蜂箱奇露尾甲均感染了蜜蜂病毒,但未调查来源于同一群的蜜蜂与蜂箱奇露尾甲是否感染了同一病毒,且在进化上有无差异.本研究调查了来源于同一群的中蜂及蜂箱奇露尾甲感染病毒情况,发现只感染了蜜蜂残翅病毒A型(DWV-A),未发现其它病毒.对DWV-A的3个片段进化分析表明,虽然DWV-A同时感染了中华蜜蜂和蜂箱奇露尾甲,但是VP3基因在不同的寄主中呈现明显不同的亲缘关系.这一结果为蜜蜂残翅病毒在蜜蜂及寄生病害中的进化提供了数据,为进一步研究病毒在不同寄主间的传播与进化奠定了基础.     

三种微小RNA在意大利蜜蜂工蜂蛹期发育过程中的表达谱及潜在功能

[目的]意大利蜜蜂Apis mellifera ligustica是生产性能优良的西方蜜蜂Apis mellifera亚种.本研究旨在利用分子生物学手段揭示3种微小RNA(microRNA,miRNA) ame-miR-13b,ame-miR-100和ame-miR-bantam调控蛹期变态发育的分子机理.[方法]通过...     

中华蜜蜂mab-21基因序列分析及表达特征

【目的】探究中华蜜蜂 Apis cerana cerana Male abnomal 21（mab-21）基因在不同发育阶段的表达特性及染螨条件下mab-21基因表达变化规律。【方法】本研究利用RT-PCR方法,克隆了中华蜜蜂mab-21的基因编码区;采用荧光定量PCR方法检测了中华蜜蜂mab-21在不同发育时期（新出房蜂、哺育蜂、守卫蜂及采集蜂）工蜂头部中mRNA的表达量以及接种大蜂螨 Varroa destructor 前后mab-21基因的表达变化。【结果】克隆获得中华蜜蜂mab-21 cDNA,命名为 Accmab 21（GenBank登录号KR000001）。序列分析显示, 该编码区开放阅读框长为1 098 bp,编码365个氨基酸,推测的编码蛋白的相对分子量和等电点分别为41.63 kD和8.53。系统发育分析表明中华蜜蜂mab-21与西方蜜蜂 Apis mellifera mab-21、小蜜蜂Apis florea mab-21和熊蜂Bombus impatiens mab-21聚成一支。该基因在中华蜜蜂工蜂的不同发育时期均有表达,其中哺育蜂阶段显著高于新出房蜂、守卫蜂和采集蜂(P<0.05)。接种大蜂螨后,哺育蜂和守卫蜂中mab-21基因的表达下降显著(P<0.05);而在新出房蜂和采集蜂中表达量变化不显著。【结论】该基因可能与中华蜜蜂抗螨行为相关。     

枣花病对意大利蜜蜂消化酶和解毒酶系活性的影响

[目的]枣花期蜜蜂易发生枣花中毒(又称枣花病),致使大量采集蜂死亡.本研究初步探讨了枣花病对蜜蜂的毒性作用机理.[方法]以正常的意大利蜜蜂Apis mellifera ligustica工蜂为对照,采用分光光度法测定山西临猗、太谷和临县患枣花病的意大利蜜蜂工蜂体内淀粉酶、蛋白酶、蔗糖酶、果胶酶、海藻糖酶、酯酶(este...     

全基因组重测序揭示青藏高原东缘及东南缘东方蜜蜂遗传多样性与适应性进化

【目的】对青藏高原东缘和东南缘及邻近地区东方蜜蜂Apis cerana样本进行群体遗传多样性与适应性进化研究,为进一步解析青藏高原东方蜜蜂的遗传资源多样性、种群扩散规律以及适应高原生境的分子进化机制提供参考。【方法】对青藏高原东缘和东南缘及邻近地区77群东方蜜蜂样本进行全基因组重测序;运用群体遗传学方法,基于群体结构、主成分分析、系统进化树、遗传分化指数、线粒体基因组单倍型以及选择信号分析对东方蜜蜂这77群及GenBank数据库下载的90群的全基因组重测序原始数据进行分析。【结果】这167群东方蜜蜂区分出来自川西高原、藏南高原、滇北高原和川北高原的4个高原型群,分属于两个进化支,平均遗传分化指数(Fst=0.1178)高于非高原型群的(Fst=0.0411)。基于种群间最小遗传距离分析,东南缘的藏南高原群、滇北高原群与滇南群具有更近的遗传关系;东缘的川西高原群、川北高原群分别与川西山地和秦巴群具有更近的遗传关系。结合线粒体基因组单倍型分析,初步推断出青藏高原东缘及东南缘高原型群体祖先单倍型及来源。选择性分析鉴定到了潜在的参与脂肪酸代谢、光转导、温度适应、卵巢发育等信号通路相关的潜在受选基因。在藏南高原和滇北高原群体中发现两个共同受选择基因ISL-1和FOXO,主要参与胰岛素分泌以应对细胞压力,暗示它们在东方蜜蜂适应青藏高原东南缘生境中发挥着重要作用。【结论】青藏高原东方蜜蜂遗传多样性丰富,东南缘的藏南高原和滇北高原群体以及东缘的川西高原和川北高原群体在遗传分化上明显区分,这4个高原群是由邻近地区非高原群扩散后的地理隔离形成了种群分化;初步筛选获得东方蜜蜂高原环境适应性潜在基因。本研究为进一步探析青藏高原东方蜜蜂适应高原生境下的分子进化机制奠定了基础。     

亚致死剂量吡虫啉胁迫对意大利蜜蜂内勤蜂与外勤蜂免疫解毒相关基因表达及免疫解毒酶系活力的影响

摘要: 【目的】吡虫啉(imidaclorprid)是广泛使用的新烟碱类杀虫剂之一。大量研究表明亚致死剂量吡虫啉影响意大利蜜蜂Apis mellifera ligustica(简称“意蜂”)幼虫的发育和成年蜜蜂的采集、学习等行为。本实验旨在探究亚致死剂量吡虫啉对意大利蜜蜂内勤蜂(1日龄成年工蜂)与外勤蜂(21日龄成年工蜂)免疫解毒相关基因表达及免疫解毒酶系活力的影响,进而为蜜蜂健康的维护提供科学依据。【方法】测定饲喂含0.1 ng/μL吡虫啉的50%蔗糖溶液不同时间后意蜂成年工蜂的存活率;利用荧光定量PCR检测饲喂含0.1 ng/μL吡虫啉的50%蔗糖溶液6 d后其体内免疫基因多酚氧化酶基因(PPOA3, GenBank登录号： GB43738), Abaecin类抗菌肽基因(ABA, GenBank登录号： GB18323),葡萄糖脱氢酶基因(GLD, GenBank登录号： GB43007)和解毒基因细胞色素P450基因(CYP450 6a2, GenBank登录号： GB49876)的表达,并采用双抗体一步夹心法酶联免疫吸附试验测定其体内细胞色素P450酶(cytochrome P450, CYP450)含量和多酚氧化酶(polyphenol oxidase,PPO)活力。【结果】1日龄和21日龄意蜂成年工蜂连续饲喂6 d含0.1 ng/μL吡虫啉的50%蔗糖溶液后,其存活率与对照组(饲喂含0.1 ng/μL丙酮的50%蔗糖溶液)无显著差异;连续饲喂9 d含0.1 ng/μL吡虫啉的50%蔗糖溶液后,1日龄意蜂成年工蜂存活率与对照组无显著差异,而21日龄意蜂成年工蜂存活率与对照组有显著差异。1日龄意蜂成年工蜂自由取食含0.1 ng/μL吡虫啉的蔗糖溶液6 d后, PPOA3, CYP450 6a2, ABA和GLD表达水平,细胞色素P450含量以及多酚氧化酶活力与对照组相比均有显著下调趋势;而21日龄意蜂成年工蜂取食该药液6 d后,CYP450 6a2, ABA和GLD表达水平及多酚氧化酶活力与对照组相比均有显著下调趋势,PPOA3表达水平和细胞色素P450含量有显著上调趋势。【结论】亚致死剂量吡虫啉影响意大利蜜蜂内勤蜂与外勤蜂免疫解毒相关基因的表达及免疫解毒酶系活力;吡虫啉短期胁迫对意大利蜜蜂内勤蜂与外勤蜂的存活无显著影响,长期胁迫则会影响意大利蜜蜂内勤蜂与外勤蜂的存活。     

蜜蜂瓦螨敏感卫生行为研究进展

蜜蜂具有很高的生态价值和经济价值,对农业生产帮助巨大。然而,狄斯瓦螨Varroa destructor寄生给西方蜜蜂Apis mellifera蜂群造成重大损失,对蜜蜂健康构成严重威胁,因此,狄斯瓦螨的防治变得尤为紧要。虽然化学防治是防治狄斯瓦螨常用且有效措施,但仍存在许多缺点,如造成蜂产品污染、导致蜂螨产生抗药性等。另一方面,培育抗螨蜂种被证明是可持续的狄斯瓦螨防治方法。瓦螨敏感卫生行为(Varroa sensitive hygiene, VSH)是蜜蜂重要的抗螨性状之一。本文从狄斯瓦螨的生活周期、对蜜蜂的危害、蜜蜂抗螨行为、瓦螨敏感卫生行为调控和遗传育种等方面进行综述,为狄斯瓦螨防治和抗螨蜂种选育提供参考。     

基于线粒体Cyt b基因部分序列的中国东方蜜蜂不同地理种群的系统发育

从线粒体DNA水平上对中国境内东方蜜蜂不同地理种群的系统发育进行研究,为保护和合理利用这一宝贵的蜂种资源提供理论基础。对我国10个省市的21群东方蜜蜂mtDNA Cyt b基因片段进行了扩增、测序,并以意大利蜂、卡尼鄂拉蜂、印尼蜂作为外群进行序列分析。采用Minimum-Spanning Network方法构建系统进化树。结果显示：扩增片段长度429bp;在得到的21条同源序列中,共检测出13个变异位点,其中转换数为9,颠换数为4,并且变异位点大都发生在密码子的第三位,无碱基插入或缺失;聚类分析结果显示西方蜜蜂、印尼蜂和东方蜜蜂各为独立的蜂种;在中国东方蜜蜂群体中,吉林、海南和云南的东方蜜蜂各为一个独立的类群,其它地区的东方蜜蜂为一个类群[动物学报54（6）：1005-1013,2008]。     

中华蜜蜂气味结合蛋白AcerOBP14的表达及配体结合特性分析

【目的】气味结合蛋白(odorant binding proteins, OBPs)在昆虫寄主定位、产卵地选择等行为中发挥着重要作用,明确中华蜜蜂Apis cerana cerana AcerOBP14与配体的结合特性有助于阐明中华蜜蜂嗅觉识别的分子机制。【方法】通过qRT-PCR测定OBP14在20日龄中华蜜蜂成年工蜂、20日龄中华蜜蜂成年雄蜂、中华蜜蜂采粉蜂和意大利蜜蜂Apis mellifera ligustica采粉蜂触角中的表达量。构建原核表达载体pET28a/AcerOBP14,表达并分离纯化重组蛋白AcerOBP14。利用荧光竞争结合实验检测AcerOBP14与37 种气味配体化合物的结合特性。【结果】qRT-PCR分析发现, OBP14在中华蜜蜂采粉蜂触角中的表达量极显著高于20日龄中华蜜蜂成年工蜂和雄蜂以及意大利蜜蜂采粉蜂中的。荧光竞争结合实验表明,AcerOBP14与蜂王信息素、告警信息素、那氏信息素及多种植物挥发物都具有结合能力,其中与β 罗勒烯的结合能力最强,解离常数 Ki=0.297 μmol/L。【结论】AcerOBP14的配体结合谱较宽,暗示其可能参与了中华蜜蜂的多种生理行为反应,且在中华蜜蜂的采粉行为中发挥着重要作用。     

中国东方蜜蜂不同地理种群的RAPD分析

本研究采用RAPD技术对来自中国境内8个省市的28群东方蜜蜂Apis cerana进行了研究,探讨了它们之间的亲缘关系.结果表明:15条随机引物共扩增产生了155条谱带,其中单个弓l物扩增的谱带数为4～16条,多态性条带为122条,多态比率为78.7%;不同种群东方蜜蜂遗传多样性丰富,遗传距离变化较大,最大为0.4197(海南-北京),最小为0.0933(江西-甘肃);聚类结果显示8个地理种群的东方蜜蜂可聚为3个聚类簇,且各地理种群的亲缘关系与空间距离正相关.