花生根際土壤細菌群落多樣性對鹽脅迫的響應

戴良香 徐 揚 張冠初 史曉龍 秦斐斐 丁 紅,* 張智猛,*

花生根際土壤細菌群落多樣性對鹽脅迫的響應

戴良香1,**徐 揚1,**張冠初1史曉龍2秦斐斐1丁 紅1,*張智猛1,*

1山東省花生研究所, 山東青島 266100;2沈陽農業大學農學院, 遼寧沈陽 110866

為明確鹽脅迫條件下花生根際土壤細菌群落結構的變化, 采用盆栽試驗設置不同鹽脅迫強度處理, 以開花下針期和成熟期花生根際土壤為研究對象, 提取其總DNA并構建細菌16S rRNA基因文庫, 利用高通量測序技術進行測序, 并對測序結果進行生物信息學分析。結果顯示, 各處理樣本根際土壤優勢菌門均為變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、Patescibacteria、酸桿菌門(Acidobacteria)、綠彎菌門(Chloroflexi)。優勢菌目分別為(Saccharimonadales)、β-變形菌目(Betaproteobacteria)、鞘脂單胞菌目(sphingomonadales)、芽單胞菌目(Gemmatimonadales)和根瘤菌目(Rhizobiales)。鹽脅迫提高了變形菌門(Proteobacteria)的相對豐度, 但降低了放線菌門(Actinobacteria)的相對豐度, 且均隨鹽脅迫強度提高其升降幅度增大。鹽脅迫下基施鈣肥使β-變形菌目(Betaproteobacteria)、芽單胞菌目(Gemmatimonadales)和鞘脂單胞菌目(Sphingomonadales)相對豐度顯著升高, 并均受鹽脅迫強度、生育時期和外源施鈣肥的影響。聚類分析結果表明, 較高濃度鹽脅迫處理樣本的優勢菌屬組成相似而聚為一組, 非鹽脅迫處理樣本屬水平豐度依生育時期相同而相近, 各聚為一組。鹽脅迫強度、生育時期對花生根際土壤微生物菌群類型影響較大, 基施鈣肥影響較小。土壤微生物功能預測分析顯示, 高鹽脅迫可明顯降低次生代謝產物、聚糖的生物合成與代謝, 以及氨基酸和脂肪酸代謝等功能基因在根際土壤中的富集?；ㄉ⑸L期、低鹽脅迫和基施鈣肥處理使得功能基因豐度大幅提高, 可能對花生生長及脅迫應答具有重要功能。根際土壤微生物菌群多樣性的研究將為通過改良土壤微生物環境來提高植物脅迫耐受性提供重要途徑。

花生; 鹽脅迫; 根際; 土壤微生物群落; 16S rRNA基因

土地鹽漬化已經成為世界性的資源環境問題之一, 是制約農業發展的巨大環境壓力, 每年造成數十億農業經濟損失[1]。我國土壤鹽堿化現象日趨嚴重, 鹽堿地面積約9.91×107hm2, 主要分布于氣候比較干旱的山東、河北、河南、新疆等省區, 大部分處于待開發狀態, 是制約該區農業可持續發展的首要因素[2]?；ㄉ俏覈匾挠土献魑锖徒洕魑? 具有耐瘠、固氮、培肥地力和中度耐鹽堿等特性, 可作為鹽堿土區種植業結構調整中較為適宜的替代作物; 同時, 限制鹽堿地花生高產的關鍵是缺苗、早衰和養分失衡。因此, 擴大鹽堿地花生種植面積, 提高產量, 加強花生抗鹽生理機制研究顯得更為迫切。鹽堿土壤性質的改善主要通過減少土壤中鹽離子含量、改善土壤結構、調節土壤酸堿度、增加土壤中可被植物利用營養元素含量等途徑, 提高土壤肥力, 通常以石膏、磷石膏、硅酸鈣和工業生產中含鈣廢棄物等為主要原料, 或配施有機物料為改良劑進行鹽堿地改良。微生物菌劑與磷石膏聯合施用改良鹽堿土能夠增強鹽堿土的改良效果[3]。目前國內外改良鹽堿土技術主要包括物理改良、水利工程改良、化學改良和生物改良等, 但費用高、見效慢、推廣難, 嚴重制約了鹽堿土區種植業結構的抗風險能力和生態環境的可持續發展。

植物根際微生境是土壤中活性最強的微生境, 也是植物獲取養分的主要區域。在此微域里, 植物—微生物—土壤—環境間相互作用, 共同維持著根際微生態系統的平衡, 影響作物生產[4]。研究表明, 鹽脅迫使土壤微生物群落結構和組成發生變化, 植物根際微生態區系失衡, 可能涉及到某些功能菌群、有益菌群和有害菌群比例的平衡失調[5]。根際微生物類型和其群落組成受鹽脅迫及程度的影響, 其數量和群落組成與環境和植物種類有關[6], 施用鈣肥可改善鹽堿土壤結構、降低Na+含量、調節pH[3]。關于鹽堿土壤微生物群落的研究主要集中在不同植物群落土壤微生物的時空動態和區系特征[7-8]、微生物量與土壤酶[9]、不同改良措施及植被類型對鹽堿土壤微生物的影響效果等[9-12], 有關不同鹽脅迫強度下配施鈣肥對花生根際微生物多樣性、種群結構及其生態關系的研究鮮見報道, 為此, 本試驗以不同鹽脅迫強度下配施鈣肥處理的花生根際土壤為研究對象, 利用高通量測序技術對土壤中細菌種類及多樣性進行系統分析, 探討鹽脅迫下不同生育時期花生根際土壤細菌群落組成和多樣性的變化, 有助于認識鹽堿地花生根際土壤微生物群落功能調節并發掘新的功能類群, 對改良鹽堿土壤和開發利用鹽堿土區資源、提高鹽堿土區農業種植結構抗風險能力及鹽堿地花生高產優質栽培均具重要意義。

