基于线粒体COⅠ和COⅡ基因的沙果小食心虫与梨小食心虫的分子鉴定

基于线粒体COⅠ和COⅡ基因的沙果小食心虫与梨小食心虫
的分子鉴定
王康;李玉婷;郑燕;段辛乐;张蒙;彭雄;陈茂华
【摘要】[目的]建立基于线粒体COⅠ和COⅡ基因序列的沙果小食心虫、梨小食心虫分子鉴定方法.[方法]田间采集试虫,利用PCR和基因测序技术,扩增2种食心虫的线粒体CO Ⅰ和COⅡ基因,将获得的CO Ⅰ和COⅡ基因序列在GenBank中进行BLAST比对,并通过GeneDoc软件分析基因序列相似性,利用Mega 5.05软件分别计算基于CO Ⅰ和COⅡ基因序列的遗传距离,并分别构建CO Ⅰ和COⅡ基因的系统发育树.[结果]在分析的2种食心虫样本中,沙果小食心虫COⅠ基因的种内相似度在99.4％以上,梨小食心虫CO Ⅰ基因的种内相似度在99.6％以上,2种食心虫CO Ⅰ基因的种间相似度为95.0％～95.6％;沙果小食心虫COⅡ基因的种内相似度在99.0％以上,梨小食心虫COⅡ基因的种内相似度在99.3％以上,2种食心虫COⅡ基因的种间相似度为94.6％～95.5％;2种食心虫CO Ⅰ基因种间存在30个稳定变异位点,COⅡ基因种间有26个稳定变异位点.2种食心虫种内遗传距离为0.001～0.006(CO Ⅰ)和0.001～0.010(COⅡ),种间遗传距离为0.046～0.052(CO Ⅰ)和0.047～0.057(COⅡ),基于CO Ⅰ和COⅡ基因的种间遗传距离均显著大于种内遗传距离.分别构建CO Ⅰ和COⅡ基因单倍型的系统发育树,结果显示,在2个基因的系统发育树上,2种食心虫的基因序列分别位于不同的进化支,且置信度均达到100％,同种食心虫的基因序列位于相同的进化枝.[结论]可以根据本研究所用的CO Ⅰ和COⅡ基因序列的差异性,进行沙果小食心虫和梨小食心虫的分子鉴定.
【期刊名称】《西北农林科技大学学报（自然科学版）》
【年(卷),期】2016(044)002
【总页数】9页(P156-164)
【关键词】沙果小食心虫;梨小食心虫;COⅠ基因;COⅡ基因;单倍型;分子鉴定
【作者】王康;李玉婷;郑燕;段辛乐;张蒙;彭雄;陈茂华
【作者单位】西北农林科技大学植物保护学院,农业部西北黄土高原作物有害生物
综合治理重点实验室,植保资源与病虫害治理教育部重点实验室,陕西杨凌712100;西北农林科技大学植物保护学院,农业部西北黄土高原作物有害生物综合治理重点
实验室,植保资源与病虫害治理教育部重点实验室,陕西杨凌712100;西北农林科技大学植物保护学院,农业部西北黄土高原作物有害生物综合治理重点实验室,植保资
源与病虫害治理教育部重点实验室,陕西杨凌712100;西北农林科技大学植物保护
学院,农业部西北黄土高原作物有害生物综合治理重点实验室,植保资源与病虫害治
理教育部重点实验室,陕西杨凌712100;西北农林科技大学植物保护学院,农业部西北黄土高原作物有害生物综合治理重点实验室,植保资源与病虫害治理教育部重点
实验室,陕西杨凌712100;西北农林科技大学植物保护学院,农业部西北黄土高原作物有害生物综合治理重点实验室,植保资源与病虫害治理教育部重点实验室,陕西杨
凌712100;西北农林科技大学植物保护学院,农业部西北黄土高原作物有害生物综
合治理重点实验室,植保资源与病虫害治理教育部重点实验室,陕西杨凌712100【正文语种】中文
【中图分类】S763.42;Q78
沙果小食心虫(Grapholita dimorpha Komai)和梨小食心虫
(Grapholita molesta Busck)是重要的蛀果害虫，二者均属于鳞翅目(Lepidoptera)卷蛾科(Tortricidae)。

沙果小食心虫寄主植物为苹果、沙果、李等，在我国部分地区报道有发生[1-2]。

梨小食心虫又名东方果蛀蛾，寄主植物主要为桃、苹果、梨、
杏、李、山楂等，是世界性的主要蛀果害虫之一，在我国分布广泛[3]。

沙果小食心虫和梨小食心虫经常混合发生，而且发生期部分重叠，二者幼虫和成虫外部形态都非常相似，而且性信息素主要成分亦相似[4-5]，因此在田间调查和监测种群动态时，准确鉴定这2种食心虫十分困难[1,6-7]。

