邱潔 侯怡鈴 徐麗麗 錢葉 丁祥

摘要：【目的】探究不同品種桑樹(shù)根際土壤細(xì)菌群落結(jié)構(gòu)及其分布規(guī)律，為促進(jìn)桑樹(shù)的生長(zhǎng)發(fā)育、提高其經(jīng)濟(jì)利用價(jià)值和培育更高品質(zhì)的桑樹(shù)品種提供參考?！痉椒ā繉?duì)3個(gè)品種桑樹(shù)（粵椹大10、嘉陵30號(hào)和紅果1號(hào)）根際土壤細(xì)菌的16S rDNA序列V3～V5高變區(qū)進(jìn)行PCR擴(kuò)增，并對(duì)PCR擴(kuò)增產(chǎn)物進(jìn)行高通量測(cè)序，分析土壤中細(xì)菌群落多樣性及分布規(guī)律?！窘Y(jié)果】3個(gè)土壤樣品中共檢測(cè)出細(xì)菌26門76綱88目149科370屬。不同品種桑樹(shù)根際土壤樣品中的細(xì)菌群落組成和結(jié)構(gòu)存在一定差異，豐富度方面表現(xiàn)為嘉陵30號(hào)>粵椹大10>紅果1號(hào);多樣性方面表現(xiàn)為紅果1號(hào)>嘉陵30號(hào)>粵椹大10。3個(gè)品種桑樹(shù)的根際土壤細(xì)菌在門綱目科屬水平上的優(yōu)勢(shì)菌群及所占比例分別為變形菌門（Proteobacteria，37.2%）、α-變形菌綱（Alphaproteobacteria，24.3%）、擬桿菌目（Bacteroidales，24.3%）、擬桿菌科（Bacteroidaceae，14.2%）、擬桿菌屬（Bacteroides，14.2%）;隨著分類的細(xì)化，不同桑樹(shù)品種對(duì)土壤細(xì)菌群落組成和分布的影響越大?！窘Y(jié)論】變形菌門在3種桑樹(shù)根際土壤中均是最優(yōu)勢(shì)細(xì)菌類群，嘉陵30號(hào)桑樹(shù)品種根際土壤細(xì)菌種類最多，紅果1號(hào)桑樹(shù)品種土壤細(xì)菌分布最均勻。根際土壤細(xì)菌群落種類和均勻度可作為桑樹(shù)差異規(guī)?；N植的指標(biāo)。

關(guān)鍵詞： 桑樹(shù);根際土壤;細(xì)菌;多樣性;高通量測(cè)序

中圖分類號(hào)： S888.5? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?   文獻(xiàn)標(biāo)志碼： A 文章編號(hào)：2095-1191（2019）03-0585-08

0 引言

