甘肅小麥禾谷孢囊線蟲(chóng)rDNA-ITS和28S rDNA-D2/D3區(qū)序列特征及ITS-RFLP分析

葉文興, 徐秉良, 彭德良, 黃文坤

(1.甘肅農(nóng)業(yè)大學(xué)草業(yè)學(xué)院,蘭州 730070;2.中國(guó)農(nóng)業(yè)科學(xué)院植物保護(hù)研究所,植物病蟲(chóng)害生物學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100193)

禾谷孢囊線蟲(chóng)(Heterodera avenae Wollenweber,1924,cereal cyst nematode,CCN)是溫帶禾谷類作物上的世界性重要病原線蟲(chóng),最早于1874年Kühn在德國(guó)發(fā)現(xiàn),隨后在英國(guó)、俄羅斯、挪威、澳大利亞、加拿大、日本、印度、美國(guó)、新西蘭、中國(guó)等40多個(gè)國(guó)家發(fā)生危害,嚴(yán)重影響禾谷類作物尤其是麥類作物的產(chǎn)量和質(zhì)量。自1989年我國(guó)首次在湖北天門(mén)縣發(fā)現(xiàn)禾谷孢囊線蟲(chóng)危害小麥后[1-2],目前已發(fā)現(xiàn)在河北、河南、北京、山西、內(nèi)蒙古、青海、湖北、安徽、山東等11省市都有CCN的分布,CCN已成為嚴(yán)重威脅我國(guó)小麥生產(chǎn)的主要病害[3-5]。

隨著分子診斷技術(shù)的快速發(fā)展,已探索出許多線蟲(chóng)種群分類鑒定的分子診斷方法,利用分子技術(shù)探索不同種群 rDNA-ITS和 rDNA 28S的D2/D3區(qū)的序列特征已成為近年來(lái)線蟲(chóng)分子診斷的研究熱點(diǎn)。28S是核糖體大亞基rRNA的轉(zhuǎn)錄基因,28S的D2/D3區(qū)是 rDNA中進(jìn)化最快的一個(gè)編碼區(qū),具有高度保守的序列將其分割,能夠在種的水平上進(jìn)行種間的分子鑒別[6]。Subbotin等[7]從28S rDNA-D2/D3區(qū)分析了墊刃目(Ty lenchida)線蟲(chóng)在亞目、科、種系統(tǒng)發(fā)育中與其他種群的親緣關(guān)系。李佳等[8]對(duì)香蕉穿孔線蟲(chóng) rRNA基因的D2/D3區(qū)序列進(jìn)行了分析,蔣立琴等[9]也成功地利用 28S rDNA-D2/D3區(qū)特征將傘滑刃線蟲(chóng)屬部分種群進(jìn)行區(qū)分,所以,D2/D 3區(qū)常常用于研究不同地理來(lái)源的不同線蟲(chóng)種間的遺傳多樣性[10]。同樣,rDNA-ITS區(qū)序列和 ITS-PCR-RLFP圖譜也能夠有效地區(qū)分和鑒別線蟲(chóng)種間的親緣關(guān)系和遺傳特征。ITS序列是一個(gè)多用途的遺傳標(biāo)記,位于18S和28S核糖體DNA基因的重復(fù)簇之間,中間被5.8S核糖體DNA基因分開(kāi)[11]。Subbotin等[12]擴(kuò)增了5種孢嚢線蟲(chóng)rDNA的28S基因的D3擴(kuò)展區(qū)和ITS1-5.8S-ITS2區(qū),結(jié)果發(fā)現(xiàn) G.rostochiensis和G.pallida D3擴(kuò)展區(qū)的核苷酸序列相同,并用RFLP分析成功區(qū)分了來(lái)自南美茄科植物上的3種孢嚢線蟲(chóng) G.rostochiensis、G.tabacum、G.pallida和1個(gè)未描述的種,同時(shí)構(gòu)建了遺傳上親緣關(guān)系的系統(tǒng)發(fā)育樹(shù)。Ferris等[13]分析來(lái)自俄國(guó)H.filip jevi群體的rDNA-ITS序列,證實(shí)俄國(guó)不同地區(qū)的 H.filip jevi基因一致,與瑞典的 H.avenae相似性達(dá)到99.7%,支持了瑞典的H.filip jevi和 H.avenae是近緣種,同時(shí)也提示形態(tài)相似,地理隔離的群體有可能是同屬種。在RFLP圖譜分析中,有時(shí)限制性片段長(zhǎng)度的總和大于沒(méi)有酶切的陽(yáng)性對(duì)照擴(kuò)增產(chǎn)物,這說(shuō)明禾谷孢囊線蟲(chóng)群體的ITS區(qū)域存在異質(zhì)性[14]。Subbotin等[15]對(duì)來(lái)自世界多個(gè)國(guó)家和地區(qū)的禾谷孢囊線蟲(chóng)群體 rDNA中ITS區(qū)的PCR-RFLP分析,發(fā)現(xiàn)了ITS有異質(zhì)性,命名為 ITS“A”型和 ITS“B”型。鄭經(jīng)武等[16]對(duì)來(lái)自山西太谷的禾谷孢囊線蟲(chóng)群體與國(guó)外上述兩種ITS型群體rDNA中ITS區(qū)的PCR-RFLP比較分析,發(fā)現(xiàn)中國(guó)禾谷孢囊群體的ITS特征系是一個(gè)新的類型“C”型,用 H in fI和 Tru9 I可將中國(guó)禾谷孢囊線蟲(chóng)與其他禾谷孢囊線蟲(chóng)分開(kāi)。彭德良等[17]用 H in fⅠ、AluⅠ和 RsaⅠ將中國(guó)禾谷孢囊線蟲(chóng)和摩洛哥禾谷孢囊線蟲(chóng)群體存在差異,AvaⅠ和H in dⅢ2種限制性內(nèi)切酶不能酶切禾谷孢囊線蟲(chóng)的ITS產(chǎn)物。

