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

研究了西方蜜蜂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个亚种的鉴别提供依据。     

意大利蜜蜂工蜂中肠SNP和InDel突变位点的挖掘及分析

本研究利用已获得的意大利蜜蜂Apis mellifera ligustica工蜂中肠的转录组数据对意蜂的单核苷酸多态性（Single nucleotide polymorphism, SNP）和插入缺失（Insertion-Deletion, InDel）突变位点进行挖掘和分析,共鉴定到232 678个SNP位点,其中发生转换和颠换的SNP位点数分别为196 087和36 591个;最丰富的突变类型是G/A,最少的突变类型为T/G;分布在内含子的SNP位点最多,其次为外显子和基因间区;密码子突变类型为同义突变的SNP位点数最多,其次是非同义突变、终止子增加和终止子减少;此外,SNP位点所在基因可注释到 50个 GO条目和351条KEGG通路。共鉴定到38 715个InDel位点,最丰富的突变类型为移码插入,其次是移码缺失;分布InDel位点数较多的基因组区域为内含子和基因间区;另外,InDel位点所在基因可注释到50个功能条目和340条通路。研究结果丰富了西方蜜蜂Apis mellifera的SNP和InDel位点信息,并为开发和利用意蜂的新型分子标记提供基础。     

亚历山大藻属微卫星标记的筛选

采用生物信息学方法,从塔玛亚历山大藻的ESTs数据库中筛选微卫星(Simple Sequence Repeat, SSR)位点,设计SSR引物,利用毛细管电泳对4株塔玛亚历山大藻(Alexandrium tamarense)、2株相关亚历山大藻(Alexandrium affine)和1株链状亚历山大藻(Alexandrium catenella)进行种内、种间遗传多样性分析,获得以下主要结果:(1)从6830条非冗余的ESTs中共查找到222个SSR,平均每18.5kb就有一个SSR,三核苷酸重复在所有的重复类型中所占比例最大,达到48.2%,重复单元CTG/CAG出现频率最高.进一步筛选出12条SSR引物,用于遗传多样性分析.(2)种内遗传多样性水平整体偏低,多态性引物比率仅为41 67%,Nei's基因多样性平均为0.2130.种间遗传多样性水平较高,多态性引物比率为75%,Nei's基因多样性平均水平为0 4643;另外,种间遗传分化系数较高,平均水平为0.7051.(3)对7个样品的聚类图分析发现,A. tamarense CCMP116 与A. affine CCMP112遗传距离极为接近;A. tamarense ATHK9301、CCMP1598和ATDH01相互聚在一起;A. catenella ACDH01与A. tamarense分属两枝.以上结果表明,SSR标记在亚历山大藻种内遗传多样性分析方面是一种非常有用的工具.     

蜜蜂表皮蛋白apidermin (apd)基因家族3个新成员的特性鉴定及昆虫APD家族序列特征的分析

Apidermin (APD)蛋白家族是一个新的昆虫结构性表皮蛋白家族。本研究结合生物信息学和RT-PCR扩增, 对意大利蜜蜂Apis mellifera ligustica（简称“意蜂”）的apd-1-like, apd-3-like和中华蜜蜂Apis cerena cerena（简称“中蜂”）的apd-2 等3个新的apd基因的结构特征和表达进行了分析, 并分析了昆虫APD蛋白家族的序列特征。结果显示, 在西方蜜蜂Apis mellifera（简称“西蜂”）中, apd基因家族的6个成员串联排列在基因组序列第4号连锁群上, 它们在A. m. ligustica雄蜂头部中的转录水平差异明显, 且其启动子序列所含顺式元件也不同。中蜂apd-2和意蜂apd-1-like都含有3个外显子和2个内含子, 而意蜂apd-3-like则由4个外显子和3个内含子组成。蛋白序列分析结果显示, 目前已知的10条APD蛋白序列N末端均具有相似的信号肽序列, 其成熟蛋白分子量为6.0~37.0 kD, pI为6.2~10.8。其中西蜂的APD1-3、APD-like和东方蜜蜂Apis cerena的APD-2等5条较短的多肽中疏水氨基酸残基达52%~67%, 且Ala含量最为丰富（占25%~34%）; 而丽蝇蛹集金小蜂Nasonia vitripennis的APD 1-3和西蜂APD-1-like, APD-3-like等另外5条APD多肽富含Gly（21%~30%）, 其序列中疏水氨基酸残基含量为35%~41%。多肽序列多重比对和系统进化分析结果显示, APD家族可划分为2个亚家族。亚家族Ⅰ含有西蜂APD 1-3和东方蜜蜂APD-2等4条较短的多肽序列, 其N末端为一个长33 aa的保守基序; 亚家族Ⅱ由另外6条相对较长的多肽序列组成, 其N末端保守基序长50 aa, C末端保守基序长16 aa。本文所描述的APD蛋白家族序列特征有助于以后从其他昆虫中鉴定新的apd基因。     

