基于種特異性COI標記的西花薊馬溫室品系和羽扇豆品系一步雙重快速鑒定技術(shù)

張蓉，王玉生，楊麗梅，萬方浩,2，張桂芬,2

基于種特異性COI標記的西花薊馬溫室品系和羽扇豆品系一步雙重快速鑒定技術(shù)

張蓉1，王玉生1，楊麗梅3，萬方浩1,2，張桂芬1,2

（1中國農(nóng)業(yè)科學(xué)院植物保護研究所植物病蟲害生物學(xué)國家重點實驗室/農(nóng)業(yè)農(nóng)村部作物有害生物綜合治理重點實驗室，北京 100193；2農(nóng)業(yè)農(nóng)村部外來入侵生物預(yù)防與控制研究中心，北京 100193；3北京市大興區(qū)榆垡鎮(zhèn)農(nóng)業(yè)技術(shù)推廣站，北京 102602）

【】西花薊馬（）是世界性農(nóng)業(yè)和園藝作物害蟲，也是我國重要的對外對內(nèi)檢疫對象。西花薊馬具有兩個品系，即溫室品系（greenhouse race，GR）和羽扇豆品系（lupin race，LR），二者在寄主范圍、抗藥性、生存環(huán)境等方面均具有明顯差異，但外部形態(tài)相似，無法準確鑒別?！尽恳晕骰ㄋE馬溫室品系和羽扇豆品系為靶標，以田間常見的其他14種薊馬為參照，采用基于線粒體DNA細胞色素C氧化酶亞基I（mitochondrial DNA cytochrome c oxidase subunit I，mtDNA COI）基因的種特異性（species-specific COI，SS-COI）PCR法，研究一步雙重品系快速檢測技術(shù)。【】采用mtDNA COI基因通用引物獲得西花薊馬溫室品系和羽扇豆品系以及其他常見薊馬的COI序列，根據(jù)測序結(jié)果以及數(shù)據(jù)庫已公開數(shù)據(jù)，設(shè)計篩選獲得可同時擴增兩個品系的特異性組合引物（TF6/GR32/LR12），其中TF6為品系通用上游引物，GR32為溫室品系特異性下游引物，LR12為羽扇豆品系特異性下游引物，其擴增片段大小分別為溫室品系362 bp、羽扇豆品系541 bp；對組合引物比例和退火溫度進行優(yōu)化；同時，對組合引物的種/品系特異性和靈敏性/檢測閾值進行檢驗?！尽慨?dāng)組合引物TF6/GR32/LR12比例為1.0/0.2/0.8、退火溫度為44℃時，擴增效果最好。種/品系特異性檢測結(jié)果證實，該檢測技術(shù)僅對西花薊馬溫室品系和羽扇豆品系具有擴增能力，對田間常見的其他14種薊馬包括花薊馬（）、禾花薊馬（）、黃胸薊馬（）、大薊馬（）、棕櫚薊馬（）、煙薊馬（）、普通大薊馬（）、豆喙薊馬（）、草木樨近絹薊馬（）、蔗腹齒薊馬（）、稻簡管薊馬（）、榕端寬管薊馬（）、榕母管薊馬（）和橫紋薊馬（）等不具有擴增能力。靈敏性檢測結(jié)果顯示，該組合引物不僅對單一品系或混合品系的成蟲具有良好的擴增效果，對單粒卵、1齡和2齡幼蟲以及預(yù)蛹和蛹均具有同樣的擴增效能；對溫室品系的最低檢測閾值為35.90 pg·μL-1，相當(dāng)于1/10 240頭雌成蟲，對羽扇豆品系的最低檢測閾值為146.95 pg·μL-1，相當(dāng)于1/2 560頭雌成蟲；同時，對來自不同地域的同一品系不同單倍型的成蟲亦具有良好的擴增效能?！尽拷⒌募夹g(shù)體系完全可以用于西花薊馬品系的快速鑒定及檢疫監(jiān)管，對有效阻截其進一步傳播擴散及靶向防控措施制定意義重大。

西花薊馬溫室品系；西花薊馬羽扇豆品系；一步雙重PCR；SS-COI標記；品系特異性組合引物；快速鑒定

0 引言

【研究意義】西花薊馬（）又稱苜蓿薊馬，屬纓翅目（Thysanoptera）薊馬科（Thripidae）花薊馬屬（），是一種世界性檢疫性害蟲，也是我國重要的對外對內(nèi)檢疫對象。西花薊馬包括兩個品系，即溫室品系（greenhouse race，GR）和羽扇豆品系（lupin race，LR），二者在寄主植物范圍[1]、生存環(huán)境[2]、抗藥性[3]等方面均有明顯差異。相較于羽扇豆品系，溫室品系具有更廣的寄主植物范圍（溫室品系60余科，500余種[4-6]；羽扇豆品系18科，20余種[1,7-11]）、更強的殺蟲劑抗性[3]，更喜歡在溫暖干燥的環(huán)境下生存[2]，以及更強的繁殖能力和更快的發(fā)育速率[12]。西花薊馬溫室品系和羽扇豆品系的快速檢測鑒定，對西花薊馬的檢疫監(jiān)管及其靶向性防控具有重要意義?！厩叭搜芯窟M展】西花薊馬原產(chǎn)于美國西北部[13]，是美國加州常見害蟲，之后隨花卉產(chǎn)業(yè)發(fā)展以及國際貿(mào)易活動迅速擴散，目前，在以色列[14]、南非[15]、馬來西亞[16]、日本[17]、澳大利亞[18]、韓國[19]、印度[20]等均有發(fā)生的報道。西花薊馬為多食性害蟲，寄主植物有60余科500多種，包括蔬菜、果樹、花卉及園林植物等，既可直接銼吸危害，亦可通過傳播植物病毒間接危害，嚴重影響產(chǎn)品的質(zhì)量和產(chǎn)量以及花卉園林植物的觀賞價值，已成為一種世界性農(nóng)業(yè)和園藝作物害蟲[21-22]。2000年中國昆明國際花卉節(jié)上，我國首次在來自緬甸的盆景植物上檢疫到了西花薊馬[23]，2003年在北京市郊的溫室大棚中首次發(fā)現(xiàn)其嚴重危害辣椒[4]，目前已在我國北京、山東[24]、貴州[25]、湖南[26]、江蘇[27]、新疆[28]、西藏[29]、四川[30]及河北、山西、內(nèi)蒙古、遼寧、吉林、黑龍江、河南、云南、陜西、甘肅、寧夏[9,31]19個?。ㄗ灾螀^(qū)、直轄市）發(fā)生危害。其中溫室品系在上述19個?。ㄗ灾螀^(qū)、直轄市）發(fā)生危害，而羽扇豆品系目前僅在北京、山東、遼寧、吉林、黑龍江、云南、陜西、寧夏、甘肅9個?。ㄗ灾螀^(qū)、直轄市）有發(fā)現(xiàn)[8-9,31]。西花薊馬個體微小，雌成蟲體長1.4—1.7 mm，雄成蟲體長1—1.15 mm，且體色多變[21,32]。有研究認為，深色型的西花薊馬更有可能是羽扇豆品系[1]。然而，鑒于煙薊馬（）的體色變化與溫度有關(guān)，溫度越低、體色越深[33]；因此，多數(shù)學(xué)者認為溫度應(yīng)該是昆蟲體色變化的主要誘因[34]。顯然，基于成蟲外部特征的傳統(tǒng)形態(tài)學(xué)鑒定法無法滿足植物檢疫和檢測監(jiān)測的基本需求，尤其當(dāng)樣本為幼體或殘體時?；趍tDNA COI[1-2,8-9,11,35]和COII基因[11]、rDNA 28S D2基因[1,8]，以及基于核蛋白編碼基因[1]等標記技術(shù)，曾用于西花薊馬溫室品系和羽扇豆品系的識別鑒定?！颈狙芯壳腥朦c】然而上述標記技術(shù)，或鑒定流程復(fù)雜、成本較高，或鑒定結(jié)果不穩(wěn)定、重復(fù)性差。如基于COI和COII基因標記的鑒定技術(shù)，需要對PCR產(chǎn)物進行堿基序列測定以及系統(tǒng)發(fā)育分析[1-2,24,29-31]；基于rDNA 28S D2基因的品系特異性檢測，雖然無需堿基序列測定和系統(tǒng)發(fā)育分析[1]，但實際應(yīng)用中卻發(fā)現(xiàn)重復(fù)性差，無法擴增出理想的品系特異性靶標片段。【擬解決的關(guān)鍵問題】以西花薊馬溫室品系和羽扇豆品系為靶標，以田間常見的其他14種薊馬為參照，采用基于線粒體DNA細胞色素C氧化酶亞基I（mitochondrial DNA cytochrome c oxidase subunit I，mtDNA COI）基因的種特異性（species-specific COI，SS-COI）PCR法，研究建立西花薊馬品系特異性一步雙重快速檢測技術(shù)體系。