本研究應(yīng)用PCR技術(shù)對(duì)采自甘肅6個(gè)市縣的7個(gè)禾谷孢囊線蟲(chóng)群體和1個(gè)河南群體、1個(gè)安徽群體的28S rDNA-D2/D3區(qū)和ITS區(qū)進(jìn)行了序列特征分析,并對(duì)甘肅禾谷孢囊線蟲(chóng) ITS區(qū)進(jìn)行了RLFP圖譜分析,以揭示甘肅省內(nèi)不同地理種群間的親緣關(guān)系。

1 材料與方法

1.1 禾谷孢囊線蟲(chóng)孢囊的分離方法

在2009年2月和 6月,從甘肅省文縣、天水、甘谷、榆中、武威、民勤以及河南省安陽(yáng)縣和安徽省蚌埠市等8個(gè)縣市(表1)采集受害的小麥根及根際土壤,取200 g土壤在小桶容器中使土塊分散,用強(qiáng)水流沖洗,用木棍攪動(dòng),沉淀片刻后,將水溶液通過(guò)20目和60目的樣品篩,重復(fù)以上步驟3次,用小水流輕輕清洗60目篩殘余物,將其60目篩中殘余物淋洗到三角瓶中,用濾紙過(guò)濾,在解剖鏡下挑取飽滿可育的孢囊冷藏備用[17]。

1.2 禾谷孢嚢線蟲(chóng)DNA的提取

挑取1個(gè)飽滿的孢囊,放入滅菌的Eppendorf管中,加入14μL滅菌的重蒸水,在液氮中速凍1min取出。用75%乙醇消毒的玻璃棒在Eppendorf管中轉(zhuǎn)動(dòng)至冰融化,使孢囊破碎研磨1 min后加入3μL的10×PCR buffer(Promega,USA)和3μL蛋白酶 K溶液(600mg/ML)(TaKaRa,Biotech),然后在-20℃下冷凍至少2 h。之后,將Eppendorf管從冰箱中取出,置65℃下溫育1.5 h;接著在95℃10min,使蛋白酶K變性,最后12 000 r/min下離心1 min,取上清DNA懸浮液于-20℃保存?zhèn)溆肹17]。

表1 禾谷孢嚢線蟲(chóng)采樣及注冊(cè)的ITS和D2/D3序列

1.3 28S rDNA D2/D3區(qū)域PCR擴(kuò)增