1 材料與方法

1.1 試驗材料

采用山東省花生研究所萊西試驗站耕層土壤(0~30 cm)作為供試土壤, 含土壤有機質含量13.23 g kg-1、全磷(P2O5) 0.84 g kg-1、全鉀(K2O) 11.53 g kg-1、全氮1.70 g kg-1、水解氮(N) 92.1 mg kg-1、速效磷(P2O5) 11.7 mg kg-1、速效鉀(K2O) 112.2 mg kg-1, 土壤pH 7.1, 土壤交換性Na+0.82 cmol kg-1, Cl–28.5 mg kg-1, 土壤含水量9.72%。選擇花生品種花育25號(HY25)。

1.2 試驗設計

試驗于山東省花生研究所試驗站防雨棚中進行,土壤風干、過篩(1 cm)后裝入同批次、同規格塑料盆中(底部直徑為36 cm, 高為26 cm), 每盆裝土量為18 kg。設置0 (Y0)、1.5 (Y1)、3.0 g kg-1(Y2) 3個鹽脅迫梯度, 采用分析純NaCl重量法控制鹽脅迫梯度。各處理均配合基施鈣肥(分析純CaO) 150 kg hm-2(C, 150 kg hm-2CaO), 以過篩后不施鈣肥的土壤為對照, 隨機排列, 6次重復。播種前將供試肥料按N∶P2O5∶K2O = 1.0∶1.5∶1.5 (15.0、22.5、22.5 kg hm-2)比例連同CaO以基肥形式均勻施入盆中。隨后各處理均澆入2 L清水沉實, 待土壤墑情適宜時再次混勻每盆中的土壤, 以使土壤含鹽量和施入的鈣肥均勻分布。選取飽滿均勻的種子, 每盆播6粒, 播深均為3 cm, 留4株長勢一致的幼苗。

分別于開花下針期(播后75 d)和收獲期(播后125 d)以處理和重復為單位采集花生植株根際土壤樣本。根際土壤樣本的采集以處理為單元進行, 采用多點混合樣本采集方法, 以盆為單位分別小心地從盆中取出植株, 去除根際附著不緊密的土壤, 然后用無菌刷收集附著緊密的土壤, 混合每盆中植株根際土壤, 每3重復根際土壤樣品混合為1個生物學樣本重復, 每處理均獲取2個生物學重復, 封入無菌袋, 置于冰盒帶回實驗室, 保存于-80℃冰箱備用。于收獲期, 各處理收獲實測計產, 室內考察花生出米率、百果質量、百仁質量和莢果產量。

開花下針期各處理土壤樣本分別表示為CK、CY0、CY1、CY2, 收獲期各處理土壤樣本分別表示為HCK、HCY0、HCY1、HCY2。由北京諾賽基因公司進行相關測試。

1.3 土壤DNA提取

收集所得土壤樣品, 利用OMEGA土壤總DNA提取試劑(OMEGA soil DNA kit)盒提取DNA。采用1.5%瓊脂糖凝膠電泳和Nanodrop 2000分光光度計檢測DNA的純度和濃度。

1.4 16srRNA文庫構建及高通量測序

以提取的土壤樣本DNA為模板, 利用引物340F: 5￠-CCTACGGGNBGCASCAG-3￠以及805R: 5￠-GACTACNVGGGTATCTAATCC-3￠, 對土壤細菌16s rRNA基因V3~V4區進行PCR擴增。擴增程序如下: 95℃預變性3 min; 95℃30 s, 50℃30 s, 72℃60 s, 30個循環; 72℃7 min。使用1.5%濃度的瓊脂糖凝膠電泳檢測PCR產物; 根據PCR產物濃度進行等濃度混樣, 充分混勻后, 使用0.5×TBE濃度1.5%的瓊脂糖膠電泳純化PCR產物, 割膠回收目標條帶。產物純化試劑盒使用QIAGEN公司的MinElute膠回收試劑盒。最后使用HiSeq2500進行250PE測序。

1.5 生物信息學分析

對花生根際土壤微生物的多樣性和豐富度進行樣本內物種多樣性即Alpha多樣性和樣本間Beta多樣性分析, 包括繪制稀釋性曲線(Rarefaction curve)和分析Alpha多樣性指數(包括Chao、Shannon、ACE、Simpson、coverage和sobs等指數)。同時進行群落結構統計分析和基于OTUs的物種組成聚類分析(heatmap分析), 最后, 進行反映樣本間的多樣性距離關系和生物群落間分化程度的Beta多樣性分析(PCA統計分析), 挖掘樣本間物種組成差異和生物功能分析。

2 結果與分析

2.1 對莢果性狀和產量的影響

由表1可知, 鹽脅迫處理明顯影響花生出米率、百果質量、百仁質量和莢果產量, 并隨鹽脅迫程度的提高作用顯著, 施用鈣肥可顯著提高無鹽脅迫處理下的百果質量和莢果產量, 對出米率和百仁質量有促進作用但不顯著。鹽脅迫處理下百仁質量較CK和CY0分別降低12.94%、20.44%和38.70%、43.98%, 較高鹽脅迫處理的莢果產量僅為CK處理的1/2, CY0處理下產量是CY2處理2.5倍。表明鹽脅迫嚴重抑制花生籽仁發育和產量形成。