【研究意義】土壤細(xì)菌是土壤生態(tài)系統(tǒng)的重要組成成分，參與養(yǎng)分元素循環(huán)、有機(jī)質(zhì)分解及能量轉(zhuǎn)換，在保持生態(tài)系統(tǒng)的功能、穩(wěn)定性及生產(chǎn)力等方面均發(fā)揮關(guān)鍵作用，是衡量土壤生產(chǎn)力和質(zhì)量的重要指標(biāo)（王衛(wèi)霞等，2013;馮丹妮等，2014;張文元等，2015;尹玉玲等，2017）。在土地覆被過(guò)程中，土壤特性、氣候因子和植被群落等均會(huì)影響土壤的細(xì)菌群落（Li et al.，2014;Tkacz and Poole，2015）。土壤有機(jī)質(zhì)含量與細(xì)菌群落組成和結(jié)構(gòu)間一般具有很好的線性關(guān)系（王傳杰等，2018）。土壤細(xì)菌不僅推動(dòng)了土壤的物質(zhì)循環(huán)和生化反應(yīng)（劉銀銀等，2013），土壤細(xì)菌群落的變化還會(huì)影響土壤養(yǎng)分的轉(zhuǎn)化和吸收，從而影響植物生長(zhǎng)代謝（孫波等，2017）。不同植物類型根系分泌物和凋落物的質(zhì)量又反作用于土壤細(xì)菌群落的生長(zhǎng)發(fā)育，造成某些細(xì)菌群落的出現(xiàn)或消亡（陳心想等，2014;潘彥碩等，2018;王美溪等，2018）。桑（Morus alba L.）為?？粕僦参?，其葉是家蠶重要的飼料，是蠶絲產(chǎn)業(yè)發(fā)展的必須物質(zhì)基礎(chǔ)（Jia et al.，2014;向仲懷等，2017），因此，了解桑樹(shù)根標(biāo)土壤細(xì)菌組成多樣性對(duì)保障桑葉生產(chǎn)和促進(jìn)蠶絲產(chǎn)業(yè)發(fā)展均具有重要意義?！厩叭搜芯窟M(jìn)展】Carvalhais等（2013）通過(guò)擬南芥在無(wú)菌土壤和非無(wú)菌土壤中生長(zhǎng)對(duì)比試驗(yàn)，證明擬南芥在整個(gè)細(xì)菌群落存在條件下比細(xì)菌貧乏時(shí)表現(xiàn)出更強(qiáng)的生長(zhǎng)力。Muehe等（2015）通過(guò)對(duì)鼠耳芥根際細(xì)菌組成與植物對(duì)鎘和鋅的吸收研究，證明土壤細(xì)菌可直接或間接影響植物對(duì)金屬元素的吸收。崔佩佩等（2018）采用Biolog微平板法研究不同施肥條件下高粱根際土壤細(xì)菌功能多樣性的變化，結(jié)果證明施肥能改善土壤細(xì)菌的功能多樣性。蔣景龍等（2018）采用Illumina MiSeq高通量測(cè)序技術(shù)對(duì)西洋參根腐病病株和健株根際土壤細(xì)菌的組成及多樣性變化進(jìn)行研究，結(jié)果顯示鞘脂菌屬（Sphingobium）和紅游動(dòng)菌屬（Rhodoplanes）可能是引起西洋參根腐病的重要菌群。對(duì)于桑樹(shù)根際土壤細(xì)菌的研究，目前主要采用傳統(tǒng)方法探究不同施肥方案、作物間作等對(duì)根際土壤細(xì)菌種類變化、重金屬含量、酶活性差異的影響（楊統(tǒng)一等，2015;張萌萌等，2015）。鄧文等（2016）采用Biolog微生態(tài)技術(shù)研究偏施氮肥桑樹(shù)根際土壤細(xì)菌對(duì)碳源利用的能力變化，證明4年偏施氮肥桑樹(shù)根際土壤細(xì)菌多樣性及活性均高于其他方法處理的土壤。樊芳玲等（2016）利用磷脂脂肪酸（PLFA）法分析不同模式的坡耕地—桑樹(shù)系統(tǒng)對(duì)旱坡地紫色土壤中細(xì)菌群落的影響，證明種植桑樹(shù)能顯著提高土壤中PLFA含量，改善細(xì)菌群落結(jié)構(gòu)?！颈狙芯壳腥朦c(diǎn)】高通量測(cè)序技術(shù)能產(chǎn)生測(cè)序覆蓋深度更高的基因數(shù)據(jù)，在細(xì)菌群落組成研究中能真實(shí)地揭示環(huán)境中細(xì)菌群落的多樣性和復(fù)雜性。目前針對(duì)不同品種桑樹(shù)根際土壤細(xì)菌多樣性的高通量測(cè)序分析研究鮮見(jiàn)報(bào)道?！緮M解決的關(guān)鍵問(wèn)題】對(duì)不同品種桑樹(shù)根際土壤細(xì)菌的16S rDNA序列V3～V5高變區(qū)進(jìn)行PCR擴(kuò)增，并對(duì)PCR擴(kuò)增產(chǎn)物進(jìn)行高通量測(cè)序，分析土壤中細(xì)菌群落多樣性及分布規(guī)律，了解不同品種桑樹(shù)根際土壤細(xì)菌的群落結(jié)構(gòu)、物種組成和差異，探究根際土壤細(xì)菌與桑樹(shù)的相互關(guān)系，為促進(jìn)桑樹(shù)的生長(zhǎng)發(fā)育、提高其經(jīng)濟(jì)利用價(jià)值及培育出更高品質(zhì)的桑樹(shù)品種提供參考。

1 材料與方法

1. 1 試驗(yàn)材料

土壤樣品采自四川省農(nóng)業(yè)科學(xué)院蠶業(yè)研究所瀠溪桑園，分別以粵椹大10（廣東省農(nóng)業(yè)科學(xué)院蠶業(yè)與農(nóng)產(chǎn)品加工研究所選育）、嘉陵30號(hào)和紅果1號(hào)（四川省農(nóng)業(yè)科學(xué)院蠶業(yè)研究所選育）3種果桑根際土壤作為采樣區(qū)，隨機(jī)在每個(gè)采樣區(qū)設(shè)4個(gè)采樣點(diǎn)，每個(gè)采樣點(diǎn)采集10份土壤樣品（0～10 cm），使用S形采樣法，剔除石塊和植物根系，混合均勻并過(guò)2 mm篩，以四分法篩取所得土樣為同一樣品，再裝入已滅菌的采樣袋，分別編號(hào)為L(zhǎng)CF20、LCF21和LCF22，于冰盒中保存。

1. 2 總DNA提取、16S rDNA序列擴(kuò)增及高通量測(cè)序

