張桂芬，劉萬(wàn)學(xué)，郭建英，呂志創(chuàng)，萬(wàn)方浩，申香菊中國(guó)農(nóng)業(yè)科學(xué)院植物保護(hù)研究所，植物病蟲(chóng)害生物學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京009;農(nóng)業(yè)部外來(lái)入侵生物預(yù)防與控制研究中心，北京0008;北京市順義區(qū)植保植檢站，北京000

美洲斑潛蠅Liriomyza sativae Blanchard，屬雙翅目潛蠅科斑潛蠅屬，是一種世界性檢疫害蟲(chóng)。歐洲和地中海植物保護(hù)組織 (European and Mediterranean Plant Protection Organization，EPPO)將其列為A2 類(lèi)檢疫性有害生物(http:∥www.eppo.org/);我國(guó)農(nóng)業(yè)部于1995年將其列為進(jìn)境地植物檢疫潛在性害蟲(chóng)，并制定了美洲斑潛蠅“九五”監(jiān)測(cè)治理規(guī)劃(陳文龍和張潤(rùn)杰，2003)。美洲斑潛蠅為多食性害蟲(chóng)，寄主廣泛，如僅在北京地區(qū)其寄主已涉及20科130多種植物(雷仲仁和王音，2005)。

美洲斑潛蠅原產(chǎn)于巴西，20世紀(jì)40年代以來(lái)陸續(xù)暴發(fā)于佛羅里達(dá)、夏威夷等地，70年代后又于大洋洲的一些島嶼和阿拉伯半島南部地區(qū)發(fā)現(xiàn)該蟲(chóng)。多年來(lái)該蟲(chóng)隨寄主植物及其他介質(zhì)在世界范圍內(nèi)傳播蔓延，現(xiàn)已在40多個(gè)國(guó)家和地區(qū)嚴(yán)重發(fā)生(EPPO/CABI，2011)。由于美洲斑潛蠅具有繁殖能力強(qiáng)，發(fā)育周期短，寄主范圍廣，能傳播多種植物病毒，可在較短時(shí)間暴發(fā)成災(zāi)等特點(diǎn)，許多國(guó)家將其列為植物檢疫一類(lèi)危險(xiǎn)性害蟲(chóng)(Deeming，1992;EPPO/CABI，2011;Zitter ＆ Tsai，1980)。我國(guó)于1994年初在海南省三亞市蔬菜上首次發(fā)現(xiàn)美洲斑潛蠅，同年底即在華南、華中等地普遍發(fā)生(問(wèn)錦曾等，1996)，1997年蔓延至全國(guó)25個(gè)省/直轄市/自治區(qū);到目前為止，該蟲(chóng)已廣泛分布于我國(guó)29個(gè)省/直轄市/自治區(qū)，對(duì)瓜果、蔬菜、煙草、棉花等經(jīng)濟(jì)作物和花卉造成嚴(yán)重危害，已成為我國(guó)農(nóng)業(yè)生產(chǎn)上的一個(gè)突出問(wèn)題(陳兵等，2002;范文中和王中武，2005;胡長(zhǎng)效和蘇新林，2003;雷仲仁和王音，2005;王軍等，1999)。因此，加強(qiáng)對(duì)美洲斑潛蠅的檢疫和監(jiān)測(cè)成為維護(hù)我國(guó)對(duì)外貿(mào)易信譽(yù)，保障人民健康生活的必要前提。美洲斑潛蠅以卵和幼蟲(chóng)藏匿在寄主植物組織內(nèi)，隨切花、切條、帶葉的瓜、果、豆、菜、其他寄主葉片以及盆景植物進(jìn)行遠(yuǎn)距離傳播;蛹隨寄主植株土壤傳播(EPPO/CABI，2011)。因此，要加強(qiáng)對(duì)美洲斑潛蠅的檢疫與監(jiān)測(cè)必需建立一套快速準(zhǔn)確的分子檢測(cè)技術(shù)。

