西瓜連作對根際微生物群落的影響

聶園軍 李瑞珍 趙佳 薄曉峰

摘 要： 為明確根際微生物群落演替變化與西瓜連作之間的關(guān)系，調(diào)查2014—2016年連續(xù)3 a（年）種植西瓜的產(chǎn)量和發(fā)病情況，并利用Illumina Miseq平臺對連作西瓜根際土壤微生物群落進行擴增子測序分析。結(jié)果表明，連續(xù)3 a種植西瓜其產(chǎn)量顯著降低而枯萎病發(fā)病率和病情指數(shù)顯著提高。不同種植年限根際微生物群落多樣性存在顯著差異，西瓜的連年種植打破了其根際微生物群落原有的平衡，使得細菌群落多樣性下降而真菌群落多樣性上升。在門水平上，Proteobacteria、Acidobacteria和Ascomycota的相對豐度升高，而Firmicutes、Actinobacteria、Gemmatimonadetes和Chytridiomycota的相對豐度下降。在屬水平上，Bacillus，Lysobacter與Pseudomonas等許多有益微生物相對豐度降低而病原菌所在的Fusarium相對豐度增加，最終表現(xiàn)為西瓜枯萎病爆發(fā)，產(chǎn)量降低。

關(guān)鍵詞： 西瓜； 根際微生物群落； 擴增子測序； 連作障礙

Abstract： Watermelon rhizosphere soil were taken as the experiment objectives， we explored the relationship between rhizosphere microbial communities and continuous cropping. The relationships between yield and disease of watermelon from 2014 to 2016 were investigated. Nine soil samples from three years were sequenced using Illumina Miseq. The continuous cropping decreased the yield， increased the disease incidence and disease index. The diversity of rhizosphere microbial communities in different years were significant difference. The balance of rhizosphere microbial communities was broken because of continuous cropping of watermelon， which resulted the reduction of bacteria community diversity and the increase of fungal community diversity. At the phyla level， the relative abundance of Proteobacteria，Acidobacteria and Ascomycota increased， while the relative abundance of Firmicutes，Actinobacteria，Gemmatimonadetes and Chytridiomycota decreased. At the genus level， some beneficial microbes like Bacillus， Lysobacter， and Pseudomonas decreased and pathogenic microbes （Fusarium） increased. Finally， the Fusarium wilt outbroke and the yield decreased.

Key words： Watermelon； Rhizosphere microbial communities； Amplifier sequencing； Continuous cropping obstacle

我國人均耕地面積較少，經(jīng)濟作物連續(xù)多年在同一地塊種植的連作模式非常普遍，而多年的連作常引起作物生長受阻、土傳病害頻發(fā)[1-3]。西瓜是世界上重要的園藝作物之一，我國栽培總面積為180.15萬hm2，年產(chǎn)量約為7 071.3萬 t[4]?？菸『κ俏鞴线B作過程中最主要病害之一，它是由尖孢鐮刀菌西瓜?；虵usarium oxysporum f. sp. niveum侵染引起的一種土傳真菌病害，在西瓜各個生育期都可能發(fā)病，尤以坐果期和膨瓜期發(fā)病最重[5-7]?？菸∈沁B作西瓜產(chǎn)量降低和品質(zhì)下降的重要原因[8-9]。因此，弄清連作西瓜根際微生物群落演替變化與西瓜生長、病害發(fā)生之間的互作關(guān)系，是解決西瓜連作障礙的關(guān)鍵。

根際土壤微生物群落是土壤生態(tài)系統(tǒng)持續(xù)發(fā)揮作用的重要媒介，是反應(yīng)土壤健康狀態(tài)的重要指標[10]。目前，國內(nèi)外對根際微生物群落的研究主要集中在細菌群落上，根際真菌群落報道較少，如Ling等[9]研究了西瓜連作條件下施用不同肥料根際細菌群落的變化情況；肖蓉等[11]研究了患炭疽病草莓與健康草莓根際土壤細菌群落的多樣性，這些研究沒有系統(tǒng)地反映根際微生物群落整體的演替變化。在根際微生物的研究中，基因測序技術(shù)已成為主流，對根際的細菌和真菌群落都可以進行系統(tǒng)的分析。擴增子高通量測序技術(shù)通過提取根際土壤微生物總DNA，擴增其中的特異性片段并進行測序分析，測序結(jié)果由于可精確定位到種屬而被科研人員廣泛應(yīng)用[12-13]。郝海婷等[14]利用Illumina MiSeq測序平臺研究了堆肥與化肥對蘭州百合根際土壤微生物群落的影響，物種注釋結(jié)果顯示細菌的門、綱、目、科、屬的數(shù)量在施用堆肥的土壤中都大于施用化肥的土壤。擴增子高通量測序技術(shù)無需對根際微生物進行培養(yǎng)，具有周期短、精度高的優(yōu)點，是研究根際微生物群落的有效手段。

