黃 振，侯有明 ，郭瓊霞 ，陳韶萍,

（1. 福建農(nóng)林大學(xué)植物保護(hù)學(xué)院，福建 福州 350002；2. 福州長(zhǎng)樂機(jī)場(chǎng)海關(guān)，福建 福州 350209；3. 福州海關(guān)技術(shù)中心，福建 福州 350001）

0 引言

【研究意義】南瓜實(shí)蠅Bactrocera tau（Walker），屬雙翅目Diptera、實(shí)蠅科Tephritidae、果實(shí)蠅屬Bactrocera[1]，是我國(guó)進(jìn)境植物檢疫性害蟲。卵、幼蟲、蛹能隨寄主如果蔬、土壤、包裝物等可進(jìn)行遠(yuǎn)距離傳播。南瓜實(shí)蠅食性雜、寄主廣，達(dá)16個(gè)科80余種[2]，主要危害南瓜、角瓜、絲瓜、黃瓜、葫蘆、苦瓜、桑椹、西瓜、菠蘿蜜、楊桃、芒果、番石榴、番茄、扁蒲、菜豆等瓜果植物，對(duì)葫蘆科植物危害尤為嚴(yán)重[3]。Walker1849年首次在中國(guó)報(bào)道[4]，目前主要報(bào)道分布在越南、緬甸、泰國(guó)、老撾、日本、柬埔寨、馬來(lái)西亞、菲律賓、新加坡、不丹、斯里蘭卡、印度、印度尼西亞等地[5?6]。國(guó)際果蔬貿(mào)易給南瓜實(shí)蠅的傳入傳播帶來(lái)了極大的潛在風(fēng)險(xiǎn)[7]，日常分類鑒定工作通常以成蟲外部形態(tài)特征作為分類鑒定依據(jù)，而日常檢疫中經(jīng)常截獲幼蟲、卵或蛹，甚至翅、足等殘?bào)w，無(wú)法快速、準(zhǔn)確鑒定，通常是在實(shí)驗(yàn)室內(nèi)飼養(yǎng)為成蟲后進(jìn)行鑒定，耗時(shí)較長(zhǎng)，嚴(yán)重影響果蔬進(jìn)出口貿(mào)易的通關(guān)速度。因此，建立快速、特異、有效的PCR技術(shù)，對(duì)快速鑒定南瓜實(shí)蠅具有重要的現(xiàn)實(shí)意義?！厩叭搜芯窟M(jìn)展】采用分子生物學(xué)的方法開展昆蟲種類鑒定，可以解決通過(guò)傳統(tǒng)形態(tài)特征不能進(jìn)行實(shí)蠅種類鑒定的問(wèn)題。越來(lái)越多的分子鑒定技術(shù)被運(yùn)用于昆蟲的鑒定[8]。Mun等[9]應(yīng)用限制性內(nèi)切酶Apa I、Nhe I和Sac I，將南瓜實(shí)蠅、桔小實(shí)蠅和地中海實(shí)蠅等4種實(shí)蠅區(qū)分開來(lái)。吳佳教等[10]用PCR-RFLP技術(shù)對(duì)我國(guó)口岸截獲頻率較高的9種檢疫性實(shí)蠅開展了研究，用限制性內(nèi)切酶MSE1和DRAI對(duì)PCR擴(kuò)增產(chǎn)物進(jìn)行酶切，用得到的酶切位點(diǎn)將供試的9種實(shí)蠅區(qū)分開來(lái)。張亮等[11]通過(guò) RAPD 技術(shù)構(gòu)建了南瓜實(shí)蠅、黑膝實(shí)蠅B. scutellaris等6種實(shí)蠅的指紋圖譜的分種鑒定。林麗莉等[12]應(yīng)用RFLP技術(shù)，對(duì)我國(guó)部分地區(qū)發(fā)生分布的蜜柑大實(shí)蠅B. tsuneonis（Miyake）和橘大實(shí)蠅B. minax（Enderlein）開展鑒別研究。Chua等[13]采用RFLP技術(shù)成功地區(qū)分了洋桃實(shí)蠅和木瓜實(shí)蠅，進(jìn)而使桔小實(shí)蠅復(fù)合種的有效識(shí)別得以實(shí)現(xiàn)。Kakouli等[14]用擴(kuò)增片段長(zhǎng)度多態(tài)性技術(shù) AFLP進(jìn)行了實(shí)蠅種類的檢測(cè)鑒定。然而，上述對(duì)實(shí)蠅的檢測(cè)鑒定技術(shù)各有優(yōu)缺點(diǎn)，如 RFLP和RAPD技術(shù)具有需求DNA的量較少、基因的多態(tài)性豐富、對(duì)基因組的覆蓋范圍比較廣等優(yōu)點(diǎn)；但RFLP技術(shù)由于操作繁瑣，需通過(guò)檢測(cè)內(nèi)切酶識(shí)別位點(diǎn)上的變異，來(lái)進(jìn)行昆蟲種間親緣關(guān)系和昆蟲種屬特異性鑒定，多態(tài)性檢出率低，而RAPD 技術(shù)的穩(wěn)定性較差，對(duì)試驗(yàn)環(huán)境因子變化較為敏感，擴(kuò)增的產(chǎn)物可重復(fù)性、穩(wěn)定性較低等，限制了RFLP和RAPD技術(shù)的廣泛應(yīng)用[15]。雖然AFLP技術(shù)可以較好地彌補(bǔ)上述兩種技術(shù)方法的不足，但AFLP技術(shù)對(duì)檢測(cè)樣本DNA的質(zhì)量要求較高，且所需設(shè)備的費(fèi)用較為昂貴。因此，上述的檢測(cè)方法難以滿足口岸進(jìn)出境果蔬貿(mào)易快速通關(guān)、準(zhǔn)確鑒定及低成本的檢疫需求。因?yàn)閙tDNA COⅠ基因相對(duì)保守，排列的順序比較穩(wěn)定緊密，容易被通用引物擴(kuò)增，又有足夠的特異性能將不同物種區(qū)分開，所以mtDNA COⅠ基因作為標(biāo)記基因，為生物分類學(xué)提供了信息化的分類標(biāo)準(zhǔn)和有效的分類手段，成為目前昆蟲分類鑒定及昆蟲分子進(jìn)化系統(tǒng)研究應(yīng)用較多的基因之一，并被應(yīng)用于實(shí)蠅近緣種的系統(tǒng)發(fā)育研究[16?18]。【本研究切入點(diǎn)】基于RAPD、AFLP、RFLP等技術(shù)在實(shí)際應(yīng)用中表現(xiàn)的不足[19]，本試驗(yàn)在前人利用COI標(biāo)記鑒定技術(shù)的基礎(chǔ)上，利用mt DNA COI基因序列，篩選并設(shè)計(jì)種特異引物[20]，采用SS-PCR方法快速鑒定南瓜實(shí)蠅?！緮M解決的關(guān)鍵問(wèn)題】應(yīng)用種特異性引物快速鑒定南瓜實(shí)蠅的方法，實(shí)現(xiàn)對(duì)口岸進(jìn)境水果攜帶并截獲的疑似南瓜實(shí)蠅的卵、幼蟲、蛹、成蟲及成蟲殘?bào)w的快速鑒定。