目前，世界上的斑潛蠅屬昆蟲(chóng)約有300種，但真正危害農(nóng)作物的僅有23種(Parrella，1987)，我國(guó)已知有14種(楊龍龍，1995)。斑潛蠅屬害蟲(chóng)體型較小(如美洲斑潛蠅成蟲(chóng)體長(zhǎng)1.3 ～2.3 mm，卵0.2～0.3 mm)，且外部形態(tài)極為相似，卵、幼蟲(chóng)和蛹更難以從形態(tài)上加以區(qū)分(附錄)。然而，檢疫工作中截獲的斑潛蠅多為卵、幼蟲(chóng)或蛹。通常斑潛蠅屬害蟲(chóng)幼期階段多因形態(tài)不穩(wěn)定或特征不明顯而難以鑒定分類(lèi)，無(wú)法及時(shí)進(jìn)行種類(lèi)識(shí)別，為檢疫工作帶來(lái)諸多不便(陳萍等，2002;Frick，1952;Sasakawa，1961)，而分子生物學(xué)技術(shù)的發(fā)展與應(yīng)用為檢疫害蟲(chóng)的及時(shí)鑒定提供了快捷而準(zhǔn)確的途徑。本研究以重要檢疫害蟲(chóng)美洲斑潛蠅為靶標(biāo)，以我國(guó)菜田常見(jiàn)潛葉蠅類(lèi)害蟲(chóng)，包括南美斑潛蠅Liriomyza huidobrensis Blanchard、三葉斑潛蠅Liriomyza trifolli Brugess、蔥斑潛蠅Liriomyza chinensis Kato、豌豆?jié)撊~蠅Phytomyza horticola Goureau等為參照，采用種特異性PCR(speciesspecific PCR，SS-PCR)進(jìn)行快速鑒定技術(shù)的研究，旨在為有效阻截美洲斑潛蠅的進(jìn)一步傳播擴(kuò)散，保障國(guó)家貿(mào)易信譽(yù)和農(nóng)業(yè)生產(chǎn)安全提供技術(shù)支撐。

1 材料與方法

1.1 供試蟲(chóng)源

美洲斑潛蠅飼養(yǎng)于中國(guó)農(nóng)業(yè)科學(xué)院植物保護(hù)研究所南區(qū)溫室，寄主植物為菜豆Phaseolus vulgaris L.;南美斑潛蠅采自云南省昆明市呈貢縣蔬菜生產(chǎn)基地，寄主植物為芹菜Apium graveolens L.;三葉斑潛蠅由華南農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院昆蟲(chóng)生態(tài)研究室提供;蔥斑潛蠅采自市售大蔥Allium fistulosum L.，產(chǎn)地為山東省章丘市;豌豆?jié)撊~蠅采自北京市密云縣河南寨鄉(xiāng)，寄主植物為二月蘭Orychophragmus violaceus(L.)Schulz。

1.2 潛葉蠅總DNA的制備

參照王莉萍等(2007)和周志湘等(2007)的方法并稍加改動(dòng)。取美洲斑潛蠅、南美斑潛蠅、三葉斑潛蠅、蔥斑潛蠅和豌豆?jié)撊~蠅各1頭(粒)，置于滴有20 μL DNA 提取緩沖液(pH 8.0，0.05 mol Tris-HCl、0.001 mol EDTA、0.02 mol NaCl、1%SDS)的Parafilm膜上，用0.2 mL的PCR管底部作為勻漿器充分研磨勻漿，勻漿液用微量移液器移入1.5 mL離心管;然后用100 μL提取緩沖液分2次沖洗勻漿器，并移入同一離心管，混勻;向管中加入8 μL(20 mg·mL-1)蛋白酶K，充分混勻后，于60℃水浴1.5 h(中途混勻1次);然后沸水浴8 min，加入220 μL氯仿/異戊醇(v∶v=24∶1)抽提液，輕柔混勻數(shù)10次后，冰浴30 min;4 ℃、12000 r·min-1離心 20 min，取上清液，加入440 μL預(yù)冷無(wú)水乙醇，輕輕混勻，待出現(xiàn)少量絮狀沉淀后于-20℃放置30 min;4℃、12000 r·min-1離心15 min，棄去上清液。加入400 μL 75%預(yù)冷乙醇洗滌，4 ℃、12000 r·min-1離心 15 min，棄去上清液;然后將離心管倒扣于潔凈濾紙上，自然干燥30 min后每管加入30 μL超純水(其中美洲斑潛蠅的蛹以及1～3齡幼蟲(chóng)以20 μL超純水復(fù)溶，卵以10 μL超純水復(fù)溶)，充分溶解后于-20℃保存?zhèn)溆谩?/p>