1 材料與方法

試驗于2017年7月至2019年2月在中國農(nóng)業(yè)科學(xué)院植物保護研究所完成。

1.1 供試蟲源

西花薊馬溫室品系和羽扇豆品系為室內(nèi)飼養(yǎng)種群，田間常見的其他14種薊馬包括花薊馬（Frankliniella intonsa）、禾花薊馬（Frankliniella tenuicornis）、黃胸薊馬（Thrips hawaiiensis）、大薊馬（Thrips major）、棕櫚薊馬（Thrips palmi）、煙薊馬、普通大薊馬（Megalurothrips usitatus）、豆喙薊馬（Mycterothrips glycines）、草木樨近絹薊馬（Sussericothrips melilotus）、蔗腹齒薊馬（Fulmekiola serrata）、稻簡管薊馬（Haplothrips aculeatus）、榕端寬管薊馬（Mesothrips jordani）、榕母管薊馬（Gynaikothrips ficorum）和橫紋薊馬（Aeolothrips fasciatus）均采自田間，具體信息詳見表1。

1.2 DNA提取

以西花薊馬溫室品系和羽扇豆品系以及田間常見的其他14種薊馬為對象，取單頭雌成蟲，參照喬瑋娜等[36]的方法提取DNA，以20 μL超純水復(fù)溶，充分溶解后以超微量分光光度計（NanoPhotometerTMP330，Implen，Munich，Germany）測定DNA濃度，-20℃保存?zhèn)溆谩?/p>

1.3 COI基因擴增與序列測定

分別以西花薊馬溫室品系、羽扇豆品系，以及上述14種其他薊馬共16種/品系薊馬DNA為模板，以DNA條形碼通用型引物L(fēng)CO1490（GGTCAACAAAT CATAAAGATATTGG），HCO2198（TAAACTTCAGG GTGACCAAAAAATCA）[37]進行PCR擴增，擴增片段大小約為710 bp。PCR反應(yīng)體系為25 μL，其中DNA模板2 μL，上游引物和下游引物各0.5 μL，10×Buffer 2.5 μL，dNTPs（2.5 mmol·L-1）0.5 μL，TaqDNA聚合酶（5 U·μL-1）0.2 μL，ddH2O補足至25 μL。擴增條件：94℃預(yù)變性5 min；94℃ 30 s，51℃ 30 s，72℃ 1 min，35個循環(huán)；72℃延伸5 min。取5 μL擴增產(chǎn)物，加1 μL上樣緩沖液，以DNA Marker為參照，在含有染色劑Gel stain（北京全式金生物技術(shù)有限公司）的1%瓊脂糖凝膠上進行電泳分離25 min（150 V），以凝膠成像系統(tǒng)（Gel DocTMXR+Imaging System Universal Hood II，Bio-Rad，USA）分析檢測結(jié)果；每種薊馬分別檢測10頭。電泳檢測質(zhì)量合格的PCR產(chǎn)物，直接送生工生物工程（上海）股份有限公司進行雙向堿基序列測定。

表1 薊馬樣本信息

1.4 一步雙重品系特異性引物設(shè)計、篩選與體系優(yōu)化

DNAMAN（DNAMAN Version 9）比對分析顯示，西花薊馬兩個品系的堿基序列相似度為96.67%，堿基差異位點19個。以本研究獲得的其他14種薊馬以及數(shù)據(jù)庫中已公開的薊馬COI基因序列（包括西花薊馬所有單倍型）為參照，在兩個品系所有單倍型堿基序列相似度比較高的區(qū)段設(shè)計通用上游引物和通用下游引物各1條；然后，在兩個品系所有單倍型堿基序列差異位點比較集中且差異位點在同一品系不同單倍型間較保守的區(qū)段，分別靶向溫室品系和羽扇豆品系，設(shè)計6條品系特異性上游引物和3條品系特異性下游引物。即品系通用引物共計2條，其中品系通用上游引物1條，品系通用下游引物1條，分別命名為TF6和TR3。品系特異性引物共計18條，其中溫室品系特異性上游引物6條、下游引物3條，分別命名為GF1-GF6和GR1-GR3；羽扇豆品系特異性上游引物6條、下游引物3條，分別命名為LF1-LF6和LR1-LR3；并據(jù)此構(gòu)建一步雙重PCR反應(yīng)體系。體系構(gòu)建時引物組合方式為，1條通用上游引物與2條特異性下游引物（溫室品系和羽扇豆品系各1條）或1條通用下游引物與2條特異性上游引物（溫室品系和羽扇豆品系各1條）組合，共計獲得組合45種。其中30種組合對兩個品系的擴增片段大小基本一致，其混合擴增產(chǎn)物無法在凝膠電泳上有效分離，不符合一步雙重PCR快速檢測的基本要求[38]，予以舍棄。對其余15種引物組合進行品系特異性預(yù)實驗驗證和引物優(yōu)化（基于特異性引物設(shè)計中GC含量需＞40%的原則，進行AT→GC堿基轉(zhuǎn)換。優(yōu)化后品系特異性下游引物，溫室品系命名為GR32，羽扇豆品系命名為LR12；優(yōu)化后品系特異性上游引物，溫室品系命名為GF62，羽扇豆品系命名為LF52）以及再驗證，對優(yōu)化后依然不能進行品系特異性擴增的引物組合予以舍棄，最終確定TF6/GR32/LR12為西花薊馬品系特異性一步雙重PCR檢測的最佳引物組合，其對溫室品系和羽扇豆品系的特異性擴增片段大小分別為362和541 bp（表2）。然后，以該組合引物為靶標進行引物比例和退火溫度優(yōu)化。

表2 西花薊馬品系特異性引物序列及擴增產(chǎn)物

1.5 一步雙重品系特異性PCR檢測體系的建立與優(yōu)化

為確定一步雙重品系特異性PCR檢測的最佳反應(yīng)體系，對組合引物的比例和退火溫度進行優(yōu)化研究。PCR反應(yīng)體系為25 μL，引物組合共計11種，即TF6/GR32/LR12分別以1.0/0/1.0、1.0/0.1/0.9、1.0/0.2/0.8、1.0/0.3/0.7、1.0/0.4/0.6、1.0/0.5/0.5、1.0/0.6/0.4、1.0/0.7/0.3、1.0/0.8/0.2、1.0/0.9/0.1、1.0/1.0/0 μL的比例進行組合（引物濃度為10 μmol·L-1），10×Buffer 2.5 μL，dNTPs（2.5 mmol·L-1）0.5 μL，Taq DNA聚合酶（5 U·μL-1）0.2 μL，DNA模板1 μL，ddH2O補足至25 μL。擴增條件：94℃預(yù)變性5 min；94℃ 30 s，梯度退火溫度（36—54℃）1 min，72℃ 1 min，35個循環(huán)；72℃延伸5 min。電泳檢測方法同1.3。

1.6 品系特異性組合引物的種/品系特異性檢測

