董艷輝,于宇鳳,溫 鑫,王亦學(xué),聶園軍,侯麗媛,李亞莉,劉 江,任 元,王育川,曹秋芬,吳慎杰,王 斌,秦永軍

(1.山西省農(nóng)業(yè)科學(xué)院 生物技術(shù)研究中心，山西 太原 030031；2.農(nóng)業(yè)部黃土高原作物基因資源與種質(zhì)創(chuàng)新重點(diǎn)實(shí)驗(yàn)室，山西 太原 030031；3.山西省農(nóng)業(yè)科學(xué)院 農(nóng)業(yè)資源與經(jīng)濟(jì)研究所，山西 太原 030031；4.山西省農(nóng)業(yè)科學(xué)院 農(nóng)作物品種資源研究所，山西 太原 030031;5. 山西省農(nóng)業(yè)科學(xué)院 農(nóng)業(yè)環(huán)境與資源研究所，山西 太原 030031)

連作是農(nóng)民常用的一種簡(jiǎn)單、長(zhǎng)期、有效地提高土地利用率獲取經(jīng)濟(jì)利益的一種做法，藜麥屬于莧科植物，作為一種引進(jìn)于南美的全營(yíng)養(yǎng)作物，營(yíng)養(yǎng)價(jià)值高，被列為全球十大健康食物之一[1]，2008年，我國(guó)山西開始引種并試種成功，山西靜樂作為藜麥產(chǎn)業(yè)的首個(gè)生產(chǎn)基地，目前種植面積已經(jīng)超過0.2萬hm2，且甘肅、內(nèi)蒙古、青海、吉林、河北、寧夏等省份種植面積也在逐漸擴(kuò)大[2]。由于藜麥作為一種有著較高營(yíng)養(yǎng)價(jià)值的新型雜糧，市場(chǎng)前景比較廣闊，相對(duì)于其他作物有較高的經(jīng)濟(jì)收入，農(nóng)民對(duì)種植藜麥的熱情比較高，在土地資源有限的情況下對(duì)藜麥進(jìn)行連年種植，然而，藜麥第2年重茬種植就會(huì)出現(xiàn)出苗率低、植株生長(zhǎng)緩慢、病蟲害加重、產(chǎn)量急劇下降等現(xiàn)象，采取土壤簡(jiǎn)單消毒、增施底肥等措施后效果也不明顯。藜麥連作障礙現(xiàn)象在藜麥大規(guī)模種植地區(qū)基本都存在，導(dǎo)致個(gè)體種植戶和藜麥種植公司有著不同程度的損失，科研人員進(jìn)行了各方面的嘗試，但因?yàn)閷?duì)藜麥各方面的研究還處于初級(jí)階段，目前還沒有很好的解決辦法。馬紅梅等[3]對(duì)連作靈芝土壤微生物研究表明，連作改變了靈芝栽培土壤中原來的微生態(tài)平衡，特別是病原真菌改變了原有的土壤微生境，但具體是哪一種菌還有待進(jìn)行進(jìn)一步的分子鑒定。薛超等[4]總結(jié)了大量前人對(duì)連作機(jī)制的研究后認(rèn)為，施用有機(jī)肥尤其是將有機(jī)肥與功能微生物相結(jié)合制成微生物有機(jī)肥后施用，對(duì)土傳病害有一定的防治作用，但效果不是太穩(wěn)定。

藜麥作為抗逆性較強(qiáng)的新型雜糧作物，我國(guó)引進(jìn)的時(shí)間較短，對(duì)其生理生化特性以及遺傳機(jī)制了解還不夠深入，面對(duì)產(chǎn)生的連作障礙問題，很難達(dá)到對(duì)癥下藥，但根際土壤是根際微生物直接作用于植物的重要場(chǎng)所，也是根際微生物受植物根系和分泌物影響最直接的區(qū)域，作物根際微生物能夠分解并轉(zhuǎn)化根際養(yǎng)分供給植物根系吸收，提高植物對(duì)生物和非生物脅迫的抵抗力[5-7]。夏圍圍等[8]比較了新一代高通量測(cè)序技術(shù)與傳統(tǒng)的變性梯度凝膠電泳(Denaturing gradient gel electrophoresis,DGGE)指紋圖譜技術(shù)，評(píng)價(jià)了這2種技術(shù)研究土壤微生物群落結(jié)構(gòu)的優(yōu)缺點(diǎn)，在不同的微生物分類水平，高通量測(cè)序草地土壤檢測(cè)到22門、54綱、60目、131科、350屬；而DGGE僅檢測(cè)到6門、9綱、8目、10科、10屬、表明DGGE顯著低估了土壤微生物的群落組成。