筆者以連作西瓜根際微生物群落為研究對象，利用擴增子高通量測序技術(shù)從整體上分析西瓜連作過程中根際微生物群落變化情況，并與西瓜產(chǎn)量和病害等因子進行關(guān)聯(lián)分析，找出影響西瓜生長的關(guān)鍵菌群，明晰根際微生物群落演替變化與西瓜連作障礙之間的關(guān)系，為克服西瓜連作障礙提供理論依據(jù)。

1 材料與方法

1.1 材料

供試植物：西瓜品種為‘早佳（8424），由山西省農(nóng)業(yè)科學(xué)院植物保護研究所提供。

試劑及儀器：美國OMEGA公司E.Z.N.A. Soil DNA Isolation kit試劑盒、瓊脂糖試劑盒、 Qubit 2.0 DNA 檢測試劑盒，上海派森諾生物科技有限公司。DYY-6C電泳儀，北京六一儀器廠；BG-gdsAUTO凝膠成像儀，北京百晶生物技術(shù)有限公司；Illumina Miseq平臺測序，上海派森諾生物科技有限公司。

1.2 方法

1.2.1 西瓜連作田間試驗及土樣采集 2014—2016年在山西省晉中市山西省農(nóng)業(yè)科學(xué)院東陽基地進行西瓜連作試驗。該試驗田從未種植過西瓜，2013年種植作物為玉米。2014年種植西瓜前，試驗田土壤有機質(zhì)含量（w，后同）16.59 g·kg-1，全氮含量1.98 g·kg-1，有效磷含量16.43 mg·kg-1，有效鉀含量232 mg·kg-1，pH 7.08。西瓜品種為‘早佳（8424），每年3月25號營養(yǎng)缽育苗，35 d左右2葉1心時移栽到試驗田，6月25號左右開始采收；分3個小區(qū)，每個小區(qū)寬5 m長10 m，小區(qū)間隔5 m，面積為50 m2，寬窄行種植，寬行行距2 m，窄行0.5 m，株距0.8 m，常規(guī)水肥管理。連續(xù)種植3 a，每年種植前每個小區(qū)施入10 kg復(fù)合肥（mN∶mP∶mK=15∶15∶15）。

每年西瓜收獲后，用土鉆在每個小區(qū)隨機選擇5處根際土壤每處100 g，采集深度為15～20 cm，均勻地混合為一個樣品，然后裝密封袋低溫帶回實驗室，去除植物根系和石塊并過2 mm篩后于-20 ℃冰箱保存，備用；每年采集3個土壤樣品，連續(xù)3 a共采集9個樣品。

1.2.2 連作西瓜病害調(diào)查及產(chǎn)量統(tǒng)計 每年西瓜收獲前7 d對試驗小區(qū)內(nèi)的所有西瓜植株進行觀察，記錄枯萎病的發(fā)生情況，計算發(fā)病率和病情指數(shù)[15-16]。發(fā)病率=發(fā)病株數(shù)/總植株數(shù)×100%。病害等級劃分標準：0級-整株無癥狀；1級-整株無變黃葉片；2級-1～2個葉片萎蔫；3級-全株>1/3葉片萎蔫；4級-全株>1/2葉片萎蔫；5級-全株>3/4葉片萎蔫或死亡。病情指數(shù)=Σ（病害等級×該級對應(yīng)植株數(shù)）/總植株數(shù)。每年收獲期，各小區(qū)西瓜分批全部采收、稱量，按小區(qū)統(tǒng)計產(chǎn)量。