分別以西花薊馬溫室品系和羽扇豆品系以及其他14種田間常見薊馬DNA為模板，以優(yōu)化的反應(yīng)體系進行PCR擴增，檢測組合引物的種/品系特異性。PCR反應(yīng)體系為25 μL，其中10×Buffer 2.5 μL，dNTPs（2.5 mmol·L-1）0.5 μL，最佳配比組合引物TF6/GR32/LR12共計2.0 μL，Taq DNA聚合酶（5 U·μL-1）0.2 μL，DNA模板1 μL，ddH2O補足至25 μL。擴增條件：94℃ 5 min；94℃ 30 s，最佳退火溫度 1 min，72℃ 1 min，35個循環(huán)；72℃延伸5 min。電泳檢測方法同1.3。每種/品系薊馬分別檢測10頭。

1.7 品系特異性組合引物對不同蟲態(tài)的擴增效果

以品系特異性組合引物為靶標，以提取的單一品系不同蟲態(tài)（卵、1齡和2齡幼蟲、預(yù)蛹和蛹）與不同性別（雄成蟲、雌成蟲）DNA，以及混合品系（1﹕1）各蟲態(tài)/性別DNA為模板，進行PCR擴增，電泳檢測方法同1.3。每種性別或蟲態(tài)分別檢測10頭。

1.8 品系特異性組合引物的檢測閾值

分別取兩個品系單頭雌成蟲DNA模板10 μL，以2倍遞減梯度稀釋至2-12倍，其中溫室品系初始DNA模板濃度為（18.38±5.40）×103pg·μL-1，最終稀釋濃度為4.49 pg·μL-1；羽扇豆品系初始DNA模板濃度為（18.81±5.71）×103pg·μL-1，最終稀釋濃度為4.59 pg·μL-1。分別取1 μL作為DNA模板進行PCR擴增，電泳檢測方法同1.3。分別檢測10頭。

1.9 不同地理種群西花薊馬的檢測

以采自我國寧夏、陜西、山東、甘肅、吉林、云南、北京、河北8個省（自治區(qū)、直轄市）10個不同地區(qū)的西花薊馬（表3）成蟲DNA（1 μL）為模板，以品系特異性組合引物TF6/GR32/LR12進行PCR擴增，電泳檢測方法同1.3。根據(jù)標本采集獲得情況，每個地理種群分別檢測2—25頭。電泳檢測質(zhì)量合格的PCR產(chǎn)物，送生工生物工程（上海）股份有限公司進行堿基序列測定。然后，用DnaSP 5.10軟件分析不同地理種群西花薊馬的單倍型組成。

2 結(jié)果

2.1 薊馬類昆蟲線粒體DNA COI基因序列擴增與分析

以西花薊馬溫室品系和羽扇豆品系及其他14種田間常見薊馬DNA為模板，以DNA條形碼通用型引物 LCO1490/HCO2198進行PCR擴增。電泳檢測結(jié)果顯示，16種/品系薊馬均能擴增出一條清晰的靶標片段（圖1），擴增片段大小為629—681 bp，與數(shù)據(jù)庫中對應(yīng)薊馬種類的堿基序列同源性為98.7%—100%。表明所有薊馬種類的mtDNA提取無誤，可以用于后續(xù)PCR擴增反應(yīng)。

表3 不同地理種群西花薊馬樣本信息

2.2 一步雙重品系特異性PCR檢測體系建立

2.2.1 組合引物比例的優(yōu)化 以篩選出的品系通用上游引物TF6以及溫室品系特異性下游引物GR32和羽扇豆品系特異性下游引物L(fēng)R12為靶標，設(shè)置不同比例的引物組合，以兩個品系DNA為模板分別進行PCR擴增。結(jié)果顯示，當(dāng)組合引物TF6/GR32/LR12以1.0/0.2/0.8的比例進行組合時，無論是溫室品系（圖2-A）還是羽扇豆品系（圖2-B）均能擴增出清晰的靶標條帶。

2.2.2 退火溫度的優(yōu)化 分別以兩個品系DNA為模板，以優(yōu)化后的品系特異性引物組合（TF6/GR32/LR12配比為1.0/0.2/0.8），比較36—54℃不同退火溫度下的擴增效果。結(jié)果顯示，當(dāng)退火溫度為44℃時，擴增效果最好，擴增出的分別靶向溫室品系（圖3-A）或羽扇豆品系（圖3-B）的特異性條帶最為清晰。

M：標準DNA分子質(zhì)量 DNA ladder；1：西花薊馬溫室品系F. occidentalis greenhouse race；2：西花薊馬羽扇豆品系F. occidentalis lupin race；3：花薊馬F. intonsa；4：禾花薊馬F. tenuicornis；5：黃胸薊馬T. hawaiiensis；6：大薊馬T. major；7：棕櫚薊馬T. palmi；8：煙薊馬T. tabaci；9：普通大薊馬M. usitatus；10：豆喙薊馬M. glycines；11：草木樨近絹薊馬S. melilotus；12：蔗腹齒薊馬F. serrata；13：稻簡管薊馬H. aculeatus；14：榕端寬管薊馬M. jordani；15：榕母管薊馬G. ficorum；16：橫紋薊馬A. fasciatus；17：陰性對照（模板為超純水）negative control (ultrapure water)

M：標準DNA分子質(zhì)量 DNA ladder；1—11：TF6/GR32/LR12比例分別為The ratios of TF6/GR32/LR12 are 1.0/0/1.0, 1.0/0.1/0.9, 1.0/0.2/0.8, 1.0/0.3/0.7, 1.0/0.4/0.6, 1.0/0.5/0.5, 1.0/0.6/0.4, 1.0/0.7/0.3, 1.0/0.8/0.2, 1.0/0.9/0.1, and 1.0/1.0/0 μL, respectively；12：陰性對照（模板為超純水）negative control (ultrapure water)

M：標準DNA分子質(zhì)量 DNA ladder；1—10：退火溫度分別為 Annealing temperatures are 36, 38, 40, 42, 44, 46, 48, 50, 52, 54℃, respectively；11：陰性對照（模板為超純水）negative control (ultrapure water)

2.3 組合引物的種/品系特異性檢驗

以西花薊馬溫室品系和羽扇豆品系以及田間常見的其他14種薊馬DNA為模板，以優(yōu)化的一步雙重PCR反應(yīng)體系，檢驗組合引物TF6/GR32/LR12的種/品系特異性。電泳檢測結(jié)果顯明，該組合引物的特異性強，當(dāng)以溫室品系為靶標時，只擴增出362 bp的溫室品系特異性條帶，當(dāng)以羽扇豆品系為靶標時，只擴增出541 bp的羽扇豆品系特異性條帶，對田間常見的其他14種薊馬沒有擴增效果（圖4），表明該組合引物為西花薊馬的品系特異性引物。

2.4 品系特異性組合引物對不同蟲態(tài)的擴增效果

分別以單一品系不同蟲態(tài)（卵、1齡和2齡幼蟲、預(yù)蛹和蛹）和不同性別（雄成蟲、雌成蟲）DNA，以及混合品系（1﹕1）各蟲態(tài)/性別DNA為模板，以品系特異性組合引物TF6/GR32/LR12進行一步雙重PCR擴增。電泳檢測結(jié)果顯示，當(dāng)以單一品系DNA為模板時，溫室品系各蟲態(tài)/性別均能穩(wěn)定擴增出清晰的362 bp特異性條帶，羽扇豆品系各蟲態(tài)/性別均能穩(wěn)定擴增出清晰的541 bp特異性條帶（圖5-A）；當(dāng)以混合品系DNA為模板時，各蟲態(tài)/性別均能穩(wěn)定地同時擴增出分別靶向溫室品系和羽扇豆品系的特異性條帶（圖5-B）。

2.5 品系特異性組合引物的檢測閾值