用SSR和AFLP技术分析花生抗青枯病种质遗传多样性的比较

由Ralstonia solanacearum E.F.Smith引起的青枯病是若干亚洲和非洲国家花生生产的重要限制因子,利用抗病品种是防治这一病害最好的措施。虽然一大批抗青枯病花生种质资源材料已被鉴定出来,但对其遗传多样性没有足够的研究,限制了在育种中的有效利用。本研究以31份对青枯病具有不同抗性的栽培种花生种质为材料,通过简单序列重复(SSR)和扩增片段长度多态性(AFLP)技术分析了它们的遗传多样性。通过78对SSR引物和126对AFLP引物的鉴定,筛选出能显示抗青枯病种质多态性的SSR引物29对和AFLP引物32对。所选用的29对多态性SSR引物共扩增91条多态性带,平均每对引物扩增3.14条多态性带;32对多态性AFLP引物共扩增72条多态性带,平均扩增2.25条多态性带。在所筛选引物中,4对SSR引物(14H06,7G02,3A8,16C6)和1对AFLP引物(P1M62)检测花生多态性的效果优于其他引物。SSR分析获得的31个花生种质的遗传距离为0.12-0.94,平均为0.53,而AFLP分析获得的遗传距离为0.06～0.57,平均为0.25,基于SSR分析的遗传距离大于基于AFLP分析的遗传距离,疏枝亚种组的遗传分化相对大于密枝亚种组。基于两种分析方法所获得的聚类结果基本一致,但SSR数据聚类结果与栽培种花生的形态分类系统更为吻合。根据分析结果,对构建青枯病抗性遗传图谱群体的核心亲本和抗性育种策略提出了建议。     

越冬西方蜜蜂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基因的表达量较高,而卡尼...     

中国落叶松-杨栅锈菌遗传多样性研究

采用SSR分子标记技术对采自全国19个地区、不同年份的落叶松-杨栅锈菌Melampsora larici-populina（简称MLP）36个单孢子堆菌系的遗传多样性和种群遗传结构进行了研究。用5对SSR多态性引物共检测到62个观察等位基因（Na）,有效等位基因数（Ne）为5.42–9.82,平均有效等位基因数为7.05。供试MLP样本在5个SSR位点上的多态性信息量（PIC）变化范围为0.64–0.89,平均值0.79。平均观察杂合度（Ho）、平均期望杂合度（He）分别为0.36和0.86,不同位点的S     

中国35个苦荞审定品种EST-SSR指纹图谱构建与遗传多样性分析

核心引物对种质资源遗传多样性分析、品种鉴定、指纹图谱构建等研究具有重要价值。本研究以35个苦荞（Fagopyrum tataricum（L.） Gaertn）审定品种为材料,从91对苦荞EST-SSR引物中筛选出50对多态性引物。综合考虑引物多态性信息量（PIC）大小、鉴别力（DP）,筛选出等位变异位点数在2~4,PIC值在0.60~0.78之间的6对引物（SSR9007、SSR6873、SSR7642、SSR2234、SSR6789、SSR68216）构建了供试品种的分子指纹图谱。遗传多样性聚类分析结果表明,供试品种的相似系数为0.50~0.99。当遗传相似系数为0.60时,可将供试品种分为4大类群,其中54.3%的供试品种被聚为一类,表明苦荞审定品种遗传组成差异较小,遗传基础狭窄。聚类结果表明各类群间没有明显的地域分布趋势,但能较好的反映供试品种间的亲缘关系。     