新一代高通量測(cè)序技術(shù)的發(fā)展擺脫了研究土壤微生物分離培養(yǎng)的瓶頸，使得根際微生物和連作障礙之間的關(guān)系再次成為研究的熱點(diǎn)。本試驗(yàn)采用新一代高通量測(cè)序技術(shù)，以靜樂縣娑婆鄉(xiāng)連作藜麥試驗(yàn)地為研究對(duì)象，分析了第1年種植藜麥與重茬種植藜麥的根際土壤樣品的細(xì)菌種群豐度和多樣性的變化，并對(duì)細(xì)菌功能進(jìn)行初步預(yù)測(cè)，試圖探索藜麥連作障礙與土壤根際細(xì)菌之間的關(guān)系，以期了解藜麥連作對(duì)根際土壤細(xì)菌菌群多樣性的影響，探究藜麥連作障礙機(jī)制，旨在從藜麥根際土壤微生物方面入手為藜麥的生產(chǎn)實(shí)踐和可持續(xù)種植提供理論基礎(chǔ)。

1 材料和方法

1.1 試驗(yàn)地概況

試驗(yàn)地選取我國(guó)藜麥之鄉(xiāng)靜樂縣娑婆鄉(xiāng)，于2015-2016年連續(xù)在同一地塊進(jìn)行試驗(yàn)。該地屬于溫帶季風(fēng)氣候，夏季暖熱且晝夜溫差大，冬季寒冷，年降雨量380～500 mm，無霜期120～135 d。

1.2 試驗(yàn)方法

用無菌毛刷將藜麥根際土壤刷進(jìn)無菌自封袋內(nèi)，去除石塊、大顆粒等雜質(zhì)，同時(shí)在根周圍用無菌鏟取土壤放入無菌袋內(nèi)，去除石塊、大顆粒等雜質(zhì)。每個(gè)重復(fù)取1個(gè)土壤樣品，每個(gè)土壤樣品取3個(gè)取樣點(diǎn)然后混合。取樣后立刻放入冰盒帶回實(shí)驗(yàn)室，根際土壤放入-80 ℃冰箱用于DNA提取。

試驗(yàn)設(shè)2個(gè)處理，處理1為種植1 a的藜麥田(2016年)，該處理設(shè)3次重復(fù)，編號(hào)分別為11，12，13，處理2為重茬種植的藜麥田(2015，2016年)，該處理設(shè)3次重復(fù)，編號(hào)分別為C11、C12、C13。

1.3 土壤細(xì)菌的16S rRNA基因測(cè)序

1.3.1 DNA的提取及16SrRNA基因V4區(qū)片段的擴(kuò)增 利用土壤DNA提取試劑盒提取6個(gè)樣本的DNA，針對(duì)16SrRNA基因的V4區(qū)域進(jìn)行PCR擴(kuò)增預(yù)試驗(yàn)，然后進(jìn)行大量PCR擴(kuò)增試驗(yàn)。上游引物為530F(5′-AYTGGGYDTAAAGNG-3′)；下游引物為805R(5′-TACNVGGGTATCTAATCC-3′)；PCR反應(yīng)體系為：Q5高保真DNA聚合酶：0.25 μL，反應(yīng)緩沖液 5 μL，高GC緩沖液 5 μL，dNTP(10 mmol/L)0.5 μL，模板 DNA 1 μL，正向引物(10 μmol/L)1 μL，反向引物(10 μmol/L)1 μL，水11.25 μL；PCR擴(kuò)增程序：98 ℃預(yù)變性2 min；98 ℃變性15 s，55 ℃退火30 s，72 ℃延伸30 s，共25個(gè)循環(huán)；最后72 ℃終延伸5 min結(jié)束。

