中国舞毒蛾六个地理种群的RAPD分析及SCAR标记构建

舞毒蛾Lymantria dispar L.是世界性农林害虫, 包含不同的亚种, 其中亚洲舞毒蛾的雌蛾具有较强的飞行能力, 已成为国际性的重要检疫性有害生物。然而, 不同舞毒蛾亚种及种群间形态难辨, 因此采用传统的手段鉴别舞毒蛾亚种种群是很困难的。本研究首先采用RAPD标记分析了中国舞毒蛾6个地理种群的遗传多态性。结果表明, 所检测的舞毒蛾种群的遗传分化系数G_st为0.7571, 由此推算出的平均有效迁移数（基因流参数）N_em为0.1604, 说明不同舞毒蛾种群间的遗传分化程度较高, 缺乏广泛的基因流动。本研究在RAPD遗传分析基础之上, 筛选出了4个舞毒蛾种群的特异性遗传位点, 然后对这些特异性位点进行了克隆测序、 序列分析和位点特异性引物设计。结果表明, 其中2个舞毒蛾种群的位点特异性引物可产生序列特征性扩增区域（SCAR）标记。经验证, 这些标记可被用来鉴别特定的舞毒蛾地理种群, 因此有助于对这些舞毒蛾地理种群的分布与扩散进行监测。     

欧洲型舞毒蛾气味受体基因OrCo的克隆及序列分析

【目的】克隆欧洲型舞毒蛾Lymantria dispar Linnaeus的OrCo气味受体基因,并分析其序列特征,为进一步研究舞毒蛾嗅觉机制提供有益参考。【方法】本研究利用RT-PCR和RACE方法,克隆获得欧洲型舞毒蛾OrCo受体基因cDNA全长序列,将该基因命名为LdisOrCo,在GenBank中的登录号为KF482409。【结果】序列分析结果显示,LdisOrCo开放阅读框全长为1 311 bp,编码436个氨基酸,序列中有7个跨膜区和高度保守的C端区域。序列联配分析表明,LdisOrCo基因的氨基酸序列与近缘种灰翅夜蛾(Spodoptera littoralis Or83b)的同源性高达89%,与鞘翅目台湾黑金龟(Holotrichia plumbea Or83b)同源性高达64%,与已经报道的其他昆虫的嗅觉受体同源性都在60%以上,特别是在C端几乎完全一致。【结论】嗅觉受体OrCo(Olfactory receptor coreceptor)在不同昆虫体内高度保守,克隆欧洲型舞毒蛾OrCo基因可以为进一步研究舞毒蛾气味受体的功能,OrCo基因的进化以及揭秘嗅觉机制奠定基础。     

舞毒蛾Ldcht10基因序列的克隆及其功能分析

【目的】研究舞毒蛾Lymantria dispar几丁质酶基因的时空表达特性及其在蜕皮发育过程中的生物学功能,筛选在舞毒蛾发育过程中致死性的几丁质酶基因,为实现基于RNAi的舞毒蛾有效控制提供重要的基础数据。【方法】设计简并引物克隆舞毒蛾几丁质酶基因Ldcht10关键序列,通过使用实时荧光定量PCR方法检测该基因在舞毒蛾不同龄期与组织中的相对表达量,选取部分序列的双链RNA(ds RNA)研究该基因功能。【结果】本研究克隆了长度为2 057 bp的舞毒蛾几丁质酶基因Ldcht10序列;各组织与不同龄期RT-PCR结果表明Ldcht10的时空表达特性存在明显差异,Ldcht10在各个龄期均表达活跃且在前肠与后肠中的表达水平最高;Ldcht10的RNA干扰试验表明:注射Ldcht10 ds RNA后Ldcht10的表达受到极大抑制,该基因被沉默后24.3%的舞毒蛾幼虫因无法完成蜕皮或化蛹而死亡。【结论】舞毒蛾中几丁质酶基因Ldcht10在各个龄期与组织中的表达存在差异,且在舞毒蛾蜕皮过程中具有十分重要的生物学功能,该基因被沉默后部分舞毒蛾因无法完成蜕皮而导致死亡。     

荧光素对舞毒蛾核型多角体病毒不同地理品系的增效与光保护作用