1.3 美洲斑潛蠅SS-PCR引物的設(shè)計(jì)

在GenBank中調(diào)用已知的美洲斑潛蠅線粒體DNA細(xì)胞色素氧化酶亞基I基因(mitochondrial DNA cytochrome oxidase I，mtDNA COⅠ)序列(GenBank 登錄號(hào)為 EU219613)(Shang et al.，2008)，運(yùn)用 Primer Premier 5軟件將美洲斑潛蠅COⅠ序列與其近緣種南美斑潛蠅、三葉斑潛蠅、蔥斑潛蠅的COⅠ序列進(jìn)行比較分析，選取美洲斑潛蠅與其他潛葉蠅類(lèi)昆蟲(chóng)同源性最低的區(qū)段，根據(jù)該區(qū)段的堿基序列設(shè)計(jì)1對(duì)美洲斑潛蠅特異性引物L(fēng)SCZE1/LSCZF1，上游引物堿基序列為T(mén)CGAGCAGAATTAGGACAT，下游引物堿基序列為ATTGAAGAAAGTGGAGGGT(上海生工生物工程技術(shù)服務(wù)有限公司協(xié)助合成)，引物擴(kuò)增片段大小為294 bp。

1.4 SS-PCR引物的種特異性檢驗(yàn)

分別以單頭南美斑潛蠅、三葉斑潛蠅、蔥斑潛蠅、豌豆?jié)撊~蠅的DNA為模板，美洲斑潛蠅為陽(yáng)性對(duì)照，檢驗(yàn)美洲斑潛蠅特異片段擴(kuò)增引物L(fēng)SCZE1/LSCZF1的種特異性。PCR反應(yīng)體系為20 μL，其中超純水 11.4 μL、10 × 緩沖液 2.0 μL、10 mmol·L-1dNTP 0.4 μL、5 U·μL-1L Taq 聚合酶 0.2 μL(美國(guó)NEB 公司)、20 ng·μL-1上游引物和下游引物各 2.0 μL、模板DNA 2.0 μL。擴(kuò)增程序:94 ℃預(yù)變性3 min;94 ℃ 60 s，68 ℃ 60 s，72 ℃ 60 s，40 個(gè)循環(huán);最后72 ℃延伸5 min。擴(kuò)增反應(yīng)在ABI-9700 PCR基因擴(kuò)增儀上運(yùn)行。取5 μL PCR擴(kuò)增產(chǎn)物在1.5%(重量)瓊脂糖(amresco)凝膠上以80 V電泳(Bio-Rad Power-Pac Basic)分離 50 min，然后以 GelDoc Universal Hood II型凝膠成像系統(tǒng)(Bio-Rad Laboratories)分析電泳結(jié)果。每個(gè)種類(lèi)分別檢測(cè)10頭，并隨機(jī)以其他2種型號(hào)(Bio-Rad iCycler和 ABI Veriti)的PCR儀進(jìn)行驗(yàn)證。

1.5 SS-PCR引物對(duì)不同蟲(chóng)態(tài)的擴(kuò)增效果以及最低檢測(cè)閾值的測(cè)定