分別以不同濃度西花薊馬溫室品系和羽扇豆品系的雌成蟲DNA為模板，進行最低檢測閾值測定。結(jié)果表明，當(dāng)以溫室品系為檢測靶標時，其最低檢測濃度為35.90 pg·μL-1，相當(dāng)于1/10 240頭雌成蟲（圖6-A）；當(dāng)以羽扇豆品系為檢測靶標時，其最低檢測濃度為146.95 pg·μL-1，相當(dāng)于1/2 560頭雌成蟲（圖6-B）。

M：標準DNA分子質(zhì)量DNA ladder；1：西花薊馬溫室品系F. occidentalis greenhouse race；2：西花薊馬羽扇豆品系F. occidentalis lupin race；3：花薊馬F. intonsa；4：禾花薊馬F. tenuicornis；5：黃胸薊馬T. hawaiiensis；6：大薊馬T. major；7：棕櫚薊馬T. palmi；8：煙薊馬T. tabaci；9：普通大薊馬M. usitatus；10：豆喙薊馬M. glycines；11：草木樨近絹薊馬S. melilotus；12：蔗腹齒薊馬F. serrata；13：稻簡管薊馬H. aculeatus；14：榕端寬管薊馬M. jordani；15：榕母管薊馬G. ficorum；16：橫紋薊馬A. fasciatus；17：陰性對照（模板為超純水）negative control (ultrapure water)

A：單一品系Single race；2—8：溫室品系Greenhouse race；9—15：羽扇豆品系lupin race；M：標準DNA分子質(zhì)量DNA ladder；2、9：卵egg；3、10：1齡幼蟲1st instar larva；4、11：2齡幼蟲2nd instar larva；5、12：預(yù)蛹pre-pupa；6、13：蛹pupa；7、14：雄成蟲male adult；8、15：雌成蟲female adult；1、16：陰性對照（模板為超純水）negative control (ultrapure water)。B：混合品系Mixed races；M：標準DNA分子質(zhì)量DNA ladder；1：卵egg；2：1齡幼蟲1st instar larva；3：2齡幼蟲2nd instar larva；4：預(yù)蛹pre-pupa；5：蛹pupa；6：雄成蟲male adult；7：雌成蟲female adult；8：陰性對照（模板為超純水）negative control (ultrapure water)

2.6 不同地理種群西花薊馬的檢測鑒定

以采自我國8個?。ㄗ灾螀^(qū)、直轄市）10個不同地區(qū)的127頭西花薊馬成蟲DNA為模板，以品系特異性組合引物TF6/GR32/LR12進行PCR擴增。檢測結(jié)果顯示，各地理種群均能擴增出或靶向溫室品系（362 bp）、或靶向羽扇豆品系（541 bp）的COI基因片段，部分種群個體（10個地理種群45個個體）電泳檢測結(jié)果如圖7所示。其中，114頭被鑒定為溫室品系，來自全部10個地理種群；13頭被鑒定為羽扇豆品系，分別來自寧夏銀川、陜西寶雞、山東煙臺、吉林四平和吉林長春。進一步的堿基序列測定和比對分析顯示，溫室品系和羽扇豆品系的堿基序列與數(shù)據(jù)庫中對應(yīng)品系堿基序列的同源性分別為98.33%—100%和98.34%—100%；品系單倍型分析結(jié)果表明，114頭溫室品系包含有3種單倍型，13頭羽扇豆品系包含有2種單倍型（表4）。

3 討論

薊馬類為小型昆蟲，其物種鑒別主要依據(jù)成蟲外部形態(tài)特征，對于非專業(yè)從事薊馬分類的人員而言難以實現(xiàn)，而近年來隨著分子生物學(xué)技術(shù)的迅速發(fā)展，這一難題迎刃而解。目前，已有多種DNA標記技術(shù)用于薊馬類昆蟲的識別與鑒定，包括隨機多態(tài)性DNA（random amplified polymorphic DNA，RAPD）標記技術(shù)[39]、基于線粒體DNA COI/COII/COIII基因的DNA序列測序法[40-42]、基于COI和ITS-1/ITS-2序列的限制片段長度多態(tài)性（restriction fragment length polymorphism，RFLP）標記技術(shù)[40,43-45]、特征序列擴增區(qū)域（sequence characterized amplified regions，SCAR）標記技術(shù)[46]、種特異性PCR（species-specific PCR，SSP）擴增技術(shù)[47-48]、實時熒光定量PCR（real-time quantitative PCR，qPCR）技術(shù)[49-51]、DNA條形碼技術(shù)[35-36,52-54]、基因芯片技術(shù)[55]等，并取得了顯著成效。其中，基于線粒體DNA COI基因的DNA條形碼技術(shù)在薊馬類以及其他昆蟲類群的識別鑒定中起到了重要作用[56-57]。

A：溫室品系 Greenhouse race；M：標準DNA分子質(zhì)量 DNA ladder；1—13：18.38×103, 9.19×103, 4.60×103, 2.30×103, 1.15×103, 574.38, 287.19, 143.59, 71.80, 35.90, 17.95, 8.97, and 4.49 pg·μL-1, respectively；14：陰性對照（模板為超純水）negative control (ultrapure water)。B：羽扇豆品系 Lupin race；M：標準DNA分子質(zhì)量 DNA ladder；1—13：18.81×103, 9.41×103, 4.70×103, 2.35×103, 1.18×103, 587.81, 293.91, 146.95, 73.48, 36.74, 18.37, 9.18, and 4.59 pg·μL-1, respectively；14：陰性對照（模板為超純水）negative control (ultrapure water)

M：標準DNA分子質(zhì)量 DNA ladder；1—5：銀川YC；6—10：寶雞BJI；11—15：煙臺YT；16—20：敦煌DH；21、22：四平SP；24—28：長春CC；29—33：臨滄LC；34—38：蘭州LZ；39—43：延慶YQ；44—46：廊坊LF；23、47：陰性對照（模板為超純水）negative control (ultrapure water)

表4 基于mtDNA COI基因的各單倍型在西花薊馬10個地理種群中的分布

線粒體DNA為多拷貝基因，具有分子量小、進化速度快、遺傳自主性和嚴格的母系遺傳等特性，基于線粒體DNA細胞色素C氧化酶亞基Ⅰ（mtDNA COI）基因的DNA條形碼技術(shù)即由此發(fā)展而來[58]，目前已廣泛應(yīng)用于個體識別鑒定、生物分類，以及分子進化、遺傳結(jié)構(gòu)和親緣關(guān)系分析等[56-57]。本研究采用基于線粒體DNA COI基因的SS-COI標記技術(shù)，研發(fā)出了靶向西花薊馬溫室品系和羽扇豆品系的一步雙重特異性組合引物TF6/GR32/LR12及其快速檢測鑒定方法。該組合引物特異性強，只對西花薊馬的溫室品系和羽扇豆品系具有擴增效果，只需一步PCR反應(yīng)，即可獲得分別靶向溫室品系（362 bp）和羽扇豆品系（541 bp）的特異性條帶，對田間常見的14種其他種類薊馬（包括兩種花薊馬屬薊馬）不具有擴增能力（圖4）。同時，該技術(shù)體系的檢測靈敏性高，無論單一品系還是混合品系，既對成蟲具有良好的擴增效果，也對單粒卵、初孵幼蟲、2齡幼蟲、預(yù)蛹及蛹具有同樣的擴增效能（圖5）；靶向溫室品系，其最低檢測閾值為35.90 pg·μL-1，相當(dāng)于1/10 240頭雌成蟲，靶向羽扇豆品系其最低檢測閾值為146.95 pg·μL-1，相當(dāng)于1/2 560頭雌成蟲（圖6）。進一步的田間樣本檢測分析顯示，該技術(shù)體系能快速鑒定出來自不同?。ㄗ灾螀^(qū)、直轄市）的西花薊馬（包括不同的單倍型）所屬的品系類型。此外，數(shù)據(jù)庫同源性比對分析結(jié)果表明，沒有任一同屬（）近緣種薊馬（包括東方花薊馬、花薊馬、首花薊馬、禾花薊馬、佛羅里達花薊馬、褐花薊馬、梳缺花薊馬、澳洲花薊馬、西印花薊馬、單色花薊馬，以及、、、、、、、、等19種）或其他種類的薊馬或昆蟲的DNA序列與西花薊馬溫室品系或羽扇豆品系特有的362 bp（同源性＜88%）或541 bp（同源性＜85%）片段完全一致，顯示該技術(shù)體系完全可用于西花薊馬品系的快速鑒定，以及國際國內(nèi)農(nóng)產(chǎn)品（包括水果、蔬菜和花卉）貿(mào)易或調(diào)運過程中的檢驗檢疫和監(jiān)管。同時，還可用于指導(dǎo)田間西花薊馬品系靶向性防控措施的制定，如靶向抗性比較強的溫室品系單優(yōu)地域種群[3]，應(yīng)少用或不用化學(xué)農(nóng)藥；而靶向品系混發(fā)地域種群，應(yīng)采用劑量控制、合理用藥（輪用、混用等）等措施，延緩抗藥性產(chǎn)生。此外，兩個品系在危害程度、傳播病毒能力等方面是否存在差異，尚不得而知，而本研究的品系特異性快速檢測技術(shù)體系的建立，為單一品系種群建立及其后續(xù)研究提供了技術(shù)支撐。