1.2.3 土壤DNA提取、PCR擴增和高通量測序 利用E.Z.N.A. Soil DNA Isolation kit 試劑盒參照其使用說明書提取土壤樣品DNA[17]。對細菌16S rDNA基因V3-V4高變區(qū)進行PCR擴增，引物序列分別為341F（CCTACGGGNGGCWGCAG） 和 805R（GACTACHVGGGTATCTAATCC）；對真菌NS1-FUNGAL區(qū)域進行PCR擴增，引物序列分別為NS1（GTAGTCATATGCTTGTCTC）和FUNGAL（ATTCCCCGTTACCCGTTG）[18]。擴增反應(yīng)程序：94℃預(yù)變性5 min；30個循環(huán)（94 ℃ 45 s；55 ℃ 30 s（16S）/58.5 ℃ 45 s（18S）；72 ℃ 1 min）；72 ℃ 5 min。對PCR產(chǎn)物進行瓊脂糖電泳，并采用瓊脂糖試劑盒回收，然后用 Qubit 2.0 DNA 檢測試劑盒對回收的DNA 精確定量，每個處理DNA取10 ng，最后按1∶1等量混合后采用Illumina Miseq平臺進行測序[19]。

1.2.4 可操作分類單元聚類與物種分類分析 利用Fast QC軟件對測序原始數(shù)據(jù)進行質(zhì)量控制，然后拼接、過濾掉低質(zhì)量的序列，而后基于這些高質(zhì)量序列進行可操作分類單元（operational taxonomic units，OTUs）聚類和物種分類分析，并將 OTU和物種注釋結(jié)合，從而得到每個樣品的OTUs 和分類譜系的基本分析結(jié)果[20]。對物種注釋后，對門和屬的分類水平上對相對豐度大于等于1%的優(yōu)勢物種進行群落結(jié)構(gòu)的統(tǒng)計分析，明晰物種構(gòu)成。

1.2.5 微生物群落Alpha多樣性分析 利用Mothur軟件進行微生物群落的Alpha多樣性分析，包括稀釋性曲線、Chao1指數(shù)和香濃指數(shù)[21-23]。稀釋曲線是從樣本中隨機抽取一定測序量的數(shù)據(jù)，統(tǒng)計它們所代表物種數(shù)目，即OTUs 數(shù)目，用數(shù)據(jù)量與物種數(shù)所構(gòu)建的曲線；當(dāng)曲線趨向平坦時表明樣品的測序數(shù)據(jù)量漸進合理，并間接反映樣品中物種的豐富程度。Chao1指數(shù)和香濃指數(shù)是衡量群落多樣性的重要指標，指數(shù)越大，群落的豐富度和多樣性越高。

1.2.6 土壤樣品與環(huán)境因子相關(guān)性分析 利用冗余分析技術(shù)（redundancy analysis，RDA）對土壤樣品與微生物群落（門水平）和西瓜病害及產(chǎn)量這些環(huán)境因子進行相關(guān)性分析[24]。通過RDA分析可將土壤樣品與環(huán)境因子數(shù)據(jù)互相擬合，并通過置換檢驗來判斷環(huán)境因子與土壤樣品、環(huán)境因子與環(huán)境因子之間的影響是否顯著。RDA分析可以將樣品和環(huán)境因子反映在同一個二維排序圖上，從圖中可以直觀地看出樣品分布和環(huán)境因子間的關(guān)系，確定微生物群落在不同年份的消長變化，明晰顯著影響西瓜病害和產(chǎn)量的菌群。

1.3 數(shù)據(jù)分析

試驗數(shù)據(jù)采用SPSS 17.0統(tǒng)計軟件進行分析，應(yīng)用鄧肯氏新復(fù)極差法檢驗（Duncans new multiple range tests）進行差異顯著性檢驗。

2 結(jié)果與分析

2.1 連作西瓜的枯萎病發(fā)病情況及產(chǎn)量