本研究采用 50μL的 PCR反應(yīng)體系,10×PCR-buffer(含 Mg2+)5.0 μL,dNTPs 4 μL,引物D2A 1.5μL,引物 D3B 1.5μL,Taq DNA 聚合酶(5 U/μL)0.4μL,模板DNA 5.0 μL,滅菌重蒸水補(bǔ)足到50μL,同時(shí)設(shè)置無(wú)DNA模板的陰性對(duì)照。試驗(yàn)所用的引物為 D2A(ACAAGTACCGTGAGGGAAAGTTG)和D3B(TCGGAAGGAACCAGCTACTA)[18-20]。在 Eppendorf PCR擴(kuò)增儀上進(jìn)行擴(kuò)增。PCR擴(kuò)增的反應(yīng)條件為：94℃預(yù)變性4Min,94℃30s,55℃45s,72℃2 min,循環(huán)35次,72℃延伸10 min后保存在4℃條件下。取5μL擴(kuò)增產(chǎn)物加1μL 6×Loading buffer在1.5%瓊脂糖凝膠上電泳,用0.5×TAE作為電泳緩沖液,100 V下電泳45 Min,用EB染色,用BIO-RAD成像儀在紫光燈下觀測(cè)并照相。

1.4 28S rDNA D2/D3區(qū)域PCR產(chǎn)物的克隆、測(cè)序和序列分析

將PCR產(chǎn)物進(jìn)行回收純化,純化后的產(chǎn)物與pGEM-T easy vector連接12 h,然后轉(zhuǎn)入感受態(tài)細(xì)胞DH 5a中,用熱擊法進(jìn)行轉(zhuǎn)化后涂在含有Ampicillin(100 Mg/L)、IPTG(24 mg/L)和 X-gal(200mg/L)的 LB培養(yǎng)基平板上,涂勻后倒置入37℃培養(yǎng)箱中12～16 h,進(jìn)行藍(lán)白斑篩選。挑取單個(gè)白斑于含有Ampicillin(100 mg/L)的液體 LB中搖培過(guò)夜,經(jīng)菌落PCR鑒定后,將含有正確片段的菌液進(jìn)行測(cè)序。所得序列用 Chromas、DNAMAN分析比對(duì),在GenBank注冊(cè),并在NCBI下載已知近緣種序列,用MEGA 4.0中的UPGMA法構(gòu)建禾谷孢囊線蟲(chóng)群體的系統(tǒng)發(fā)育樹(shù),進(jìn)行親緣關(guān)系分析。

1.5 rDNA-ITS的PCR擴(kuò)增

本研究采用 50μL的PCR反應(yīng)體系,10×PCR-buffer(含 Mg2+)5.0 μL,dNTPs 4 μL,引物TW 81 0.3μL,引物 AB28 0.3μL,Taq DNA 聚合酶(5 U/μL)0.4μL,模板DNA 3.0 μL,滅菌重蒸水補(bǔ)足到50μL,同時(shí)設(shè)置無(wú)DNA模板的陰性對(duì)照。試驗(yàn)所用的引物為 TW 81(GTTTCCGTAGGTGAACCTGC)和AB28(A TATGCTTAAGTTCAGCGGGT)[21]。在Eppendorf PCR擴(kuò)增儀上進(jìn)行擴(kuò)增。PCR擴(kuò)增的反應(yīng)條件為：94℃預(yù)變性5Min,94℃45s,55℃1 min,72℃2Min,循環(huán)35次,72℃延伸10 min后保存在4℃條件下。取5μL擴(kuò)增產(chǎn)物加1μL 6×Loading bu ffer在1.5%瓊脂糖凝膠上電泳,用 0.5×TAE作為電泳緩沖液,100 V下電泳45Min,用EB染色,用BIO-RAD成像儀在紫光燈下觀測(cè)并照相。

1.6 ITS的PCR擴(kuò)增產(chǎn)物的克隆、測(cè)序和序列分析

PCR產(chǎn)物的回收純化、連接、轉(zhuǎn)化、測(cè)序等方法與上述方法相同,所得序列在GenBank注冊(cè),并在NCBI下載已知近緣種序列,用 MEGA 4.0進(jìn)行UPGMA聚類分析。

1.7 ITS-RFLP分析

分別取5μL純化后的PCR產(chǎn)物用9種限制性內(nèi)切酶 AluⅠ 、AvaⅠ 、Bsh1236Ⅰ 、Bsu RⅠ 、C foⅠ 、Hin fⅠ、RsaⅠ、TaqⅠ和 MvaⅠ酶切。AluⅠ 、AvaⅠ、Bsh1236 Ⅰ、Bsu R Ⅰ、C fo Ⅰ 、H in fⅠ 、Rsa Ⅰ、MvaⅠ在37℃水浴中16 h,Taq I在65℃水浴16 h。酶切后的DNA加1μL 10×Loading buffer在1.5%瓊脂糖凝膠上電泳,用0.5×TAE作為電泳緩沖液,90 V下電泳1.5 h,用EB染色,用BIORAD成像儀在紫外光下觀測(cè)和照相。