一步多重PCR（multiplex PCR）又稱多重引物PCR或復(fù)合PCR，是指在同一PCR反應(yīng)體系中添加兩對以上引物，同時擴增出多個核酸片段的PCR反應(yīng)，且每對引物只對應(yīng)一個靶標片段[59]，與常用的RAPD技術(shù)（一對引物對應(yīng)多個擴增產(chǎn)物）[38]全然不同。多重PCR技術(shù)可用于病原微生物、昆蟲、軟體動物等的檢測鑒定，所用的靶標基因主要有COI基因、16S rDNA基因[60-64]和ITS基因[65]等。如基于16S rDNA[60-61]、EF1[66]等基因的一步雙重/多重PCR技術(shù)常用于細菌鑒定；基于病毒保守基因（DNA/RNA）的一步多重PCR技術(shù)常用于5種病毒的快速檢測[67-68]。同時，多重PCR技術(shù)在外來入侵種的檢測監(jiān)測中亦發(fā)揮了重要作用[69]。如Zhang等[70]基于SCAR標記技術(shù)，研究構(gòu)建了靶向煙粉虱（）B隱種和西花薊馬的一步雙重快速檢測鑒定技術(shù)體系，只需一步PCR反應(yīng)，即可將兩種同域發(fā)生的重大入侵種區(qū)分開，無需堿基序列測定和比對分析；Kurata等[61]以COI基因為靶標構(gòu)建一步多重PCR擴增體系，只需一步擴增反應(yīng)即可有效區(qū)分入侵種煙粉虱（MEAM1，MED Q1，MED Q2隱種）和本地種煙粉虱（Asia II 6隱種）；Nakamura等[63]以COI基因為靶標構(gòu)建了一步多重PCR引物組合體系，用以快速鑒定5種斑潛蠅屬（）害蟲（包括、蔥斑潛蠅、南美斑潛蠅、美洲斑潛蠅和三葉草斑潛蠅等），并且針對美洲斑潛蠅存在隱存種或進化支的現(xiàn)象，分別設(shè)計了靶向支系A(chǔ)和L的特異性上游引物，與潛蠅屬通用下游引物組合成新的多重PCR引物組合，用于鑒定美洲斑潛蠅的3個隱存種，取得了預(yù)期效果；Cooke等[62]以COI基因為靶標構(gòu)建多重PCR技術(shù)體系，無需測序即可快速鑒別入侵種福壽螺（）與當(dāng)?shù)貎煞N同屬近緣種蝸牛，即淡水蘋果蝸牛（）和。基于SS-PCR技術(shù)的種特異性一步雙重/多重PCR檢測鑒定技術(shù)，只需一步PCR擴增反應(yīng)，即可根據(jù)靶標片段的大小快速地鑒定物種，無需堿基序列測定與比對分析，具有操作簡便，經(jīng)濟、高效、系統(tǒng)性強等特性，適用于區(qū)分近緣種、隱存種等[38]。

本研究研發(fā)的一步雙重品系特異性SS-COI組合引物TF6/GR32/LR12及其快速檢測技術(shù)體系，提高了西花薊馬溫室品系和羽扇豆品系檢測鑒定的準確性和靈敏度，節(jié)約了檢測時間，在西花薊馬品系檢測監(jiān)測及其靶向性防控中具有重要應(yīng)用價值。將本研究獲得的靶標序列與數(shù)據(jù)庫中已有序列進行同源性比對，結(jié)果顯示，本研究溫室品系的靶標片段與數(shù)據(jù)庫[71]中溫室品系其他單倍型（105種）的同源性為98.33%—100%，與羽扇豆品系（包括不同的單倍型）的同源性為95.03%—96.13%，而與其他近緣種的同源性＜87.68%；本研究羽扇豆品系的靶標片段與數(shù)據(jù)庫中羽扇豆品系其他單倍型（104種）的同源性為98.34%—100%，與溫室品系（包括不同的單倍型）的同源性為96.12%—96.30%，而與其他近緣種的同源性＜84.73%。表明該檢測技術(shù)完全可以用于西花薊馬的溫室品系或羽扇豆品系鑒定，并且在同一品系不同單倍型間具有較好的通用性，在不同品系之間具有良好的特異性（圖7）。不過薊馬類昆蟲種類比較多，且同一種類（如西花薊馬）包含多種單倍型，有一些常見的種類/種群/單倍型本研究沒有涉及，因此該組合引物的特異性，尤其靶向同一品系不同單倍型的通用性以及靶向不同品系不同單倍型的特異性，還需進一步驗證；而有關(guān)該組合引物對溫室品系和羽扇豆品系最低檢測閾值的差異，可能與其擴增的兩個品系靶標片段的大小不同有關(guān)。

4 結(jié)論

采用種特異性SS-COI標記法建立了西花薊馬品系一步雙重特異性快速檢測鑒定技術(shù)體系。該檢測鑒定技術(shù)體系特異性強，僅對西花薊馬的溫室品系和/或羽扇豆品系具有擴增能力，對田間常見的其他14種薊馬不具有擴增效能。靈敏度高，無論單一品系還是混合品系，或同一品系不同單倍型，不僅對成蟲、預(yù)蛹和蛹具有良好的擴增效果，對單粒卵、1齡和2齡幼蟲亦具有同等的擴增能力；對溫室品系的最低檢測閾值為35.90 pg·μL-1，相當(dāng)于1/10 240頭雌成蟲，對羽扇豆品系的最低檢測閾值為146.95 pg·μL-1，相當(dāng)于1/2 560頭雌成蟲。該檢測鑒定技術(shù)體系的建立對西花薊馬的檢疫監(jiān)管、有效阻截西花薊馬兩個品系的進一步傳播擴散及其靶向性防控措施制定具有重要意義。

[1] Rugman-Jones P F, Hoddle M S, Stouthamer R. Nuclear-mitochondrial barcoding exposes the global pest western flower thrips (Thysanoptera: Thripidae) as two sympatric cryptic species in its native California., 2010, 103(3): 877-886.

[2] Brunner P C, Frey J E. Habitat-specific population structure in native western flower thrips(Insecta, Thysanoptera)., 2010, 23(4): 797-804.

[3] Martin N A, Workman P J. Confirmation of a pesticide- resistant strain of western flower thrips in New Zealand//, 1994: 144-148.

[4] 張友軍, 吳青君, 徐寶云, 朱國仁. 危險性外來入侵生物——西花薊馬在北京發(fā)生危害. 植物保護, 2003, 29(4): 58-59.