根據(jù)Fast DNA SPIN Kit for Soil（MP Biomedical，USA）的使用說(shuō)明，加入溶菌酶（上海源葉生物科技有限公司）60 mg/mL輔助提取土壤細(xì)菌的總DNA。用1%瓊脂糖凝膠對(duì)總DNA進(jìn)行電泳檢測(cè)，電泳條件為120 V、45 min，總DNA條帶單一清晰，無(wú)嚴(yán)重拖尾現(xiàn)象，可用于后續(xù)體外擴(kuò)增。用引物515F/907R擴(kuò)增細(xì)菌16S rDNA序列的V3～V5區(qū)域。反應(yīng)體系25.0 μL：Dream Taq PCR Master Mix（2×）12.5 μL，DNA模板1.0 μL，515F/907R引物各0.5 μL，去離子水補(bǔ)足至25.0 μL。擴(kuò)增程序：94 ℃預(yù)變性5 min;94 ℃ 30 s，56 ℃ 30 s，72 ℃ 1 min，進(jìn)行10個(gè)循環(huán)，每個(gè)循環(huán)降低0.5 ℃;94 ℃ 30 s，56 ℃ 30 s，72 ℃ 1 min，進(jìn)行25個(gè)循環(huán);72 ℃延伸10 min;4 ℃保存。擴(kuò)增產(chǎn)物用1%瓊脂糖凝膠進(jìn)行電泳檢測(cè)，將單一清晰沒(méi)有拖尾現(xiàn)象條帶對(duì)應(yīng)的樣本送至北京金諾銳杰基因科技有限公司測(cè)序，并采用MiSeq測(cè)序儀完成測(cè)序工作（Lu et al.，2015;Navarrete et al.，2015;Lossius et al.，2016）。

1. 3 高通量測(cè)序數(shù)據(jù)處理

原始數(shù)據(jù)經(jīng)過(guò)低質(zhì)量過(guò)濾、去接頭等后得到高質(zhì)量的基因序列Reads，然后合并正反向Reads，獲得16S rDNA序列V3～V5區(qū)擴(kuò)增子。利用Flash 8.0合并正方向Reads，合并后序列輸出，并使用FastQC進(jìn)行質(zhì)控分析。采用Usearch61檢查嵌合體，De novo方式去除嵌合體序列，再對(duì)小片段文庫(kù)進(jìn)行組裝，獲得該物種的基因組序列圖譜。對(duì)所有分離的高質(zhì)量樣本序列進(jìn)行聚類分析，即一個(gè)OTU（同源性為97%），再劃分操作分類單元（OTUs），并構(gòu)建稀釋曲線及進(jìn)行聚類分析。

Chao1指數(shù)計(jì)算公式：Schao1=Sobs+n1（n1-1）/2（n2+1），其中，Schao1為估計(jì)的OUT數(shù)，Sobs為觀測(cè)到的OTU數(shù)，n1為只有一條序列的OUT數(shù)目，n2為只有兩條序列的OUT數(shù)目。

Shannon指數(shù)計(jì)算公式：H'=-∑（Pi）（lnPi），其中，Pi為此細(xì)菌個(gè)體數(shù)占總個(gè)體數(shù)比例。

相對(duì)豐度計(jì)算公式：D=（S-1）/lnN，其中，S為群落中的特定細(xì)菌總數(shù)目，N為群落中所有細(xì)菌的總數(shù)。

2 結(jié)果與分析

2. 1 不同品種桑樹(shù)根際土壤細(xì)菌多樣性分析結(jié)果

用Chao1指數(shù)表示細(xì)菌群落的豐富度，其值越高反映細(xì)菌群落物種豐富度越高，細(xì)菌種類越多。顯示3個(gè)樣品的細(xì)菌群落豐富度均隨測(cè)序深度的增加而增加，隨后逐漸趨于平緩，此時(shí)嘉陵30號(hào)的細(xì)菌群落豐富度最高，粵椹大10次之，紅果1號(hào)最低。用Shannon指數(shù)表示樣品的細(xì)菌多樣性程度，其值越高反映細(xì)菌群落的物種多樣性越高，分布也越均勻。當(dāng)3個(gè)樣品的Shannon指數(shù)曲線最終趨于平緩時(shí)，紅果1號(hào)的細(xì)菌多樣性程度最高，粵椹大10最低。

2. 2 在門水平上的細(xì)菌群落豐度分析結(jié)果

從門水平來(lái)看，細(xì)菌種類最多的品種是嘉陵30號(hào)（26門），最少的是粵椹大10（24門）。變形菌門（Proteobacteria）、擬桿菌門（Bacteroidetes）和厚壁菌門（Firmicutes）的相對(duì)豐度較高，分別為37.2%、25.4%和14.1%，在3個(gè)土壤混合樣品中3種細(xì)菌門的相對(duì)豐度之和占土壤細(xì)菌總量的76.7%。將群落相對(duì)豐度水平位于前10的細(xì)菌（相對(duì)豐度>0.4%）按照樣品和分類進(jìn)行聚類后繪制熱圖，結(jié)果顯示，3個(gè)樣品中的細(xì)菌大致分為兩簇，梭桿菌門（Fusobacteria）、芽單胞菌門（Gemmatimona-detes）和厚壁菌門等親緣關(guān)系較近的聚為一簇，硝化螺旋菌門（Nitrospirae）、綠彎菌門（Chloroflexi）和酸桿菌門（Acidobacteria）等聚為另一簇。不同品種桑樹(shù)根際土壤中細(xì)菌的相對(duì)豐度有一定差異，其中，酸桿菌門在紅果1號(hào)和嘉陵30號(hào)中的相對(duì)豐度較高，在粵椹大10中的相對(duì)豐度較低;硝化螺旋菌門和綠彎菌門在紅果1號(hào)中的相對(duì)豐度較高，在嘉陵30號(hào)和粵椹大10中的相對(duì)豐度較低。