表1 施鈣對鹽脅迫下花生產量及構成因素的影響

CK、CY0、CY1、CY2分別表示開花下針期對照和各鹽脅迫梯度下施鈣肥各處理根際樣本; HCK、HCY0、HCY1、HCY2分別表示收獲期相應鹽脅迫下各處理土壤樣本。同列不同小寫字母表示處理間差異顯著(< 0.05)。

CK, CY0, CY1, and CY2 represent the rhizosphere samples of each treatment with calcium fertilizer under each salt stress gradient at the needle stage at the flowering and needle stage, respectively; HCK, HCY0, HCY1, and HCY2 represent the soil samples of each treatment under the corresponding salt stress at the harvest stage, respectively. Values within a column followed by different lowercase letters mean significant differences among the treatments at the 0.05 probability level.

2.2 花生根際微生物群落測序數據分析

供試的16個根際土壤樣本共檢測出細菌30個門、73個綱、204個目、351個科、650個屬和1279個種。獲得有效序列1,187,782.5條, 其中HCY2樣本有效序列條數最小為127,455.5條, CY0樣本最大為174,409.0條, 多數樣本序列平均長度在445.0~452.0 bp之間, 平均長度482.0 bp (表2)。通過維恩圖分析, 各處理根際土壤樣本共有2446個OTUs, 存在于CY1、CK、CY2中的OTUs數分別為29、25和23個, 而HCK、HCY0和HCY2處理樣本中均不足10個, 表明各處理樣本微生物組成相似, 鹽脅迫強度、外源施鈣肥和生育時期對根際土壤細菌群落結構均無顯著影響(圖1)。

表2 花生根際微生物群落測序質量

處理同表1。Treatments are the same as those given in Table 1.

處理同表1。Treatments are the same as those given in Table 1.

2.3 Alpha多樣性分析

2.3.1 稀釋曲線 由圖2可知, 隨測序數據量的增加, 物種豐富度呈現前期迅速增加的趨勢, 各處理樣本在測序量達到30,000條以上時, 稀釋曲線趨近平緩, 各樣本微生物群落測序數據達到飽和, 測序量能夠覆蓋花生根際微生物組群落的絕大部分物種。各處理樣本重復間取其均值進行比較發現, 對照、鹽脅迫強度和生育時期及外源施用鈣肥各處理OTUs數量基本一致, 相互間均無顯著差異。

2.3.2 多樣性指數 由表3可知, 供試樣本的coverage測序深度指數均在98.90%以上, 樣本微生物群落中OTUs的chao和ACE菌群豐富度指數分別在3522.40~3762.58和3448.57~3685.00之間, Shannon多樣性指數為6.57~7.11, Simpson和sobs指數分別在0.001,946~0.005,682和3113.00~3428.50之間。ACE、chao和sobs指數均以CK處理最小, ACE、chao豐富度指數以CY2處理最大, 而Shannon和sobs指數則均以CY0最高。可見, 配施鈣肥可提高鹽脅迫下花生根際微生物群落豐富度和多樣性, 且對開花下針期的影響作用優于收獲期, 無鹽脅迫下施用鈣肥處理根際土壤細菌多樣性更為豐富。

2.4 物種群落組成分析

2.4.1 全樣本門水平菌落結構分析 所有樣本細菌群落組成的30個門中, 其Top10 (others合并< 0.025)相對豐度之和達93.40%。優勢菌門均為變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、Patescibacteria、酸桿菌門(Acidobacteria)、綠彎菌門(Chloroflexi)、疣微菌門(Verrucomicrobia)、芽單胞菌門(Gemmatimonadetes)和擬桿菌門(Bacteroidetes)等8種, 其相對豐度均在5.17%以上, 其中以變形菌門(Proteobacteria)平均相對豐度最高(31.70%), 放線菌門(Actinobacteria)次之(16.12%), 綠彎菌門(Chloroflexi)和酸桿菌門(Acidobacteria)分別為8.23%、8.65%。芽單胞菌門(Gemmatimonadetes)僅為變形菌門(Proteobacteria)的16.32% (圖3-A)。鹽脅迫提高根際變形菌門(Proteobacteria)的相對豐度, 降低放線菌門(Actinobacteria)豐度且均隨鹽脅迫強度提高其升降幅度增大, 尤以開花下針期作用明顯, 并均受外源配施鈣肥的正向調節, 開花下針期和收獲期不施鈣肥和施鈣肥高鹽脅迫處理的變形菌門(Proteobacteria)相對豐度較其對照分別提高27.00%、10.46%和32.44%%、50.09%, 而放線菌門(Actinobacteria)豐度較其對照分別降低55.58%、48.83%和40.90%、28.87%。較高鹽脅迫處理的擬桿菌門豐度明顯升高, 并以收獲期較為明顯。

處理同表1。Treatments are the same as those given in Table 1.

處理同表1。Treatments are the same as those given in Table 1.

對各處理根際樣本進行優勢門(Top10) (圖3-B)組間差異顯著性檢驗分析表明, 各處理間均呈顯著或極顯著差異。生育時期、施用鈣肥和鹽脅迫各因子均使得花生根際微生物群落分異明顯。