基于PCR技术的分子鉴定可以检测不同物种间及相同物种内不同个体间的细微差异，这为物种鉴定提供了新的手段[8-10]。

线粒体DNA具有基因序列相对保守、无重组和单拷贝、进化速率适中、易于检测等优点，因此常被应用于昆虫的分子鉴定与系统进化研究之中，其中线粒体细胞色素氧化酶亚基Ⅰ(COⅠ)和线粒体细胞色素氧化酶亚基Ⅱ(COⅡ)等基因应用较为广泛[11-12]。

常虹等[13]基于线粒体COⅠ基因构建了齿小蠹属(Ips)昆虫的快速鉴定方法，很好地解决了植物检疫中截获的齿小蠹幼虫或残缺成虫的鉴定难问题；游中华等[14]报道了基于线粒体COⅠ基因的9种蓟马的分子鉴定方法；通过COⅠ基因序列的分析发现，弄蝶(Astraptes fulgerator)是由多个物种组成的混合种[15]；王戎疆等[16]利用线粒体COⅡ基因构建蛱蝶属(Polyura) 5种蝴蝶的系统发育树并进行聚类关系分析，结果和形态学结果相一致；杨飞龙等[17]通过比较大量中华按蚊(Anopheles sinensis)标本的COⅠ、COⅡ基因及几种核糖体基因序列，建立了该虫基于COⅡ等基因的分子鉴定标准。

本研究从田间采集沙果小食心虫和梨小食心虫试虫，提取试虫的基因组DNA，PCR扩增试虫的线粒体COⅠ和COⅡ基因序列，并对扩增的序列进行测序，以此建立沙果小食心虫和梨小食心虫分子鉴定方法，为这2种害虫的准确鉴定及预测预报提供了重要依据。

利用沙果小食心虫(Grapholita dimorpha Komai)和梨小食心虫(Grapholita molesta Busck)信息素诱捕器(诱捕器为三角胶粘式，购自北京中捷四方生物科技公司)，在2种害虫的发生期，于陕西省的5个苹果园，分别诱集2
种食心虫成虫样本(表1)。

诱芯的载体中空，由硅橡胶制成，每个质量0.3～0.5 g，每个诱芯的性信息素含量不低于0.001 g，纯度90%～97%。

选择较为完整的试虫，在实验室利用体视显微镜进行准确的形态学鉴定，鉴定后的样本保存于无水乙醇中，置于-20 ℃冰箱中保存备用。

在采集的样本中，随机选取沙果小食心虫和梨小食心虫样本，利用EasyPureTM Genomic DNA Kit(北京全式金生物技术有限公司)提取单头试虫基
因组DNA，基因组DNA的提取方法按照试剂盒说明书步骤进行，用1%琼脂糖
凝胶电泳检测提取的DNA质量，提取的基因组DNA置于-20 ℃冰箱保存备用。

参照已有的报道，选用扩增2种食心虫COⅠ和COⅡ基因的引物[18]，稍作修改。

扩增2种食心虫COⅠ基因的上、下游引物序列分别为：5′-CAA-CATTTATTTTGATTTTTTGG-3′，5′-TCCAATGCACTAATCTGCCATATTA-3′；扩增2种食心虫COⅡ基因的上、下游引物序列分别为：5′-ATG-GCAGATTATATGTAATGG-3′，5′-GTTTAAG-AGACCAGTACTTG-3′。

PCR扩增反应总体系为25 μL，其中含2×PCR master mix 12.5 μL，上、下游引物各3 μL，DNA模板2 μL(含20～50 ng DNA)，加ddH2O至25 μL。

PCR扩
增程序为：94 ℃预变性2 min；94 ℃变性30 s，(CO Ⅰ基因：52 ℃/COⅡ基因：57 ℃)退火延伸30 s，72 ℃延伸30 s，共35个循环；循环结束后72 ℃延
伸5 min。

利用1%琼脂糖凝胶电泳检测PCR扩增产物。

将扩增的PCR产物直接进行测序，测序由生工生物工程(上海)有限公司完成。

将所测2种食心虫COⅠ和COⅡ基因在GenBank中进行BLAST搜索，并将获得的COⅠ和COⅡ基因登录到GenBank。

通过GeneDoc软件分析基因序列相似性，利用Mega 5.05软件中的Kimura2-parameter(K2P)双参数模型计算遗传距离，并利用邻接法(Neighbor-Joining，NJ)构建系统发育树[19]。

在所有沙果小食心虫和梨小食心虫试虫中都成功扩增到线粒体COⅠ和COⅡ基因
序列，其中COⅠ基因序列长度为720 bp，COⅡ基因序列长度为682 bp，所有
序列均无碱基插入和缺失。