2. 3 在綱水平上的細(xì)菌群落豐度分析結(jié)果

3個(gè)土壤樣品中的細(xì)菌在綱水平上群落組成最多的是粵椹大10和嘉陵30號(hào)（76綱），最少的是紅果1號(hào)（74綱）。α-變形菌綱（Alphaproteobacteria）、β-變形菌綱（Betaproteobacteria）、γ-變形菌綱（Gammaproteobacteria）、擬桿菌綱（Bacteroidia）和梭狀芽胞桿菌綱（Clostridia）等5個(gè)細(xì)菌綱的相對(duì)豐度較高，分別為24.3%、16.0%、12.3%、11.2%和7.7%，在3個(gè)土壤混合樣品中5個(gè)細(xì)菌綱的相對(duì)豐度之和占土壤細(xì)菌總豐度的71.5%。將群落相對(duì)豐度水平位于前20的細(xì)菌（相對(duì)豐度>0.4%）按照樣品和分類進(jìn)行聚類后繪制熱圖，結(jié)果顯示，嘉陵30號(hào)與粵椹大10細(xì)菌種類的組成相比紅果1號(hào)而言更相似;3個(gè)樣品中的細(xì)菌大致分為兩簇，α-變形菌綱和擬桿菌綱等7個(gè)綱親緣關(guān)系較的近聚為一簇，γ-變形菌綱和β-變形菌綱等13個(gè)綱聚為另一簇。

2. 4 在目水平上的細(xì)菌群落豐度分析結(jié)果

3個(gè)土壤樣品中的細(xì)菌在目水平上的群落組成分別是紅果1號(hào)88目、嘉陵30號(hào)71目和粵椹大10 78目。擬桿菌目（Bacteroidales）、腸桿菌目（Enterobacteriales）、梭菌目（Clostridiales）、酸桿菌目（Acidobacteriales）和伯克氏菌目（Burkholderiales）等5個(gè)細(xì)菌目的相對(duì)豐度較高，分別為24.3%、12.3%、10.6%、4.1%和3.8%，在3個(gè)土壤混合樣品中5個(gè)細(xì)菌目的相對(duì)豐度之和占土壤細(xì)菌總豐度的55.1%。將相對(duì)豐度位于前30的細(xì)菌（相對(duì)豐度>0.4%）按照樣品和分類進(jìn)行聚類后繪制熱圖，3個(gè)樣品中的細(xì)菌大致分為兩簇，擬桿菌目和梭菌目等14個(gè)目親緣關(guān)系較近的聚為一簇，酸桿菌目和伯克氏菌目等16個(gè)目聚為另一簇。在目水平上，3個(gè)樣品中的細(xì)菌組成存在差異，如黃單胞桿菌目（Xanthomonadales）和紅螺菌目（Rhodospirales）在紅果1號(hào)中的相對(duì)豐度較高，在嘉陵30號(hào)和粵椹大10中的相對(duì)豐度較低。

2. 5 在科水平上的細(xì)菌群落豐度分析結(jié)果

3個(gè)土壤樣品中的細(xì)菌在科水平上的群落組成最多的是紅果1號(hào)（149科），嘉陵30號(hào)（116科）和粵椹大10（118科）相對(duì)較少。擬桿菌科（Bacteroidaceae，14.2%）、腸桿菌科（Enterobacteriaceae，10.6%）、普雷沃氏菌科（Prevotellaceae，8.3%）、瘤胃菌科（Ruminococcaceae，5.6%）、酸桿菌科（Acidobacteriaceae，3.2%）、毛螺菌科（Lachnospiraceae，2.6%）和韋榮氏菌科（Veillonellaceae，3.2%）等7個(gè)細(xì)菌的相對(duì)豐度較高，在3個(gè)土壤混合樣品中7個(gè)細(xì)菌科的相對(duì)豐度之和占土壤細(xì)菌總豐度的47.7%。將相對(duì)豐度位于前40的細(xì)菌（相對(duì)豐度>0.2%）按照樣品和分類進(jìn)行聚類后繪制熱圖，3個(gè)樣品中的細(xì)菌大致分為兩簇，梭桿菌科（Fusobacteriaceae）和黃單胞菌科（Xanthomonadaceae）等26個(gè)科聚為一簇，瘤胃菌科和韋榮氏菌科等14個(gè)科聚為另一簇。在科水平上，隨著分類的細(xì)化，桑樹(shù)品種對(duì)細(xì)菌群落組成的影響越大。如鞘酯菌科（Sphingomonadaceae）在嘉陵30號(hào)中的相對(duì)豐度較高，在粵椹大10和紅果1號(hào)中的相對(duì)豐度較低。