分別以提取的不同性別(雌性成蟲(chóng)、雄性成蟲(chóng))與蟲(chóng)態(tài)(卵、1～3齡幼蟲(chóng)和蛹)的美洲斑潛蠅DNA為模板，進(jìn)行靶標(biāo)片段擴(kuò)增效果檢驗(yàn)，每種性別或蟲(chóng)態(tài)分別檢測(cè)10頭。此外，取單頭雌性成蟲(chóng)原模板DNA溶液2 μL(即1/15頭成蟲(chóng))以2倍濃度遞減梯度稀釋至1/30720頭成蟲(chóng)DNA，然后以不同濃度的成蟲(chóng) DNA 模板(即 1/15、1/30、1/60、1/120、1/240、1/480、1/960、1/1920、1/3840、1/7680、1/15360、1/30720頭成蟲(chóng))進(jìn)行最低檢出閾值的測(cè)定，每個(gè)濃度共計(jì)檢測(cè)10頭。

2 結(jié)果與分析

2.1 美洲斑潛蠅SS-PCR引物的種特異性檢驗(yàn)

以美洲斑潛蠅及其近緣種屬潛葉蠅類(lèi)昆蟲(chóng)的DNA為模板，以所設(shè)計(jì)的特異性引物L(fēng)SCZE1/LSCZF1進(jìn)行PCR擴(kuò)增。檢測(cè)結(jié)果顯示，該引物只對(duì)美洲斑潛蠅具有擴(kuò)增能力，對(duì)其他種類(lèi)的潛葉蠅沒(méi)有擴(kuò)增能力，表明該引物為美洲斑潛蠅的特異性引物(圖1)。

圖1 SS-PCR引物L(fēng)SCZE1/LSCZF1對(duì)美洲斑潛蠅及其近緣種屬潛葉蠅的PCR擴(kuò)增結(jié)果Fig.1 Amplification pattern of mitochondrial DNA of L.sativae and its related leafminer species using SS-PCR primers LSCZE1/LSCZF1

2.2 SS-PCR引物對(duì)美洲斑潛蠅不同蟲(chóng)態(tài)的擴(kuò)增效果

以不同性別和蟲(chóng)態(tài)的美洲斑潛蠅DNA為模板，以SS-PCR引物L(fēng)SCZE1/LSCZF1進(jìn)行擴(kuò)增，結(jié)果顯示，不同性別與蟲(chóng)態(tài)的美洲斑潛蠅均能穩(wěn)定擴(kuò)增出294 bp的特異性片段(圖2)。

2.3 SS-PCR引物對(duì)美洲斑潛蠅成蟲(chóng)的最低檢測(cè)閾值

當(dāng)將單頭雌性成蟲(chóng)的DNA模板進(jìn)行遞減梯度稀釋至1/3840頭時(shí)，所有的個(gè)體均能擴(kuò)增出特異性的靶標(biāo)片段;當(dāng)稀釋為1/7680頭時(shí)，仍有60%的個(gè)體可擴(kuò)增出清晰的條帶(圖3)。

圖2 SS-PCR引物L(fēng)SCZE1/LSCZF1對(duì)美洲斑潛蠅不同蟲(chóng)態(tài)和性別的PCR擴(kuò)增結(jié)果Fig.2 Amplification pattern of mitochondrial DNA from different developmental stages and sexes of L.sativae using SS-PCR primers LSCZE1/LSCZF1

圖3 SS-PCR引物L(fēng)SCZE1/LSCZF1對(duì)美洲斑潛蠅的最低檢出閾值Fig.3 The detection efficiency for L.sativae，using a dilution series of mitochondrial DNA isolated from a single adult，using the specifically designed SS-PCR primers LSCZE1/LSCZF1

3 討論

