利用DNA條形碼對(duì)10種蒼耳屬雜草的鑒定

袁俊杰+馬新華+龍+陽(yáng)+魏霜+張娜+陳文+楊卓瑜

摘要：使用ITS、ITS2、psbA-trnH、rbcL和matK等序列對(duì)截獲的10種蒼耳屬（Xanthium）雜草進(jìn)行擴(kuò)增測(cè)序，比較各序列的擴(kuò)增和測(cè)序效率，采用最近距離法和相似性搜索算法評(píng)價(jià)不同序列的鑒定效率。采用最佳鑒定序列計(jì)算種間K2P距離并構(gòu)建系統(tǒng)進(jìn)化樹。結(jié)果表明，ITS、psbA-trnH 2個(gè)序列的擴(kuò)增效率和鑒定效率最高；ITS序列的種間變異最大，其次是psbA-trnH序列，matK序列的種間變異最小；ITS序列的Neighbor-Joining（NJ）樹可明顯地將不同蒼耳屬雜草種分開；psbA-trnH序列可明顯地將國(guó)內(nèi)蒼耳種和檢疫性蒼耳種分開。說(shuō)明利用DNA條形碼能夠準(zhǔn)確地鑒別蒼耳屬雜草，ITS和psbA-trnH序列是鑒別蒼耳屬雜草較理想的條形碼組合。該研究結(jié)果為蒼耳屬雜草的分子鑒別提供了科學(xué)依據(jù)與新的思路。

關(guān)鍵詞：psbA-trnH；ITS；DNA條形碼；蒼耳屬；雜草鑒定

中圖分類號(hào)：S451 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1003-935X（2016）03-0011-06

Identification of Ten Xanthium weeds Using DNA Barcoding Sequences

YUAN Junjie1，MA Xinhua1，LONG Yang1，WEI Shuang2，ZHANG Na1，CHEN Wen1，YANG Zhuoyu1

（1.Zhanjiang Entry-exit Inspection and Quarantine Bureau，Zhanjiang 524022，China；

2.Shantou Entry-exit Inspection and Quarantine Bureau，Shantou 515000，China）

Abstract：DNA from ten weed species of the genus Xanthium was amplified and sequenced using five candidate sequences （ITS＼ITS2＼psbA-trnH＼rbcL and matK）. The PCR amplification and sequencing efficiencies were compared and the identification efficiency was assessed using BLAST1 and Distance methods. The best identification of sequence was analyzed and the NJ tree was constructed. The amplification and sequencing efficiencies of ITS，psbA-trnH were the best. ITS showed the largest interspecific variation，followed by psbA-trnH sequence； matK was the worst. The ten species of Xanthium weeds could be distinguished by ITS2 region using inter-specific K2P genetic distance and NJ tree. Quarantine and common Xanthium plants could be distinguished by psbA-trnH region. Using DNA barcode can accurately identify ten species of Xanthium weeds. ITS2 and psbA-trnH could be used as a potential DNA barcode combination to correctly identify the Xanthium weeds. The research provides scientific basis and new ideas for the molecular identification of Xanthium weeds.

Key words：psbA-trnH；ITS；DNA barcode；Xanthium；weed identification

蒼耳屬（Xanthium）為雙子葉綱菊科植物，全球約有25種[1]，蒼耳屬（非中國(guó)種）雜草為我國(guó)禁止進(jìn)境的植物檢疫性有害生物[2]，其適應(yīng)性強(qiáng)、競(jìng)爭(zhēng)力強(qiáng)、繁殖力強(qiáng)，與農(nóng)作物爭(zhēng)奪生存空間，幼苗具有毒性，且總苞具鉤狀的硬刺，易劃破動(dòng)物的皮膚，威脅牲畜健康。蒼耳屬雜草危害較大，一直是世界各國(guó)重點(diǎn)關(guān)注的惡性雜草，在口岸防止蒼耳屬雜草的入侵具有重要意義。然而在進(jìn)出口貨物中，蒼耳屬雜草一般以籽實(shí)的形態(tài)混雜于糧谷類產(chǎn)品及其包裝物中，種類繁多，近似種相似性較高，鑒定困難；且蒼耳屬雜草籽在運(yùn)輸過(guò)程中總苞表面堅(jiān)硬的刺以及密布的柔毛等重要鑒定特征會(huì)發(fā)生較大磨損，這就給形態(tài)學(xué)鑒定帶來(lái)很大的難度[3-4]。因此，建立穩(wěn)定且準(zhǔn)確的分子生物學(xué)分類方法作為形態(tài)學(xué)鑒定的補(bǔ)充方法具有重要意義。