由表1可知，隨著種植年限的增加，西瓜的發(fā)病率和病情指數(shù)逐年上升而產(chǎn)量逐年下降，且不同年份間差異顯著。2014年，西瓜的發(fā)病率極低僅為2%，產(chǎn)量表現(xiàn)正常。而連續(xù)種植3 a后，西瓜枯萎病爆發(fā)，發(fā)病率與病情指數(shù)急速上升，分別達到61.3% 和1.87，而小區(qū)產(chǎn)量嚴重下降僅有102.6 kg，為第1年產(chǎn)量的29%。

2.2 微生物群落分類分析

在門水平，2014—2016年連作西瓜根際土壤微生物群落的優(yōu)勢物種出現(xiàn)消長變化。從細菌群落看，變形菌門（Proteobacteria）、厚壁菌門（Firmicutes）、浮霉菌門（Planctomycetes）、放線菌門（Actinobacteria）、擬桿菌門（Bacteroidetes）、芽單胞菌門（Gemmatimonadetes）、酸桿菌門（Acidobacteria）、綠彎菌門（Chloroflexi）、疣微菌門（Verrucomicrobia）和硝化螺旋菌門（Nitrospirae）是所有樣品中相對豐度≥1%的10個優(yōu)勢門，變形菌門與酸桿菌門的相對豐度隨著連作年限的增加而逐漸升高，而厚壁菌門、放線菌門與芽單胞菌門的相對豐度隨著連作年限的增加而逐漸下降（圖版A，見彩色插頁第12頁）。從真菌群落看，壺菌門（Chytridiomycota）、接合菌門（Zygomycota）、子囊菌門（Ascomycota）、纖毛菌門（Ciliophora）、擔(dān)子菌門（Basidiomycota）、絲足菌門（Cercozoa）、球囊菌門（Glomeromycota）與羅茲菌門（Rozellomycota）是所有樣品中相對豐度≥1%的8個優(yōu)勢門，其中子囊菌門的相對豐度在3 a中均在70%以上，超過其它微生物群落之和，處于支配地位。隨著連作年限的增加，子囊菌門的相對豐度逐漸升高，而壺菌門的相對豐度逐漸降低（圖版B，見彩色插頁第12頁）。

在細菌屬水平，相對豐度處于前10位的優(yōu)勢屬中，在連作3 a后芽孢桿菌屬（Bacillus）、紅游動菌屬（Rhodoplanes）、溶桿菌屬（Lysobacter）、土壤紅桿菌科（Solirubrobacter）、德沃斯氏菌屬（Devosia）與假單胞菌屬（Pseudomonas）的相對豐度處于最低值；而浮霉菌屬（Planctomyces）和豐祐菌屬（Opitutus）的相對豐度處于最高值（圖版C，見彩色插頁第12頁）。在真菌屬水平，相對豐度處于前10位的優(yōu)勢屬中，在連作3 a后尖孢鐮刀菌屬（Fusarium）、被孢霉屬（Mortierella）與疣孢漆斑菌（Myrothecium）的相對豐度處于最高值；而紅根菌屬（Rhizophlyctis）的相對豐度處于最低值。最重要的是，西瓜枯萎病病原菌所在的尖孢鐮刀菌屬，連作3 a后相對豐度最高（圖版D，見彩色插頁第12頁）。

2.3 連作西瓜土壤微生物群落的多樣性分析

2.3.1 稀釋曲線分析 9個樣品中，細菌16S rRNA序列總數(shù)為436 187（45 603-64 320）；真菌NS1–FUNGAL序列總數(shù)為310 734（25 341-83 499）。土壤樣本細菌和真菌隨著序列數(shù)的增加的稀釋曲線均趨于平緩，說明測序量已經(jīng)基本覆蓋到樣品中所有的物種。隨著種植年限的增加西瓜根際土壤中細菌群落的OTUs顯著降低（圖1-A）；與之相反真菌群落的OTUs隨著種植年限的增加而升高（圖1-B）。

2.3.2 Chao1指數(shù)和香濃指數(shù)分析 連作西瓜根際土壤樣品中細菌的Chao1指數(shù)和香濃指數(shù)總體上均高于真菌，細菌在土壤微生物中占主體地位，微生物群落Chao1指數(shù)和香濃指數(shù)因種植年限不同而產(chǎn)生顯著差異；細菌群落的Chao1指數(shù)和香濃指數(shù)隨種植年限的增加而顯著降低，在連作3 a后達到最低值（4 254.817和8.857）；而真菌群落的Chao1指數(shù)和香濃指數(shù)隨種植年限的增加而顯著升高，在連作3 a后達到最高值（641.462和6.246）。西瓜根際土壤中的細菌和真菌群落多樣性隨連作年限的增加呈現(xiàn)此消彼長的趨勢（表2）。