2.4.2 全樣本目水平菌落結構分析 所有樣本共包含204個菌目, 其中Top16 (others合并<0.03)相對豐度之和達60.40%, 相對豐度在5.0%以上的優勢菌目分別為Saccharimonadales (7.94%)、β-變形菌目(Betaproteobacteria, 7.42%)、鞘脂單胞菌目(Sphingomonadales, 6.41%)、芽單胞菌目(Gemmatimonadales, 5.08%)和根瘤菌目(Rhizobiales, 5.00%)。除優勢菌目β-變形菌目(Betaproteobacteria)和外, 其他優勢菌目開花下針期樣品的相對豐度均明顯低于收獲期。鹽脅迫下基施鈣肥使β-變形菌目(Betaproteobacteria)、芽單胞菌目(Gemmatimonadales)和鞘脂單胞菌目(Sphingomonadales)豐度顯著升高, 而Saccharimonadales、微球菌目()菌目相對豐度顯著降低, 并隨鹽脅迫強度升高其相對豐度呈降低趨勢, 且收獲期降幅較大(圖4-A)。

處理同表1。*、**、***分別表示處理間達0.05、0.01和0.001顯著性水平。

Treatments are the same as those given in Table 1. *, **, and *** indicate significant differencesatthe 0.05, 0.01, and 0.001 levels, respectively.

處理同表1。*、**、***分別表示處理間達0.05、0.01和0.001顯著性水平。

Treatments are the same as those given in Table 1. *, **, and *** indicatesignificant differences at the 0.05, 0.01, and 0.001levels, respectively.

各樣本細菌群落目水平物種組成heatmap圖顯示, 不同處理根際樣本細菌目水平豐度高低有異, CY1處理的β-變形菌目、土圈菌目(Pedosphaerales)和unclassified_k__norank_d__Bacteria菌目豐度均顯著升高, CY2、HCY2兩處理均使β-變形菌目和鞘脂單胞菌目(Sphingomonadales)豐度明顯升高, 而HCY0和HCK則使Saccharimonadales菌目豐度升幅極為明顯。聚類樹分析結果表明, 同一生育時期根際樣本優勢菌目組成較為相似聚為一類, 鹽脅迫處理與非鹽脅迫處理根際樣本中的優勢菌目組成存在差異各聚為一類(圖4-B)。對各處理根際樣本進行優勢目(Top10) (圖4-C)組間差異顯著性檢驗分析表明, 除根瘤菌目(Rhizobiales)外各處理間均呈顯著或極顯著差異。

2.4.3 全樣本屬水平菌落結構分析 各樣本共包含650個菌屬, 其中、鞘氨醇單胞菌屬)、、、芽單胞桿菌屬()、等6種為優勢菌屬。各處理樣本中的和兩優勢菌屬豐度隨鹽脅迫強度升高而升高,、、和非優勢菌屬菌屬則均隨脅迫強度升高而降低,和兩生育時期分別較其相應對照降低72.27%、67.08%和81.20%、51.75%, 非優勢菌屬羅河桿菌屬()和豐祐菌屬()相對豐度則均隨鹽脅迫強度增加顯著升高, 分別是相應對照的1.31、4.23倍和3.02、2.55倍。鹽脅迫顯著影響不同生育時期根際土壤細菌群落組成, 尤以收獲期優勢菌群豐度降低顯著, 而非優勢菌群豐度往往升幅明顯(圖5-A)。

各樣本細菌群落屬水平物種組成heatmap圖和聚類樹分析結果表明, 較高濃度鹽脅迫處理樣本優勢菌屬組成相似聚為一組, 非鹽脅迫處理樣本屬豐度依生育時期相同而相近, 各聚為一組, 1.5 g kg-1較低鹽脅迫處理根際土壤細菌優勢菌屬較為豐富(圖5-B)。

2.5 樣本間Beta多樣性分析

通過PCA分析發現, 各處理樣本分布于不同象限, 其細菌群落組成具有明顯差異, PC1和PC2對結果的解釋度分別為25.88%和15.87%, 二者分別由生育時期因子和鹽脅迫因子貢獻。HCK和HCY0與第1、第2主成分均呈正相關關系, CY2和CY1則與第1、第2主成分均呈負相關關系, HCY1和HCY2與第1主成分為正相關關系而與第2主成分為負相關關系, 而CY0和CK則與之相反(圖6-A)。說明鹽脅迫、生育時期和施鈣肥與否對根際土壤微生物菌群類型影響較大。層級聚類分析表明, 依生育時期聚為兩大類的同時, 又以鹽脅迫及其強度各聚一簇, 開花下針期受鹽脅迫強度影響較小, 成熟期根際土壤微生物菌群類型影響較大(圖6-B)。對種水平上菌群分型分析結果亦表明, 鹽脅迫強度、生育時期對花生根際土壤微生物菌群類型影響較大, 基施鈣肥影響較小(圖6-C)。

處理同表1。Treatments are the same as those given in Table 1.

處理同表1。Treatments are the same as those given in Table 1.

2.6 16S功能預測分析