種類(lèi)鑒定是害蟲(chóng)檢疫的首要工作。目前，斑潛蠅屬昆蟲(chóng)的識(shí)別主要依據(jù)成蟲(chóng)的形態(tài)特征(Sasakawa，1961;Shiao et al.，1991;Spencer，1986);成蟲(chóng)前各時(shí)期(卵、幼蟲(chóng)、蛹)的鑒定由于形態(tài)不穩(wěn)定或特征不明顯而無(wú)法實(shí)施(陳萍等，2002;Frick，1952;Sasakawa，1961)，通常的做法是將其進(jìn)行室內(nèi)飼養(yǎng)，待羽化后再做鑒定，不僅耗時(shí)長(zhǎng)，而且常因標(biāo)本數(shù)量少而使鑒定工作無(wú)法完成(邱一中等，2000)。因此，研究開(kāi)發(fā)快速、簡(jiǎn)便、準(zhǔn)確且適用于昆蟲(chóng)不同發(fā)育時(shí)期的害蟲(chóng)鑒定新技術(shù)，是有效阻止危險(xiǎn)性害蟲(chóng)進(jìn)一步傳播擴(kuò)散，增強(qiáng)檢驗(yàn)檢疫和檢測(cè)監(jiān)測(cè)能力的必要前提。

目前，以分子生物學(xué)技術(shù)進(jìn)行檢疫性害蟲(chóng)種類(lèi)鑒定識(shí)別以及檢測(cè)監(jiān)測(cè)的研究時(shí)有報(bào)道并逐漸得以推廣應(yīng)用，其中RAPD(random amplification of polymorphic DNA)(邱一中等，2000;??l et al.，2006)、RFLP(restriction fragment length polymorphism)(陳萍等，2002;Scheffer et al.，2001)、多重 PCR(mutiplex PCR)(Miura et al.，2004)以及DNA條形編碼(DNA barcoding)(Scheffer et al.，2006)等技術(shù)均已用于斑潛蠅屬害蟲(chóng)的鑒定研究，并取得良好成效。但是，RAPD標(biāo)記技術(shù)對(duì)環(huán)境因子的變化較為敏感，擴(kuò)增產(chǎn)物的可重復(fù)性相對(duì)較差，結(jié)果的可靠性低，因此目前已較少使用;RFLP技術(shù)由于所需DNA量大，步驟多，周期長(zhǎng)，而且制備探針及檢測(cè)中要用到放射性同位素，雖可利用非放射性同位素標(biāo)記方法代替，但成本高，成功率低，因此影響了其推廣使用。mtDNA COⅠ基因?yàn)槎嗫截悾赶颠z傳基因，既具有相對(duì)的保守性又有足夠的變異性，進(jìn)化速率比核DNA快，很少存在插入和缺失，檢測(cè)靈敏度高，適用于種間鑒別(Hoelzel，1991;Hoy，1994)?；?COⅠ基因的DNA條形編碼技術(shù)盡管需要經(jīng)過(guò)產(chǎn)物回收、純化、測(cè)序以及序列比對(duì)等步驟，但目前已成為物種鑒定的新方法。SS-PCR技術(shù)由COⅠ標(biāo)記技術(shù)轉(zhuǎn)化而來(lái)，具有重現(xiàn)性強(qiáng)、譜帶單一等優(yōu)點(diǎn)，適宜于大規(guī)?？焖贆z測(cè)分析。