意大利蜜蜂工蜂中肠响应Nosema ceranae胁迫的高表达基因分析

东方蜜蜂微孢子虫Nosema ceranae是专性寄生蜜蜂中肠上皮细胞的单细胞真菌,对意大利蜜蜂(意蜂)具有较强的侵染性。本研究利用RNA-seq技术对正常10 d (Apis mellifera ligustica control group,Am CK)和N. ceranae胁迫10 d的意蜂工蜂中肠(Apis mellifera ligustica treatment group,Am T)进行测序,共得到160 847 237 100条原始读段,过滤后得到1 066 955 298条有效读段。主成分分析结果显示Am CK与Am T测序样品的组内重复性较好,组间的基因表达模式差异明显。GO分类结果显示,Am CK与Am T的前100位高表达基因(HEGs)分别分布于32和33个GO terms,基因富集数最多的均为代谢进程。KEGG代谢通路(pathway)富集分析结果显示,Am CK的前100位HEGs富集在21个pathways,基因富集数最多的是内吞作用、信号通路和嘌呤代谢; Am T的前100位HEGs富集在26个pathways,基因富集数最多的是内吞作用、RNA转运和蛋白质的内质网加工。前100位HEGs的Venn分析结果显示Am CK与Am T的共有HEGs为87个,特有HEGs分别均为13个。研究结果揭示了N. ceranae胁迫意蜂工蜂中肠过程中的基因表达谱信息,也为解析N. ceranae致病的分子机理提供了有益信息和线索。     

SSR标记在毛木耳种质资源遗传多样性研究中的应用

本研究利用基于毛木耳全基因组开发的SSR标记对27份毛木耳菌株（野生14株、栽培13株）的遗传多样性进行分析。首先随机选取3个菌株（2个野生菌株、1个栽培菌株）的DNA为模板,从144对SSR引物中筛选出扩增条带清晰、稳定性强、多态性丰富的引物24对。24对SSR引物共检测到116个多态性SSR片段,每对引物的多态性片段有3-7个,引物平均检测效率为4.83个,Shannon’s遗传多样性指数范围是0.866-1.885,多态性位点比率100%。供试菌株遗传相似系数范围是0.618-0.971,说明毛木耳种质资源具有丰富的遗传多样性。野生菌株与栽培菌株间平均遗传相似系数分别为0.746、0.779,说明毛木耳野生菌株遗传多样性更为丰富。经聚类分析,在遗传相似系数为0.680时,可将供试菌株分为无色（白色）类群Ⅰ和有色（浅黄色到红褐色）类群Ⅱ。遗传相似系数为0.704时,可将供试菌株中栽培菌株和野生菌株明显区分（14株野生菌株均在类群Ⅱ-2中,13株栽培菌株分别在类群Ⅰ和Ⅱ-1中）。本研究表明基于全基因组的SSR标记能从分子水平上揭示各菌株间的遗传差异,丰富毛木耳遗传多样性的研究手段,并为进一步进行毛木耳的品种选育、遗传学研究等提供有力手段。     

亚致死剂量吡虫啉胁迫对意大利蜜蜂工蜂嗅觉学习行为的影响及其脑部转录组分析