2.4 土壤樣品與環(huán)境因子相關(guān)性分析

將樣品、環(huán)境因子進行關(guān)聯(lián)分析（圖2），不論細菌類群還是真菌類群，不同連作年限的樣品各自分開，產(chǎn)量與2014年的3個樣品呈正相關(guān)，而與2016年的3個樣品呈負相關(guān)；發(fā)病率和病情指數(shù)的結(jié)果與之相反。從細菌類群看，厚壁菌門、Actinobacteria和Gemmatimonadetes及產(chǎn)量因子和2014年的3個樣品呈正相關(guān)，而與發(fā)病率、病情指數(shù)、2016年的3個樣品呈負相關(guān)；Proteobacteria、Acidobacteria和Planctomycetes的結(jié)果與之相反。從真菌類群看，Chytridiomycota和Glomeromycota及產(chǎn)量因素和2014年的3個樣品呈正相關(guān)，而與發(fā)病率、病情指數(shù)、2016年的3個樣品呈負相關(guān)；而Ascomycota、Ciliophora、Basidiomycota和Cercozoa的結(jié)果與之相反。

3 討 論

根際土壤微生物群落是反應(yīng)土壤健康狀況的重要指標之一，不同的農(nóng)業(yè)耕作模式會影響根際土壤微生物多樣性甚至作物的生長[25-26]。西瓜連作對土壤的影響同樣反應(yīng)在根際微生物群落的變化上。從西瓜根際微生物門、屬水平物種分類上看，所有土壤樣品中的物種類群大體相似但在相對豐度上差別很大，Li等[8]在番茄的研究中也有類似報道。這說明種植的作物或農(nóng)藝條件不同對微生物的物種類群的影響不顯著，但對物種相對豐度的影響尤為顯著。如：Bacillus、Lysobacter與Pseudomonas的相對豐度隨著連作年限的增加而迅速降低，而西瓜枯萎病病原菌所在的屬Fusarium隨著連作年限的增加而迅速累積。Bacillus能分泌多抗菌物質(zhì)來拮抗土傳病原菌，從而保護植株不受侵害[27-28]。Lysobacter也被證明有拮抗病原菌和促進植株生長的作用[29]。Rumberger等[30]報道，Pseudomonas與植株生長密切相關(guān)。因此，可以確認Bacillus、Lysobacter與Pseudomonas這些有益微生物相對豐度的降低和Fusarium相對豐度的升高，與西瓜枯萎病害的發(fā)生呈正相關(guān)。

通過對稀釋性曲線、Chao1指數(shù)和香濃指數(shù)的分析發(fā)現(xiàn)，西瓜根際微生物群落中細菌群落占優(yōu)勢地位，但隨著連作年限的增加呈現(xiàn)細菌群落多樣性降低而真菌群落多樣性升高的現(xiàn)象。這說明微生物群落結(jié)構(gòu)有從細菌性向真菌性轉(zhuǎn)變的趨勢。這一現(xiàn)象在很多研究中得到確認，如Luo和Shen等[31-32]發(fā)現(xiàn)細菌群落多樣性與連作年限成反比，而真菌群落多樣性與連作年限成正比。本試驗亦證明了低豐度的細菌群落多樣性和高豐度的真菌群落多樣性與西瓜枯萎病害呈正相關(guān)。