2 結(jié)果與分析

2.1 28S rDNA D2/D3區(qū)域和rDNA-ITS的PCR擴(kuò)增結(jié)果

經(jīng)PCR擴(kuò)增的禾谷孢嚢線蟲(chóng)28S rDNA的D2/D3區(qū)和ITS區(qū)序列,擴(kuò)增片段長(zhǎng)度分別約為780 bp(圖1)和1 040 bp(圖2)。序列測(cè)定結(jié)果表明,GsWx、GsMq、GsTs1、HnAy1、GsYz、AhBb、GsTs2 7個(gè)禾谷孢囊線蟲(chóng)群體的D2/D3區(qū)的序列長(zhǎng)度為782 bp,G sGg、G sWw 2個(gè)禾谷孢囊線蟲(chóng)群體D2/D3區(qū)序列長(zhǎng)度為784 bp;這9個(gè)禾谷孢囊線蟲(chóng)的ITS區(qū)的序列長(zhǎng)度都是1 045 bp,將所測(cè)的禾谷孢囊線蟲(chóng)D2/D3區(qū)和 ITS區(qū)序列在GenBank上注冊(cè),得到了9個(gè)序列注冊(cè)號(hào)(表1)。

圖1 9個(gè)禾谷孢囊線蟲(chóng)(Heteroderaavenae)28S rDNA-D2/D3區(qū)擴(kuò)增產(chǎn)物

圖2 9個(gè)禾谷孢囊線蟲(chóng)(Heteroderaavenae)rDNA-ITS擴(kuò)增產(chǎn)物

2.2 序列分析結(jié)果

2.2.1 28S rDNA-D2/D3區(qū)序列分析和聚類分析

用DNAMAN軟件對(duì)比這9個(gè)禾谷孢囊線蟲(chóng)的D2/D3區(qū)的序列,其序列的一致性達(dá)99.77%,其中GsGg與GsWw種群的序列沒(méi)有差異,比其他7個(gè)群體的序列在第136、137兩個(gè)位點(diǎn)插入2個(gè)堿基(T和C);GsTs1群體在第284位點(diǎn)、GsYz群體的第331位點(diǎn)有堿基的替換 ;G sW x、GsWq、GsTs2、A h-Bb1、HnAy1 5個(gè)群體的序列無(wú)差異。這表明在甘肅禾谷孢囊線蟲(chóng)群體間D2/D3區(qū)的一致性很高、具有很高的保守性。

通過(guò)NCBI中的BLAST比對(duì)分析,將本研究中所涉及的9個(gè)禾谷孢囊線蟲(chóng)群體的D2/D3序列與GenBank中下載的孢囊屬(Heterodera)、球形孢囊屬(Globodera)、棘皮屬(Cactodera)的13個(gè)群體的D2/D3序列進(jìn)行比對(duì)分析,用UPGMA方法構(gòu)建了系統(tǒng)發(fā)育樹(shù)(圖3)。從圖3中可以看出,在遺傳距離為0.02的水平上,可分為3類：第1類有GsW x、GsMq、G sTs1、H nAy1、GsYz、AhBb1 、G sTs2、Gs-Gg、GsWw 、新西蘭的 H.auck landica(DQ328688)和 H.litoralis(DQ328691)、敘利亞的H.latipons(DQ328687)、印度的 H.sorghi(DQ328689)、比利時(shí)的 H.salixophila(DQ328690)、美國(guó)的 H.glycines(DQ328692)15個(gè)群體,其中7個(gè)甘肅禾谷孢囊線蟲(chóng)群體與 HnAy1、A hBb1、新西蘭的 H.aucklandica(DQ328688)為同一組群,遺傳距離很小,并且與敘利亞的H.latipons(DQ328687)有著96%的置信度;第2類有德國(guó)的C.cacti(DQ328702)、波蘭的G.pallida(EU855119)和G.rostochiensis(EU-855120)、荷蘭的 G.pa llida(AY 592992)、美國(guó)的G.tabacum(GQ294492)5個(gè)群體,棘皮屬和球形孢囊屬的遺傳距離比較接近;第3類有比利時(shí)的H.urticae(DQ328696)、德國(guó)的 H.goetting iana(DQ-328697)2個(gè)群體,這2個(gè)群體之間的遺傳距離小于0.01,置信度在100%。D2D3區(qū)的系統(tǒng)發(fā)育樹(shù)表明孢囊屬、球形孢囊屬、棘皮屬之間有著較近的遺傳關(guān)系,在孢囊屬有兩個(gè)大的分支,這些遺傳特征之間的差異與寄主有著密切的關(guān)系,孢囊屬在禾谷類寄主和豆科類寄主就存在明顯差異,不同寄主間有小種的分化,但在不同屬之間又有著近緣關(guān)系。