ZHANG Y J, WU Q J, XU B Y, ZHU G R. Dangerous alien invasive species-western flower thrips occur and damage in Beijing., 2003, 29(4): 58-59. (in Chinese)

[5] LOOMANS A J M, VAN LENTEREN J C. Biological control of thrips pests: a review on thrips parasitoids//, 1995, 95(1): 89-201.

[6] MORITZ G. The biology of thrips is not the biology of their adults: a developmental view//:Australian National Insect Collection, Canberra, 2002: 259-267.

[7] 武曉云, 程曉非, 張仲凱, 桂富榮, 李正躍. 西花薊馬() rDNA ITS2和COI基因5末端序列的克隆與比較分析. 浙江大學(xué)學(xué)報(農(nóng)業(yè)與生命科學(xué)版), 2009, 35(4): 355-364.

WU X Y, CHENG X F, ZHANG Z K, GUI F R, LI Z Y. Cloning and comparative analysis of rDNA ITS2 and 5 terminal sequence of COI gene of., 2009, 35(4): 355-364. (in Chinese)

[8] Yang X M, Sun J T, Xue X F, Li J B, Hong X Y. Invasion genetics of the western flower thrips in China: evidence for genetic bottleneck, hybridization and bridgehead effect., 2012, 7(4): e34567.

[9] 張桂芬, 喬瑋娜, 古君伶, 閔亮, 萬方浩. 我國西花薊馬線粒體DNA-COI基因變異及群體遺傳結(jié)構(gòu)分析. 生物安全學(xué)報, 2014, 23(3): 196-209.

ZHANG G F, QIAO W N, GU J L, MIN L, WAN F H. Genetic variability of mtDNA COI and population structure of(Pergande) in China., 2014, 23(3): 196-209. (in Chinese)

[10] 沈登榮, 宋文菲, 袁盛勇, 田學(xué)軍, 和紹禹, 張宏瑞. 基于mtDNA-COI的云南西花薊馬的遺傳分析. 植物保護, 2014, 40(5): 75-79.

SHEN D R, SONG W F, YUAN S Y, TIAN X J, HE S Y, ZHANG H R. Genetic analysis ofbased on mtDNA-COI in Yunnan, China., 2014, 40(5): 75-79. (in Chinese)

[11] 田虎, 張蓉, 王玉生, 萬方浩, 張桂芬. 基于線粒體I和II基因的5種不同寄主植物西花薊馬種群遺傳多樣性研究. 植物保護, 2018, 44(1): 27-36.

TIAN H, ZHANG R, WANG Y S, WAN F H, ZHANG G F. Genetic diversity among five different host plant populations of(Pergande) based onI andII gene sequences., 2018, 44(1): 27-36. (in Chinese)

[12] Nielsen M C, Teulon D A, Chapman R B, Butler R C, Drayton G M, Philipsen H. Comparison of life history parameters of two(Thysanoptera: Thripidae) strains in New Zealand., 2010, 39(2): 303-311.

[13] Bryan D E, Smith R F. The(Pergande) complex in California., 1956, 10(6): 359-410.

[14] Argaman Q, Klein Z, Ben-Dov Y, Mendel Z.(Thysanoptera: Thripidae), an injurious intruder., 1989, 69(7): 1268-1269.

[15] Giliomee J H. First record of western flower thrips,(Pergande) (Thysanoptera: Thripidae) from South Africa., 1989, 52(1): 179-182.

[16] Fauziah I, Saharan H A.. Asian Vegetable Research and Development Center Publication, 1991: 29-33.

[17] Ohta I. Host plant resistance in Japanese chrysanthemums against(Thysanoptera: Thripidae) during the non-flowering stage., 2002, 37(2): 271-277.

[18] Malipatil M B, Postle A C, Osmelak J A, Hill M, Moran J. First record of(Pergande) in Australia (Thysanoptera: Thripidae)., 1993, 32: 378.

[19] Chung B K, Kang S W, Kwon J H. Chemical control system of(Thysanoptera: Thripidae) in greenhouse eggplant., 2000, 3(1): 1-9.

[20] Tyagi K, Kumar V. First report of western flower thrips,(Pergande) (Thripidae: Thysanoptera) from India- a potential Havoc to Indian agriculture., 2015, 6: 1-3.

[21] Kirk W D J, Terry L I. The spread of the western flower thrips(Pergande)., 2003, 5(4): 301-310.

[22] Wan F H, Yang N W. Invasion and management of agricultural alien insects in China., 2016, 61: 77-98.

[23] 蔣小龍, 白松, 肖樞, 楊碧. 為中國昆明國際花卉節(jié)把關(guān)服務(wù). 植物檢疫, 2001, 15(2): 115-117.

JIANG X L, BAI S, XIAO S, YANG B. Inspection and quarantine service for Kunming International Flower Festival, China., 2001, 15(2): 115-117. (in Chinese)

[24] 鄭長英, 劉云虹, 張乃芹, 趙希麗. 山東省發(fā)現(xiàn)外來入侵有害生物—西花薊馬. 青島農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)版), 2007, 24(3): 172-174.

ZHENG C Y, LIU Y H, ZHANG N Q, ZHAO X L. Invaded insect pest-first reported in Shandong province., 2007, 24(3): 172-174. (in Chinese)

[25] 袁成明, 郅軍銳, 李景柱, 張勇. 貴州省蔬菜薊馬的種類、分布及綜合防治. 湖北農(nóng)業(yè)科學(xué), 2008, 47(12): 1442-1444.

YUAN C M, ZHI J R, LI J Z, ZHANG Y. The distribution and integrated management of vegetables thrips in Guizhou Province., 2008, 47(12): 1442-1444. (in Chinese)

[26] 劉佳, 張林, 盧焰梅, 張宏瑞. 湖南外來入侵害蟲西花薊馬初步調(diào)查. 安徽農(nóng)業(yè)科學(xué), 2010, 38(25): 13800-13801, 13804.

LIU J, ZHANG L, LU Y M, ZHANG H R. Preliminary investigation of invasive pest ofin Hunan Province., 2010, 38(25): 13800-13801, 13804. (in Chinese)

[27] 嚴丹侃, 湯云霞, 賀子義, 孫雷, 王鳴華, 薛曉峰, 范加勤. 南京地區(qū)西花薊馬發(fā)生調(diào)查及其分子檢測. 南京農(nóng)業(yè)大學(xué)學(xué)報, 2010, 33(4): 59-63.

YAN D K, TANG Y X, HE Z Y, SUN L, WANG M H, XUE X F, FAN J Q. Survey in Nanjing and the PCR diagnosis of., 2010, 33(4): 59-63. (in Chinese)

[28] 楊華, 崔元玗, 張升, 孫曉軍. 外來入侵害蟲——西花薊馬在新疆的發(fā)生危害. 新疆農(nóng)業(yè)科學(xué), 2010, 47(11): 2252-2253.

YANG H, CUI Y Y, ZHANG S, SUN X J. The occurrence and damage of the exotic invasive pest: western flower thrip () in Xinjiang., 2010, 47(11): 2252-2253. (in Chinese)

[29] 王海鴻, 雷仲仁, 李雪, 代安國, 陳翰秋. 西藏發(fā)現(xiàn)重要外來入侵害蟲——西花薊馬. 植物保護, 2013, 39(1): 181-183.

WANG H H, LEI Z R, LI X, DAI A G, CHEN H Q. An important invasive pest,, inspected in Tibet., 2013, 39(1): 181-183. (in Chinese)

[30] 董華芳, 周正東, 許延波, 張曉云, 張旭東, 劉永碧. 不同顏色粘蟲板對西昌市草莓西花薊馬的集誘效果比較. 安徽農(nóng)業(yè)科學(xué), 2018, 46(19): 156-158.

DONG H F, ZHOU Z D, XU Y B, ZHANG X Y, ZHANG X D, LIU Y B. Comparison of trapping effect of different color sticky board on strawberryin Xichang., 2018, 46(19): 156-158. (in Chinese)