2. 6 在屬水平上的細(xì)菌群落豐度分析結(jié)果

3個(gè)土壤樣品中共有細(xì)菌屬370個(gè)，其中單個(gè)樣品最多的是嘉陵30號(hào)（370屬），紅果1號(hào)（354屬）和粵椹大10（351屬）相對(duì)較少。擬桿菌屬（Bacteroides，14.2%）、普氏菌屬（Prevotella，8.3%）、棲糞桿菌屬（Faecalibacterium，3.1%）、奈瑟氏球菌屬（Neisseria，2.1%）、鏈球菌屬（Streptococcus，1.2%）、薩特氏菌屬（Sutterella，1.2%）、考拉桿菌屬（Phascolarctobacterium，1.1%）、不動(dòng)桿菌屬（Acinetobacter，1.1%）和巨單胞菌屬（Megamonas，1.0%）等9個(gè)細(xì)菌屬相對(duì)豐度較高，在3個(gè)土壤混合樣品中9個(gè)細(xì)菌屬的相對(duì)豐度之和占土壤細(xì)菌總豐度的33.3%。將相對(duì)豐度位于前50的細(xì)菌（相對(duì)豐度>0.1%）分別按照樣品和分類進(jìn)行聚類后繪制熱圖，3個(gè)樣品中的細(xì)菌大致分為兩簇，羅氏菌屬（Roseburla）和薩特氏菌屬等37個(gè)屬聚為一簇，硝化螺菌屬（Nitrospira）和黃桿菌屬（Flavobacterium）等13個(gè)屬聚為另一簇。各樣品屬級(jí)分類聚類分析結(jié)果顯示，在屬的水平上，莫拉氏菌（Kaistobacter）在嘉陵30號(hào)中的相對(duì)豐度較高，在粵椹大10和紅果1號(hào)中則相對(duì)較低;艾克曼菌屬（Akkermansia）在紅果1號(hào)中未檢測(cè)出;分枝桿菌屬（Ramlibacter）在粵椹大10中未檢測(cè)出。

3 討論

土壤作為細(xì)菌的天然培養(yǎng)基，含有大量且豐富的細(xì)菌種類。相比于傳統(tǒng)測(cè)序技術(shù)，使用高通量測(cè)序技術(shù)對(duì)土壤細(xì)菌多樣性的研究鮮見(jiàn)報(bào)道，且高通量測(cè)序產(chǎn)出的數(shù)據(jù)量是傳統(tǒng)測(cè)序技術(shù)無(wú)法比擬（明磊等，2018）。本研究采用高通量測(cè)序技術(shù)對(duì)3個(gè)不同品種的桑樹(shù)根際土壤細(xì)菌組成和結(jié)構(gòu)進(jìn)行研究，分析了3個(gè)不同品種桑樹(shù)根際土壤細(xì)菌在門、綱、目、科、屬不同分類水平上的優(yōu)勢(shì)類群，共檢測(cè)出包括擬桿菌門、厚壁菌門和變形菌門等優(yōu)勢(shì)菌門在內(nèi)的26門，其中變形菌門是最豐富的細(xì)菌類群，與Constancias等（2015）、Docherty等（2015）的研究結(jié)果一致。但Heiko等（2016）發(fā)現(xiàn)酸桿菌門是歐洲山毛櫸和挪威云杉根際土壤中含量最豐富的細(xì)菌類群;而Kristin等（2016）在溫帶草原和森林土壤中發(fā)現(xiàn)放線菌門（Actinobacteria）是占主導(dǎo)地位的細(xì)菌類群，說(shuō)明不同植物根際土壤中存在不同的優(yōu)勢(shì)細(xì)菌類群。

本研究共測(cè)出包括α-變形菌綱和β-變形菌綱等優(yōu)勢(shì)菌綱在內(nèi)的76綱;擬桿菌目和腸桿菌目等優(yōu)勢(shì)菌目在內(nèi)的88目，擬桿菌科和腸桿菌科等優(yōu)勢(shì)菌科在內(nèi)的149科，擬桿菌屬和普氏菌屬等優(yōu)勢(shì)菌屬在內(nèi)的370屬。3種不同品種桑樹(shù)根際間土壤細(xì)菌在門綱目科屬水平上的優(yōu)勢(shì)菌群及所占總比例分別為變形菌門（37.2%）、α-變形菌綱（24.3%）、擬桿菌目（24.3%）、擬桿菌科（14.2%）、擬桿菌屬（14.2%）。但隨著分類的細(xì)化，不同的桑樹(shù)品種對(duì)細(xì)菌群落組成和分布的影響越大，在目科屬水平上差異更明顯。如黃單胞桿菌目和紅螺菌目在紅果1號(hào)中的相對(duì)豐度較高，鞘酯菌科在嘉陵30號(hào)中的相對(duì)豐度較高，而分枝桿菌屬在粵椹大10中未檢測(cè)出。土壤細(xì)菌群落多樣性不僅與植物種類有關(guān)，還與土壤落葉層、腐殖質(zhì)及理化性質(zhì)有關(guān);同時(shí)，桑樹(shù)根系的生理活動(dòng)能改善土壤理化性質(zhì)，使之適應(yīng)土壤中細(xì)菌的生長(zhǎng)（Tkacz et al.，2015）。實(shí)際應(yīng)用中可根據(jù)桑樹(shù)品種的不同，選擇合適的栽培土壤，改善土壤理化性質(zhì)并施用高效細(xì)菌肥料，以促進(jìn)桑樹(shù)根際間土壤優(yōu)勢(shì)細(xì)菌群落的形成，從而促進(jìn)桑樹(shù)生長(zhǎng)。