圖3 用UPGM方法構(gòu)建的28S rDNA-D2/D3區(qū)的系統(tǒng)發(fā)育樹(shù)

2.2.2 rDNA-ITS區(qū)序列分析和聚類分析

用DNAMAN比對(duì)后的結(jié)果為：禾谷孢囊線蟲(chóng)rDNA-ITS區(qū)僅在第620和第897位點(diǎn)有堿基替換,序列的一致性達(dá)99.9%。在NCBI上用BLAST比對(duì)分析后,在GenBank下載與本研究9個(gè)禾谷孢囊線蟲(chóng)有近緣關(guān)系的16個(gè)群體,包括Avenae組的有效種有澳大利亞的 H.australis(AY148395)、中國(guó)(AY 148382)、法國(guó)(AY 148374)、德國(guó)(AY148360)、以色列(AY 148366)、印度(AF274397)的 H.avenae、德國(guó)的 H.Mani(AY148377)和 H.p ratensis(AY 148383)、美國(guó)的 H.filip jevi(GU 079654)、英國(guó)的H.arenaria(AF274396);Schachtii組有伊朗的 H.glycines(AF498387)、比利時(shí)的 H.schachtii(AY 166438);Goettingiana組的有德國(guó)的H.goettingiana(AF274414)、荷 蘭的 H.cruciferae(AF274411)、意大利的 H.carotae(AY 347917);還有比利時(shí)的H.urticae(AF274412)。用MEGA 4.0的UPGMA方法構(gòu)建ITS區(qū)的系統(tǒng)發(fā)育樹(shù),其遺傳距離為0.05(圖4)。在圖4中,本研究的9個(gè)禾谷孢囊線蟲(chóng)群體和在GenBank下載的16個(gè)孢囊線蟲(chóng)群體可分為3類：本研究中的9個(gè)禾谷孢囊線蟲(chóng)群體與上述Avenae組的10個(gè)群體為一類,這19個(gè)群體間的遺傳距離很小,除美國(guó)的 H.filip jevi(GU 079654)群體外,其余18個(gè)群體很接近,其中澳大利亞的H.australis(AY148395)群體與中國(guó)通州的H.avenae(AY148382)和本研究的9個(gè)群體可分為同一組群;很明顯,伊朗的 H.glycines(AF498387)、比利時(shí)的H.schachtii(AY166438)為一類,其置信度為100;Goettingiana組3個(gè)群體和比利時(shí)的 H.urticae(AF27-4412)為一類。這3類群體中,Avenae組與其余2類群體存在明顯的異質(zhì)性,在同寄主群體間存在地域和寄主對(duì)群體分化的影響。

2.3 ITS-RFLP分析

用9種限制性內(nèi)切酶酶切7個(gè)甘肅禾谷孢囊線蟲(chóng)群體的rDNA-ITS擴(kuò)增產(chǎn)物后,共產(chǎn)生22個(gè)清晰可見(jiàn)的酶切片段(圖5)。用 AluⅠ酶切7個(gè)甘肅禾谷孢囊線蟲(chóng)群體的rDNA-ITS擴(kuò)增產(chǎn)物,均產(chǎn)生2個(gè)相同的560 bp和480 bp RFLP片段(圖5a);C foⅠ酶切后,7個(gè)群體均產(chǎn)生相同的3個(gè) 750、160、110 bp的 RFLP 片段(圖 5b);Bsh1236Ⅰ酶切后產(chǎn)生900 bp和140 bp的2個(gè)片段分布都相同(圖5c);AvaⅠ不能酶切這7個(gè)甘肅禾谷孢囊線蟲(chóng)的rDNA-ITS產(chǎn)物(圖5d);RsaⅠ酶酶切后產(chǎn)生相同的2個(gè)720 bp和320 bp的片段(圖5e);Bsu RⅠ酶切片段產(chǎn)生的3個(gè)RFLP分布圖均相同,片段為 480、380 bp和180 bp(圖 5f);H in fⅠ酶切產(chǎn)生 2個(gè)相同的860 bp和180 bp的片段,由于酶切效率低,有未切完全的殘留片段(1 040 bp)(圖5g);MvaⅠ酶切后產(chǎn)生的3個(gè)片段也相同,分別是420、330 bp和290 bp(圖5h);TaqⅠ產(chǎn)生4個(gè)均相同的400、280、140 bp和 130 bp的 RFLP片段(圖 5i)。