通過KEGG代謝途徑的預測差異分析發現, CK、HCK、HCY2處理樣本中微生物群落的功能豐度均明顯降低, 尤以次生代謝產物的代謝與合成(丙酸酯代謝、2-氧代羧酸代謝、乙醛酸和二羧酸酯代謝、氨?；?tRNA的生物合成)、聚糖的生物合成與代謝(糖酵解/糖異生、氨基糖和核苷酸糖代謝)、氨基酸代謝(半胱氨酸和蛋氨酸的代謝、甘氨酸, 絲氨酸和蘇氨酸的代謝)、脂肪酸代謝等功能基因的豐度降低顯著。CY0、CY1、CY2、HCY0和HCY1樣本中各功能基因豐度均相對升高, 尤以CY1和HCY0升高明顯(圖7)。表明無鹽脅迫和全生育期高強度鹽脅迫可明顯降低次生代謝產物的代謝與合成、聚糖的生物合成與代謝、氨基酸代謝和脂肪酸代謝等功能基因在根際土壤中的富集, 花生旺盛生長期、低鹽脅迫和基施鈣肥處理使得功能基因豐度大幅提高。

處理同表1。Treatments are the same as those given in Table 1.

3 討論

根際微生物是土壤—根系間養分轉化和轉運的調節器, 土壤微生物群落多樣性反映了群落總體的動態變化。植被對土壤微生物群落的數量具有顯著正向的影響, 鹽脅迫下, 除土壤微生物生長受到抑制, 植物本身的生理機能也迅速遭到破壞, 根際土壤細菌群落的豐富度顯著降低[13-14], 隨鹽脅迫強度增加微生物種群數量、群落結構和功能穩定性明顯降低[15-17]。耐鹽植物根際土壤微生物量隨種植時間的延長而增加并顯著高于非根際土壤[18-19]。但不同鹽堿地類型花生根際土壤微生物種類、優勢種群數量和群落功能多樣性存在差異, 含鹽量較高土壤更為豐富。黃河三角洲濱海鹽堿土和不同含鹽量的鹽堿土花生根際土壤微生物優勢菌門均以變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、綠彎菌門(Chloroflexi)和酸桿菌門(Acidobacteria)等4種菌群為優勢菌群[20]。鹽堿土花生根際芽單胞菌門(Gemmatimonadetes)和擬桿菌門(Bacteroidetes)的豐度顯著高于非鹽堿土壤約2倍[21]。本試驗條件下, 鹽脅迫強度和花生生長期不影響花生根際細菌群落多樣性和豐富度, 各處理微生物種群數量和優勢菌群無差異, 但對細菌群落結構組成有影響。各處理優勢菌門均為變形菌門、放線菌門、Patescibacteria、酸桿菌門、綠彎菌門、疣微菌門(Verrucomicrobia)、芽單胞菌門和擬桿菌門。優勢菌目為Saccharimonadales、β-變形菌目、鞘脂單胞菌目、芽單胞菌目和根瘤菌目。鹽脅迫使根際微生物優勢菌群數量明顯增加, 旺盛生長的開花下針期較為明顯。

耐鹽植物根際環境中細菌群落組成與耐鹽性有關, 鹽生植物根際細菌群落受土壤鹽分的影響較小[22]。鹽脅迫顯著提高可加快碳氮循環的藍藻菌門的含量, 而益于花生抵御逆境脅迫[20]。中度耐鹽植物根際環境中耐鹽性較高的變形菌門和厚壁菌門豐度較高, 輕度耐鹽植物根際環境中酸桿菌門和芽單胞菌門較豐富, 說明酸桿菌門和芽單胞菌門不具耐鹽性或更易受高鹽影響。本試驗條件下, 花生根際土壤優勢菌群除變形菌門和酸桿菌門外, 還包括放線菌門、綠彎菌門、擬桿菌門和疣微菌門。酸桿菌門在土壤生態過程中起著重要作用, 變形菌門的增多可更有效固定氮源。放線菌門具有共生固氮和解磷作用[23-26], 綠彎菌存在于根際土壤樣本中, 通過光合作用產生能量但不能固氮, 但具有較好的生物解磷作用[27-28]。擬桿菌門具有非常強的營養物質代謝能力, 如復雜有機物、蛋白質和脂類等[29]。解磷菌懸液可顯著增加鹽堿土壤酶活性和土壤微生物數量, 促進與土壤營養元素循環相關的菌屬酸桿菌(Acidobacteria)、綠彎菌(Chloroflexi)和浮霉菌門(Planctomycetes)等微生物群落富集, 全面改善鹽堿土性質, 促進植株生長[3]。本試驗結果表明, 鹽脅迫顯著影響不同生育時期根際土壤細菌群落組成, 尤以收獲期優勢菌屬豐度降低顯著, 而非優勢菌屬豐度往往升幅明顯。開花下針期β-變形菌目(Betaproteobacteria)和芽單胞菌目(Gemmatimonadales)明顯高于收獲期。表明花生不同生長期和外源施鈣對鹽脅迫下花生根際微生物群落具有適應性和調節功能, 根際土壤較高的變形菌門、放線菌門、擬桿菌門和酸桿菌門可能有助于維持根系養分吸收和微環境的平衡, 以改善土壤環境, 增強抵御鹽脅迫的能力, 是鹽堿土區生物改良鹽脅迫的寶貴資源。同時, 鹽脅迫可大幅提高變形菌門、擬桿菌門的相對豐度, 且配施鈣肥的調節作用明顯增強并保持花生全生育期, 從而改善了花生根際氮磷及營養物質代謝的動態平衡而增強了對鹽脅迫的抵御作用。

