蓋 凱，唐 葉，權永剛，陳愛民，陳全家，吳家和，胡保民，

（1.新疆農(nóng)業(yè)大學農(nóng)學院農(nóng)業(yè)技術重點實驗室棉花教育部工程研究中心，烏魯木齊830052；2.農(nóng)業(yè)部西北棉花品種創(chuàng)新重點實驗室九圣禾種業(yè)股份有限公司，昌吉831100；3.中國科學院微生物研究所植物基因組國家重點實驗室，北京100101）

棉花是一種重要的經(jīng)濟作物［1］。棉花纖維是紡織 工業(yè)的重要原料［1］。由大麗輪枝菌（Verticillium dahli?ae）引起的棉花黃萎病是影響當前棉花生產(chǎn)的最重要因素之一，因此如何防治黃萎病發(fā)生成為當前棉花生產(chǎn)的一個主要問題。

棉花黃萎病俗稱棉花的“癌癥”，在生產(chǎn)上難以防治［2］。發(fā)生面積涉及我國各大種植棉區(qū)。然而，在不同年份、不同地區(qū)甚至不同棉花單株黃萎病發(fā)生程度各不相同；棉花黃萎病發(fā)病程度的不均勻性是其主要特征，因此，往往采取病指來對抗感材料進行分類［3?5］。一般認為棉花黃萎病發(fā)病不均勻原因包括：土壤中病原菌含量不同、單株之間的生長狀況不同、農(nóng)田小氣候差異等［6?12］。

土壤中黃萎病菌含量應該包括其孢子、菌絲和微菌核等［13?14］。測定土壤中大麗輪枝菌的含量可能對棉花抗黃萎病生產(chǎn)有著重要的意義。魏鋒等［15］利用土樣水篩提純微菌核和實時熒光定量檢測對土壤中黃萎病微菌核量；羅舒文等［10］利用培養(yǎng)土壤中微菌核進行定量檢測，探明棉花黃萎病菌的分布區(qū)域及在土壤中的動態(tài)變化過程；李艷［16］通過菌絲生長速率法，對棉花黃萎病菌的菌絲及微菌核進行統(tǒng)計。但是，對于棉田棉株發(fā)病程度差異尚鮮見報道。

本研究利用盆栽棉花作為實驗模型，探討不同時間點病原菌含量的變化與植株發(fā)病程度直接關系，并結(jié)合一種抑制大麗輪枝菌生長的抑制菌——枯草芽孢桿菌和常年發(fā)病嚴重的病圃地土壤為正負對照，探索棉花黃萎病發(fā)病程度與土壤中大麗輪枝菌含量的相關性，以期為棉花抗病生產(chǎn)提供一定的理論依據(jù)和實踐保證。

1 材料與方法

1.1 材料

1.1.1 供試材料和棉花栽培

本實驗所用大麗輪枝菌（V.dahliae）菌株“V991”以及枯草芽孢桿菌（Bacillus subtilis）于本實驗室保存。陸地棉（Gosspium hirsutum）新陸早57 由九圣禾種業(yè)股份公司提供，盆栽供試土壤樣品采集于新疆昌吉市九圣禾產(chǎn)業(yè)園區(qū)內(nèi)，沒有大麗輪枝菌污染；病圃地土壤來自新疆瑪納斯縣六戶地鎮(zhèn)黃萎病抗性鑒定試驗地。

將接菌和對照土壤裝入花盆（15 cm × 15 cm）中，然后將盆埋入土中10 cm；在盆里直接播種，最后定苗為每盆1 株。澆水、施肥等農(nóng)田管理和日常農(nóng)田管理一致；實驗重復3次。

1.1.2 相關試劑盒與引物

DNA 提取試劑盒DNeasy powerSoil Kit 購自德國QIAGEN 公司；DNA 聚合酶2×Super Pfx Master Mix 購自聚合美生物公司；其余試劑盒購自北京擎科生物有限公司。測序及引物合成均由北京擎科公司完成。

1.2 方法

1.2.1 黃萎病菌和枯草芽孢桿菌的活化培養(yǎng)及接種

將用于本實驗的黃萎病菌菌株“V991”和枯草芽孢桿菌活化，分別接種至Czapek 和LB 培養(yǎng)液中，使大麗輪枝菌孢子液濃度為1×108mL-1，枯草芽孢桿菌菌液濃度OD600為1.39。將收集到的孢子液及菌液接種土壤，按10 kg土壤50 mL菌量接種。

1.2.2 土壤黃萎病菌基因組DNA提取

將盆栽和溫室中接種了黃萎病菌和枯草芽孢桿菌的土壤及棉花病圃地土壤樣品分別稱取0.1 g，根據(jù)土壤樣品基因組提取試劑盒DNeasy powerSoil Kit 說明書提取樣品DNA進行qPCR分析。

1.2.3 實時熒光定量PCR 定量分析土壤中大麗輪枝菌含量