DNA條形碼（DNA barcoding）由加拿大動(dòng)物學(xué)家Hebert 等首次提出[5]，具有不受個(gè)體形態(tài)特征限制、不受個(gè)體發(fā)育階段影響、從基因水平上客觀區(qū)分物種、操作方法相對(duì)簡(jiǎn)便等特點(diǎn)[6-10]，是近十幾年來(lái)研究進(jìn)展最迅速的學(xué)科之一[11]，其通用序列的篩選一直是該領(lǐng)域研究的熱點(diǎn)問(wèn)題。2009年，生命條形碼聯(lián)盟植物工作組（the plant working group of the consortium for the barcode of life，CBOL）分析了550個(gè)物種907個(gè)樣品后，推薦以復(fù)合序列rbcL+matK作為整個(gè)植物界的DNA條形碼序列，但其在物種水平的鑒定成功率僅為72%[12]，所以仍在進(jìn)行驗(yàn)證和新序列的研究工作。

目前DNA條形碼在植物類群中的研究和應(yīng)用還處于起步階段。理想的條形碼序列在單一科屬內(nèi)大量近緣物種鑒別時(shí)應(yīng)具備準(zhǔn)確性和高效性，同時(shí)在大量科屬間進(jìn)行物種鑒別時(shí)應(yīng)該具有通用性。盡管已有學(xué)者在近緣屬、種中開展植物條形碼的研究工作[13-16]，但在被考察序列受關(guān)注程度、樣本的數(shù)目、取樣密度及代表性、候選序列的評(píng)價(jià)指標(biāo)等方面存在一定的局限性。學(xué)者們推薦了幾種植物的DNA條形碼候選序列及其組合[17-19]，但這些序列是否適合每一個(gè)具體的植物科或?qū)?，尚需進(jìn)行針對(duì)性的探索研究。在此背景下，本研究選取近年來(lái)幾個(gè)重點(diǎn)推薦的序列，如ITS、ITS2、rbcL、matK以及psbA-trnH，比較其對(duì)蒼耳屬雜草種類的鑒別能力，考察各候選序列作為該屬雜草DNA條形碼通用序列的適用性，以期為DNA條形碼技術(shù)在蒼耳屬雜草籽鑒定中的應(yīng)用提供依據(jù)。

1 材料與方法

1.1 材料

試驗(yàn)所用蒼耳屬雜草籽從近5年的進(jìn)境糧谷（美國(guó)黃大豆、巴西黃大豆、烏克蘭玉米、澳大利亞油菜籽、加拿大油菜籽、美國(guó)高粱等）中截獲（表1）。經(jīng)中國(guó)檢驗(yàn)檢疫科學(xué)研究院鑒定，憑證標(biāo)本保存于湛江出入境檢驗(yàn)檢疫局。

1.2 方法

1.2.1 DNA提取 采用改良后的DNeasy Plant Mini Kit試劑盒法提取DNA，參照試劑盒說(shuō)明書進(jìn)行操作。