4 結(jié)論

變形菌門在供試3種桑樹(shù)品種根際土壤中均是最優(yōu)勢(shì)細(xì)菌類群，嘉陵30號(hào)桑樹(shù)品種根際土壤細(xì)菌種類最多，紅果1號(hào)桑樹(shù)品種土壤細(xì)菌分布最均勻。根際土壤細(xì)菌群落種類和均勻度可作為桑樹(shù)差異規(guī)模化種植的指標(biāo)。

參考文獻(xiàn)：

陳心想，耿增超，王森，趙宏飛. 2014. 施用生物炭后塿土土壤微生物及酶活性變化特征[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，33（4）： 751-758. [Chen X X，Geng Z C，Wang S，Zhao H F. 2014. Effects of biochar amendment on microbial biomass and enzyme activities in loess soil[J]. Journal of Agro-Environmental Science，33（4）： 751-758.]

崔佩佩，武愛(ài)蓮，王勁松，董二偉，南江寬，白文斌，焦曉燕. 2018. 不同施肥處理對(duì)高粱根際土壤微生物功能多樣性的影響[J]. 華北農(nóng)學(xué)報(bào)，33（5）： 195-202. [Cui P P，Wu A L，Wang J S，Dong E W，Nan J K，Bai W B，Jiao X Y. 2018. Effect of different fertilization treatments on soil microbial function diversity in rhizosphere of sorghum[J]. Acta Agriculturae Boreali-Sinica，33（5）： 195-202.]

鄧文，熊雙偉，于翠，李勇，莫榮利，熊超，胡興明. 2016. 長(zhǎng)期偏施氮肥對(duì)桑樹(shù)根際土壤微生物碳源利用能力的影響[J]. 湖北農(nóng)業(yè)科學(xué)，55（23）： 6108-6112. [Deng W，Xiong S W，Yu C，Li Y，Mo R L，Xiong C，Hu X M. 2016. E-ffects of long-term partial application of nitrogen fertilizer on soil carbon source utilization ability in mulberry rhizosphere[J]. Hubei Agricultural Sciences，55（23）： 6108-6112.]

樊芳玲，謝德體，魏朝富，倪九派，唐振亞. 2016. 坡耕地—桑樹(shù)系統(tǒng)土壤微生物群落結(jié)構(gòu)的PLFA分析[J]. 生態(tài)學(xué)報(bào)，36（12）： 3812-3823. [Fan F L，Xie D T，Wei C F，Ni J P，Tang Z Y. 2016. PLFA analysis showed that the confi-guration modes of sloping crop-mulberry systems affect the soil microbial community structure[J]. Acta Ecologica Sinica，36（12）： 3812-3823.]

馮丹妮，伍鈞，楊剛，張璘瑋，張乙涵，王靜雯. 2014. 連續(xù)定位施用沼液對(duì)水旱輪作耕層土壤微生物區(qū)系及酶活性的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，33（8）： 1644-1651. [Feng D N， Wu J， Yang G，Zhang L W，Zhang Y H，Wang J W. 2014. Influence of long-term applications of biogas slurry on microbial community composition and enzymatic activities in surface soil under rice-rape rotation[J]. Journal of Agro-Environmental Science，33（8）：1644-1651.]

蔣景龍，余妙，李麗，焦成瑾，徐皓，任緒明，楊玲娟. 2018. 西洋參根腐病發(fā)生與根際土壤細(xì)菌群落結(jié)構(gòu)變化關(guān)系研究[J]. 中草藥，49（18）： 4399-4407. [Jiang J L，Yu M，Li L，Jiao C J，Xu H，Ren X M，Yang L J. 2018. Relationship between occurrence of root-rot and changes of bacterial community structure in rhizosphere soil of Panax quinquefolius[J]. Chinese Traditional and Herbal Drugs，49（18）： 4399-4407.]

劉銀銀，李峰，孫慶業(yè)，謝永宏. 2013. 濕地生態(tài)系統(tǒng)土壤微生物研究進(jìn)展[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào)，19（3）：547-552. [Liu Y Y，Li F，Sun Q Y，Xie Y H. 2013. Review on the study of soil microorganisms in wetland ecosystems[J]. Chinese Journal of Applied and Environmental Biology，19（3）： 547-552.]

明磊，戴習(xí)林，江宗冰，丁福江. 2018. 高通量測(cè)序在羅氏沼蝦抗病力比較中的應(yīng)用[J]. 河南農(nóng)業(yè)科學(xué)，47（8）：134-141. [Ming L，Dai X L，Jiang Z B，Ding F J. 2018. Application of high-throughput sequencing in the comparison of disease resistance of Macrobrachium rosenbergii[J]. Journal of Henan Agricultural Sciences，47（8）：134-141.]

潘彥碩，朱清禾，李聰，王幽靜，吳坤，李烜楨，張世敏，吳宇澄. 2018. 纖維素、秸稈和木屑對(duì)農(nóng)田土壤硝化作用及微生物的影響[J]. 河南農(nóng)業(yè)大學(xué)學(xué)報(bào)，52（5）：785-792.[Pan Y S，Zhu Q H，Li C，Wang Y J，Wu K，Li X Z，Zhang S M，Wu Y C. 2018. Effects of cellulose，straw and sawdust on soil nitrification and microorganisms in farmland[J]. Journal of Henan Agricultural University，52（5）：785-792.]

孫波，王曉玥，呂新華. 2017. 我國(guó)60年來(lái)土壤養(yǎng)分循環(huán)微生物機(jī)制的研究歷程——基于文獻(xiàn)計(jì)量學(xué)和大數(shù)據(jù)可視化分析[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào)， 23（6）： 1590-1601. [Sun B， Wang X Y， Lü X H. 2017. The historical venation in research on microbial mechanisms of soil nu-trient cycling in the past 60 years—Based on bibliometric analysis and big data visualization[J]. Journal of Plant Nutrition and Fertilizer，23（6）： 1590-1601.]

王傳杰，王齊齊，徐虎，高洪軍，朱平，徐明崗，張文菊. 2018. 長(zhǎng)期施肥下農(nóng)田土壤—有機(jī)質(zhì)—微生物的碳氮磷化學(xué)計(jì)量學(xué)特征[J]. 生態(tài)學(xué)報(bào)，38（11）： 3848-3858. [Wang C J，Wang Q Q，Xu H，Gao H J，Zhu P，Xu M G，Zhang W J. 2018. Carbon，nitrogen，and phosphorus stoichiometry characteristics of bulk soil，organic matter，and soil microbial biomass under longterm fertilization in cropland[J].Acta Ecologica Sinica，38（11）： 3848-3858.]

王美溪，劉珂藝，邢亞娟. 2018. 氣候變化背景下土壤微生物與植物物種多樣性關(guān)聯(lián)分析[J]. 中國(guó)農(nóng)學(xué)通報(bào)，34（20）：111-117. [Wang M X，Liu K Y，Xing Y J. 2018. Association analysis of soil microorganism and plant species diversity under climate change[J]. Chinese Agricultural Science Bulletin，34（20）：111-117.]

王衛(wèi)霞，史作民，羅達(dá)，劉世榮，盧立華. 2013. 南亞熱帶3種人工林土壤微生物生物量和微生物群落結(jié)構(gòu)特征[J]. 應(yīng)用生態(tài)學(xué)報(bào)，24（7）： 1784-1792. [Wang W X，Shi Z M，Luo D，Liu S R，Lu L H. 2013. Characteristics of soil microbial biomass and community composition in three types of plantations in southern subtropical area of China[J]. Chinese Journal of Applied Ecology，24（7）： 1784-1792.]

向仲懷，何寧佳，黃先智. 2017. 桑與畜牧業(yè)[J]. 草業(yè)學(xué)報(bào)，26（2）：1-9. [Xiang Z H，He N J，Huang X Z. 2017. Mulberry and animal husbandry[J]. Acta Pratacultration Sinica，26（2）： 1-9.]

楊統(tǒng)一，曲靜怡，趙衛(wèi)國(guó)，李龍. 2015. 不同抗病性桑樹(shù)品種根際土壤微生態(tài)環(huán)境特征分析[J]. 蠶業(yè)科學(xué)，41（2）： 211-217. [Yang T Y，Qu J Y，Zhao W G，Li L. 2015. Analysis on micro-ecological environmental characterisitics in rhizosphere soil of different resistant mulberry varieties[J]. Science of Sericulture，41（2）： 211-217.]

尹玉玲，湯泳萍，謝啟鑫，羅紹春，趙萍. 2017. 豆蔻酸對(duì)茄子根際土壤微生物生理類群和土壤酶活性的影響[J]. 江蘇農(nóng)業(yè)學(xué)報(bào)，33（1）：181-184. [Yin Y L，Tang Y P，Xie Q X，Luo S C，Zhao P. 2017. Effect of myristic acid on soil microbial population and enzyme activities in rhizosphere of the eggplants[J]. Jiangsu Journal of Agricultu-ral Sciences，33（1）：181-184.]

張萌萌，敖紅，李鑫，張景云，王寧，鞠成梅，王佳，蔡敦江，孫廣玉. 2015. 桑樹(shù)/苜蓿間作對(duì)根際土壤酶活性和微生物群落多樣性的影響[J]. 草地學(xué)報(bào)，23（2）： 302-309. [Zhang M M，Ao H，Li X，Zhang J Y，Wang N，Ju C M，Wang J，Cai D J，Sun G Y. 2015. Effects of intercropping between mulberry and alfalfa on soil enzyme activity and microbial community diversity in rhizosphere[J]. Acta Agrestia Sinica，23（2）： 302-309.]

張文元，盛可銀，范承芳，劉順，陳淼淼，鐘石泉，溫衛(wèi)華，吳珍花，涂淑萍，郭曉敏，胡冬南. 2015. 贛南毛竹根際土壤微生物群落結(jié)構(gòu)的PLFA分析[J]. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào)，37（3）： 475-483. [Zhang W Y，Sheng K Y，F(xiàn)an C F，Liu S，Chen M M，Zhong S Q，Wen W H，Wu Z H，Tu S P，Guo X M，Hu D N. 2015. Phospholipid fatty acid analysis of microbial community structure in rhizosphere soil of Phyllostachys edulis in Gannan[J]. Acta Agriculturae Universitatis Jiangxiensis，37（3）： 475-483.]