qPCR 特異引物參照中國農(nóng)業(yè)大學張力群教授提供的大麗輪枝菌glyceraldehyde?3?phosphate dehydroge?nase（G3PD）基因序列。通過Primer premier 5軟件設計引物qG3PD?F/R，該對引物在NCBI 上進行全物種Blast，結(jié)果顯示為大麗輪枝菌特異序列。實時熒光定量PCR 按照擎科SYBR Green qPCR Master mix 試劑盒說明書進行操作，以qG3PD?F/R 為引物擴增，其中qG3PD?F：5'?cacggcgtcttcaagggt?3'，qG3PD?R：5' cagtg?gactcgacgacgtac?3'，反應體系為20 μL：SYBR Green qP?CR Master mix 10 μL，稀釋模板DNA 20 ng／μL，引物qG3PD?F/R 各1.0 μL，加ddH2O 補齊至20 μL。擴增條件為：95 ℃預變性1 min；95 ℃變性10 s，58 ℃退火5s，72 ℃延伸15 s，40 個循環(huán)。將已知濃度的GFP 質(zhì)粒作為熒光定量的質(zhì)粒作為內(nèi)參，進行實時熒光定量PCR擴增［17］。

1.2.4 水洗涂布菌落計數(shù)和PCR鑒定大麗輪枝菌落

稱取土樣5 g，根據(jù)Yang 等［18］的實驗方法，對土壤中的黃萎病菌進行水篩，并涂布于PDA 培養(yǎng)基上，待菌落形成時，利用PCR對大麗輪枝菌菌落進行鑒定，然后計數(shù)，實驗重復3次。

常規(guī)PCR反應所用大麗輪枝菌特異引物序列由中國農(nóng)業(yè)大學張力群教授提供，引物擴增基因為G3PD。其中G3PD?F：5'?gtctatcatgtgctggggttctcc?3'，G3PD?R：5'?gagatgatgaccttcttggcgcc?3'，擴增體系為25.0 μL：模板DNA 20 ng／μL，引物G3PD?F/R 各0.5 μL，2×Super Pfx Master Mix 10 μL，加ddH2O 補齊至25 μL。擴增條件：98 ℃預變性30 s；98 ℃變性10 s，56 ℃退火30 s，72 ℃延伸30 s，32個循環(huán)；72 ℃再延伸10 min。PCR 產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測。

1.2.5 發(fā)病程度分析

參照農(nóng)業(yè)部棉花抗性鑒定方法與評價標準進行病級劃分，發(fā)病等級分為0、1、2、3 和4，統(tǒng)計5 個病級的棉花植株，根據(jù)植株的發(fā)病等級計算接菌植株的發(fā)病指數(shù)（Disease index，DI），病指計算參照Wang 等［4］并進行適當修改。

將浸染后的植株截取子葉節(jié)以下相同部位的莖稈進行徒手斜剖，觀察其木質(zhì)部褐化程度，參照Fradin等［2］方法，將徒手橫剖的植株莖稈進行表面消毒，于PDA培養(yǎng)基上恢復培養(yǎng)。

1.2.6 數(shù)據(jù)統(tǒng)計分析

通過Excel 2010、GraphPad 5和SPSS 21.00等軟件進行數(shù)據(jù)統(tǒng)計分析。

2 結(jié)果與分析

2.1 棉花植株之間發(fā)病程度差異分析

棉花植株發(fā)病程度差異較大，幼苗接種等量的黃萎病菌時，植株發(fā)病程度存在明顯差異（圖1），發(fā)病程度5 個級別植株分別為26%、11%、17%、6%和40%，發(fā)病指數(shù)為51%。統(tǒng)計分析表明，不同發(fā)病程度植株的株高、始節(jié)高度、果枝數(shù)、有效果指枝數(shù)、結(jié)鈴數(shù)和鈴重等農(nóng)藝性狀呈現(xiàn)顯著性差異（表1）。

圖1 棉花個體之間的發(fā)病程度差異性Figure 1 The difference of disease degree among cotton individuals

表1 不同盆栽植株農(nóng)藝性狀Table 1 Agronomic character of different potted plants

2.2 土壤大麗輪枝菌含量分析

利用試劑盒提取盆栽土壤中的總DNA，電泳分析表明這些DNA 質(zhì)量較好（圖2），滿足qPCR 分析要求。qPCR 分析結(jié)果表明：不同取樣時間土壤中大麗輪枝菌的相對含量存在顯著差異；同一時間點，不同盆栽土壤病菌量差異也為顯著水平（圖3）。

圖2 土壤提取DNA凝膠分析Figure 2 Gel analysis of DNA extracted from soil

圖3 土壤中大麗輪枝菌相對含量Figure 3 The relative contents of V.dahliae in soil