【目的】分析吡虫啉处理对意大利蜜蜂Apis mellifera ligustica嗅觉学习行为及脑部基因转录的影响,为新烟碱类杀虫剂对蜜蜂的负面影响提供依据。【方法】实验室条件下一次性饲喂意大利蜜蜂成年工蜂含有4 ng吡虫啉的50%蔗糖溶液,以饲喂不含吡虫啉的50%蔗糖溶液为对照,通过伸吻反应(proboscis extension response, PER)行为实验测定其对意大利蜜蜂成年工蜂嗅觉学习行为的影响;收集上述测试的意大利蜜蜂工蜂脑部提取RNA,进行RNA-Seq测序和生物信息学分析,利用实时荧光定量PCR挑选6个差异表达基因(DEGs)检测其在脑部的表达量以验证RNA-Seq测序结果。【结果】一次性饲喂意大利蜜蜂成年工蜂含4 ng吡虫啉的50%蔗糖溶液后,意大利蜜蜂的嗅觉学习能力与对照组(饲喂50%蔗糖溶液)相比显著降低。RNA-Seq测序结果显示饲喂意大利蜜蜂含4 ng吡虫啉蔗糖溶液后,处理组与对照组之间共有123个DEGs［校正后的P值(padj)<0.05)］,包括82个下调DEGs和41个上调DEGs。GO聚类分析发现下调DEGs主要富集在S-腺苷甲硫氨酸依赖性甲基转移酶活性、酸性磷酸酶活性、氧化还原酶活性、蛋白质异二聚化活性等,上调DEGs主要富集在跨膜受体活性、分子传感器活性、神经学系统过程等功能条目。KEGG富集分析显示下调DEGs主要富集在核糖体和溶酶体等细胞器、碳代谢和色氨酸代谢等代谢通路及Toll和Imd信号通路,上调DEGs未富集在KEGG代谢通路。荧光定量PCR结果显示测试的6个DEGs相对表达量变化趋势与RNA-Seq测序结果FPKM(fragments per kilobase million)值变化趋势一致,证明RNA-Seq测序结果可靠。【结论】亚致死剂量吡虫啉胁迫意大利蜜蜂成年工蜂显著降低意大利蜜蜂的嗅觉学习能力,影响意大利蜜蜂脑部免疫解毒等相关基因的表达和酶活性及氧化还原等生物代谢过程;亚致死剂量吡虫啉短时胁迫会刺激意大利蜜蜂的嗅觉感觉过程及神经信号传导过程。     

球囊菌侵染对意大利蜜蜂幼虫肠道免疫与解毒相关基因表达及蛋白质、脂质和糖类含量的影响

[目的]球囊菌Ascosphaera apis是一种真菌性病原体,可以侵染意大利蜜蜂Apis mellifera ligustica幼虫,导致蜂群患白垩病.本研究通过球囊菌侵染意大利蜜蜂幼虫,探究对意大利蜜蜂幼虫肠道免疫解毒相关基因表达及蛋白质、脂质和糖类含量的影响,分析球囊菌侵染胁迫下意大利蜜蜂免疫应答反应与营养物质...     

意大利蜜蜂染色体DNA非编码区抗白垩病相关的SNP的筛选与验证

【目的】筛选验证意大利蜜蜂Apis mellifera ligustica染色体DNA非编码区与抗白垩病相关的SNP。【方法】本研究将蜜蜂球囊菌Ascosphaera apis孢子接种于人工饲养的意大利蜜蜂3日龄幼虫,根据是否存在白垩病症状进而筛选出抗病个体和易感个体。基于前期重测序结果中意大利蜜蜂第2和11号染色体DNA非编区与抗白垩病相关的SNP信息,利用PCR测序的方法筛选并验证意大利蜜蜂幼虫第2和11号染色体DNA非编码区与幼虫抗白垩病相关的55个SNP。【结果】发现位于意大利蜜蜂第11号染色体LOC100578413基因5′端的非编码区的SNP(T14570310C)在抗病个体中T等位基因频率高于C等位基因频率,且抗病个体中的T等位基因频率显著高于易感幼虫中的T等位基因频率,表明该SNP位点与抗白垩病相关。该分子标记对抗性个体和易感个体的判断结果与前期筛选的编码区SNP(C2587245T)分子标记的结果一致。【结论】筛选并验证意大利蜜蜂第11号染色体DNA非编码区的SNP(T14570310C)与抗白垩病相关。该位点为抗白垩病分子辅助选育提供新的分子标记,在意大利蜜蜂白垩病早期检测和培育白垩病抗性的蜂种方面具有重要意义。     

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

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