1.3.2 PCR產(chǎn)物純化與定量 對(duì)通過PCR聚合酶鏈?zhǔn)椒磻?yīng)獲得不同土壤樣品的PCR產(chǎn)物后，進(jìn)行瓊脂糖凝膠電泳分析，針對(duì)目標(biāo)條帶進(jìn)行割膠回收，獲得純化的PCR產(chǎn)物。

1.3.3 文庫(kù)驗(yàn)證及測(cè)序 利用熒光分光光度計(jì)方法定量DNA，通過Agilent 2100對(duì)富集片段進(jìn)行電泳完成質(zhì)量控制，驗(yàn)證DNA文庫(kù)片段大小及分布規(guī)律(Agilent 2100 bioanalyzer，Agilent 2100;Agilent High Sensitivity DNA Kit,Agilent, 5067-4626)，然后送至上海派森諾生物科技有限公司，在Illumina-MiSeq平臺(tái)上進(jìn)行高通量測(cè)序。

1.4 數(shù)據(jù)處理與分析

為獲得更為精準(zhǔn)、高質(zhì)量的生物學(xué)信息，首先對(duì)原始數(shù)據(jù)進(jìn)行質(zhì)量控制，獲得最終用于分析的序列，然后應(yīng)用QIIME軟件，根據(jù)序列97%的相似度，將序列歸并并劃分為多個(gè)OUTs，并利用軟件Mothur計(jì)算豐富度指數(shù)Chao1和ACE以及多樣性指數(shù)Simpson和Shannon，并進(jìn)行Alpha多樣性分析，同時(shí)進(jìn)行細(xì)菌群落分布、聚類分析和細(xì)菌功能預(yù)測(cè)分析。

2 結(jié)果與分析

2.1 不同處理土壤樣本測(cè)序結(jié)果及質(zhì)量控制