将所获得的序列在GenBank中进行BLAST分析，结
果显示，本研究获得的沙果小食心虫COⅠ基因和COⅡ基因序列与GenBank中
登录的沙果小食心虫COⅠ基因(登录号：AB603522)和COⅡ基因(登录号：
JQ404438)对应序列的相似性分别为99.17%～99.31%和 99.19%～100%；测序获得的梨小食心虫COⅠ基因和COⅡ基因序列与GenBank中登录的梨小食心虫
线粒体全基因组(登录号：HQ116416)上对应的COⅠ基因和COⅡ基因序列的相
似性分别为 99.72%～100%和99.56%～100%。

沙果小食心虫和梨小食心虫COⅠ和COⅡ基因序列的比对结果分别见图1和图2。

供试样本共获得3个沙果小食心虫COⅠ基因单倍型(GenBank登录号：
KC923197～KC923199)，这3个单倍型在4个碱基位点上具有多态性，3个单倍型之间的相似度在99.4%以上。

供试样本共获得6个梨小食心虫COⅠ基因单倍型(GenBank登录号：KC923200～KC923205)，6个单倍型在5个碱基位点上具有多态性，6个单倍型之间的相似度在99.6%以上。

沙果小食心虫和梨小食心虫
COⅠ基因种间的相似度在95.0%～95.6%，COⅠ基因在2种食心虫之间存在30
个变异位点(图1，表2)。

供试样本共获得4个沙果小食心虫COⅡ基因单倍型(GenBank登录号：
KC923206～KC923209)，这4个单倍型在7个碱基位点上具有多态性，4个单倍型之间的相度在99.0%以上(表3)。

供试样本共获得8个梨小食心虫COⅡ基因单倍型(GenBank登录号：
KC923210～KC923217)，这8个单倍型在8个碱基位点上具有多态性，8个单倍型之间的相似度在99.3%以上。

沙果小食心虫和梨小食心虫COⅡ基因种间的相似度在94.6%～95.5%，种间存在26个变异位点(图2，表3)。

由表2可以看出，沙果小食心虫COⅠ基因序列间的遗传距离在0.001～0.006，
梨小食心虫COⅠ基因序列间的遗传距离在0.001～0.004，2种食心虫C OⅠ基因种间遗传距离为0.046～0.052，种间遗传距离显著高于种内遗传距离
(df=29,t=2.05，P<0.001)。

由表3可以看出，COⅡ基因中，沙果小食心虫序列间遗传距离在0.001～0.010，梨小食心虫序列间遗传距离在0.001～0.007，2种食心虫COⅡ基因种间序列遗传距离在0.047～0.057，种间遗传距离显著高于种内遗传距离(df=58,t=2.00，
P<0.001)。

以苹果蠹蛾(Cydia pomonella)COⅠ基因序列(GenBank登录号：JQ928797)作外类群，将本研究获得的沙果小食心虫和梨小食心虫COⅠ基因单倍型序列和1个已登录的沙果小食心虫COⅠ基因序列(GenBank登录号：AB603522)，以及3个已登录的梨小食心虫COⅠ基因序列(GenBank登录号：AB603521、FJ449660、HQ700338)构建系统发育树，结果(图3)显示，沙果小食心虫和梨小食心虫的COⅠ基因单倍型各自聚在一支。

以苹果蠹蛾(C.pomonella)线粒体全基因(GenBank登录号：JX407107)中COⅡ基因序列作外类群，将本研究获得的2种食心虫COⅡ基因单倍型序列和已登录的相应沙果小食心虫COⅡ基因序列(GenBank登录号：JQ404438)和梨小食心虫COⅡ基因序列(GenBank登录号：EU671060)构建系统发育树，结果(图4)表明，所有的沙果小食心虫COⅡ基因单倍型聚为一支，所有的梨小食心虫单倍型聚为另外一支。

线粒体COⅠ和COⅡ基因序列目前已广泛地应用到昆虫近缘种的鉴定之中[20]。

在利用线粒体基因进行分子鉴定时，不同物种间相同的线粒体基因核苷酸序列差异是最重要的标准[21-22]。

Hebert等[23]通过分析比较200种鳞翅目昆虫COⅠ基因序列后发现，COⅠ基因核苷酸序列3%的差异可以作为鉴定不同物种的标准，即不同物种COⅠ基因差异大于3%，同一物种COⅠ基因序列差异小于3%，目前
大量的相关研究证实了这一标准的准确性[24-25]。