土壤類型和花生根系分布密集程度對土壤微生物菌群類型影響較大, 鹽堿土壤花生根際微生物菌群類型依據土壤含鹽量高低和根系分布深度不同而聚為不同簇, 非鹽堿土壤、鹽堿土壤根系分布不同密集程度層0~20 cm、20~40 cm各歸為一類。非根際土壤與根際微生物群落分離明顯, 非根際土壤群落之間, 以及根際土壤群落之間各具有較高的相似性[20,30-31]。本試驗條件下層級聚類分析表明, 生育時期、鹽脅迫及強度各聚為一簇, 花生根際土壤微生物菌群類型受鹽脅迫強度和生育時期影響較大, 基施鈣肥影響較小, 并以花生生育末期根際微生物菌群類型受影響較大。

16S功能預測顯示, 鹽脅迫影響花生根際微生物群落的功能豐度譜, 尤以某些代謝與合成相關功能如次生代謝產物的代謝與合成、聚糖生物合成與代謝、氨基酸代謝等功能基因豐度變化明顯。在非鹽脅迫和全生育期高強度鹽脅迫處理下, 這些次生代謝產物的代謝與合成、聚糖生物合成與代謝、氨基酸代謝等功能基因豐度顯著降低, 而在花生旺盛生長期、低鹽脅迫和基施鈣肥處理均使得這些功能基因豐度大幅提高。說明花生旺盛生長時期、較低的鹽脅迫強度和外源鈣肥均可提高花生根際微生物與代謝與合成相關功能基因豐度, 以抵御鹽脅迫逆境。然而, 花生根際微生物的富集與代謝和合成相關的功能基因不足以抵抗全生育期內高強度鹽脅迫逆境, 致使花生產量顯著降低。

4 結論

花生根際土壤微生物優勢菌群組成結構不受鹽脅迫和生育時期影響, 但其相對豐度與鹽脅迫強度有關, 隨鹽脅迫強度的提高, 變形菌門(Proteobacteria)的相對豐度和放線菌門(Actinobacteria)豐度升降幅度增大。鹽脅迫強度、生育時期對花生根際土壤微生物菌群類型影響程度大于基施鈣肥處理?；ㄉ⑸L時期和基施鈣肥使得鹽脅迫下花生根際微生物菌群功能基因豐度大幅提高, 高鹽脅迫明顯降低次生代謝產物、聚糖的生物合成與代謝, 以及氨基酸和脂肪酸代謝等功能基因在根際土壤中的富集。通過改良根際土壤微生物環境提高植物耐鹽脅迫能力是生物改良鹽堿脅迫的可行途徑。

[1] Jamil A, Riaz S, Ashraf M, Foolad M R. Gene expression profi-ling of plants under salts tress., 2011, 30: 435–458.

[2] 賈敬敦, 張富. 依靠科技創新推進我國鹽堿地資源可持續利用. 中國農業科技導報, 2014, 16(5): 1–7. Jia J D, Zhang F. Sustainable utilization of saline-alkali land resources through scientific and technological innovation in China., 2014, 16(5): 1–7 (in Chinese with English abstract).

[3] 姜煥煥. 耐鹽堿解磷菌與磷石膏聯用改良鹽堿土的效果與機制. 哈爾濱工業大學博士學位論文, 黑龍江哈爾濱, 2019. Jiang H H. Saline-alkali Soil Remediation by the Combined Application of Halotolerant Phosphate Solubilizing Microorganism and Rock Phosphate. PhD Dissertation of Harbin Institute of Technology, Harbin, Heilongjiang, China, 2019 (in Chinese with English abstract).

[4] 張旭龍, 馬淼, 吳振振, 張志政, 高睿, 石靈玉. 不同油葵品種對鹽堿地根際土壤酶活性及微生物群落功能多樣性的影響. 生態學報, 2017, 37: 1659–1666. Zhang X L, Ma M, Wu Z Z, Zhang Z Z, Gao R, Shi L Y. Effects of Helianthus annuus varieties on rhizosphere soil enzyme activities and microbial community functional diversity of saline-alkali land in Xinjiang., 2017, 37: 1659–1666 (in Chinese with English abstract)．

[5] Zhang Q, Yu X. Allelopathy in replant problem in forest soil., 2001, 8: 51–64.

[6] Kielak A, Pijl A S, Veen J A, Kowalchuk G A. Differences in vegetation composition and plant species identity lead to only minor changes in soil-borne microbial communities in a former arable field., 2008, 63: 372–382.

[7] Nelson E B. Microbial dynamics and interactions in the spermosphere., 2004, 42: 271–309.

[8] Buyer J S, Roberts D P, Russek-Cohen E. Microbial community structure and function in the spermosphere as affected by soil and seed type., 1999, 45: 138–144.

[9] 葉淑紅, 王艷, 萬惠萍, 吉云秀, 林學政, 丁德文. 遼東灣濕地微生物量與土壤酶的研究. 土壤通報, 2006, 37: 897–900． Ye S H, Wang Y, Wan H P, Ji Y X, Lin X Z, Ding D W. Wetland microbial biomass and soil enzyme activities in Liaodong bay., 2006, 37: 897–900 (in Chinese with English abstract).

[10] 李鳳霞, 王學琴, 郭永忠, 許興, 楊建國, 季艷清. 寧夏不同類型鹽漬化土壤微生物區系及多樣性. 水土保持學報, 2011, 25(5): 107–111． Li F X, Wang X Q, Guo Y Z, Xu X, Yang J G, Ji Y Q. Microbial flora and diversity in different types of saline-alkali soilin Ningxia., 2011, 25(5): 107–111 (in Chinese with English abstract).