舞毒蛾Lymantria dispar L.是源于欧亚大陆的多食性叶部害虫,取食300多种乔灌木,现已分布于北美、北非,成为世界性危险害虫之一,给林业生产带来了巨大损失。舞毒蛾核型多角体病毒(LdMNPV)是控制舞毒蛾种群动态的重要生物因素,可引起舞毒蛾种群急剧下降。在室内采用青杨枝条饲养的方法,测定了来自中国、北美和日本的3个LdMNPV品系(分别为LdMNPV-H,LdMNPV-D和LdMNPV-J)对危害青杨的亚洲型舞毒蛾幼虫的毒力,并测定了荧光素Tinopal LPW对它们的增效和光保护作用。结果表明Tinopal LPW对LdMNPV 3个地理品系均有增效和光保护作用,而且随着Tinopal LPW浓度的增加,增效作用增强,1%Tinopal LPW的增效作用最好。添加1%Tinopal LPW的LdMNPV-D品系、LdMNPV-H品系和LdMNPV-J品系对取食青杨的舞毒蛾幼虫的致死中浓度(LC50)分别为1.0、1.6、17.6 OBs/μL,不添加1%Tinopal LPW时,它们的LC50分别为32.9、39.0、1076.4 OBs/μL,分别降低了33、24、61倍。不添加1%Tinopal LPW时,D、H和J品系对舞毒蛾二龄幼虫的致LC95分别是2125.5、1275.8、303540.0 OBs/μL,添加1%Tinopal LPW后LC95分别为73.0、285.4、2360.8OBs/μL,分别降低了26、4.5、128.6倍。此外,1%Tinopal LPW的荧光素使3个品系的致死中时间(LT50)分别缩短了2.9d、5.3d、1.2d。LdMNPV-D和LdMNPV-H品系对亚洲型舞毒蛾表现出低致死中浓度、较短的致死中时间和较大的斜率,二者的毒力较LdMNPV-J品系高,在生产实践中应选择LdMNPV-D添加1%Tinopal LPW。Tinopal LPW对LdMNPV-D、LdMNPV-H和LdMNPV-J 3个品系均有光保护作用,添加1%Tinopal LPW后在距离30W紫外灯40cm下照射16h后,它们毒力保持系数比未添加Tinopal LPW分别高1.8、2.6、1.8倍。     

绿原酸提高舞毒蛾核型多角体病毒(LdNPV)的致病力

【目的】明确绿原酸(chlorogenic acid,CA)对舞毒蛾核型多角体病毒Lymantria dispar necleopolyhydrovirus(Ld NPV)致病力的影响,为舞毒蛾的防治提供参考依据。【方法】采用食料给毒法进行生物测定,测定舞毒蛾2龄幼虫对单独Ld NPV及添加绿原酸的Ld NPV(CA+Ld NPV)的剂量及致死时间响应。【结果】CA+Ld NPV与单独Ld NPV对舞毒蛾2龄幼虫的剂量及时间响应间有显著差异,二者对舞毒蛾2龄幼虫的致死中浓度(LC50)分别为161.8 OBs/μL(95%的置信区间为105.6~235.3 OBs/μL)和264.4 OBs/μL(95%的置信区间为178.6~384.0 OBs/μL),前者对舞毒蛾的致病力较后者强。CA+Ld NPV对舞毒蛾2龄幼虫的致死中时间(LT_(50))较Ld NPV的短,当Ld NPV浓度为590 OBs/μL,CA+Ld NPV及Ld NPV对舞毒蛾2龄幼虫的LT_(50)分别为9.9 d和12.3 d;当Ld NPV浓度为5 900 OBs/μL时,则LT_(50)分别为6.9 d和8.0 d。绿原酸降低了Ld NPV对舞毒蛾幼虫的致死中浓度,缩短了致死中时间。【结论】绿原酸可提高Ld NPV对舞毒蛾2龄幼虫的致病力,其机理有待进一步研究。     

温度胁迫对舞毒蛾营养成分和总抗氧化力的影响