高通量測(cè)序結(jié)果顯示，6個(gè)樣本共獲得有效序列276 400條(Sequences)，運(yùn)用QIIME軟件(Quantitative Insights Into Microbial Ecology，v1.8.0，http://qiime.org/)剔除序列長(zhǎng)度<150 bp的序列和 5′端引物錯(cuò)配堿基數(shù)> 1的序列以及含有連續(xù)相同堿基數(shù)>8的序列，然后調(diào)用USEARCH軟件檢查并剔除嵌合體序列。按照優(yōu)化標(biāo)準(zhǔn)，去除不合格序列后共得到高質(zhì)量序列251 346條，每個(gè)樣本高質(zhì)量序列占有效序列的比例都在89%以上(表1)。高質(zhì)量序列用于樣本間微生物豐富度和多樣性的評(píng)估。

表1 各樣本序列數(shù)據(jù)統(tǒng)計(jì)Tab.1 Statistics of sample sequences datas

稀疏曲線(Dilution curve)揭示了每個(gè)樣品的取樣深度，可以用來評(píng)價(jià)當(dāng)前測(cè)序深度是否足以反映該群落樣本所包含的微生物多樣性。從圖1可以看出，當(dāng)測(cè)序量超過40 000讀長(zhǎng)時(shí)，仍然會(huì)有新的OUTs出現(xiàn)，但是由于曲線已經(jīng)趨于平緩，表明該樣品取樣基本合理，在真實(shí)的環(huán)境中細(xì)菌群落結(jié)構(gòu)的置信度比較高，能夠比較真實(shí)地反映土壤中的細(xì)菌群落多樣性，即當(dāng)前測(cè)序深度足以反映該群落樣本所包含的細(xì)菌群落多樣性。

圖1 各土壤樣品的稀疏曲線(97%)Fig.1 Dilution curve of each soil sample(97%)

2.2 不同處理土壤細(xì)菌群落豐富度和多樣性分析

利用Alpha多樣性分析土壤中細(xì)菌群落的豐富度和多樣性，各土壤樣本的群落豐富度指數(shù)Chao1和ACE以及兼顧群落均勻度指數(shù)Shannon和Simpson如表2所示。

表2 土壤樣本中細(xì)菌豐富度和多樣性指數(shù)Tab.2 Index of the bacteria OUTs′ abundance and diversity in soil samples

就物種豐富度而言，種植1 a藜麥的土壤樣本較種植2 a的土壤樣本Chao1指數(shù)平均值增加16.4%，ACE指數(shù)平均值增加22.9%，表明種植2 a的土壤中細(xì)菌的物種豐富度大量減少，這與大部分研究結(jié)果一致。說明細(xì)菌物種豐富度減少可能是導(dǎo)致藜麥連作障礙的原因之一。 Simpson指數(shù)顯示，種植1 a的土壤較種植2 a的土壤Simpson指數(shù)平均值增加0.4%，變化不明顯，表明二者土壤中細(xì)菌群落的均勻度和優(yōu)勢(shì)OUTs數(shù)變化不明顯；而Shannon指數(shù)顯示，種植1 a的土壤較種植2 a的土壤Shannon指數(shù)平均值增加2.3%，表明種植1 a藜麥的土壤細(xì)菌群落多樣性和稀有OUTs較種植2 a藜麥的土壤要高，說明重茬種植后不僅導(dǎo)致土壤細(xì)菌多樣性下降，還導(dǎo)致部分稀有細(xì)菌種群多樣性減少；土壤細(xì)菌在土壤養(yǎng)分循環(huán)的各個(gè)環(huán)節(jié)占有重要地位，其數(shù)量的下降和群落結(jié)構(gòu)的變化均會(huì)引起微生物功能的失調(diào)，進(jìn)而導(dǎo)致土壤養(yǎng)分和肥力的下降，這與徐雪雪[9]的研究結(jié)果基本一致。

2.3 不同處理土壤菌群的分類學(xué)組成分析

將各分類水平的群落組成數(shù)據(jù)根據(jù)分類單元的豐度分布或樣本間的相似程度加以聚類分析并繪制熱圖(圖2)，種植1 a藜麥的土壤樣品和種植2 a藜麥的土壤樣品分別聚類在不同的類別，說明樣品一致性較好，有一定的代表性。聚類結(jié)果表明，重茬種植后倫茨氏菌屬(Lentzea)、溶桿菌屬(Lysobacter)、中慢生根瘤菌屬(Mesorhizobium)、極地單胞菌屬(Polaromonas)多樣性增多；分枝桿菌屬(Mycobacterium)、藤黃單胞菌屬(Luteimonas)、芽單胞菌屬(Gemmatimonas)、浮霉菌屬(Planctomyces)等細(xì)菌多樣性減少。增多和減少的細(xì)菌種群分屬不同的類別。其中，藤黃單胞菌屬細(xì)菌為革蘭氏陰性、好氧的棒狀桿菌，目前共有10個(gè)種。李廣寧等[10]分離出1株藤黃單胞菌屬細(xì)菌 HF-1127，通過研究該菌株的形態(tài)、培養(yǎng)特征、生理生化特征和遺傳特性，表明該細(xì)菌具有抗菌活性，其代謝產(chǎn)物對(duì)大腸桿菌和金黃色葡萄球菌有抗性；范曉陽[11]從海水中分離出一種深淵藤黃單胞菌，研究表明，其可以降解明膠、淀粉、吐溫-20、吐溫-40以及吐溫-80等物質(zhì)；黃佩蓓等[12]對(duì)浮霉?fàn)罹芯空J(rèn)為，浮霉?fàn)罹鷮偌?xì)菌是一類具有重要生態(tài)作用的環(huán)境微生物，參與環(huán)境中的碳和氮的循環(huán)。

圖2 結(jié)合聚類分析的屬水平群落組成熱圖Fig.2 Heat map of the community composition in genus level combined with the cluster analysis