[31] 田虎. 基于mtDNA、和rDNA多基因條形碼的西花薊馬鑒定及遺傳結(jié)構(gòu)研究[D]. 北京: 中國農(nóng)業(yè)科學(xué)院, 2017.

TIAN H. Molecular identification and genetic structure of(Pergande) based on multiple gene barcoding, mtDNA,and rDNAmarkers[D]. Beijing: Chinese Academy of Agricultural Sciences, 2017. (in Chinese)

[32] 雷仲仁, 問錦曾, 王音. 危險性外來入侵害蟲——西花薊馬的鑒別、危害及防治. 植物保護, 2004, 30(3): 63-66.

LEI Z R, WEN J Z, WANG Y. The identification, prevention and harm of dangerous invasive alien pests—western flower thrips., 2004, 30 (3): 63-66. (in Chinese)

[33] Murai T, Toda S. Variation ofin colour and size//:, 2002: 377-378.

[34] Protas M E, Patel N H. Evolution of coloration patterns., 2008, 24: 425-446.

[35] 田虎, 張蓉, 張金良, 王玉生, 萬方浩, 張桂芬. 入侵種西花薊馬與本地近緣種花薊馬的雙基因鑒定. 中國生物防治學(xué)報, 2017, 33(5): 612-622.

TIAN H, ZHANG R, ZHANG J L, WANG Y S, WAN F H, ZHANG G F. Identification of invasive speciesand native speciesbased on double gene markers., 2017, 33(5): 612-622. (in Chinese)

[36] 喬瑋娜, 萬方浩, 張愛兵, 閔亮, 張桂芬. DNA條形碼技術(shù)在田間常見薊馬種類識別中的應(yīng)用. 昆蟲學(xué)報, 2012, 55(3): 344-356.

QIAO W N, WAN F H, ZHANG A B, MIN L, ZHANG G F. Application of DNA barcoding technology for species identification of common thrips (Insecta: Thysanoptera) in China., 2012, 55(3): 344-356. (in Chinese)

[37] Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochromeoxidase subunit I from diverse metazoan invertebrates., 1994, 3(5): 294-299.

[38] Edwards M C, Gibbs R A. Multiplex PCR: advantages, development, and applications., 1994, 3: S65-S75.

[39] Bayar K, T?rjék O, Kiss E, Gyulai G, Heszky L. Intra- and interspecific molecular polymorphism of thrips species., 2002, 53(3): 317-324.

[40] Brunner P C, Fleming C, Frey J E. A molecular identification key for economically important thrips species (Thysanoptera: Thripidae) using direct sequencing and a PCR-RFLP-based approach., 2002, 4(2): 127-136.

[41] 游中華, 路虹, 張憲省, 馮紀年, 石寶才, 宮亞軍, 黃大衛(wèi). 入侵害蟲西花薊馬及其他8種常見薊馬的分子鑒定. 昆蟲學(xué)報, 2007, 50(7): 720-726.

YOU Z H, LU H, ZHANG X S, FENG J N, SHI B C, GONG Y J, HUANG D W. Molecular identification of the introduced western flower thrips,(Pergande) and other eight common thrips species (Thysanoptera: Thripidae)., 2007, 50(7): 720-726. (in Chinese)

[42] 魏書軍, 馬吉德, 石寶才, 宮亞軍, 劉靜, 康總江, 陳學(xué)新, 路虹. 我國新入侵外來害蟲美洲棘薊馬的外部形態(tài)和分子鑒定. 昆蟲學(xué)報, 2010, 53(6): 715-720.

WEI S J, MA J D, SHI B C, GONG Y J, LIU J, KANG Z J, CHEN X X, LU H. External morphology and molecular identification of the newly found invasive pest(Thysanoptera: Thripidae) in China., 2010, 53(6): 715-720. (in Chinese)

[43] Moritz G, Delker C, Paulsen M, Mound L A, Burgermeister W. Modern methods for identification of Thysanoptera./, 2000, 30(3/4): 591-593.

[44] Toda S, Komazaki S. Identification of thrips species (Thysanoptera: Thripidae) on Japanese fruit trees by polymerase chain reaction and restriction fragment length polymorphism of the ribosomal ITS2 region., 2002, 92(4): 359-363.

[45] Rugman-Jones P F, Hoddle M S, Mound L A, Stouthamer R. Molecular identification key for pest species of(Thysanoptera: Thripidae)., 2006, 99(5): 1813-1819.

[46] 孟祥欽, 閔亮, 萬方浩, 周忠實, 王文凱, 張桂芬, 西花薊馬的SCAR分子檢測技術(shù). 昆蟲學(xué)報, 2010, 53(3): 323-330.

MENG X Q, MIN L, WAN F H, ZHOU Z S, WANG W K, ZHANG G F. SCAR marker for rapid identification of the western flower thrips,(Pergande) (Thysanoptera: Thripidae)., 2010, 53(3): 323-330. (in Chinese)

[47] Asokan R, Krishna Kumar N K, Kumar V, Ranganath H R. Molecular differences in the mitochondrial cytochrome oxidase I (mtCOI) gene and development of a species-specific marker for onion thrips,Lindeman, and melon thrips,Karny (Thysanoptera: Thripidae), vectors of tospoviruses (Bunyaviridae)., 2007, 97(5): 461-470.

[48] 周力兵, 劉忠善, 李春艷, 丁元明. PCR法鑒定西花薊馬. 植物檢疫, 2007, 21(2): 78-81.

ZHOU L B, LIU Z S, LI C Y, DING Y M. Method of the polymerase chain reaction for identification of., 2007, 21(2): 78-81. (in Chinese)

[49] Walsh K, Boonham N, Barker I, Collins D W. Development of a sequence-specific real-time PCR to the melon thrips(Thysan., Thripidae)., 2005, 129(5): 272-279.

[50] Huang K S, Lee S E, Yeh Y, Shen G S, Mei E, Chang C M. Taqman real-time quantitative PCR for identification of western flower thrips () for plant quarantine., 2010, 6(4): 555-557.

[51] 吳霞, 張桂芬, 萬方浩. 基于TaqMan實時熒光定量PCR技術(shù)的西花薊馬快速檢測. 應(yīng)用昆蟲學(xué)報, 2011, 48(3): 497-503.

WU X, ZHANG G F, WAN F H. TaqMan real-time fluorescent quantitative PCR for identification of., 2011, 48(3): 497-503. (in Chinese)

[52] Glover R H, Collins D W, Walsh K, Boonham N. Assessment of loci for DNA barcoding in the genus(Thysanoptera: Thripidae)., 2010, 10(1): 51-59.

[53] Karimi J, Hassani-Kakhki M, Awal M M. Identifying thrips (Insecta: Thysanoptera) using DNA barcodes., 2010, 2(1): 35-41.

[54] Rebijith K B, Asokan R, Krishna V, Ranjitha H H, Krishna Kumar N K, Remamurthy V V. DNA barcoding and elucidation of cryptic diversity in thrips (Thysanoptera)., 2014, 97(4): 1328-1347.

[55] 馮毅, 王莉, 白云峰, 王潔, 馮紀年. 基于COI序列快速鑒定花薊馬的DNA條形碼芯片初探. 生物技術(shù)通報, 2009(8): 169-173.

FENG Y, WANG L, BAI Y F, WANG J, FENG J N. Molecular identification ofbased on COI sequences by DNA barcoding chip., 2009(8): 169-173. (in Chinese)

[56] 趙廣宇, 李虎, 楊海林, 彩萬志. DNA條形碼技術(shù)在昆蟲學(xué)中的應(yīng)用. 植物保護學(xué)報, 2014, 41(2): 129-141.

ZHAO G Y, LI H, YANG H L, CAI W Z. Application of DNA barcoding in entomology: a review., 2014, 41(2): 129-141. (in Chinese)

[57] Taylor H R, Harris W E. An emergent science on the brink of irrelevance: a review of the past 8 years of DNA barcoding., 2012, 12(3): 377-388.