本研究根據(jù)已知的美洲斑潛蠅mtDNA COⅠ基因序列，設(shè)計(jì)美洲斑潛蠅特異性引物1對(duì)，并通過(guò)SS-PCR技術(shù)證實(shí)該引物只對(duì)美洲斑潛蠅mtDNA具有擴(kuò)增能力，對(duì)同域發(fā)生的其他種類(lèi)的潛葉蠅類(lèi)害蟲(chóng)，如南美斑潛蠅、三葉斑潛蠅、蔥斑潛蠅、豌豆?jié)撊~蠅等沒(méi)有擴(kuò)增能力;該技術(shù)體系不僅對(duì)美洲斑潛蠅的成蟲(chóng)具有良好的擴(kuò)增效果，而且對(duì)單頭1～3齡幼蟲(chóng)、蛹以及單粒卵也具有同樣的擴(kuò)增效果;同時(shí)，該技術(shù)體系還具有較強(qiáng)的靈敏性，其最低檢出閾值為1/3840頭成蟲(chóng);對(duì)靶標(biāo)序列的測(cè)序及Gen-Bank比對(duì)結(jié)果顯示，該片段與美洲斑潛蠅的同源性為 100%(Shang et al.，2008;Yang et al.，2011)，而與其他雙翅目昆蟲(chóng)的同源性均低于96%(http:∥blast.ncbi.nlm.nih.gov/Blast.cgi)。此外，SS-PCR 分析技術(shù)操作簡(jiǎn)便快捷，縮短了通關(guān)時(shí)間，提高了工作效率。因此，該技術(shù)體系可以用于美洲斑潛蠅的檢測(cè)鑒定工作，并在害蟲(chóng)檢疫檢驗(yàn)與檢測(cè)監(jiān)測(cè)中具有重要應(yīng)用價(jià)值。然而，由于我國(guó)菜田常見(jiàn)的潛葉蠅類(lèi)害蟲(chóng)種類(lèi)較多，本研究中SS-PCR引物的特異性仍有待進(jìn)一步驗(yàn)證。

致謝:北京城市學(xué)院生物技術(shù)專(zhuān)業(yè)的袁波同學(xué)和長(zhǎng)江大學(xué)農(nóng)學(xué)院植物保護(hù)專(zhuān)業(yè)的陳婷婷同學(xué)協(xié)助完成部分試驗(yàn)，特表感謝!

陳兵，趙云鮮，康樂(lè).2002.外來(lái)斑潛蠅入侵和適應(yīng)機(jī)理及管理對(duì)策.動(dòng)物學(xué)研究，23(2):155-160.

陳萍，溫碩洋，曾玲.2002.PCR-RFLP技術(shù)應(yīng)用于鑒定美洲斑潛蠅和番茄斑潛蠅的初步研究.武夷科學(xué)，18(1):60-64.

陳文龍，張潤(rùn)杰.2003.美洲斑潛蠅生物學(xué)特性研究.中山大學(xué)學(xué)報(bào):自然科學(xué)版，42(6):91-92.

范文中，王中武.2005.美洲斑潛蠅的發(fā)生及防治.吉林農(nóng)業(yè)科技學(xué)院學(xué)報(bào)，14(4):11-13.

胡長(zhǎng)效，蘇新林.2003.我國(guó)美洲斑潛蠅發(fā)生及防治研究進(jìn)展.江西農(nóng)業(yè)學(xué)報(bào)，15(1):48-54.

雷仲仁，王音.2005.美洲斑潛蠅∥萬(wàn)方浩，鄭小波，郭建英.重要農(nóng)林外來(lái)入侵物種的生物學(xué)與控制.北京:科學(xué)出版社，177-205.

邱一中，吳文哲，蕭旭峰，石正人.2000.RAPD-PCR在六種斑潛蠅(Liriomyza spp.)(雙翅目:潛蠅科)快速鑒定技術(shù)之應(yīng)用.中華昆蟲(chóng)，20(4):293-309.

王軍，石寶才，宮亞軍，廖度，宋婧，路虹.1999.美洲斑潛蠅寄主植物調(diào)查名錄.北京農(nóng)業(yè)科學(xué)，17(1):37-39.

王莉萍，杜予州，何婭婷，陸亞娟，陸自強(qiáng).2007.不同地理種群美洲斑潛蠅及近緣種的rDNA-ITS1序列分析和比較.昆蟲(chóng)學(xué)報(bào)，50(6):597-603.