【目的】为了探究温度胁迫对舞毒蛾Lymantira dispar幼虫营养成分和总抗氧化力的影响。【方法】测定了不同温度(15、20、25、30、35、40℃)下处理3 h的舞毒蛾幼虫可溶性蛋白、碳水化合物、海藻糖、脂质含量以及总抗氧化力,以饲养温度25℃的舞毒蛾幼虫为对照。【结果】舞毒蛾幼虫在不同温度下可溶性蛋白含量变化无显著差异;碳水化合物与海藻糖含量均随着温度的升高,呈现逐渐下降趋势;脂质含量无论处理温度低于还是高于对照温度均会增加;而舞毒蛾幼虫总抗氧化力在经不同温度处理后均显著降低。【结论】温度变化可以显著影响舞毒蛾幼虫碳水化合物、海藻糖、脂质3种营养成分含量以及舞毒蛾幼虫的总抗氧化力的能力,而对舞毒蛾幼虫可溶性蛋白含量无显著影响。本研究对舞毒蛾幼虫的耐热性及其在不同温度下的防治具有指导意义。     

舞毒蛾保幼激素结合蛋白家族基因克隆及对CO_2胁迫的响应

保幼激素结合蛋白(JHBP)是存在于血淋巴和细胞内的一类载体蛋白,与保幼激素结合运送到靶标组织。本文研究舞毒蛾JHBP基因特性及其对高浓度CO_2胁迫的响应,为明确全球气候变化下舞毒蛾适应性机制提供理论依据。通过舞毒蛾转录组文库分析结合RT-PCR克隆鉴定出7个JHBP基因,并进行基因特性和发育阶段特异性分析,同时利用密闭式CO_2人工气候箱在不同CO_2浓度(397μL/L、550μL/L和750μL/L)下将舞毒蛾卵饲养至3龄幼虫,利用qRT-PCR技术测定3龄幼虫JHBP基因表达量变化。结果表明,舞毒蛾JHBP家族7个基因全长开放阅读框大小为714～756 bp,编码237～251个氨基酸,分子质量为28.22～28.54 kDa,理论等电点为5.33～8.47。7个基因在不同发育阶段的表达存在差异,幼虫期LdJHBP1、2、5、6基因表达量较高,而LdJHBP3、4和7在蛹期和成虫期高表达。进化树分析表明舞毒蛾LdJHBP1、LdJHBP2、LdJHBP6分别与棉铃虫Helicoverpa armigera、家蚕Bombyx mori、冬尺蛾Operophtera brumata的JHBP亲缘关系较近。高浓度CO_2下舞毒蛾3龄幼虫的LdJHBP表达量下降。JHBP基因表达可能影响JH结合与运输从而调节生长发育。     

Cd胁迫对舞毒蛾幼虫食物利用的影响及其Cd的外泌机制

为探明植食性昆虫对受重金属胁迫的寄主植物的生理生态响应机制,本研究用Cd胁迫下银中杨的叶片饲养舞毒蛾幼虫,分析舞毒蛾幼虫对食物的利用情况以及其对Cd的排毒代谢机制.结果表明: 取食Cd胁迫下银中杨的叶片后,舞毒蛾3～6龄幼虫体内的Cd浓度和Cd含量均显著高于对照,但随着幼虫龄期增长,其体内Cd浓度显著降低,而Cd含量有不同程度的提高;舞毒蛾幼虫粪便和虫蜕中的Cd浓度均显著高于对照;舞毒蛾3～5龄幼虫的食物消耗率显著高于对照,而转化率显著低于对照;3～4龄幼虫的食物利用率均与对照差异不显著,但在5龄时显著低于对照.说明在Cd胁迫下,舞毒蛾幼虫能通过有效的排毒代谢途径将体内富集的部分Cd排出体外,且高龄幼虫的排毒代谢能力强于低龄幼虫;舞毒蛾幼虫体质量的增加会对体内的Cd浓度形成一种稀释效应;舞毒蛾幼虫能通过调整食物消耗率和转化率之间的比例,来维持其正常生长发育所需的食物利用率,但超过一定限度后仍会造成食物利用率降低.     

舞毒蛾末龄幼虫部分口器的扫描电镜观察

利用扫描电子显微镜对亚洲型舞毒蛾末龄幼虫口器的下颚和下唇的形态结构特征进行了观察。结果表明：在下颚和下唇须上共存在5种感器和1种角质齿,即毛形感器、刺形感器、锥形感器、栓锥形感器、钟形感器和刺状角质齿,其中毛形感器最发达。     

丁香酸可增强LdNPV对舞毒蛾致死及亚致死作用