1.2.2 PCR擴(kuò)增和電泳檢測(cè) PCR反應(yīng)使用寶生物工程（大連）有限公司rTaq酶進(jìn)行擴(kuò)增，反應(yīng)采用50 μL體系：3 μL DNA模板，5 μL 10×PCR Buffer，4 μL dNTPs，0.25 μL rTaq，各1 μ上、下游引物（表2），補(bǔ)充ddH2O至50 μL。ITS、ITS2、rbcL采用相同的PCR條件：95 ℃、4 min；94 ℃、30 s，52 ℃、45 s，72 ℃、45 s，35個(gè)循環(huán)；72 ℃、10 min。matK的PCR條件：95 ℃、4 min；94 ℃、45 s，52 ℃、1.5 min，72 ℃、1.5 min，35個(gè)循環(huán)；72 ℃、10 min。psbA-trnH的PCR條件：95 ℃、4 min；94 ℃、30 s，52 ℃、1 min，72 ℃、1 min，35個(gè)循環(huán)；72 ℃、10 min。擴(kuò)增產(chǎn)物用1%瓊脂糖凝膠電泳檢測(cè)，用凝膠成像分析系統(tǒng)進(jìn)行分析。

使用ABI3730XL測(cè)序儀（Applied Biosystems Co.，USA）直接對(duì)PCR產(chǎn)物進(jìn)行雙向測(cè)序。利用序列拼接軟件CodonCode Aligner V3.7.1對(duì)測(cè)序峰圖文件進(jìn)行拼接，去除引物區(qū)及低質(zhì)量序列。使用Clustal（MEGA 5.0）軟件對(duì)拼接后的序列進(jìn)行比對(duì)分析。運(yùn)用相似性搜索算法（BLAST1）和最近距離法（K2P nearest distance）考察各序列的鑒定成功率。整合處理后結(jié)合樣品序列，并運(yùn)用MEGA 5.0建立Neighbor Joining（NJ）樹，分析蒼耳屬各樣本的親緣關(guān)系。

2 結(jié)果與討論

2.1 序列信息及鑒定成功率

從表3可以看出，5條片段的 PCR 擴(kuò)增效率和測(cè)序效率均較高，除matK的擴(kuò)增效率為56.0%，測(cè)序效率為88%外，其他片段的擴(kuò)增和測(cè)序效率均高于90%，其中ITS、ITS2、psbA-trnH擴(kuò)增和測(cè)序效率均高于95%。片段長(zhǎng)度方面，ITS2 的長(zhǎng)度最短，為 494～498 bp；其次是 psbA-trnH，為 521～543 bp；rbcL 和 ITS 較長(zhǎng)，分別為 720 bp 和 720～727 bp；matK最長(zhǎng)，為844 bp。各片段的 GC 量以psbA-trnH的最低，為 26.5%～26.9%；ITS2的最高，為 548%～56.4%。psbA-trnH與ITS序列均有53個(gè)變異位點(diǎn)，其次為ITS2，包含26個(gè)變異位點(diǎn)，matK和 rbcL 分別有7、3個(gè)變異位點(diǎn)。psbA-trnH有22 bp的插入/缺失，為5個(gè)序列中最多的。

2.2 鑒定效率分析

采用Distance 和BLAST1方法分析不同DNA條形碼序列對(duì)10種蒼耳屬雜草籽的鑒定效果，結(jié)果如表4所示，應(yīng)用ITS序列2種方法的鑒定成功率均最高，分別為94%、99%；其次是psbA-trnH，鑒定成功率分別為 92%、89%；ITS2稍差于 psbA-trnH，鑒定成功率為 91%、86%；而rbcL和matK的鑒定效率較低，2種方法的成功率均不足40%。

2.3 10種檢疫性蒼耳屬雜草的psbA-trnH和ITS序列分析

2.3.1 K2P遺傳距離 同一地區(qū)樣本個(gè)體之間序列完全相同，種內(nèi)K2P遺傳距離為0。利用 MEGA 7.0 軟件，采用 Kimura 2-parameter model模型計(jì)算種間遺傳距離（表4），蒼耳屬植物ITS序列的遺傳距離最高值為 0.068 0，是蒼耳與沃式蒼耳的遺傳距離；最小值為 0.000 0，是河岸蒼耳與賓州蒼耳，意大利蒼耳與球果蒼耳，西方蒼耳與意大利蒼耳、球果蒼耳的遺傳距離。這說(shuō)明國(guó)內(nèi)蒼耳與各檢疫蒼耳之間存在較大的遺傳距離，存在較大的遺傳差異；河岸蒼耳與賓州蒼耳，意大利蒼耳與球果蒼耳，西方蒼耳與意大利蒼耳、球果蒼耳之間遺傳距離為0，不存在遺傳差異。