圖4 用UPGM方法構(gòu)建的rDNA-ITS區(qū)的系統(tǒng)發(fā)育樹(shù)

圖5 用9種酶酶切7個(gè)甘肅禾谷孢囊線蟲(chóng)(Heterodera avenae)群體

3 結(jié)論與討論

本研究通過(guò)對(duì)7個(gè)甘肅禾谷孢囊線蟲(chóng)群體和HnAy1、AhBb1群體的28S rDNA-D2/D3區(qū)的序列特征分析發(fā)現(xiàn),7個(gè)甘肅禾谷孢囊線蟲(chóng)群體和HnAy1、AhBb1群體的28S rDNA-D2/D3區(qū)都具有相同的保守性,在序列特征上只有2個(gè)堿基的插入(GsGg、GsWw)和 2個(gè)位點(diǎn)的堿基替換,序列差異很小,序列的一致性達(dá)99.8%。構(gòu)建的28S rDNAD2/D3區(qū)系統(tǒng)發(fā)育樹(shù)將 GsWx、GsMq、GsTs1、HnAy1 、GsYz、AhBb1 、GsTs2 、GsGg 、GsWw 9 個(gè)群體歸為同一分支,同新西蘭的 H.auck landica(DQ328688)有著很近的親緣關(guān)系。

分析rDNA-ITS區(qū)的序列和RFLP圖譜,對(duì)孢囊線蟲(chóng)種的鑒定和近緣種辨別有著重要的意義。Subbotin等[22]研究了70個(gè)H.avenae復(fù)合種群的ITS序列和RFLP圖譜,構(gòu)建了H.avenae復(fù)合種群的系統(tǒng)關(guān)系,揭示出H.avenae是一個(gè)并系分類群,中國(guó)禾谷孢囊線蟲(chóng)不同于其他的 H.avenae群體。本研究中的 GsW x、GsMq、GsTs1、HnAy1、GsYz、AhBb1、GsTs2、GsGg、GsWw 9 個(gè)群體與澳大利亞的 H.australis(AY 148395)、中國(guó)的 H.avenae(AY 148382)為同一分支(圖4),與歐洲和印度的H.avenae或A venae復(fù)合種群的其他群體相比較,其遺傳距離最為接近。所以,可以在分子水平上診斷出7個(gè)甘肅禾谷孢囊線蟲(chóng)群體都在A venae組,且與美國(guó)的H.filip jevi(GU079654)有明顯的差異。同時(shí),在種群的分化上存在地域差異的影響,ITS-RFLP圖譜分析可以直觀地鑒別種內(nèi)和種間的異質(zhì)性,有利于群體間的遺傳多態(tài)性的研究。Subbotin等[15]從禾谷孢囊線蟲(chóng)群體中的ITS檢測(cè)到的兩種不同的單元型,稱之為“A 型”和“B型”,“A 型”不能被 AluⅠ和RsaⅠ這兩種酶酶切,而“B型”能被A luⅠ酶切成2個(gè)片段(560 bp和500 bp),被RsaⅠ酶切產(chǎn)生720、320 bp 2個(gè)片段,因此,AluⅠ和RsaⅠ是區(qū)分禾谷孢囊線蟲(chóng)“A型”和“B型”這兩單元型的主要限制性內(nèi)切酶[23],并且在歐洲禾谷孢囊線蟲(chóng)群體中以ITS“A型”為典型,而在印度禾谷孢囊線蟲(chóng)中是以“B型”為典型的。Zheng等[16]發(fā)現(xiàn) H in fⅠ酶切中國(guó)禾谷孢囊線蟲(chóng)群體后產(chǎn)生的圖譜與“A型”和“B型”存在異質(zhì)性,將其劃分為“C型”,AvaⅠ則不能酶切禾谷孢囊線蟲(chóng)群體。用Bsu RⅠ、TaqⅠ可把 H.avenae和H.filip jevi群體分離出[15],C foⅠ可以將H.ciceri與其他群體分開(kāi),MvaⅠ和Bsh1236Ⅰ可以把 H.hordecalis分離出來(lái)。本研究的ITS-RFLP結(jié)果表明,7個(gè)甘肅禾谷孢囊線蟲(chóng)群體均能被 A luⅠ、RsaⅠ和H in fⅠ酶切,并且酶切結(jié)果與H.avenae“C型”完全一致。本文中的9種酶切結(jié)果與所知文獻(xiàn)完全相符[15-23]。綜上所述,甘肅的這7個(gè)禾谷孢囊線蟲(chóng)群體應(yīng)為禾谷孢囊線蟲(chóng)(H.avenae),ITS-RFLP的圖譜特征應(yīng)為“C型”。

