陈冬平，陈志廷，徐汉虹，林 菲

(天然农药与化学生物学教育部重点实验室/华南农业大学 农学院，广东 广州 510642)

草地贪夜蛾Spodoptera frugiperda是一种原产于美洲的热带和亚热带地区的杂食性害虫，严重威胁农业生产。其寄主范围广，能够为害玉米、水稻、棉花、大豆等80 多种作物[1]。不同的草地贪夜蛾群体通常对特定作物具有取食趋性[2]。田间草地贪夜蛾主要分为玉米型(Corn strain)和水稻型(Rice strain)2 个亚群[3]。玉米型主要取食玉米、高粱、棉花，而水稻型的取食性相对较为广谱，主要取食水稻、苜蓿、牧草、小米等作物[4]。除了取食偏好外，2 个品系的草地贪夜蛾在生理特征、交配行为和性外激素的组成方面有很大的不同。但是，由于2 种寄主型形态学和行为方面的相似性，在田间难以被区分开[5]。因此，对2 种寄主型的草地贪夜蛾开展分子鉴定，将有助于了解该入侵物种的遗传结构，掌握害虫的入侵和扩散规律以及制定有效的防控策略[6-7]。

草地贪夜蛾亚群的分类主要根据田间采集的幼虫寄主范围，并结合分子标记的分化进行界定。基于线粒体细胞色素C 氧化酶亚基I 基因(COI)的分子标记不仅能够区分玉米型和水稻型2 个亚群的，而且具备鉴别其亚型的能力，是目前最常用的有效分子条形码[7-10]。另外，有研究表明，根据Z 染色体上编码磷酸丙糖异构酶基因(Tpi)开发的分子标记，能够更好地鉴定到草地贪夜蛾雌性“水稻型”和雄性“玉米型”的杂交型[7,11]。2 种标记有机结合能够提高鉴定的准确率[12-13]。

2018 年底草地贪夜蛾通过中缅边境进入我国境内，目前已在云南、湖南、广西、广东等地发现其对玉米的为害[14]。广东省是水稻主要种植区，严密监测草地贪夜蛾是否存在为害水稻的风险显得格外重要[15]。另外，广东省为热带和亚热带季风气候区，草地贪夜蛾可在该省越冬，是可阻止其北迁的桥头堡[16]。本研究通过基于COI和Tpi2 个基因片段的分子标记，比较它们对入侵广州的草地贪夜蛾的鉴别能力，为监测虫源扩散规律以及预测预警虫情提供参考。

1 材料与方法

1.1 供试昆虫采集

供试草地贪夜蛾于2019 年5 月12 日在广州花都花山镇儒林村玉米地采集。寄主为喇叭口期的春玉米，所有样本均为4～6 龄幼虫，6 个样本分别标记为1#、2#、3#、4#、5#和6#。单头活体幼虫采集后放入含有新鲜玉米叶片的50 mL 离心管，带回实验室进行形态学鉴定后，用液氮速冻储存于−80 ℃用于DNA 提取。

1.2 DNA 提取和扩增

采用DNA/RNA 共提取试剂盒(TIANGEN，中国)分别提取单头幼虫DNA 和RNA。虫体放入含有1.5mL DNA 抽提液的2mL 离心管内，加入钢珠，放入样品研磨机(Bullet Blender®，美国)处理后，按照试剂盒说明书操作。DNA 使用NanoDrop分光光度计(2000，赛默飞世尔)测定质量浓度，稀释到10 ng/μL 后，用作PCR 扩增模板。

用于PCR 扩增COI和Tpi部分基因片段的引物如表1 所示，选用的PCR 引物在COI基因的相对位置如图1 所示。采用PrimeSTAR®Max DNA Polymerase PCR(TaKaRa，日本)进行PCR 扩增，反应体系包括PrimeSTAR®Max DNA Polymerase PCR(2×)10μL，上游引物和下游引物各0.4μL，模板1μL，加ddH2O 至20μL；PCR 反应为：98℃预变性2min；98℃变性10s，60℃退火15s，72℃延伸30s，35个循环；72℃延伸5min。

表1 扩增COI 和Tpi 基因片段所用引物Table 1 Prim ers used for am p lifying COI and Tpi gene fragments

图1 用于单倍型分析的COI 基因3 个片段Fig.1 Three fragments of COI gene used for haplotype analysis

使用15 g/L 的琼脂糖凝胶电泳检测扩增片段特异性后，PCR 产物由广州天一辉远基因科技有限公司进行切胶测序。对于测序结果为双峰或者低质量的样品，则将其PCR 产物克隆到pEASY®-Blunt Cloning Vector 后，每个样品各挑取10 个克隆进行测序。

1.3 序列分析

利用Lastergene 软件包对正反向测序序列进行人工拼接、峰型校正和序列比对。从GenBank 数据库下载已经发布的草地贪夜蛾COI和Tpi基因序列。针对已报道的基于COI和Tpi基因序列的“水稻型”和“玉米型”单倍型特点，分别进行差异位点分析，确定本研究中草地贪夜蛾的具体亚型。

2 结果与分析