同一地區(qū)樣本個(gè)體之間序列完全相同，種內(nèi)K2P遺傳距離為0。利用 MEGA 7.0 軟件，采用Kimura 2-parameter model模型計(jì)算種間遺傳距離（表5），蒼耳屬植物 psbA-trnH 序列的遺傳距離最高值為0.093 5，是刺蒼耳與柱果蒼耳、意大利蒼耳、河岸蒼耳的遺傳距離；最小值為 0.000 0，是北美蒼耳與球果蒼耳、賓州蒼耳、西方蒼耳、沃式蒼耳之間及柱果蒼耳與意大利蒼耳、河岸蒼耳之間的遺傳距離。這說(shuō)明北美蒼耳與球果蒼耳、西方蒼耳、賓州蒼耳、沃式蒼耳之間不存在遺傳差異。

2.3.2 Neighbor-Joining（NJ）樹 利用 MEGA 70 軟件，采用 Neighbor-Joining 方法構(gòu)建系統(tǒng)發(fā)育樹?；?ITS 序列（圖1），蒼耳、刺蒼耳被分為一支，并且存在遺傳距離，各自成一支；賓州蒼耳、河岸蒼耳被分為一支；柱果蒼耳獨(dú)成一支；意大利蒼耳、球果蒼耳、西方蒼耳被分為一大支；北美蒼耳、沃式蒼耳各自獨(dú)成一支；從 NJ 樹圖可以看出，10種蒼耳屬雜草均可以互相區(qū)分，可直觀地展現(xiàn)不同序列的物種鑒別情況。

基于psbA-trnH序列（圖2），蒼耳、刺蒼耳被分為一支；柱果蒼耳、意大利蒼耳、河岸蒼耳被分為一支；北美蒼耳、球果蒼耳、賓州蒼耳、沃式蒼耳、西方蒼耳被分為一支；對(duì)其進(jìn)行 Bootstrap檢驗(yàn)，其分組支持率為100%，說(shuō)明這3支蒼耳屬雜草在系統(tǒng)進(jìn)化發(fā)育中是相互獨(dú)立的。柱果蒼耳、意大利蒼耳、河岸蒼耳又被劃分為不同級(jí)別的3個(gè)分支，但其Bootstrap檢驗(yàn)支持率較低，說(shuō)明這3種蒼耳屬雜草獨(dú)立性相對(duì)較弱。北美蒼耳、球果蒼耳、賓州蒼耳、沃式蒼耳、西方蒼耳也被分為各自獨(dú)立的分支，但其Bootstrap檢驗(yàn)支持率較低，說(shuō)明這5種蒼耳屬植物獨(dú)立性相對(duì)較弱。

由以上數(shù)據(jù)可以得出，ITS序列構(gòu)建的NJ樹

可以對(duì)10種蒼耳屬雜草籽進(jìn)行區(qū)分，相對(duì)于psbA-trnH 序列的區(qū)分度更高，更加適用于蒼耳屬雜草的鑒別。

3 討論

本研究比較分析了ITS、ITS2、psbA-trnH、rbcL和matK等5條DNA條形碼序列對(duì)10種蒼耳屬雜草籽的鑒別情況。綜合比較可以看出，ITS、ITS2、psbA-trnH序列的鑒定效率較高，遺傳距離較大，鑒定成功率較高，較適合作為鑒別蒼耳屬雜草籽的條形碼序列。