【目的】研究丁香酸(syringic acid,SY)对舞毒蛾核型多角体病毒(LdNPV)致死及亚致死作用的影响,为舞毒蛾Lymantria dispar L.的防治提供科学依据。【方法】采用食料给毒法进行生物测定,评估LdNPV及SY与LdNPV共同对舞毒蛾2龄幼虫的致死和亚致死作用。【结果】取食添加丁香酸人工饲料和正常人工饲料的舞毒蛾幼虫LdNPV致死中浓度(LC50)分别为91.1 OBs/μL和194.8 OBs/μL。当LdNPV浓度为2 400 OBs/μL,取食添加SY的饲料及正常人工饲料的舞毒蛾幼虫LdNPV致死中时间(LT_(50))分别为9.3 d和10.1 d,当LdNPV浓度为24 000 OBs/μL时,其LT_(50)分别为7.8 d和8.1 d,LdNPV的浓度越高致死中时间越短。丁香酸增加了LdNPV对舞毒蛾幼虫的致病力,降低了LdNPV对舞毒蛾幼虫的致死中浓度,缩短了致死中时间。LdNPV对取食正常饲料的雌、雄幼虫历期和蛹期及雌蛾寿命影响不显著,但显著缩短了雄虫寿命,降低了平均每雌产卵量。LdNPV对取食添加丁香酸饲料的幼虫产生明显的亚致死效应,并且LdNPV浓度高于240OBs/μL时雌、雄幼虫历期延长、雄蛾寿命缩短,雄性羽化率降低,雌成虫产卵量降低,产卵量较对照减少一半左右。【结论】结果提示,丁香酸可提高LdNPV对舞毒蛾幼虫的致病力,增加LdNPV对舞毒蛾的亚致死作用。     

北京松山国家级自然保护区兰科植物多样性及其保护评价

选取北京松山国家级自然保护区开展兰科植物多样性调查和分析,评估其保护现状。松山共有兰科植物12属18种,其中北方盔花兰为北京新记录种。在生活型方面,以地生兰为主,计14种（77．78％）,腐生兰4种（22．22％）;在分布区类型方面,属的区系成分以北温带分布为主（6属,5000％）,兼有旧世界温带分布（3属,25．00％）、世界分布（2属,16．67％）和旧世界热带分布（1属,8．33％）;种的区系成分以东亚分布为主（9种,50．00％）,兼有旧世界温带分布（5种,27．78％）和北温带分布（4种,22．22％）。松山兰科植物集中分布于北沟和小海坨山。通过综合分析松山兰科植物保护现状,提出生境评估、种群动态监测、建立迁地保护资源圃和加大宣传等都是切实加强兰科植物保护工作的可行措施。     

异色瓢虫与隐斑瓢虫的区别及其色斑型和横脊的频率

异色瓢虫Harmonia axyridis(Pallas)和隐斑瓢虫H.yedoensis(Takizawa)是亲缘关系很近的姐妹种,在中国(北京到南岭,西至甘肃南部和西藏东部)、日本(北海道以南)和南韩等地同域分布,常常在松树等植物上共存。长期以来,隐斑瓢虫被认为是异色瓢虫的一个异名。有时从外形上很难区分这2种瓢虫,因此对于依据野外数据,分析异色瓢虫色斑型和鞘翅横脊的发生频率及小进化会产生一些错误。本文从没有隐斑瓢虫分布的东北地区的材料及其他数据,报道异色瓢虫的色斑型及鞘翅横脊发生频率。异色瓢虫的花斑型在我国东部地区(从东北至广东北部,西至甘肃和云南)发生率很低,而在新疆西北部发生率较高。鞘翅横脊的发生率从东北佳木斯的98.85%降低到云南大理的78.26%。本文列出了区分这2种瓢虫的形态特征及自然分布。如果鞘翅具横脊,则属于异色瓢虫,但如果标本来自2种瓢虫的共存区,鞘翅没有横脊,鞘翅的斑纹呈花斑型、四窗型或二窗型,则很难从外部形态上对它们进行鉴定。但这2种瓢虫的幼虫很容易区分。本文还提供了2种瓢虫的雄性外生殖器形态图、幼虫和成虫。     

利用CLIMEX预测纳塔尔实蝇在中国的潜在地理分布