2.1 COI 标记鉴定草地贪夜蛾单倍型

利用标记COI-A，COI-B 和COI-C 对所采集的6 个样品进行了寄主型分析。结果（图2）显示，有5 个样品扩增出了预期大小的单一特异性条带，将其PCR 产物用于直接测序。

对COI-A 的序列分析表明，来自4 个样品(1#、2#、4#、6#) 的序列与已报道的草地贪夜蛾水稻型(HM 136593) 仅存在一处碱基差异，与玉米型(HM 136586)存在8 处单碱基差异，说明根据COIA 标记鉴定的寄主型与实际样品来源寄主型并不一致(图3)；对样品3#的COI-A 产物测序分析发现其为双峰杂合，克隆重新测序后，分别获得了样品3#-1 和3#-2 两种单倍型。其中，3#-1 与上述4 个样品检测到的单倍型完全一致，3#-2 则与Mythimna separate(KX863042) 对应的片段序列上存在1 个碱基的差异，说明3#样品可能是草地贪夜蛾和M.separate的杂合(图3)。

图2 6 头草地贪夜蛾CO I-A、CO I-B 和CO I-C 标记的PCR 扩增Fig.2 PCR amplification of CO I-A, COI-B and COI-C markers from six Spodoptera frugiperda

图3 基于COI-A 标记的单倍型序列比对Fig.3 Hap lotype sequence slignments based on COI-A markers

进一步用COI-B 和COI-C 标记对样品进行寄主型鉴定，结果表明，除了样品5#无法扩增出条带之外，其余的5 个样品均扩增出单一的预期条带。PCR 产物测序结果显示所测样品序列完全一致且为水稻型(图2)。

将COI-A、COI-B 和COI-C 等3 个标记的结果相比较，发现COI-A 能够区分出草地贪夜蛾和M.separate的2 种单倍型，而COI-B 和COI-C 仅能区分出草地贪夜蛾1 种单倍型，说明虽然基于COI基因片段序列鉴定到的草地贪夜蛾寄主型与实际采集寄主种类不相符，但是，COI-A 标记具有较好的物种鉴别能力，能够区分出杂合类型。

2.2 TPi 标记鉴定草地贪夜蛾单倍型

基于Tpi基因片段对所采集的6 个样品进行了寄主型分析。除了样品5#无法扩增出条带之外，其余的5 个样品均扩增出了预期大小的单一特异性条带(图4)。

图4 6 头草地贪夜蛾Tpi 基因片段的PCR 扩增Fig.4 PCR amplification of Tpigene fragment from six Spodoptera frugiperda

对PCR 产物进行双向测序并分析其序列差异，结果显示，所检测的序列与已报道玉米型一致，说明基于Tpi基因片段的鉴定结果与样品的寄主来源是一致的。样品2#和样品3#，在第126 和127 碱基处出现GA/AT 的杂合双峰(图5)，说明入侵广东的草地贪夜蛾种类并不单一。

图5 基于Tpi 基因标记的单倍型序列比对Fig.5 Haplotype sequence alignments based on Tpi gene marker

3 结论与讨论

本研究采集了入侵广州市的草地贪夜蛾样本，利用基于COI和Tpi基因片段的分子标记进行了寄主型的鉴定。基于COI-A、COI-B 和COI-C 标记的鉴定结果均显示入侵广州的草地贪夜蛾寄主单倍型为水稻型，与实际采集寄主不相符，可能原因是该标记不适合于检测入侵广州的草地贪夜蛾群体。但基于COI-A 的标记能够区分草地贪夜蛾与其他近源种的基因型，为研究入侵草地贪夜蛾的田间行为和扩散提供丰富的信息。而COI-B 和COIC 只鉴定到1 种单倍型。COI-A 和COI-B 标记事实上扩增的是COI基因的同一片段，差别在于COI-A 标记扩增的片段较小，且使用的是简并引物。因此，COI-A 标记具有比COI-B 标记和COIC 标记更强的鉴别能力。

Tpi基因标记的鉴定结果显示入侵广州的草地贪夜蛾寄主型为玉米型，与实际采集的寄主一致。基于Tpi基因标记，入侵广州的草地贪夜蛾群体被检测到具有2 种单倍亚型，在该片段第126 和第127 位碱基具有RW 杂合。先前报道的入侵云南的群体只有一种单倍型[12]，说明草地贪夜蛾可能在扩散的过程中出现遗传杂交现象。

草地贪夜蛾是联合国粮农组织全球预警的跨国界迁飞性重大害虫。目前，入侵我国的草地贪夜蛾主要为害玉米，未见在水稻上定殖为害。华南地区是我国水稻主产区，应慎重评估草地贪夜蛾为害水稻的风险。对于“玉米型”和“水稻型”的区分，仅依靠个别基因的片段进行评估是不全面的[6]。基于微卫星标记的分子指纹分析可以较全面地了解草地贪夜蛾的遗传结构动态变化规律[17]。获得与寄主取食相关基因的信息进而开发基因功能性分子标记，是精准鉴定草地贪夜蛾寄主单倍型的前提[18]。