1 材料與方法

1.1 供試?yán)ハx

供試蟲樣：南瓜實(shí)蠅B. tau（Walker）、洋桃實(shí)蠅B. carambolae Drew & Hancock、瓜實(shí)蠅B. cucurbitae（Coquillett）、辣 椒 實(shí) 蠅B. latifrons（Hendel）、腿端黑實(shí)蠅B. atrifemur Drew & Hancock、番石榴果實(shí)蠅 B. correcta（Bezzi）、二 顏 帶 實(shí) 蠅 B. cilifera（Hendel）、桔小實(shí)蠅B. dorsalis（Hendel）、銹紅果實(shí)蠅B. rubigina（Wang & Zhao）、瘤脛實(shí)蠅B. tuberculata（Bezzi）、具條實(shí)蠅B. scutellata（Hendel）、何氏華實(shí)蠅B. hochii（Zia）、近黑顏實(shí)蠅B. parater（Zhao &Lin）、黑顏實(shí)蠅B. diaphora Coquillett、黑膝實(shí)蠅B.scutellaris（Bezzi）、滇寡鬃實(shí)蠅B. modica（Hardy）、瑞麗果實(shí)蠅B. ruiliensis Wang, Long et Zhang, sp. nov.、五指山實(shí)蠅B. wuzhishana Lin et Yang, sp. nov.、棗實(shí)蠅Carpomya vesuviana Costa、瓜棍腹實(shí)蠅Dacus longicornis Wiedemann等20種實(shí)蠅，其中除Bactrocera屬18種外，另有 Carpomya、Dacus2個(gè)屬的種類各1種。蟲樣來(lái)自口岸進(jìn)境果蔬檢疫截獲、中國(guó)邊境監(jiān)測(cè)、調(diào)查和誘捕的檢疫性實(shí)蠅種類。用于分子生物學(xué)鑒定試驗(yàn)的實(shí)蠅標(biāo)本，放置于4 ℃冰箱保存。