Carvalhais L C，F(xiàn)rederico M，Tan C H，Jin H C，Schenk P M. 2013. Plant growth in Arabidopsisis assisted by compost soil-derived microbial communities[J]. Front Plant Scien-ce，4（1）： 235-245.

Constancias F，Saby N P，Terrat S，Dequeit S，Horrigue W，Nowak V，Guillemen J P，Luc B D，Nicolas C P，Ranjard L. 2015. Contrasting spatial patterns and ecological attributes of soil bacterial and archaeal taxa across a landscape[J]. Microbiology Open，4（3）： 518-531.

Docherty K M，Borton H M，Espinosa N，Gebhardt M，Juliana G L，Jessica L M，Maes P W，Mott B M，Parnell J J，Purdy G，Rodrigues P A，Stanish L F，Walser O N，Gallery R E. 2015. Key edaphic properties largely explain temporal and geographic variation in soil microbial communities across four biomes[J]. PLoS One，10（11）： e0135352.

Heiko N，Kezia G，Ingo S，Birgit P，Kristin K，Genis A C，Marion S，F(xiàn)ran?ois B，Rolf D，Tesfaye W. 2016. Fine spatial scale variation of soil microbial communities under european beech and norway spruce[J]. Frontiers in Microbiology，7： 2067-2081.