纳塔尔实蝇Ceratitis rosa Karsch属双翅目实蝇科腊实蝇属,为害30余种经济植物,被我国列为进境植物检疫性有害生物。本研究运用CLIMEX 3.0及ArcGIS 9.3对纳塔尔实蝇在我国目前及未来的潜在地理分布进行了预测。结果显示：在目前的气候条件下,纳塔尔实蝇在我国的潜在地理分布区为18.250°N-30.250°N,其中,华南和东南地区为高度潜在地理分布区,包括四川、重庆、云南、福建、广东、广西、海南以及台湾等地。在未来的气候条件下,2020和2050年潜在地理分布区的北界分别移至31.250°N和32.250°N,且中度潜在地理分布区北移明显;2100年潜在地理分布区的北界移至33.750°N,且高度适生区北移明显。因此建议目前应加强纳塔尔实蝇检疫措施,完善监测体系,监测网点主要设在我国的南方地区,尤其应对云南、广西、广东、福建和海南等地进行长期监测;同时,应对湖南、贵州、江西、湖北、江苏、安徽等省进行定期监测,严防该虫入侵。     

基于最大熵生态位元模型的入侵杂草春飞蓬在中国潜在分布区的预测

春飞蓬是源自北美的外来入侵杂草,已入侵我国的上海、江苏、安徽等地,并对当地生态系统造成了危害,表现出巨大的入侵潜力.本文将最大熵生态位元(MaxEnt)模型与地理信息系统相结合,综合分析了影响春飞蓬分布的环境因素,直观和定量地预测了春飞蓬在我国的潜在分布区,并采用接受者操作特性曲线( ROC)分析法对预测结果进行验证.结果表明:春飞蓬在我国有广泛的潜在分布区,主要适生区包括上海、江苏、浙江、安徽、河南、湖北、湖南、江西等地.其目前的实际分布远未达到最大潜在分布范围,仍有可能继续扩散.ROC曲线检验结果表明,MaxEnt模型对春飞蓬在中国潜在分布区的预测精度较高,结果可信.春飞蓬潜在入侵区分布可能主要受温度和降水的影响.应高度重视春飞蓬的危害性,采取灭除已经发生的种群、严格监控向最适分布区入侵等措施阻止春飞蓬的扩散蔓延.     

椰子织蛾传入中国及其海南省的风险性分析

【背景】椰子织蛾原产于印度和斯里兰卡,现主要分布于南亚和东南亚,是椰子、中东海枣、蒲葵等多种棕榈科植物的重要害虫。我国目前尚无椰子织蛾发生为害的报道;但从适生性的角度看,该虫的地理分布区与我国南方棕榈植物分布区具有相似的生物气候,理论上我国热带、亚热带的广大区域是椰子织蛾的潜在地理分布区。鉴于其危害严重性,一旦入侵,其危险程度将超过已在我国发生的重要入侵害虫椰心叶甲。【方法】根据国际植物检疫措施标准（ISPM）规定的有害生物风险性分析（PRA）程序,利用相关风险性分析模型,从国内外分布状况、潜在的危害性、受害寄主的经济重要性、传播扩散的可能性及危险性管理难度5个方面,对椰子织蛾在中国及其海南省的风险性进行定性、定量分析。【结果】椰子织蛾在中国及其海南省的风险值（R）分别为2．20和2．30,在我国属于高度危险性有害生物。【结论与意义】本研究可为制定椰子织蛾的检疫管理对策提供依据。     

新疆喀纳斯自然保护区藓类植物物种多样性与区系调查

从1995～2009年,作者对新疆喀纳斯自然保护区藓类植物进行了全面的调查、采集、鉴定,分析研究。结果表明,该自然保护区藓类植物共有37科113属253种（含1变种）,优势科主要体现寒、温带地区分布类型,该区藓类植物生境的多样性可以从物种多样性,如曲尾藓科、木灵藓科、泥炭藓科、紫萼藓科、丛藓科等的广泛分布反映出来。通过与该区相邻的6个地区的物种丰富度系数比较,该自然保护区综合系数S i排第三位,说明此区藓类植物物种多样性丰富,这与本区的水分条件有很大关系。藓类植物区系成分以北温带成分为绝对优势,占总种数的83.33%;东亚成分占总种数的5.10%;东亚-北美成分占总种数的4.17%,各区系成分交汇。主成分分析和聚类分析表明,喀纳斯自然保护区藓类植物区系与新疆东部天山、内蒙古大青山、吉林长白山、河北木兰围场藓类植物区系关系较接近,与贵州梵净山、云南鸡足山藓类植物区系关系疏远,体现了我国北方藓类植物分布特征。     

中国珍稀药用植物白及的潜在分布与其气候特征