2.4 不同處理土壤樣品序列的菌群代謝功能預(yù)測(cè)

大量研究表明，土壤中的群落微生物功能組成與環(huán)境因子的相關(guān)度要大于物種組成與環(huán)境因子的相關(guān)度，在相似的環(huán)境條件下，雖然微生物群落的功能相似，但是行使功能的微生物組成差異反而較大，因此，在揭示土壤環(huán)境中的微生物群落組成的基礎(chǔ)上，最終要揭示微生物群落的功能。運(yùn)用PICRUSt軟件通過將現(xiàn)有的群落組成數(shù)據(jù)與KEGG數(shù)據(jù)庫(kù)中微生物代謝功能的類別對(duì)比進(jìn)行群落樣本功能預(yù)測(cè)[13]。由圖3可知，重茬種植后，行使環(huán)境信息處理功能類別中編碼細(xì)胞膜運(yùn)輸?shù)墓δ芑虼罅繙p少；行使遺傳信息處理功能類別中編碼翻譯、復(fù)制和修復(fù)的功能基因大量減少；行使代謝功能類別中編碼異生素生物降解和代謝、聚糖生物合成和代謝的功能基因大量減少，編碼核苷酸代謝、萜類化合物和聚酮化合物的代謝、輔因子和維生素的代謝等功能基因少量減少；編碼信號(hào)轉(zhuǎn)導(dǎo)和脂質(zhì)代謝的功能基因有少量增加。

圖3 群落樣本的代謝功能預(yù)測(cè)圖Fig.3 Prediction map of the metabolic function in community samples

3 結(jié)論與討論

作物連作障礙產(chǎn)生的原因多且比較復(fù)雜，其中，土壤微生物豐富度和多樣性的變化是導(dǎo)致連作的主要原因之一。在農(nóng)業(yè)生態(tài)系統(tǒng)中，土壤微生物多樣性對(duì)維持生態(tài)系統(tǒng)的平衡具有重要作用。由于技術(shù)的限制，早期對(duì)土壤微生物研究需要依賴于人工分離培養(yǎng)，而大部分土壤微生物無法在實(shí)驗(yàn)室條件生存，因此阻礙了人們對(duì)微生物的了解，但是隨著技術(shù)的進(jìn)步，尤其是下一代高通量測(cè)序技術(shù)的發(fā)展，不再依賴人工培養(yǎng)，使研究者們能夠同時(shí)對(duì)多種微生物基因組進(jìn)行測(cè)序，通過對(duì)土壤根際微生物的豐富度和多樣性的研究有助于更好地研究微生物與作物生長(zhǎng)發(fā)育的關(guān)系，為克服作物連作障礙提供數(shù)據(jù)支撐[14]。

藜麥重茬種植會(huì)導(dǎo)致根際土壤的細(xì)菌種群數(shù)量和多樣性都有不同程度的降低，同時(shí)也導(dǎo)致致病菌數(shù)量的增多，使根際土壤細(xì)菌種群朝著不利于植物生長(zhǎng)的方向發(fā)展，這與大量作物連作障礙研究結(jié)果基本一致[15]。Bever等[16]根據(jù)多年的研究結(jié)果提出了土壤-植物反饋機(jī)制，該機(jī)制根據(jù)最終的結(jié)果又分為正反饋、負(fù)反饋和中性反饋。連作障礙的發(fā)生就是由于土壤與植物連續(xù)負(fù)反饋的結(jié)果。同時(shí)藜麥重茬種植后細(xì)菌種群結(jié)構(gòu)也發(fā)生了變化，倫茨氏菌屬、溶桿菌屬、中慢生根瘤菌屬、叢毛單菌屬等細(xì)菌數(shù)量增多；分枝桿菌屬、藤黃單胞菌屬、芽單胞菌屬浮霉菌屬等細(xì)菌數(shù)量減少。據(jù)報(bào)道，溶桿菌屬細(xì)菌對(duì)線蟲和植物病害(革蘭氏陰性菌、革蘭氏陽性細(xì)菌、真菌等)具有生防作用，同樣具有固氮作用的中慢生根瘤菌屬細(xì)菌數(shù)量也有增加[17]。但是由于相關(guān)藜麥連作障礙方面的研究尚未見報(bào)道，因此可能是由于藜麥根際微環(huán)境比較適合溶桿菌屬細(xì)菌以及中慢生根瘤菌屬細(xì)菌的繁殖，導(dǎo)致了其數(shù)量的增加；但是參與土壤中碳和氮循環(huán)的浮霉菌屬細(xì)菌數(shù)量卻降低了，藜麥連作打破了藜麥根際土壤細(xì)菌之間的平衡，這可能導(dǎo)致藜麥根際土壤細(xì)菌功能性紊亂，從而負(fù)反饋于植株體本身，導(dǎo)致了藜麥的連作障礙，但具體的作用機(jī)制和相互之間的關(guān)系還需要進(jìn)一步研究。