蒼耳屬是菊科的一個(gè)大屬，種類繁多，其非我國(guó)種均為我國(guó)禁止進(jìn)境的植物檢疫性有害生物，危害性極大。由于同屬雜草種外部形態(tài)極其相近，在口岸檢疫時(shí)易誤鑒定為國(guó)內(nèi)種而引入。DNA 條形碼是依靠 DNA 片段上承載的序列信息的不同，對(duì)物種進(jìn)行快速、準(zhǔn)確鑒定的新方法，在本研究中國(guó)際DNA條形碼鑒定聯(lián)盟組織CBOL推薦的rbcL、matK序列均較保守，序列種內(nèi)種間的重疊較大，兩者均不適用于10種蒼耳屬雜草籽的鑒定。psbA-trnH是進(jìn)化速率較快的基因間隔區(qū)之一，兩端具有較為保守的psbA和trnH，該序列的顯著特點(diǎn)是存在較多的插入/缺失，用psbA-trnH序列上的信息位點(diǎn)成功鑒定了國(guó)內(nèi)蒼耳與國(guó)外蒼耳籽，但無(wú)法有效區(qū)分國(guó)外蒼耳籽。本研究中北美蒼耳、球果蒼耳、賓州蒼耳、沃式蒼耳、西方蒼耳有22個(gè)堿基的缺失，其他蒼耳屬植物均無(wú)。在此情況下，采用ITS、ITS2序列進(jìn)行補(bǔ)充研究，ITS和ITS2 是位于nrITS 序列上5.8S～26S 核糖體 DNA 的內(nèi)轉(zhuǎn)錄間隔區(qū)，具備了 5.8S和26S兩端的保守區(qū)設(shè)計(jì)通用引物，該序列最早被應(yīng)用于系統(tǒng)發(fā)育和物種的分類研究中[20-22]，包含豐富的變異位點(diǎn)和信息位點(diǎn)，可以用于對(duì)物種進(jìn)行鑒定。研究發(fā)現(xiàn)，10種蒼耳屬雜草籽psbA-trnH 的種間變異、鑒定成功率均與ITS2序列持平，但均小于ITS序列。本研究中ITS在10種蒼耳屬雜草籽的PCR擴(kuò)增效率和測(cè)序成功率均高于95%，基于BLAST1和 Distance 方法的鑒定效率達(dá)到90%以上，通過(guò)分析10種蒼耳屬雜草籽ITS種間 K2P 遺傳距離和NJ樹可以看出，ITS序列均可準(zhǔn)確鑒定蒼耳等10種蒼耳屬雜草籽。

結(jié)果表明，采用ITS序列鑒定中，國(guó)內(nèi)蒼耳僅與刺蒼耳遺傳距離較近，區(qū)分度較低，但國(guó)內(nèi)蒼耳與其他檢疫性蒼耳均具有顯著區(qū)別，可以通過(guò)ITS序列進(jìn)行鑒定區(qū)分，且區(qū)分度較高；NJ樹分析均可將各物種區(qū)分開來(lái)，且差異顯著。psbA-trnH序列也能同時(shí)單獨(dú)鑒定出10種蒼耳屬雜草籽，能明顯區(qū)分出國(guó)內(nèi)蒼耳與檢疫性蒼耳，但各檢疫蒼耳區(qū)分度較低。兩者在鑒定能力上可互為補(bǔ)充，組合鑒定可以得出一個(gè)更為準(zhǔn)確的鑒定結(jié)論。綜上所述，采用ITS和psbA-trnH可以作為組合DNA條形碼應(yīng)用于蒼耳屬雜草籽的檢疫工作當(dāng)中。

參考文獻(xiàn)：

[1]呂益濤，侯海宮，蘇耀海，等. 蒼耳屬植物的鑒別研究[J]. 中國(guó)中藥雜志，2001，26（1）：17-20.

[2]林春華，唐賽春，韋春強(qiáng)，等. 廣西來(lái)賓市外來(lái)入侵植物的調(diào)查研究[J]. 雜草科學(xué)，2015，33（1）：38-44.