使用水洗涂布和PCR 鑒定相結(jié)合進一步測定土壤中大麗輪枝菌含量，盆栽土壤浸提液稀釋后平板真菌菌落生長情況見圖4（a），對不同類型菌落進行PCR擴增鑒定［圖4（b）］，并進行菌落計數(shù)。統(tǒng)計分析表明不同時間點和盆栽土壤病原菌含量呈顯著差異，并且隨著時間的后移差異越大（表2）。

圖4 病圃土壤中真菌菌落分析和鑒定Figure 4 Analysis and identification of fungal colonies in disease nursery soil

表2 盆栽棉株土壤不同時間點大麗輪枝菌含量的方差分析Table 2 Variance analysis of the content of V.dahliae in potted cotton plant soil at different time points

2.3 植株發(fā)病程度與土壤大麗輪枝菌含量的相關性

對不同植株的發(fā)病程度與其土壤含菌量qPCR 結(jié)果進行協(xié)方差分析，以土壤樣品的發(fā)病級數(shù)為縱軸（y），Ct值為橫軸（x），建立的曲線方程為y = -0.159x +8.289,該方程的相關系數(shù)為0.869，說明不同盆栽苗的發(fā)病程度與不同時間點土壤里病原菌含量存在相關性，隨著時間的推移相關性增大（圖5）。

圖5 植株發(fā)病程度與土壤大麗輪枝菌含量的相關性分析Figure 5 The correlation analysis between the disease degree of individual plant and the V.dahliae content in soil

2.4 枯草桿菌抑制大麗輪枝菌生長影響植株發(fā)病程度

為進一步證明土壤大麗輪枝菌含量與植株發(fā)病程度的相關性，一個能抑制大麗輪枝菌生長的菌株（枯草芽孢桿菌）被應用。在使用枯草芽孢桿菌的盆里，棉花植株發(fā)病很輕或不發(fā)病（圖6）。用涂布平板和PCR 分析，枯草芽孢桿菌處理盆栽土壤大麗輪枝菌含量顯著低（圖6）。這些數(shù)據(jù)再次證明了土壤中大麗輪枝菌含量低，其對應的植株發(fā)病程度輕或不發(fā)病。

圖6 枯草芽孢桿菌抑制黃萎病的發(fā)病程度分析Figure 6 Analysis of Verticillium wilt in habited by Bacillus subtilis

為從另一個方向證明土壤大麗輪枝菌含量與植株發(fā)病程度的相關性，選取新疆常年發(fā)病病圃土壤作為對照，這些病圃中棉花植株發(fā)病均為4 級程度［圖7（a）］。每個病圃地的微菌核量可以達到每克土壤50個以上，qPCR 分析結(jié)果表明這些土壤中的大麗輪枝菌相對含量達到1.25×105以上，為盆栽土壤含量的100倍（圖7），0～4 級植株發(fā)病程度分析分別為3.85%、6.25%、4.52%、8.94%和76.44%，發(fā)病指數(shù)為85.92%，這些數(shù)據(jù)表明土壤中大麗輪枝菌含量直接決定棉花植株發(fā)病程度。

圖7 病圃土壤大麗輪枝菌含量分析Figure 7 Analysis on the content of V.dahliae in the soil of the disease nursery

3 討論與結(jié)論

研究假設土壤里大麗輪枝菌含量存在局部不均勻是造成地面上棉花單株發(fā)病程度不同的主要原因，然而在培養(yǎng)箱、溫室和病圃中接種等量病原菌時，也會出現(xiàn)單株發(fā)病程度不同的現(xiàn)象［4］，利用盆栽棉花接種等量大麗輪枝菌分析不同時間點每盆土壤菌量的變化和其棉株發(fā)病程度的關系，結(jié)果表明隨著時間的推移，不同盆之間的大麗輪枝菌含量呈現(xiàn)顯著差異，這些差異可能是造成不同盆間棉株發(fā)病程度差異的主要原因。與前人［19?20］研究得出的大麗輪枝菌在土壤里含量的變化受多種因素的影響，比如多年種植棉花、氣候和耕作處理等結(jié)果一致。因此，研究結(jié)果證明土壤中病原菌含量是不斷變化的，甚至不同棉株下的局部土壤菌量也是變化的，這可能是影響其上棉株發(fā)病程度不同的主要原因。研究利用枯草芽孢桿菌抑制大麗輪枝菌生長作為實驗的正對照，枯草芽孢桿菌能夠抑制土壤大麗輪枝菌的生長，其機制目前還不清楚［21?22］。研究結(jié)果顯示枯草芽孢桿菌降低土壤中大麗輪枝菌的含量，減輕了棉花黃萎病的發(fā)病程度，也說明了棉花植株的發(fā)病程度與土壤里的病原菌含量呈相關性。同時研究也利用多年發(fā)病嚴重的病圃地土壤為負對照研究菌含量與發(fā)病程度的關系，結(jié)果發(fā)現(xiàn)病圃地土壤含有高的病原菌含量，這與其他報道［23?24］結(jié)果一致。也再次說明土壤里大麗輪枝菌的含量是影響棉花單株發(fā)病程度的一個重要因素。