[11] 田家怡, 李甲亮, 陳印平, 李建慶, 于祥. 黃河三角洲外來入侵物種米草對底泥微生物群落的影響. 海洋湖沼通報, 2009, (4): 157–162． Tian J Y, Li J L, Chen Y P, Li J Q, Yu X. Effect on the sludge microbial communities by alien invasive species Spartina spp. in yellow river delta., 2009, (4): 157–162 (in Chinese with English abstract).

[12] 王震宇, 辛遠征, 李鋒民, 高冬梅. 黃河三角洲退化濕地微生物特性的研究. 中國海洋大學學報, 2009, 39: 1005–1012． Wang Z Y, Xin Y Z, Li F M, Gao D M. Microbial community characteristics in a degraded wetland of the Yellow River Delta., 2009, 39: 1005–1012 (in Chinese with English abstract).

[13] 王文銅. 花生根際土壤細菌群落結構分析. 山東農業大學博士學位論文, 山東泰安, 2012. Wang W T. Analysis of Bacterial Community Structure in Rhizosphere Soil of Peanut. PhD Dissertation of Shandong Agricultural University, Tai’an, Shandong, China, 2012 (in Chinese with English abstract).

[14] 吳鳳芝, 包靜, 劉淑芹. 鹽脅迫對黃瓜根際土壤細菌群落結構和生長發育的影響. 園藝學報, 2010, 37: 741–748. Wu F Z, Bao J, Liu S Q. Effects of salt stress on rhizospheric soil bacterial community structure and cucumber yield., 2010, 37: 741–748 (in Chinese with English abstract).

[15] 朱泓, 王小敏, 黃濤, 吳文龍, 李維林. NaCl脅迫對濱梅根際細菌群落多樣性及優勢菌群的影響. 南京林業大學學報(自然科學版), 2017, 41(4): 49–54． Zhu H, Wang X M, Huang T, Wu W L, Li W L. Effect of NaCl stress on bacterial community diversity and core microbiome in rhizosphere and bulk soil of beach plum (Marshall).(Nat Sci Edn), 2017, 41(4): 49–54 (in Chinese with English abstract)．

[16] Baumann K, Dignac M F, Rumpel C, Bardoux G, Sarr A, Steffens M, Maron P A. Soil microbial diversity affects soil organic matter decomposition in a silty grass land soil., 2013, 114: 201–212.

[17] Philippot L, Spor A, Hénaul T C, Bru D, Bizouard F, Jones C M, Sarr A, Maron P. Loss in microbial diversity affects nitrogen cycling in soil., 2013, 7: 1609–1619.

[18] 馬曉梅, 尹林克, 陳理. 塔里木河干流胡楊和檉柳根際土壤微生物及其垂直分布. 干旱區研究, 2008, 25(2): 183–189． Ma X M, Yin L K, Chen L. Study on vertical distribution of microorganiss in rhizosphere ofandsp. in the lower reaches of the Tarim river, Xinjiang., 2008, 25(2): 183–189 (in Chinese with English abstract)．

[19] 林學政, 陳靠山, 何培青, 沈繼紅, 黃曉航. 種植鹽地堿蓬改良濱海鹽漬土對土壤微生物區系的影響. 生態學報, 2006, 26: 801–807. Lin X Z, Chen K S, He P Q, Shen J H, Huang X H. The effects of Suaeda salsa L. planting on the soil microflora in coastal saline soil., 2006, 26: 801–807 (in Chinese with English abstract).

[20] 張智猛, 慈敦偉, 張冠初, 丁紅, 楊吉順, 戴良香, 張岱. 山東地區鹽堿土花生種子際土壤微生物群落結構的研究. 微生物學報, 2017, 57: 582–596. Zhang Z M, Ci D W, Zhang G C, Ding H, Yang J S, Dai L X, Zhang D. Diversity of microbial community structure in the spermosphere of saline-alkali soil in Shandong area., 2017, 57: 582–596 (in Chinese with English abstract)

[21] 戴良香, 康濤, 慈敦偉, 丁紅, 徐揚, 張智猛, 張岱, 李文金. 黃河三角洲鹽堿地花生根層土壤菌群結構多樣性. 生態學報, 2019, 39: 7169–7178． Dai L X, Kang T, Ci D W, Ding H, Xu Y, Zhang Z M, Zhang D, Li W J. Comparison of the microbial community in the rhizosphere of peanuts between salinealkali and non-saline soil at different soil depths and intercropping cultivation in the Yellow River Delta., 2019, 39: 7169–7178．

[22] Yongge Y, Brune C, Kleunen M, Li J M, Jin Z X. Salinity-induced changes in the rhizosphere microbiome improve salt tolerance of Hibiscus hamabo., 2019, 443: 525–537.

[23] 趙嬌, 謝慧君, 張建. 黃河三角洲鹽堿土根際微環境的微生物多樣性及理化性質分析. 環境科學, 2020, 41: 1449–1455. Zhao J, Xie H J, Zhang J. Microbial diversity and physicochemical properties of rhizosphere microenvironment in saline-alkali soils of the Yellow River Delta., 2020, 41: 1449–1455 (in Chinese with English abstract).

[24] 鄭微. 天津不同鹽度土壤微生物群落結構特征及其適應機理研究. 天津師范大學博士學位論文, 天津, 2018. Zheng W. Study on the Structural Characteristics and Adaptation Mechanism of Microbial Community in Different Salinity Soils in Tianjin. PhD Dissertation of Tianjin Normal University, Tianjin, China, 2018 (in Chinese with English abstract).