[3]熊 穎，劉啟德，宓穗卿. 蒼耳子化學(xué)研究進(jìn)展[J]. 廣東藥學(xué)，2005，15（6）：65-68.

[4]高銳紅，龔大坤，張星明. 蒼耳子及其偽品單喙蒼耳子的比較鑒別[J]. 中藥材，2001，24（6）：410-411.

[5]Hebert P D N，Cywinska A，Ball S L. Biological identifications through DNA barcodes[J]. Proceedings of the Royal Society of London B：Biological Sciences，2003，270（1512）：313-321.

[6]陳士林，宋經(jīng)元，姚 輝，等. 藥用植物DNA條形碼鑒定策略及關(guān)鍵技術(shù)分析[J]. 中國(guó)天然藥物，2009（5）：322-327.

[7]寧淑萍，顏海飛，郝 剛，等. 植物DNA條形碼研究進(jìn)展[J]. 生物多樣性，2008，16（5）：417-425.

[8]丁 潛，邢光東，胡肄農(nóng)，等. 識(shí)別杜洛克豬個(gè)體身份的DNA條形碼[J]. 江蘇農(nóng)業(yè)學(xué)報(bào)，2014，30（5）：1058-1063.

[9]徐 晗，孫旻旻，吳品珊，等. DNA條形碼、形態(tài)學(xué)與地理分布相結(jié)合的植物果實(shí)（種子）鑒定方法——以口岸截獲的外來(lái)入侵植物犬薔薇（Rosa canina L.）為例[J]. 雜草科學(xué)，2015，33（2）：26-31.

[10]任保青，陳之端. 植物DNA條形碼技術(shù)[J]. 植物學(xué)報(bào)，2010，45（1）：1-12.

[11]閆化學(xué)，于 杰. DNA條形碼技術(shù)在植物中的研究現(xiàn)狀[J]. 植物學(xué)報(bào)，2010，45（1）：102-108.

[12]Group C P W，Hollingsworth P M，F(xiàn)orrest L L，et al. A DNA barcode for land plants[J]. Proceedings of the National Academy of Sciences，2009，106（31）：12794-12797.

[13]Song J，Yao H，Li Y，et al. Authentication of the family Polygonaceae in Chinese pharmacopoeia by DNA barcoding technique[J]. Journal of Ethnopharmacology，2009，124（3）：434-439.

[14]Lahaye R，van der Bank M，Bogarin D，et al. DNA barcoding the floras of biodiversity hotspots[J]. Proceedings of the National Academy of Sciences，2008，105（8）：2923-2928.

[15]Sass C，Little D P，Stevenson D W，et al. DNA barcoding in the cycadales：testing the potential of proposed barcoding markers for species identification of cycads[J]. PloS One，2007，2（11）：e1154.

[16]朱英杰，陳士林，姚 輝，等. 重樓屬藥用植物DNA條形碼鑒定研究[J]. 藥學(xué)學(xué)報(bào)，2010（3）：376-382.

[17]Chase M W，Cowan R S，Hollingsworth P M，et al. A proposal for a standardised protocol to barcode all land plants[J]. Taxon，2007，56（2）：295-299.

[18]Kress W J，Erickson D L. A two-locus global DNA barcode for land plants：the coding rbcL gene complements the non-coding trnH-psbA spacer region[J]. PLoS One，2007，2（6）：e508.

[19]Barcodes P P. Wanted：a barcode for plants[J]. Science，2007，318（5848）：190.

[20]Keller A，Schleicher T，Schultz J，et al. 5.8S-28S rRNA interaction and HMM-based ITS2 annotation[J]. Gene，2009，430（1）：50-57.

[21]Müller T，Philippi N，Dandekar T，et al. Distinguishing species[J]. RNA，2007，13（9）：1469-1472.

[22]Miao M，Warren A，Song W，et al. Analysis of the internal transcribed spacer 2（ITS2） region of scuticociliates and related taxa（Ciliophora，Oligohymenophorea） to infer their evolution and phylogeny[J]. Protist，2008，159（4）：519-533.