冗余分析說明西瓜根際微生物類群的演替變化與連作年限、病害發(fā)生和產(chǎn)量情況息息相關(guān)。Proteobacteria、Acidobacteria和Ascomycota相對豐度的增加會抑制西瓜生長，使得枯萎病病情逐年惡化；與之相反Firmicutes、Actinobacteria、Gemmatimonadetes和Chytridiomycota這些菌群的存在會促進西瓜生長，抑制病害發(fā)生。這再次證明了在連續(xù)種植西瓜這一外界刺激影響下，根際微生物群落發(fā)生了定向改變，某些有益微生物減少而土傳病原微生物增加了，而這種變化也影響到了西瓜的農(nóng)藝性狀，使得連作狀態(tài)下西瓜病害爆發(fā)、產(chǎn)量下降，而連作狀態(tài)下的黑椒和桃等作物都有相似的報道[33-34]。以上結(jié)果提示我們，增加土壤中有拮抗作用或促生作用的有益微生物含量，抑制病原菌的生長，使微生物群落結(jié)構(gòu)由“真菌型”轉(zhuǎn)變回“細菌型”，是減少枯萎病危害的有效途徑。

本研究明確了西瓜根際促生抑病功能菌群的物種定位，為篩選功能菌株提供了思路。在今后的研究中，應(yīng)進一步縮短采樣時間間隔，更加精確的分析西瓜連作過程中微生物群落的動態(tài)變化，從而為克服西瓜連作障礙提供理論依據(jù)。

參考文獻

[1] ZHANG M K，HE Z L.Long-term changes in organic carbon and nutrients of an ultisol under rice cropping in southeast China[J].Geoderma，2004，118（3/4）：167-179.

[2] WU X，LI Z G，LIU H J，et al.The effect of long-term continuous cropping of black pepper on soil bacterial communities as determined by 454 pyrosequencing[J].PLoS One，2015，10（8）：1-13.

[3] 王森，段亞楠，孫申義，等.蘋果連作土壤加入酒糟可減輕連作障礙[J].園藝學(xué)報，2017，44（6）：1157-1166.

[4] 黃春艷，卜元卿，單正軍，等.西瓜連作病害機理及生物防治研究進展[J].生態(tài)學(xué)雜志，2016，35（6）：1670-1676.

[5] LING N，XUE C，HUANG Q W，et al.Development of a mode of application of bioorganic fertilizer for improving the biocontrol efficacy to Fusarium wilt[J].BioControl，2010，55（5）：673-683.

[6] ZHAO S，LIU D Y，LING N，et al.Bio-organic fertilizer application significantly reduces the Fusarium oxysporum population and alters the composition of fungi communities of watermelon Fusarium wilt rhizosphere soil[J].Biology and Fertility of Soils，2014，50（5）：765-774.

[7] FAHEEM M，RAZA W，ZHONG W，et al.Evaluation of the biocontrol potential of Streptomyces goshikiensis YCXU against Fusarium oxysporum f.sp.niveum[J].Biological Control，2015，81：101-110.

[8] LI J G，REN G D，JIA Z J，et al.Composition and activity of rhizosphere microbial communities associated with healthy and diseased greenhouse tomatoes[J].Plant and Soil，2014，380（1/2）：337-347.

[9] LING N，DENG K Y，SONG Y，et al.Variation of rhizosphere bacterial community in watermelon continuous mono-cropping soil by long-term application of a novel bioorganic fertilizer[J].Microbiological Research，2014，169（7/8）：570-578.

[10] PHILIPPOT L，RAAIJMAKERS J M，LEMANCEAU P，et al.Going back to the roots：the microbial ecology of the rhizosphere[J].Nature Reviews Microbiology，2013，11（11）：789-799.

[11] 肖蓉，曹秋芬，聶園軍，等.基于高通量測序患炭疽病草莓根際與健康草莓根際細菌群落的比較研究[J].中國農(nóng)學(xué)通報，2017，33（11）：14-20.

[12] LI C G，LI X M，KONG W D，et al.Effect of monoculture soybean on soil microbial community in the Northeast China[J].Plant Soil，2010，330（1/2）：423-433.

[13] LIU X J，ZHANG S T，JIANG Q P，et al.Using community analysis to explore bacterial indicators for disease suppression of tobacco bacterial wilt[J].Scientific Reports，2016，6：36773.

[14] 郝海婷，王若愚，趙霞，等.基于高通量測序技術(shù)的堆肥對蘭州百合根際微生物多樣性的影響[J].西北農(nóng)業(yè)學(xué)報，2017，26（3）：437-447.

[15] WU H S，YANG X N，F(xiàn)AN J N，et al.Suppression of Fusarium wilt of watermelon by a bio-organic fertilizer containing combinations of antagonistic microorganisms[J].BioControl，2009，54（2）：287-300.