土壤根際微生物群落結(jié)構(gòu)的差異會(huì)表現(xiàn)在基因表達(dá)方面的差異，藜麥重茬根際土壤中編輯復(fù)制和修復(fù)以及外源性物質(zhì)降解和代謝的功能基因大量減少，受環(huán)境條件的影響越來越大，導(dǎo)致細(xì)菌對(duì)外來入侵的抵抗能力大大減小，原來有益的細(xì)菌群落豐度顯著下降，病原菌細(xì)菌群落大量增殖，因此，導(dǎo)致重茬種植藜麥時(shí)生長(zhǎng)受阻，土傳病害增加，產(chǎn)量急劇下降。也有研究表明[18]，連作不僅僅會(huì)導(dǎo)致有益細(xì)菌群落豐度降低，有益真菌的數(shù)量也會(huì)大幅度降低，而病原菌真菌的數(shù)量則會(huì)大幅度增加。由于藜麥連作障礙土壤根際細(xì)菌多樣性方面的研究尚未見報(bào)道，并且本研究所采用的高通量測(cè)序技術(shù)的第2代測(cè)序技術(shù)，對(duì)細(xì)菌的分類僅限于屬的水平，無法定位到某一種細(xì)菌，且同一屬的不同細(xì)菌種在土壤中可能會(huì)行使不同的功能，而不同屬的細(xì)菌也可能在土壤中行使同樣的功能，研究所得到的結(jié)果可能會(huì)有偏差。而隨著第3代測(cè)序技術(shù)(單分子測(cè)序技術(shù))的進(jìn)一步完善和發(fā)展，打開了一扇精準(zhǔn)研究土壤微生物的大門，該技術(shù)已經(jīng)能夠?qū)⑼寥牢⑸锒ㄎ坏椒N的水平，對(duì)根際微生物的研究將會(huì)更加全面和細(xì)化，現(xiàn)在無法解決的問題可能會(huì)在不久的將來得到解決。

本研究中，利用高通量測(cè)序技術(shù)對(duì)2個(gè)處理6個(gè)樣本的土壤微生物進(jìn)行了16SrRNA基因的V4區(qū)測(cè)序，然后對(duì)測(cè)序數(shù)據(jù)進(jìn)行了處理和分析，該項(xiàng)研究的樣本量和測(cè)序深度在藜麥的相關(guān)研究中尚屬首次，能夠更加全面地揭示許多利用傳統(tǒng)分子生物學(xué)技術(shù)研究土壤微生物無法解決的問題，結(jié)果中獲得的土壤微生物多樣性變化基本上能夠揭示藜麥連作后土壤細(xì)菌種群多樣性的變化，并且預(yù)測(cè)了土壤細(xì)菌種群代謝功能的變化，為今后研究藜麥連作障礙機(jī)制提供參考，但由于本研究中測(cè)序深度即稀釋曲線并沒有達(dá)到完全飽和，處理中還有部分細(xì)菌種群沒有被發(fā)現(xiàn)。因此，本研究結(jié)果僅代表各個(gè)土壤樣品中的大部分細(xì)菌種群，今后仍有待開展更深層次的研究。