1.2 DNA提取和質(zhì)量檢測(cè)

將20種實(shí)蠅蟲樣放置于2 mL的離心管中，加入適量液氮充分研磨，根據(jù)用OMEGA E.Z.N.A.TMInsect DNA Kit試劑盒的操作方法，對(duì)蟲樣進(jìn)行基因組DNA的提取，并放置于?20℃冰箱保存?zhèn)溆谩?/p>

DNA質(zhì)量檢測(cè)采用上游引物L(fēng)CO1490（堿基序列為：5′-GGTCAACAAATCATAAAGATATTG-3′）和下游引物HCO2198（堿基序列為：5′-TAAACTTC AGGGTGACCAAAAAATCA-3′）的1對(duì)通用引物對(duì)20種實(shí)蠅蟲樣提取的基因組DNA進(jìn)行質(zhì)量檢測(cè)，在定量梯度PCR儀上進(jìn)行PCR反應(yīng)，電泳，凝膠成像，拍攝記錄質(zhì)量檢測(cè)結(jié)果。

1.3 引物設(shè)計(jì)

通過(guò)數(shù)據(jù)庫(kù)（NCBI）查找已公布的南瓜實(shí)蠅的登錄序列，分析比較并選定登錄號(hào)為JN542420的線粒體mt DNA COⅠ基因序列，利用數(shù)據(jù)庫(kù)（NCBI）中提供的BLAST程序檢查同源序列，最終選定8種實(shí)蠅（表1）。

表 1 設(shè)計(jì)南瓜實(shí)蠅特異性引物所用實(shí)蠅種類及登錄號(hào)Table 1 Species and accession numbers of fruit flies used for designing B. tau-specific primers

應(yīng)用數(shù)據(jù)庫(kù)下載FASTA格式，對(duì)表1中的8種實(shí)蠅序列進(jìn)行序列比對(duì)，再根據(jù)SNP位點(diǎn)，利用Primer-Premier 5.0軟件進(jìn)行人工設(shè)計(jì)引物，采用Oligo 6.44對(duì)設(shè)置的引物進(jìn)行評(píng)定，最后應(yīng)用數(shù)據(jù)庫(kù)（NCBI）中 的Primer-BLAST（BLAST：Basic Local Alignment Search Tool）DNA數(shù)據(jù)庫(kù)進(jìn)行相似性比較分 析，檢查同源序列。

1.4 引物種特異性測(cè)試

測(cè)試南瓜實(shí)蠅引物的種特異性時(shí)，選取南瓜實(shí)蠅為陽(yáng)性對(duì)照，其余的20種實(shí)蠅蟲樣為陰性對(duì)照，檢測(cè)驗(yàn)證設(shè)計(jì)引物種特異性測(cè)試的反應(yīng)體系和條件，電泳檢測(cè)，凝膠成像，查看有否擴(kuò)增出預(yù)期大小 的目標(biāo)片段，拍攝并記錄結(jié)果。

1.5 引物靈敏度測(cè)試

將提取南瓜實(shí)蠅的基因組DNA模板，按照10?1、10?2、10?3倍3種不同濃度稀釋，測(cè)試提取的南瓜實(shí)蠅DNA的稀釋濃度后的靈敏度。用設(shè)計(jì)的種特異性引物NF404和NR610進(jìn)行PCR擴(kuò)增，對(duì)DNA模板進(jìn)行最低檢出閾值的測(cè)定，驗(yàn)證本檢測(cè)方法的靈 敏度。

1.6 SS-PCR方法的應(yīng)用與驗(yàn)證

將福州長(zhǎng)樂機(jī)場(chǎng)進(jìn)境旅客攜帶果蔬中截獲的南瓜實(shí)蠅并飼養(yǎng)到成蟲，經(jīng)過(guò)專家鑒定復(fù)核，收集幼蟲、蛹、成蟲、足、翅等17個(gè)蟲樣，按照本研究的SS-PCR方法，將提取的基因組DNA為模板，應(yīng)用種特異性引物NF404和NR610進(jìn)行PCR擴(kuò)增，來(lái)驗(yàn)證 該方法的穩(wěn)定性與準(zhǔn)確性。