基于白及（Bletillastriata）地理分布信息和影响白及分布的主导气象因子,利用最大熵模型和地理信息系统预测了中国白及的潜在分布与其气候特征。结果表明：所建模型经ROC曲线验证,预测效果非常好（AUC〉0．9）：中国白及潜在分布区主要位于秦岭、淮河以南大部分地区,白及主要适生省份为云南、湖北、四川、湖南、江西、浙江;刀切法测试表明,4月和10月最低气温、年温度变化范围、11月平均降水量为影响白及潜在分布的最主要气象因子。研究结果可为白及资源的调查、保护与人工种植产业的发展提供理论指导。     

Phylogenetic relationships of Loxosceles and Sicarius spiders are consistent with Western Gondwanan vicariance

The modern geographic distribution of the spider family Sicariidae is consistent with an evolutionary origin on Western Gondwana. Both sicariid genera, Loxosceles and Sicarius are diverse in Africa and South/Central America. Loxosceles are also diverse in North America and the West Indies, and have species described from Mediterranean Europe and China. We tested vicariance hypotheses using molecular phylogenetics and molecular dating analyses of 28S, COI, 16S, and NADHI sequences. We recover reciprocal monophyly of African and South American Sicarius, paraphyletic Southern African Loxosceles and monophyletic New World Loxosceles within which an Old World species group that includes L. rufescens is derived. These patterns are consistent with a sicariid common ancestor on Western Gondwana. North American Loxosceles are monophyletic, sister to Caribbean taxa, and resolved in a larger clade with South American Loxosceles. With fossil data this pattern is consistent with colonization of North America via a land bridge predating the modern Isthmus of Panama.     

Phylogeographic structuring of Plasmodium lineages across the North American range of the house finch (Carpodacus Mexicanus)

The determinants of the geographic distribution of avian hematozoa are poorly understood. Sampling parasites from one avian host species across a wide geographic range is an accepted approach to separate the potential influence of host species distribution from geographic effects not directly related to host species biology. We used polymerase chain reaction to screen samples for hematozoan infection from 490 house finches (Carpodacus mexicanus) collected at 8 sites spanning continental North America. To explore geographic patterns of parasite lineage distributions, we sequenced a portion of the mitochondrial cytochrome b gene of Plasmodium species infecting 77 house finches. We identified 5 distinct Plasmodium haplotypes representing 3 lineages that likely represent 3 species. One lineage was common at all sites where we detected Plasmodium species. The second lineage contained 3 haplotypes that showed phylogeographic structuring on a continent-wide scale, with 1 haplotype common in eastern North America and 2 common in western North America. The third divergent lineage was recovered from 1 individual host. Considered together, the partial phylogeographic structuring of Plasmodium cytochrome b lineages over the range of the house finch suggests that parasite lineage distribution is not solely dependent on host species distribution, and other factors such as arthropod vector competence and distribution may be important.     

Genetic structure,admixture and invasion success in a Holarctic defoliator,the gypsy moth (Lymantria dispar,Lepidoptera: Erebidae)

Characterizing the current population structure of potentially invasive species provides a critical context for identifying source populations and for understanding why invasions are successful. Non‐native populations inevitably lose genetic diversity during initial colonization events, but subsequent admixture among independently introduced lineages may increase both genetic variation and adaptive potential. Here we characterize the population structure of the gypsy moth (Lymantria dispar Linnaeus), one of the world's most destructive forest pests. Native to Eurasia and recently introduced to North America, the current distribution of gypsy moth includes forests throughout the temperate region of the northern hemisphere. Analyses of microsatellite loci and mitochondrial DNA sequences for 1738 individuals identified four genetic clusters within L. dispar. Three of these clusters correspond to the three named subspecies; North American populations represent a distinct fourth cluster, presumably a consequence of the population bottleneck and allele frequency change that accompanied introduction. We find no evidence that admixture has been an important catalyst of the successful invasion and range expansion in North America. However, we do find evidence of ongoing hybridization between subspecies and increased genetic variation in gypsy moth populations from Eastern Asia, populations that now pose a threat of further human‐mediated introductions. Finally, we show that current patterns of variation can be explained in terms of climate and habitat changes during the Pleistocene, a time when temperate forests expanded and contracted. Deeply diverged matrilines in Europe imply that gypsy moths have been there for a long time and are not recent arrivals from Asia.