[16] QIU M H，ZHANG R F，XUE C，et al.Application of bio-organic fertilizer can control Fusarium wilt of cucumber plants by regulating microbial community of rhizosphere soil[J].Biology and Fertility of Soils，2012，48（7）：807-816.

[17] ZHAO Q Y，DONG C X，YANG X M，et al.Biocontrol of Fusarium wilt disease for Cucumis melo melon using bio-organic fertilizer[J].Applied Soil Ecology，2011，47（1）：67-75.

[18] EDWARDS J，JOHNSON C，SANTOS-MEDELL?NA C，et al.Structure，variation，and assembly of the root-associated microbiomes of rice[J].Proceedings of the National Academy of Sciences of the United States of America，2015，112（8）：911-920.

[19] KINKEL L L，BAKKER M G，SCHLATTER D C.A coevolutionary framework for managing disease-suppressive soils[J].Annual Review of Phytopathology，2011，49：47-67.

[20] XIONG W，ZHAO Q Y，ZHAO J，et al.Different continuous cropping spans significantly affect microbial community membership and structure in a vanilla-grown soil as revealed by deep pyrosequencing[J].Microbial Ecology，2015，70（1）：209-218.

[21] LOZUPONE C，KNIGHT R.UniFrac：a new phylogenetic method for comparing microbial communities[J].Applied and Environmental Microbiology，2005，71（12）：8228-8235.

[22] ZHOU J H，WU L Y，DENG Y，et al.Reproducibility and quantitation of amplicon sequencing-based detection[J].ISME Journal，2011，5（8）：1303-1313.

[23] CHEN L，XUN W B，SUN L，et al.Effect of different long-term fertilization regimes on the viral community in an agricultural soil of Southern China[J].European Journal of Soil Biology，2014，62：121-126.

[24] SCHLOSS P D，GEVERS D，WESTCOTT S L.Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies[J].PloS One，2011，6（12）：e27310.

[25] JANVIER C，VILLENEUVE F，ALABOUVETTE C，et al.Soil health through soil disease suppression：which strategy from descriptors to indicators[J].Soil Biology and Biochemistry，2007，39（1）：1-23.

[26] FANG S，LIU D，TIAN Y，et al.Tree species composition influences enzyme activities and microbial biomass in the rhizosphere：a rhizobox approach[J].PloS One，2013，8（4）：e61461.

[27] CAPORASO J G，LAUBER C L，WALTERS W A，et al.Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample[J].Proceedings of the National Academy of Sciences of the United States of America，2011，108（S1）：4516-4522.

[28] 宋順華，吳萍，邢寶田，等.多粘類芽胞桿菌WY110對西瓜枯萎病的控制作用[J].植物保護學(xué)報，2011，38（6）：571-572.

[29] IDRIS H A，LABUSCHAGNE N，KORSTEN L.Screening rhizobacteria for biological control of Fusarium root and crown rot of sorghum in Ethiopia[J].Biological Control，2007，40（1）：97-106.

[30] RUMBERGER A，MERWIN I A，THIES J E.Microbial community development in the rhizosphere of apple trees at a replant disease site[J].Soil Biology and Biochemistry，2007，39（7）：1645-1654.

[31] LUO J，RAN W，HU J，et al.Application of bio-organic fertilizer significantly affected fungal diversity of soils[J].Soil Science Society of America Journal，2010，74（6）：2039-2048.

[32] SHEN Z Z，RUAN Y Z，CHAO X，et al.Rhizosphere microbial community manipulated by 2 years of consecutive biofertilizer application associated with banana Fusarium wilt disease suppression[J].Biology and Fertility of Soils，2015，51（5）：553-562.

[33] YANG J，RUEGGER P M，MCKENRY M V，et al.Correlations between root-associated microorganisms and peach replant disease symptoms in a California soil[J].PLoS One，2012，7（10）：e46420.

[34] ZU C，LI Z G，YANG J F，et al.Acid soil is associated with reduced yield，root growth and nutrient uptake in black pepper （Piper nigrum L.）[J].Agricultural Sciences，2014，5（5）：466-473.