Jia L，Zhang D Y，Qi X W，Ma B，Xiang Z H，He N J. 2014. Identification of the conserved and novel miRNAs in mulberry by high-throughput sequencing[J]. PLoS One，9（8）： e104409.

Kristin K，Bernd W，Vera K，F(xiàn)ranziska W，Heiko N，Ingo S，Marion S，Rolf D. 2016. Driving forces of soil bacterial community structure，diversity，and function in temperate grasslands and forests[J]. Scientific Reports，6： 33696-33702.

Li B Y，Yao Q，Zhu H H. 2014. Approach to analyze the diversity of myxobacteria in soil by semi-nested PCR-denaturing gradient gel electrophoresis（DGGE） based on taxon-specific gene[J]. PLoS One，9（10）： e108877.

Lossius P A，Johansen J，Vartdal F，Trygve H. 2016. High-throughput sequencing of immune repertoires in multiple sclerosis[J]. Annals of Clinical and Translational Neuro-logy，3（4）： 295-306.

Lu X，Zhang X X，Wang Z，Huang K L，Wang Y，Liang W G，Tan Y F，Li B，Tan J Y. 2015. Bacterial pathogens and community composition in advanced sewage treatment systems revealed by metagenomics analysis based on high-throughput sequencing[J]. PLoS One，10（5）： e0125549.

Muehe E M，Weigold P，Adaktylou I J，Britta P F，Kreamer U，Kappler A，Behrens S. 2015. Rhizosphere microbial community composition affects cadmium and zinc uptake by the metal-hyperaccumulating plant Arabidopsis halleri[J]. Applied and Environmental Microbiology，81（6）： 2173-2181.

Navarrete A A，Diniz T R，Braga L P，Silva G G，F(xiàn)ranchini J C，Rossetto R，Edwards R A，Tsai S M. 2015. Multi-analytical approach reveals potential microbial indicators in soil for sugarcane model systems[J]. PLoS One，10（6）： e0129765.

Tkacz A，Cheema J，Chandra G，Grant A，Poole P S. 2015. Stability and succession of the rhizosphere microbiota depends upon plant type and soil composition[J]. The ISME Journal，9（11）： 2349-2359.

Tkacz A，Poole P. 2015. Role of root microbiota in plant productivity[J]. Journal of Experimental Botany，66（8）： 2167-2175.

（責(zé)任編輯 麻小燕）