随着分子生物学发展，应用分子系统分析的方法研究昆虫种群分化和遗传多样性，从分子水平探讨寄主、环境因子等对昆虫适应进化的影响，以明确昆虫适应进化的分子机制。蚧虫作为一类营相对固定生活的刺吸类害虫，呈聚集型空间分布，扩散能力弱，其发生与寄主种群以及地理环境有着密切的关系^{[1, 2, 3, 4]}，在长期的系统发育中形成了相对稳定的区系组成以适应不同环境，不同的寄主种群、生态地理种群间差异明显。高宝嘉等对皱大球蚧 ( Eulecanium kuwana i)和瘤坚大球蚧(E. gigan tean)地理种群的 RAPD遗传分化的研究发现，已发生种内遗传变异，但尚未达到亚种水平^[5]。Burban等通过对松干蚧的研究指出松干蚧种群变异与寄主植物的遗传结构有关^[6]。本研究选取河北省不同草履蚧Drosicha corpulenta(kawana)种群，对其遗传分化进行RAPD分析，揭示寄主植物、地理环境等生态条件对蚧虫遗传分化的影响，为深入探讨植物昆虫协同进化关系和昆虫适应潜力提供依据。

2011年4—5月于保定、石家庄、邯郸不同寄主上采集雌成虫(表1)。不同寄主相距50m以上，每种寄主东、南、西、北方向的上、中、下部枝条分别取样。将采集的蚧虫饥饿24h以上，使其消化道内的物质充分消化，以防外源DNA干扰实验结果;-20℃冰箱中保存。

本实验所用的TaqDNA聚合酶、dNTP、蛋白酶K、RAPD引物、DNA标记(DL2000)均由上海生物工程有限公司生产。

参照Boyce等的方法^[7]，并做适当修改。取1只蚧虫于灭菌的研钵(每个种群取28个样本)，加入液氮迅速研磨成粉末，放入1.5mL的离心管中，加入700μL预冷的2×CTAB(0.1mol /L Tris2-HCl，1.4mol /L NaCl ，0.2mol/L EDTA，0.2%β巯基乙醇，0.05mol /L CT AB，pH8.3)，加入蛋白酶K至终浓度为100μg/mL，65℃下水浴1h，用12000r/min离心5min，取上清液，加入等体积得酚∶氯仿∶异戊醇(25∶24∶1)轻轻混匀，反复抽提两次。再经等体积的异丙醇沉淀 1h，4℃下10000r/min离心10min，去上清，70%的乙醇清洗1—2次，干燥后加入适量TE溶液溶解 DNA，加入RNaesA至终浓度为2μg/ L后，37℃下水浴1h，-20℃ 保存待用。

经预试验筛选出4条可用引物，其引物名称及序列见表2。参照Williams的方法^[8]，反应体系25μL，其中模板50ng，10×Buffer2.5μL，Mg²⁺114—210 mmol/L (因座位而异)，dNTP200μmol/L，引物(0.8μmol/L)，Taq酶1.5U反应条件为：94℃预变性1min；92℃变性1 min，37℃退火1 min(因座位而异)，72℃延伸2 min，45个循环；72℃再延伸10min；4℃终止反应。

扩增产物用1.5%琼脂糖凝胶(含1%核酸显色剂)分离，电泳缓冲液为0.5×TBE，上样量为10μL。电泳结束后，用凝胶成像仪拍照。

RAPD为显性标记，同一引物的扩增产物在琼脂糖凝胶电泳中迁移律相同的条带被认为是同源性的，属于同一位点的条带按清晰可见的强带或反复出现的弱带记为1，无带记为0，形成二元数据。应用Popstream32和NTSYSpc2.1软件，计算种群遗传多样性参数并根据遗传相似系数进行UPGMA种群聚类分析。

引物扩增结果见图1，不同种群遗传多样性统计见表3。供试引物对16个种群448个个体进行扩增，共产生41条清晰、稳定的片段，平均每条引物扩增10.25条片段，各片断分子量大小在200—2000bp之间，种群内和种群间并未产生各自特征性带谱。实验共检测出41个多态位点，多态位点比率为100%，说明草履蚧种群遗传多样性丰富。但不同种群的扩增产物多态性结果不同，各种群多态位点比率在31.71%—63.41%之间，其中邯郸紫薇寄主种群多态位点比率最低(31.71%)，邯郸紫叶李寄主种群(63.41%)和邯郸桑树寄主种群(60.98%)多态位点最高，遗传多样性较其它种群丰富。

图 1 引物H9对16个草履蚧种群的RAPD-PCR结果 Fig.1 RAPD-PCR results of primer H9 for 16 populations of D. corpulenta 1—4：HL、HF、HZ、HS； 5—8：HB、HW、HC、HM； 9—12：SC、SX、SZ、SB； 13—16：BY；BB；BZ；BBT； 17：DL2000

图选项