2 結(jié)果與分析

2.1 DNA 質(zhì)量檢測(cè)結(jié)果

采用1對(duì)LCO1490和HCO2198通用引物，對(duì)提取的20種實(shí)蠅蟲樣的基因組DNA模板進(jìn)行PCR擴(kuò)增，結(jié)果顯示所有的實(shí)蠅蟲樣的DNA均能擴(kuò)增出一條單一整齊且清晰透明的，長(zhǎng)度約700 bp的目標(biāo)條帶（圖1）。

圖 1 利用LCO1490和HC02198 一對(duì)通用型引物對(duì)提取的實(shí)蠅DNA模板進(jìn)行質(zhì)量檢測(cè)的結(jié)果Fig. 1 Utilization of universal primers, LCO1490 and HC02198,for testing extracted DNA templates from fruit flies

2.2 特異性引物選擇

應(yīng)用Bioedit對(duì)已下載的JN542420.1、HM561566.1、JN542417.1、AY530891.1、AF423105.1、JX855917.1、KF998674.1、KF318598.1等8種序列進(jìn)行ClustalW多重比對(duì)，通過(guò)篩選最終選擇的種特異性引物NF′404和NR′610（圖2）。

圖 2 8種果實(shí)蠅的ClustalW多重比對(duì)的部分結(jié)果（種特異性引物NF’404和NR’610序列位置用橫線劃出）Fig. 2 Partial result of ClustalW multiple comparisons on 8 fruit flies

NF′404：5′-TGCCTCGACGATATTCTGACT-3′；

NR′610：5′-AACTGTGTTCAGCAGGTGGT-3′。

將NF′404和NR′610利用NCBI數(shù)據(jù)庫(kù)提供的Primer-BLAST程序檢查同源序列，最終將正向引物NF′404第6個(gè)堿基C轉(zhuǎn)換成T，第7個(gè)堿基G轉(zhuǎn)換成A，第8、19個(gè)堿基A轉(zhuǎn)換成G，確定出引物NF404和NR610。

篩選出NF404和NR610種特異性引物，由上海英濰捷基貿(mào)易有限公司合成。特異性引物NF404和NR610的序列分別為：

NF404：5′-TGCCTTAGCGATATTCTGGCT-3′；

NR610：5′-AACTGTGTTCAGCAGGTGGT-3′。

2.3 引物種特異性

選用篩選獲得的種特異性引物NF404和NR610進(jìn)行擴(kuò)增，結(jié)果表明：僅南瓜實(shí)蠅擴(kuò)增出一條清晰透明且單一的約207 bp的條帶，其他20種實(shí)蠅種類樣本均未顯現(xiàn)目標(biāo)條帶（圖3）。表明該試驗(yàn)設(shè)計(jì)的NF404和NR610種特異性引物具有較強(qiáng)的特異性和穩(wěn)定性。

上海英濰捷基貿(mào)易有限公司對(duì)本試驗(yàn)將所得到的PCR產(chǎn)物進(jìn)行測(cè)序，將序列提交NCBI數(shù)據(jù)庫(kù)進(jìn)行相似性分析（BLAST），表明該段序列與數(shù)據(jù)庫(kù)中 的南瓜實(shí)蠅序列完全一致。

2.4 引物的靈敏度

應(yīng)用核酸蛋白分析儀對(duì)提取的南瓜實(shí)蠅DNA模板的質(zhì)量濃度進(jìn)行測(cè)試，質(zhì)量濃度為56.95 ng·μL?1。取3種不同濃度的DNA模版測(cè)定最低檢出閾值，結(jié)果表明：隨著模版質(zhì)量濃度的降低，擴(kuò)增出的條帶逐漸變淡，在質(zhì)量濃度稀釋100倍后，仍可見較為明 顯條帶（圖4），說(shuō)明本試驗(yàn)方法靈敏度較高。

2.5 應(yīng)用與驗(yàn)證

采用SS-PCR快速鑒定技術(shù)對(duì)截獲的南瓜實(shí)蠅標(biāo)

圖 3 引物NF404和NR610的種特異性驗(yàn)證Fig. 3 Verification on species-specificity of primers, NF404 and NR610

圖 4 引物NF404和NR610的靈敏度驗(yàn)證Fig. 4 Verification on sensitivity of primers, NF404 and NR610