問(wèn)錦曾，王音，雷仲仁.1996.美洲斑潛蠅中國(guó)新紀(jì)錄.昆蟲(chóng)分類(lèi)學(xué)報(bào)，18(4):311-312.

楊龍龍.1995.對(duì)斑潛蠅屬中檢疫性害蟲(chóng)的研究.植物檢疫，9(1):1-5.

周志湘，萬(wàn)方浩，張桂芬，陳斌.2007.煙粉虱基因組DNA快速提取方法.植物保護(hù)，33(5):131-133.

??l B，Tongu? A，?zgül O，Civelek H S and Kaya B.2006.The use of RAPD-PCR analysis in characterization of Liriomyza trifolii(Burgess，1880)，Liriomyza congesta(Becker 1903)，Agromyza apfelbecki Strobl，1902 and Chromatomyia horticola(Goureau，1851)species collected from Turkey.Turkiye Entomoloji Dergisi，30(4):243-253.

Deeming J C.1992.Liriomyza sativae Blanchard(Diptera:Agromyzidae)established in the old world.Tropical Pest Management，38(2):218-219.

EPPO/CABI.2011.Liriomyza sativae.Data Sheets on Quarantine Pests.http:∥www.eppo.org/QUARANTINE/listA2.htm(2011-10-10).

Frick K E.1952.A generic revision of the family Agromyzidae(Diptera)with a catalogue of New World species.University of California Publications in Entomology，8:339-452.

Hoelzel A R.1991.Analysis of regional mitochondrial DNA variation in the killer whale:implications for conservation∥Hoelzel A R.Genetic Ecology of Whales and Dolphins.Report to the International Whaling Commission.UK:Cambridge，225-233.

Hoy M A.1994.Insect Molecular Genetics:An Introduction to Principles and Applications.San Diego，California:Academic Press.

Miura K，Tagami Y，Ohtaishi M and Iwasaki A.2004.Application of molecular techniques to distinguish Liriomyza trifolii from L.sativae(Diptera:Agromyzidae)on tomato cultivation in Japan.Journal of Economic Entomology，97(3):964-969.

Parrella M P.1987.Biology of Liriomyza.Annual Review of Entomology，32:201-224.

Sasakawa M.1961.A study of the Japanese Agromyzidae(Diptera)，Part 2.Pacific Insects，3:307-472.

Scheffer S J，Lewis M L and Joshi R C.2006.DNA barcoding applied to invasive leafminers(Diptera:Agromyzidae)in the Philippines.Annals of the Entomological Society of America，99(2):204-210.

Scheffer S J，Wijesekara A，Visser D and Hallett R H.2001.Polymerase chain reaction-restriction fragment-length polymorphism method to distinguish Liriomyza huidobrensis from L.langei(Diptera:Agromyzidae)applied to three recent leafminer invasions.Journal of Economic Entomology，94(5):1177-1182.

Shang H，Zhao Y，Cui X，Li H，Zhou J and Huang T.2008.Polymerase chain reaction-restriction fragment length polymorphism assays to distinguish Liriomyza trifolii(Diptera:Agromyzidae)from associated species in China.http:∥www.ncbi.nlm.nih.gov/protein/163867295(2011-11-02).

Shiao S F，Lin F J and Wu W J.1991.Redescription of four Liriomyza species(Diptera:Agromyzidae)from Taiwan.Chinese Journal of Entomology，11:65-74.

Spencer K A.1986.Agromyzidae(Diptera)from Thailand:new species，revisionary notes and new records.Proceedings of the Indian Academy of Sciences，95:487-507.

Yang F，Du Y Z，Wang L P，Cao J M and Yu W W.2011.The complete mitochondrial genome of the leafminer Liriomyza sativae(Diptera:Agromyzidae):great difference in the A+T-rich region compared to Liriomyza trifolii.Gene，485(1):7-15.

Zitter T A and Tsai J H.1980.Flies∥Harris K F and Maramorosch K.Vectors of Plant Pathogens.New York，USA:Academic Press，165-176.