Nei指数估计的不同种群间遗传距离和遗传一致度结果见表4。由表4可知，16个种群Shannon 信息指数的变化范围在0.1708—0.2967之间，平均值为0.2306。邯郸桑树种群的Shannon 信息指数最高，邯郸紫薇种群最低。16个种群Nei′s基因多样性指数变化范围在0.1101—0.1991之间，平均值为0.1511。邯郸桑树种群的基因多样性最高，保定紫叶碧桃种群最低，这与Shannon信息指数估计的遗传多样性稍有差异。说明同一物种不同种群遗传多样性存在差异。

由表4可知，16个草履蚧种群的遗传相似度在0.6466—0.9290之间，平均值为0.7873；遗传距离在0.0701—0.4360之间，平均值为0.2395，说明不同种群间遗传差异较大。其中以邯郸枫杨和以邯郸垂丝海棠为寄主的草履蚧种群最相似，遗传距离最小(0.0701)；以石家庄紫叶李和以邯郸木槿为寄主的种群遗传距离最大(0.4360)，说明二者的遗传差异最大。

紫叶李邯郸种群与石家庄种群遗传距离为0.3330，与保定种群遗传距离为0.3426；紫叶碧桃邯郸种群与石家庄种群遗传距离为0.3170，与保定种群遗传距离为0.3816；即同种寄主不同地理种群随地理距离的增加遗传距离增加。

用NTSYSpc2.1对草履蚧种群的遗传相似系数进行UPGMA聚类分析，进一步分析种群间遗传关系(图2)。结果显示：同一地理寄主种群绝大多数个体都能聚为一类，说明同一种群中个体间的遗传特性相似程度高，遗传稳定性；草履蚧相同寄主种群首先聚类，遗传相似系数在0.67水平，说明寄主因素对草履蚧的种群分化产生明显的影响；保定、石家庄、邯郸的紫叶碧桃种群各自聚类，且遗传距离较远；同样以上3地的紫叶李种群也各自聚类，说明地理因素对种群的分化产生影响；草履蚧不同寄主的三大地理种群并未产生明显地理聚类，说明相对于地理因素，寄主因素对草履蚧种群分化的影响更加明显。同一地区不同寄主种群聚类见图3、图4、图5，可见同一地区同一寄主的个体聚为一类，而同一地区不同寄主间聚类规律不明显。

图 2 16个草履蚧种群聚类分析图 Fig.2 The UPGMA dendrogram cluster of 16 Drosicha corpul populations种群代号同表1 (See table 1 for explanation of population code); 种群代号同表1 (See table 1 for explanation of population code);为了方便比较查看，本图中每个种群采用5个样本的聚类分析; 其聚类分析趋势与全部个体聚类图相一致，具有代表性

图选项

图 3 石家庄草履蚧种群聚类分析图 Fig.3 The UPGMA dendrogram cluster of Drosicha corpul populations of Shijiazhuang

图选项

图 4 保定草履蚧种群聚类分析图 Fig.4 The UPGMA dendrogram cluster of Drosicha corpul populations of Baoding

图选项

图 5 邯郸草履蚧种群聚类分析图 Fig.5 The UPGMA dendrogram cluster of Drosicha corpul populations of Handan

图选项

本文通过以上对草履蚧不同寄主和地理种群遗传分化的RAPD分析，表明： Shannon信息指数和Nei′s基因多样性指数显示草履蚧不同地理种群以及不同寄主种群遗传多样性丰富程度不同，草履蚧不同寄主种群之间存在不同程度的遗传差异，并已产生种群分化。草履蚧相同寄主不同地理种群之间存在遗传差异，地理因素对草履蚧种群分化产生影响。聚类分析结果显示，草履蚧种群分化同时受到寄主因素和地理因素的双重影响，相对于地理因素，寄主对草履蚧种群分化的影响更加明显。

草履蚧寄主植物广泛，主动扩散距离有限，因此其发生与寄主植物关系十分密切。本研究表明草履蚧不同寄主种群之间已经产生了明显的种群分化，这与草履蚧寄主生态条件、环境胁迫有密切关系。寄主植物自身的物理、化学等因素可影响植食性昆虫的生物学习性和遗传特性，促进昆虫种群分化，而昆虫对寄主植物具有选择性，对于不同寄主植物喜好程度也不尽相同^{[9, 10]}。昆虫学家认为植食性昆虫可因嗜食的寄主不同而造成遗传隔离，由此形成不同的物种。Takada、Blackman、Weber等研究表明烟蚜在不同寄主植物上存在着种群分化^{[11, 12, 13]}。本实验中寄主植物分类科属不尽相同，其生物学特性和化学成分各异，我们推测寄主植物种类不同是造成草履蚧种群分化的重要原因。

此外，大气污染以及昆虫群落结构^{[14, 15]} 等环境胁迫均能影响寄主植物与昆虫的关系，本实验寄主植物除受到草履蚧危害外，同时受到其他不同种昆虫取食。不同寄主植物上昆虫的种类、数量以及植物所处环境的不同可能都对草履蚧种群分化期作用，下一步有待从多因素分析探讨草履蚧种群分化原因。