本進(jìn)行檢測(cè)驗(yàn)證，取經(jīng)復(fù)核鑒定的幼蟲、蛹、成蟲及成蟲的足、翅或部分殘?bào)w等17個(gè)樣品經(jīng)過(guò)檢驗(yàn)，均檢測(cè)有目標(biāo)條帶（圖5），說(shuō)明分子生物學(xué)鑒定結(jié)果與形態(tài)學(xué)鑒定結(jié)果一致，即基于mt DNA COⅠ基因建立的種特異性COⅠ引物具有很強(qiáng)的特異性，SS-PCR快速鑒定南瓜實(shí)蠅的方法可以應(yīng)用于實(shí)際鑒定。

圖 5 SS-PCR快速鑒定南瓜實(shí)蠅的應(yīng)用與驗(yàn)證Fig. 5 Application and verification of SS-PCR for rapid identification of B. tau

3 討論

3.1 引物的特異性是快速準(zhǔn)確鑒定的決定性因素

被選作分子標(biāo)記的基因很多，但是COI（cytochrome oxidase I或coxI）基因被選為DNA條形碼的靶標(biāo)基因，相對(duì)保守，種間變異大，容易被通用引物擴(kuò)增，同時(shí)又有足夠的變異能夠?qū)⒔壏N、近似種等分開。本試驗(yàn)選用mt DNA COI作為標(biāo)記基因，在基因庫(kù)中篩選查找南瓜實(shí)蠅已公布序列進(jìn)行比對(duì)分析，設(shè)計(jì)引物NF404和NR610，用南瓜實(shí)蠅作為陽(yáng)性對(duì)照，用其形態(tài)相近的3個(gè)果實(shí)蠅屬4個(gè)亞屬19種實(shí)蠅作為陰性對(duì)照，目標(biāo)種能特異性并穩(wěn)定地?cái)U(kuò)增出長(zhǎng)度207 bp清晰且單一的目標(biāo)條帶，其他實(shí)蠅種類均無(wú)條帶出現(xiàn)，表明設(shè)計(jì)引物的種特異性強(qiáng)，能夠鑒定南瓜實(shí)蠅。

3.2 SS-PCR檢測(cè)鑒定方法的實(shí)用性

采用SS-PCR檢測(cè)鑒定方法，不但可快速準(zhǔn)確鑒定特異種類南瓜實(shí)蠅，而且可檢測(cè)到的DNA模板質(zhì)量濃度的最低檢出閾值為5.695×10?3ng·μL?1。在口岸進(jìn)境的多批果蔬的檢疫檢測(cè)時(shí)，對(duì)發(fā)現(xiàn)的不同蟲態(tài)目標(biāo)實(shí)蠅（飼養(yǎng)到成蟲及經(jīng)過(guò)專家形態(tài)鑒定復(fù)核），采用SS-PCR檢測(cè)鑒定方法，可在8 h之內(nèi)快速準(zhǔn)確鑒定特異種類南瓜實(shí)蠅。同時(shí)，本試驗(yàn)建立種特異性引物PCR方法檢測(cè)南瓜實(shí)蠅的不同蟲態(tài)用時(shí)短、DNA用量少、靈敏度高、儀器設(shè)備簡(jiǎn)單、鑒定準(zhǔn)確率高，能夠滿足口岸檢驗(yàn)檢疫高通量快速通關(guān)和快速準(zhǔn)確的要求，而且本研究建立的方法對(duì)南瓜實(shí)蠅的卵或蟲體部分殘?bào)w均可檢測(cè)，不受實(shí)驗(yàn)材料存活狀態(tài)的影響，檢測(cè)靈敏度高。因此，種特異引物 PCR方法在昆蟲的快速鑒定和近緣種區(qū)分上有著 廣泛的應(yīng)用前景。

3.3 對(duì)線粒體基因組的研究有待于進(jìn)一步加強(qiáng)

目前，線粒體基因組在實(shí)蠅科害蟲中的應(yīng)用較少，僅有3屬 14種，在已公開發(fā)表的大多文獻(xiàn)中只描述了新測(cè)定物種的線粒體基因組序列及其特點(diǎn)。對(duì)于利用線粒體基因組在實(shí)蠅科害蟲中進(jìn)行物種鑒定、種群溯源、分子進(jìn)化等方面的研究尚未見報(bào)道。隨著測(cè)序技術(shù)的發(fā)展，獲得具有重要經(jīng)濟(jì)意義的實(shí)蠅種類的線粒體基因組序列研究顯得極為迫切，所以對(duì)實(shí)蠅科昆蟲線粒體基因組的研究和應(yīng)用需要進(jìn)一步加強(qiáng)。