[25] 涂勇, 楊文鈺, 劉衛國, 雍太文, 江連強, 王小春. 大豆與烤煙不同套作年限對根際土壤微生物數量的影響. 作物學報, 2015, 41: 733–742. Tu Y, Yang W Y, Liu W G, Yong T W, Jiang L Q, Wang X C. Effects of relay strip intercropping years between flue-cured tobacco and soybean on rhizospheric microbes quantities., 2015, 41: 733–742 (in Chinese with English abstract).

[26] Cui J, Li Y, Wang C, Kim K S, Wang T, Liu S. Characteristics of the rhizosphere bacterial community across different cultivation years in saline-alkaline paddy soils of Songnen Plain of China., 2018, 64: 925–936.

[27] Yu H, Si P, Shao W, Qiao X, Yang X, Gao D, Wang Z. Response of enzyme activities and microbial communities to soil amendment with sugar alcohols., 2016, 5: 604–615.

[28] Kragelund C, Levantesi C, Borger A. Identity abundance and ecophysiology of filamentous Chloroflexi species present in activated sludge treatment plants., 2007, 59: 671–682．

[29] Hill V R, Kahler A M, Jothikumar N, Johnson T B, Hahn D, Cromeans T L. Multistate evaluation of an ultrafiltration-based procedure for simultaneous recovery of enteric microbes in 100-liter tap water sample., 2007, 73: 4218–4225.

[30] 葛應蘭, 孫廷. 馬鈴薯根際與非根際土壤微生物群落結構及多樣性特征. 生態環境學報, 2020, 29: 141–148. Ge Y L, Sun T. Soil microbial community structure and diversity of potato in rhizosphere and non-rhizosphere soil., 2020, 29: 141–148 (in Chinese with English abstract).

[31] 覃瀟敏, 鄭毅, 湯利, 龍光強. 玉米與馬鈴薯間作對根際微生物群落結構和多樣性的影響. 作物學報, 2015, 41: 919–928. Qin X M, Zheng Y, Tang L, Long G Q. Effects of maize and potato intercropping on rhizosphere microbial community structure and diversity., 2015, 41: 919-928 (in Chinese with English abstract).

Response of rhizosphere bacterial community diversity to salt stress in peanut

DAI Liang-Xiang1,**, XU Yang1,**, ZHANG Guan-Chu1, SHI Xiao-Long2, QIN Fei-Fei1, DING Hong1,*, and ZHANG Zhi-Meng1,*

1Shandong Peanut Research Institute, Qingdao 266100, Shandong, China;2College of Agronomy, Shenyang Agricultural University, Shenyang 110866, Liaoning, China

To characterize the peanut rhizosphere bacteria community in response to salt stress, a pot experiment was performed with different salt concentrations. The peanut rhizosphere soils at flowering and mature stages were sampled to extract DNA for constructing bacterial 16S rRNA gene library, and then high-throughput sequencing was performed for sequencing and bioinformatics analysis. The results showed thatProteobacteria, Actinobacteria, Patescibacteria, Acidobacteria, and Chloroflexi were the dominant phyla, and the orders Saccharimonadales, Betaproteobacteria, Sphingomonadales, Gemmatimonadales, and Rhizobiales were dominated in the peanut rhizosphere soils. Comparisons of the bacterial community structure of peanuts revealed that the relative abundance of Proteobacteria dramatically increased, while that of Actinobacteria decreased in salt-treated soils, and the fluctuation increased with the increase of the salt concentration. Moreover, applying calcium fertilizer under salt stress increased the abundance of Betaproteobacteria, Gemmatimonadales, and Sphingomonadales, which were affected by salt stress, growth stages, and exogenous calcium application. Cluster analysis revealed that the dominant bacteria of soil groups with high salt concentration were similar and clustered together, while the soil samples of the same growth period were similar and clustered together according to the bacterial structure at the genus level under non-salt stress conditions. Bacterial community structure differed in the growth stages and soil salt concentrations, whereas the differences of soil groups with or without calcium application were relatively small. Function prediction analysis indicated that the sequences related to secondary metabolites, glycan biosynthesis and metabolism, and amino acid and lipid metabolism were enriched in high salt-treated soils. The functional groups increased significantly during the fast-growth period, low salt stress, and basal calcium fertilizer treatments, which may play an important role on the growth and stress response in peanut. This study of microbial communities could lay the foundation for future improvement of stress tolerance of peanuts via modification of the soil microbes.

peanut (); salt stress; rhizosphere; soil microbial community; 16S rRNA gene

10.3724/SP.J.1006.2021.04160

本研究由國家自然科學基金項目(31971856,31971854, 31901574), 山東省現代農業產業技術體系創新團隊(花生)項目(SDAIT-04-06), 山東省農業科學院創新工程項目(CXGC2018B05)資助。

This study was supported by the National Natural Science Foundation of China (31971856, 31971854, 31901574), the Modern Agricultural Industry Technical System of Shandong Province (SDAIT-04-06), and the Agricultural Scientific and Technological Innovation Project of Shandong Academy of Agricultural Sciences (CXGC2018B05).

張智猛, E-mail: qinhdao@126.com; 丁紅, E-mail: dingpeanut@163.com

**同等貢獻(Contributed equally to this work)

戴良香, E-mail: liangxiangd@163.com

2020-07-17;

2021-01-13;

2021-02-18.

URL: https://kns.cnki.net/kcms/detail/11.1809.S.20210216.1438.004.html