本研究中，2种食心虫COⅠ和COⅡ基因种内核苷酸差异小于1.0%，而种间核苷酸差异大于4.4%，因此本研究的COⅠ和COⅡ基因序列可以准确地应用于沙果小食心虫和梨小食心虫的分子鉴定中。

Aliabadian等[26]研究指出，当种间遗传距离和种内遗传距离不存在重叠区域，且分离程度较高时，目标片段能准确鉴定物种。

本研究结果显示，基于COⅠ和COⅡ基因的沙果小食心虫和梨小食心虫的种内遗传距离均小于0.010，而种间遗传距离大于0.046，基于COⅠ和COⅡ基因的种内遗传距离和种间遗传距离不存在重叠区域，并且种内遗传距离与种间遗传距离差异显著(P<0.001)。

另外，有分析认为，种间平均遗传距离一般为种内平均遗传距离的10倍以上[27]，本研究基于COⅠ基因的种间平均遗传距离(0.050 0)为种内平均遗传距离(0.002 7)的18.2倍，基于COⅡ基因的种间平均遗传距离(0.050 6)为种内平均遗传距离(0.003 8)的13.2倍。

因此，基于COⅠ和COⅡ基因的沙果小食心虫和的梨小食心虫种间和种内遗传距离分析结果表明，COⅠ和COⅡ基因可以用于这2种食心虫的分子鉴定。

另外，稳定的种间差异位点是确定种的重要手段。

朱振华等[28]通过比较6种果实蝇mtCytb基因序列，发现种间存在30个比较稳定的变异位点，并以此作为鉴定不同种的标准。

本试验发现，2种食心虫在COⅠ基因序列上有30个稳定变异位点，COⅡ基因序列上有26个稳定差异位点，这些位点的碱基及其变化可以作为梨小食心虫和沙果小食心虫分子鉴定时的重要参考。

本研究结果表明，沙果小食心虫、梨小食心虫的COⅠ和COⅡ基因种内的相似度在99.0%以上，种间相似度在95.6%以下，且2个基因的种间遗传距离显著高于种内遗传距离，2种食心虫的COⅠ基因单倍型和COⅡ基因单倍型分别在系统发育树上明显各自聚为一支，因此，本研究采用的沙果小食心虫、梨小食心虫的COⅠ和COⅡ基因可以用于这2种害虫的分子鉴定之中。

Yan S C，Liu Y Q，
Li M W．Grapholita dimorpha:A new record pest damaging fruit trees of C hina [J].Forest Pest and Disease，1999(4)：15-16．(in Chinese)
Chen M X，Luo Y Q，Zhao C J，
et al．Research advance on Grapholita molesta Busck [J].Northern Horticul ture，2009(8)：144-147．(in Chinese)
Xiao J H，Xiao H，Huang D W．DNA barcoding：
New approach of biological taxonomy [J].Acta Zoologica Sinica，2004，
50(5)：852-855．(in Chinese)
Mo B H，Qu L，Han S，et al．DNA barcoding identification：Ⅰ．Research progress and applied perspective of DNA barcoding [J].Sichu an Journal of Zoology，2008，27(2)：303-306．(in Chinese)
Niu Y D，Li M，Wei F W，
et al．Reliability of mtDNA as molecular marker and its perspective [J].Her editas，2001，23(6)：593-598．(in Chinese)
Chang H，Hao D J，Xiao R T，
et al．DNA barcoding based on the mitochondrial COⅠ gene sequences f
or Ips species ( Coleoptera: Scolytidae) [J].Acta Entomologica Sinica，2012，55(9)：1075-1081．(in Chinese)
You Z H，Lu H，Zhang X S，
et al．Molecular identification of the introduced western flower thrips,Fran kliniella occidentals (Pergande) and other eight common thrips species (Th ysanoptera: Thripidae) [J].Acta Entomologica Sinica，2007，50(7)：720-726．(in Chinese)
Wang R J，Wan H，Lei G C，
et al．Phylogenetic analysis of Polyura in China inferred from mitochondria l COⅡ sequences (Lepidoptera:Nymphalidae) [J].Acta Entomologica Sinica，2004，47(2)：243-247．(in Chinese)
Yang F L，Li X D，Yan Z T，
et al．The molecular identification markers of Anopheles sinensis [J].Journ
al of Chongqing Normal University：Natural Science Edition，2014，31(4)：40-44．(in Chinese)
Li Q Q，Li D Y，Duan Y Q，
et al．Application of DNA barcoding in lepidopteran insects [J].Chinese Bul letin of Life Sciences，2010(4)：307-312．(in Chinese)
Zhu Z H，Ye H，
Zhang Z Y．Molecular identification of six Bactrocera species (Diptera：Tephritidae) based on mtDNA [J].Acta Entomologica Sinica，2005，48(3)：386-390．(in Chinese)。