[58] HEBERT P D N, CYWINSKA A, BALL S L, DEWAARD J R. Biological identifications through DNA barcodes.:,2003, 270(1512): 313-321.

[59] Chamberlain J S, Gibbs R A, Ranier J E, Nguyen P N, Thomas CASKEY C. Deletion screening of the Duchenne muscular dystrophy locus via multiplex DNA amplification., 1988, 16(23): 11141-11156.

[60] 郝少東, 陳昱圻, 王進忠, 王合, 陶萬強, 張志勇, 石小玉, 周賽. 多重PCR法區(qū)分棗園兩種菱紋葉蟬及檢測其體內(nèi)棗瘋病植原體. 昆蟲學(xué)報, 2015, 58(3): 264-270.

HAO S D, CHEN Y Q, WANG J Z, WANG H, TAO W Q, ZHANG Z Y, SHI X Y, ZHOU S. Multiplex-PCR for identification of twospecies (Hemiptera: Cicadellidae) in jujube orchards and detection of jujube witches’ broom (JWB) phytoplasma in their bodies., 2015, 58(3): 264-270. (in Chinese)

[61] Kurata A, Fujiwara A, Haruyama N, Tsuchida T. Multiplex PCR method for rapid identification of genetic group and symbiont infection status in(Hemiptera: Aleyrodidae)., 2016, 51(1): 167-172.

[62] Cooke G M, King A G, Miller L, Johnson R N. A rapid molecular method to detect the invasive golden apple snail(Lamarck, 1822)., 2012, 4(3): 591-593.

[63] Nakamura S, Masuda T, Mochizuki A, Konishi K, Tokumaru S, Ueno K, Yamaguchi T. Primer design for identifying economically importantspecies (Diptera: Agromyzidae) by multiplex PCR., 2013, 13(1): 96-102.

[64] Laruson A J, Craig S F, Messer K J, Mackie J A. Rapid and reliable inference of mitochondrial phylogroups amongspecies, an invasive group of ship-fouling species (Bryozoa, Cheilostomata)., 2012, 4(3): 617-619.

[65] Gariepy T D, Kuhlmann U, Haye T, Gillott C, Erlandson M. A single-step multiplex PCR assay for the detection of Europeanspp., parasitoids ofspp., 2005, 15(5): 481-495.

[66] 何宛芹, 付瑤, 魯雯璐, 常小麗, 楊文鈺. 大豆根腐病致病鐮孢菌的多重PCR檢測技術(shù). 植物保護學(xué)報, 2017, 44(4): 609-616.

HE W Q, FU Y, LU W L, CHANG X L, YANG W Y. A multiplex PCR detection technique for the pathogenicspecies causing soybean root rot., 2017, 44(4): 609-616. (in Chinese)

[67] 王爽, 田雨婷, 喬奇, 秦艷紅, 張德勝, 張振臣. 侵染甘薯的菜豆金色花葉病毒屬病毒和甘薯褪綠矮化病毒多重PCR檢測方法的建立與應(yīng)用. 植物保護學(xué)報, 2018, 45(6): 1427-1428.

WANG S, TIAN Y T, QIAO Q, QIN Y H, ZHANG D S, ZHANG Z C. Development and application of multiplex PCR method for detection of sweepoviruses andin sweet potato., 2018, 45(6): 1427-1428. (in Chinese)

[68] 張曼, 韓飛. AIV、NDV和安卡拉病毒多重PCR檢測方法的建立及應(yīng)用. 西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版), 2018, 46(2): 1-6, 14.

ZHANG M, HAN F. Establishment and application of multiplex PCR assay for detecting avian influenza virus, new castle disease virus and angara disease virus., 2018, 46(2): 1-6, 14. (in Chinese)

[69] Kamenova S, Bartley T J, Bohan D A, Boutain J R, Colautti R I, Domaizon I, Fontaine C, Lemainque A, Le Viol I, Mollot G, Perga M E, Ravigne V, Massol F. Invasions toolkit: current methods for tracking the spread and impact of invasive species//, 2017, 56: 85-182.

[70] Zhang G F, Wu X, Zhou Z X, Meng X Q, Wan F H. A one-step, single tube, duplex PCR to detect predation by native predators on invasiveMEAM1 and., 2014, 150(1): 66-73.

[71] Cao L J, Wang Z H, Gong Y J, Zhu L, Hoffmann A A, Wei S J. Low genetic diversity but strong population structure reflects multiple introductions of western flower thrips (Thysanoptera: Thripidae) into China followed by human-mediated spread., 2017, 10(4): 391-401.

One-step duplex rapid identification technique ofgreenhouse and lupin races based on species-specific COI marker

ZHANG Rong1, WANG YuSheng1, YANG LiMei3, WAN FangHao1,2, ZHANG GuiFen1,2

(1State Key Laboratory for Biology of Plant Diseases and Insect Pests/Key Laboratory of Integrated Pest Management of Crop, Ministry of Agriculture and Rural Affairs, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193;2Center for Management of Invasive Alien Species, Ministry of Agriculture and Rural Affairs, Beijing 100193;3Agricultural Technology Extension Station of Yufa, Daxing, Beijing 102602)

【Background】is one of the most important pests infesting agricultural and horticultural crops worldwide, also it is a national- and international-wide important quarantine pest.has two races, i.e., greenhouse race (GR) and lupin race (LR). There are significant differences between the two races in host plant species, pesticide resistance, living environments, etc. However, it is difficult to identify these two races because of their high similarity in morphological characteristics.【】In this study, a one-step duplex PCR procedure for rapid identification ofgreenhouse race and lupin race was developed by using species-specific COI (SS-COI) marker based on mitochondrial DNA cytochrome c oxidase subunit I (mtDNA COI) gene, and other 14 common thrips species in the field were used as the control.【】The COI sequences ofgreenhouse and lupin races, and other 14 thrips species were amplified using mtDNA COI gene universal primers. The race-specific primer combination system, TF6/GR32/LR12 (TF6, the universal forward primer; GR32, the greenhouse race-specific reverse primer; LR12, the lupin race-specific reverse primer), for one-step duplex amplification of the two races was established. The fragments amplified by the combined primers were 362 bp for the greenhouse race and 541 bp for the lupin race, respectively. The system and conditions (e.g., annealing temperature) were optimized. The species/race specificity and sensitivity/ detection threshold of the combined primers were evaluated.【】When the ratio of the combined primers (TF6/GR32/LR12) was 1.0/0.2/0.8 and the annealing temperature was 44℃, the amplification effect was the best. Species/race-specific tests performed with this set of combined primers indicated that allgreenhouse race and lupin race specimens were detected positively and no cross reaction with other 14 thrips species, including,,,,,,,,,,,,,and, was identified. The system was tested on individual male and female adults, individual prepupae and pupae, individual 1st and 2nd instar larvae, and even single egg, and demonstrated to be applicable for all these stages of greenhouse and lupin races of, either in single race or two races mixed together. The detection threshold for greenhouse race was 35.90 pg·μL-1, which was equivalent to 1/10 240 of a whole female adult. The detection threshold for lupin race was 146.95 pg·μL-1, which was equivalent to 1/2 560 of a whole female adult. Meanwhile, the system worked well in testing individuals from different geographic populations with different haplotypes within greenhouse race or lupin race of the thrips. 【】The one-step duplex PCR system developed here for rapid identification ofgreenhouse and lupin races could be useful in port quarantine and interception in national and international transportation of agricultural products, and should be significant in blocking further spreading and developing targeted prevention and control measures of.

greenhouse race;lupin race; one-step duplex PCR; SS-COI marker; race-specific combined primer set; rapid identification

2019-03-27；

2019-06-06

國家自然科學(xué)基金（31572067）、國家重點研發(fā)計劃（2017YFC1200600）

張蓉，E-mail：zhangrong0120@163.com。

張桂芬，E-mail：guifenzhang3@163.com

（責(zé)任編輯 岳梅）