意大利蜜蜂(Apis mellifera ligustica)与中华蜜蜂(Apis cerana ceraca)的生态位比较

为了明确中华蜜蜂和意大利蜜蜂在皖南山区生态适应性,研究两种蜜蜂的食物生态位、时间生态位和空间生态位及其差异,结果是中蜂与意蜂食物(蜜源植物)资源生态位宽度分别是 0.923、0.765,中蜂对蜜源植物采集喜好性差异小,而意蜂差异大,中蜂对意蜂生态位重迭为0.160,中蜂对意蜂生态位相似性为0.755;油菜花期,中蜂与意蜂时间资源生态位宽度分别是0.879、0.801,枇杷花期,分别是0.760、0.677,中蜂与意蜂的空间资源生态位宽度分别是0.797、0.670.中蜂3种生态位宽度均大于意蜂,中蜂三维生态位值是意蜂的1.61倍和1.57倍.表明中蜂在皖南山区生态适应性比意蜂强.     

Hybrid origins of honeybees from italy (Apis mellifera ligustica) and sicily (A. m. sicula)

The genetic variability of honeybee populations Apis mellifera ligustica, in continental Italy, and of A. m. sicula, in Sicily, was investigated using nuclear (microsatellite) and mitochondrial markers. Six populations (236 individual bees) and 17 populations (664 colonies) were, respectively, analysed using eight microsatellite loci and DraI restriction fragment length polymorphism (RFLP) of the cytochrome oxidase I (COI)-cytochrome oxidase II (COII) region. Microsatellite loci globally confirmed the southeastern European heritage of both subspecies (evolutionary branch C). However, A. m. ligustica mitochondrial DNA (mtDNA) appeared to be a composite of the two European (M and C) lineages over most of the Italian peninsula, and only mitotypes from the African (A) lineage were found in A. m. sicula samples. This demonstrates a hybrid origin for both subspecies. For A. m. ligustica, the most widely exported subspecies, this hybrid origin has long been obscured by the fact that in the main area of queen production (from which most of the previous ligustica bee samples originated) the M mitochondrial lineage is absent, whereas it is present almost everywhere else in Italy. This presents a new view of the evolutionary history of European honeybees. For instance, the Iberian peninsula was considered as the unique refuge for the M branch during the quaternary ice periods. Our results show that the Apennine peninsula played a similar role. The differential distribution of nuclear and mitochondrial markers observed in Italy seems to be a general feature of introgressed honeybee populations. Presumably, it stems from the social nature of the species in which both genome compartments are differentially affected by the two (individual and colonial) reproduction levels.     

意大利蜜蜂和小峰熊蜂在温室桃园的传粉行为及其影响因素

意大利蜜蜂Apis mellifera ligustica和小峰熊蜂Bombus hypocrita是我国设施农业中的两种主要授粉昆虫, 二者的传粉方式和传粉效率因不同作物而各有不同。为了选用最理想的蜂种为温室作物授粉, 提高作物的授粉效率, 我们于2008-2010年连续3年在北京的温室桃园进行了意大利蜜蜂和小峰熊蜂的传粉行为及其影响因素的研究。结果表明, 两种蜂的日活动规律不同, 和意大利蜜蜂相比, 小峰熊蜂的活动起点温度低、 日工作时间长、 单花访问时间长, 采集蜂在温室内距离蜂巢不同距离之间扩散比较均匀。而意大利蜜蜂采集蜂的数量随授粉半径的增大而明显减少。温室内环境因子对两种蜂传粉活动的影响基本一致, 温度对两种蜂的传粉活动影响最大, 其次是湿度、 单花花蜜浓度和光照强度, 而单花花蜜体积对蜂活动无明显影响。本研究认为对于温室桃授粉应优先选用授粉性能稳定的小峰熊蜂, 并且适当调节温室内环境条件, 以提高其授粉效率。     

意大利蜜蜂幼虫肠道内球囊菌及其纯培养的高表达基因差异分析