通過(guò)對(duì)9種禾谷孢囊線蟲(chóng)的28S rDNA-D2/D3區(qū)和ITS區(qū)序列以及ITS-RFLP分布型的比較分析,反映了禾谷孢囊線蟲(chóng)種群間和種群內(nèi)的遺傳學(xué)關(guān)系。甘肅禾谷孢囊線蟲(chóng)種群內(nèi)的遺傳距離很小,雖然在個(gè)別群體(GsGg、GsWw)上有幾個(gè)堿基的插入或替換,這并不能影響種群內(nèi)部歸類分析,因?yàn)楦鶕?jù)rDNA基因建立的系統(tǒng)發(fā)育樹(shù)僅僅是一個(gè)基因樹(shù),和物種樹(shù)可能還存在差別[9]?；虻姆只瘯r(shí)間可能早于物種的分化時(shí)間,基因樹(shù)的拓?fù)浣Y(jié)構(gòu)可能不同于物種樹(shù)。因此,把傳統(tǒng)的形態(tài)分類特征、28S rDNA-D2/D3區(qū)和ITS區(qū)以及 RFLP圖譜特征相結(jié)合,對(duì)于禾谷孢囊線蟲(chóng)的分類鑒定會(huì)更加有效。利用28S rDNA-D2/D3區(qū)的保守性和ITS區(qū)的穩(wěn)定性,在相關(guān)的區(qū)段設(shè)計(jì)特異性引物進(jìn)行PCR擴(kuò)增,可用于禾谷孢囊線蟲(chóng)的檢疫、分類鑒定和系統(tǒng)學(xué)研究中的特異性分子檢測(cè)。

[1]陳品三,王明祖,彭德良.我國(guó)禾谷孢囊線蟲(chóng)發(fā)現(xiàn)與鑒定初報(bào)[J].中國(guó)農(nóng)業(yè)科學(xué),1991,24(5)：29.

[2]王明祖,彭德良,吳學(xué)勤.小麥孢囊線蟲(chóng)病的研究I：病原鑒定[J].華中農(nóng)業(yè)大學(xué)學(xué)報(bào),1991,10(4)：352-356.

[3]彭德良,張東升,齊淑華,等.我國(guó)小麥禾谷孢囊線蟲(chóng) Heterodera avenae發(fā)生分布區(qū)域和防治初步研究[C]∥植物保護(hù)研究進(jìn)展.北京：中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社,1995：53-56.

[4]彭德良,李惠霞,王錫鋒 ,等.我國(guó)小麥禾谷孢囊線蟲(chóng)的新發(fā)生分布地區(qū)[M]∥中國(guó)線蟲(chóng)學(xué)研究.(第 2卷)北京：中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版,2008：344-345.

[5]Peng Deliang,Nicol J,Li H M,et al.Cu rren t know ledge of cereal cy st neMatode(Heterodera avenae)on w heat in China[M]∥Riley I T,Nicol JM,Dababat A A.Cereal cyst nematodes：status,research and ou tlook.Ankara：CIMMYT,2009：29-34.

[6]Hajieghrari B,Torabi-Giglou M,Waeyenberge L.Comparative d2/d3 LSU-rDNA sequen ce study of some Iranian Pra tylenchus loosi populations[J].A frican Jou rnalof Biotechnology,2007,6(21)：2458-2466.

[7]Subbotin S A,Sturhan D,Chizhov V N,et al.Phylogenetic analysis of Tylenchida Thorne,1949 as inferred fromD2 and D 3 expansion fragments of the 28S rRNA gene sequences[J].Nematology,2006,8(3)：455-474.

[8]李佳,彭德良,黃文坤.香蕉穿孔線蟲(chóng)28S rRNA基因的D2/D3區(qū)序列分析[J].華中農(nóng)業(yè)大學(xué)學(xué)報(bào),2008,27(5)：1-5.

[9]蔣立琴,鄭經(jīng)武.傘滑刃線蟲(chóng)屬部分種群 28S RNA基因中D 2-D 3區(qū)的特征[J].植物病理學(xué)報(bào),2007,37(6)：588-594.

[10]Hancock JM,Dover G A.Molecular coevolu tion aMong cryp tically siMple expansion segMents of eukaryotic 26S/28S rRNAs[J].Molecular Biology and Evolu tion,1988,5：377-391.

[11]Subbotin S A,Vierstraete A,Ley P D,et al.Phylogenetic relationships within the cyst-forming nematodes(Nematoda,Heteroderidae)based on analysis of sequences fromthe ITS regions of ribosoMal DNA[J].Molecular Phylogenetics and Evolution,2001,21(1)：1-16.

[12]Subbotin SA,Halford PD,Warry A,etal.Variations in ribosomal DNA sequences and phylogeny of G lobodera parasitising solanaceous plants[J].Nematology,2000,2(6)：591-604.

[13]Ferris V R,Subbotin S A,IreholMA,et al.RibosoMal DNA sequence analysis of Heterodera f ilip jevi and H.1atipons fromRussia and comparisons with other nematode isolates[J].Russia Jou rnal of Nematology,1999,7(2)：121-125.

[14]歐師琪,彭德良,李玉.禾谷孢囊線蟲(chóng)(Heteroderaavenae)的分子鑒定及抗性基因研究進(jìn)展[J].植物保護(hù),2008,34(6)：7-11.

[15]Subbotin S A,W aeyenberge L,Irina A,et al.Identification of Heterodera avenae group speciesbyMorphometricsand rDNARFLPs[J].Nematology,1999,1(2)：195-207.

[16]Zheng JW,Subbotin SA,Waeyenberge L,etal.Molecular characterization of Chinese Heteroderag lycines and H.avenae populations based on RFLPs and sequences of rDNA-ITS regions[J].Russian Journalof Nematology,2000,8(2)：109-l13.

[17]彭德良,Subbotin S,MoensM.小麥禾谷孢嚢線蟲(chóng)(Heteroderaavenae)的核糖體基因(rDNA)限制性片段長(zhǎng)度多態(tài)性研究[J].植物病理學(xué)報(bào),2003,33(4)：323-329.

[18]Courtrigh t E M,W all D H,Virginia R A,et al.Nuclear and Mitocond rial DNA sequence diversity in Antarctic nematode Scottnema lindsayane[J].JNeMatol,2000,32：143-153.

[19]Kaplan D T,Tomas W K,Frisse LM,et al.Phylogenetic analysis of geographically diverse Rhadopholus siMilis via rDNA sequence reveals a monomorphicmotif[J].J Nematol,2000,32(2)：134-142.

[20]Tenente G CMV,De Ley P,De Ley IT,et al.Sequenceanalysis of the D2/D3 region of the large subunit rDNA fromdifferent Meloidogyne isolates[J].Nematropica,2004,34(1)：1-11.

[21]Joyce S A,Reid A,Driver F,et al.Application of polymerase chain reaction(PCR)methods to identification of entomopathogenic nematodes[C]∥Bumell A K M,Ehlers R U,Masson JP.Cost 812 Biotechnology：Geneticsof Entomopathogenic Nematode-bacteriuMComplexes.Belgium：European ComMissi,DGXⅡ,1994：178-187.

[22]Subbotin S A,Sturhan D,Rumpenhorst H J,et al.Molecular and Morphological characterisation of the Heterodera avenae species complex(Tylenchida：Heteroderidae)[J].NeMatology,2003,5(4)：515-538.

[23]Subbofin S A,W aeyenberge L,Moens M.Identification of cyst forMing nematodes of the genus Heterodera(Nematoda：H eteroderidae)based on the ribosomalDNA-RFLPs[J].Nematology,2000,2(2)：153-16.