【背景】球囊菌是一种典型的蜜蜂真菌性病原,特异性地侵染蜜蜂幼虫。目前,有关球囊菌在侵染过程中的基因表达规律的信息极为有限。【目的】通过深入分析胁迫意大利蜜蜂(意蜂)幼虫肠道的球囊菌及其纯培养的高表达基因(HEGs)差异,探索球囊菌在胁迫意蜂幼虫肠道后期与胁迫前的基因表达规律。【方法】利用RNA-Seq技术对球囊菌胁迫的意蜂6日龄幼虫肠道(Aam T)及球囊菌纯培养(AaCK)进行深度测序,根据FPKM值筛选得到球囊菌的HEGs,进而通过GO(Gene ontology)及KEGG(Kyoto Encyclopedia of Genes and Genomes)富集分析、Venn分析对上述HEGs进行功能预测和生物学意义挖掘。【结果】AaCK与Aam T的转录组测序共得到105 447 578条原始读段(Raw reads),经过滤得到88 466 344条有效读段(Clean reads),两端平均Q20为97.50%,平均Q30为93.81%。GO富集分析结果显示,AaCK的HEGs富集于26个GO terms,基因富集数最多的是细胞(22 Unigenes),其次是细胞组件(22 Unigenes)和代谢过程(21 Unigenes);Aam T的HEGs富集于22个GO terms,基因富集数最多的是催化活性(23 Unigenes),其次是细胞进程(18 Unigenes)和代谢过程(18 Unigenes)。KEGG代谢通路(Pathway)富集分析显示,AaCK的HEGs富集在109个Pathways上,基因富集数最多的是核糖体(179 Unigenes),其次是氨基酸生物合成(70Unigenes)和碳代谢(62 Unigenes);Aam T的HEGs富集在114个Pathways上,基因富集数量最多的是核糖体(178 Unigenes),其次是碳代谢(116 Unigenes)和氧化磷酸化(112 Unigenes)。Venn分析结果显示,AaCK与Aam T共有的HEGs有260个,二者特有的HEGs分别为2 161个和4 445个。【结论】提供了胁迫意蜂6日龄幼虫肠道的球囊菌及其纯培养的HEGs表达谱,揭示了球囊菌在胁迫前和胁迫后期的基因表达规律,也为阐明球囊菌致病的分子机理提供了有益的信息和线索。     

DNA RFLPs at a highly polymorphic locus distinguish European and African subspecies of the honey bee Apis mellifera L. and suggest geographical origins of New World honey bees

A highly polymorphic locus in the honey bee, Apis mellifera L., was detected with genomic probe pB178. Eighty-five alleles, consisting of Msp I and Dde I RFLPs, were found among the Old and New World bees tested. Forty-one Msp I and 43 Dde I restriction fragment patterns, or variants, were identified. Variants and alleles were discontinuously distributed in Old World European and African subspecies. Principal coordinate analysis of the genetic distances between the alleles resulted in the identification of three distinct groups corresponding to three groups of honey bee races with historically different geographical distributions: east European A. m. ligustica and A. m. caucasica ; west European A. m. mellifera ; and South African A. m. scutellata . The clustering of alleles into these groups is consistent with previous honey bee phylogeographic studies, employing other nuclear and mitochondrial DNA markers, which in part support the evolutionary history of the honey bee hypothesized by Ruttner based on morphometric and allozyme data. The majority of alleles in bees from the USA grouped with those found in east European bees, while other alleles grouped with alleles found in A. m. mellifera . While the majority of the alleles in neotropical bees grouped with or were identical to African alleles, other alleles grouped with alleles found in A. m. mellifera, A. m. ligustica , and A. m. caucasica . Clues to the ancestry of neotropical bees may be found in the identification of alleles that were identical or more similar to alleles found in South African and west European bees; evidence for west European ancestry has been suggested using other taxonomic characters that were not unique to west European bees. Both west European and African alleles were found in individual neotropical colonies, which may indicate that honey bee subspecies which evolved allopatrically have hybridized in the human-assisted extension of their original geographical ranges.