馬少蘭 馬彩霞 徐鵬鑫 鄭國琦 納小凡
(寧夏大學(xué)生命科學(xué)學(xué)院,銀川 750021)
寧夏枸杞(Lycium barbarumL.)是茄科枸杞屬多年生落葉灌木,廣泛分布于中國、韓國、日本、歐洲、北美和地中海等地區(qū)的干旱和半干旱區(qū)域[1]。由于近期研究證實(shí)寧夏枸杞果實(shí)(枸杞子)具有抗衰老、保護(hù)神經(jīng)、消炎、促進(jìn)代謝、控制血糖、免疫調(diào)節(jié)及抗腫瘤等多種生物學(xué)活性[1-2],枸杞子已作為食品和營養(yǎng)保健品在亞洲及歐美等多國零售市場中銷售[3-5]。受國際市場追捧,枸杞種植業(yè)已成為我國西北省份,尤其是寧夏出口創(chuàng)匯的優(yōu)勢富民產(chǎn)業(yè)之一,為地區(qū)經(jīng)濟(jì)發(fā)展做出巨大貢獻(xiàn)[6]。
然而,受枸杞自身生長特點(diǎn)及傳統(tǒng)種植模式影響,枸杞連作障礙現(xiàn)象已成為限制枸杞產(chǎn)業(yè)可持續(xù)發(fā)展的主要因素之一,嚴(yán)重影響了枸杞園土地生產(chǎn)力和枸杞生產(chǎn)效益,導(dǎo)致農(nóng)民收入下降。通過多年的田間調(diào)查和實(shí)驗(yàn),納小凡等[6]研究發(fā)現(xiàn)種植超過十年的枸杞園開始出現(xiàn)枸杞子產(chǎn)量和品質(zhì)降低、病蟲害加重等現(xiàn)象。對連作地再植枸杞生長的分析也表明,連作能夠顯著抑制再植枸杞生長,證實(shí)了枸杞連作障礙現(xiàn)象的存在[7]。張俊華和鄭國琦[8]對不同種植年限枸杞園土壤線蟲群落的分析發(fā)現(xiàn),樹齡超過9年后枸杞園土壤線蟲群落的多樣性下降,且植物寄生性線蟲的相對比例明顯增加。通過對連作地再植枸杞根際土壤細(xì)菌群落的分析也表明,連作顯著抑制再植枸杞根際土壤細(xì)菌群落多樣性并導(dǎo)致細(xì)菌群落結(jié)構(gòu)發(fā)生顯著變化[7]。這些研究均表明長期傳統(tǒng)種植模式下,枸杞園土壤理化性質(zhì)及土壤生物環(huán)境的變化可能是導(dǎo)致枸杞連作障礙的主要原因,但其作用機(jī)制目前仍不清楚。
作為土壤微生物區(qū)系的重要成員,土壤真菌與其他微生物一起參與土壤生態(tài)系統(tǒng)中物質(zhì)循環(huán)和能量流動[9-10]。研究發(fā)現(xiàn),真菌群落能夠通過互惠共生、病害及營養(yǎng)物質(zhì)吸收和循環(huán)等方式影響陸生植物群落的多樣性和組成[11]。除此以外,土壤真菌還能夠通過固氮、植物激素分泌、病原菌防護(hù)和提高植物抗旱等生物學(xué)過程影響植物生理狀況[12-14]。因此研究土壤真菌群落多樣性及其演替規(guī)律,對于探索土壤可持續(xù)利用措施及保持土壤健康至關(guān)重要[14-15]。本研究基于前期對連作地再植枸杞根際土壤細(xì)菌群落的研究,進(jìn)一步利用實(shí)時(shí)熒光定量PCR 和高通量測序技術(shù),分析連作對再植枸杞根際/非根際土壤真菌群落豐度及多樣性的影響,為深入理解枸杞連作障礙機(jī)制提供新的理論信息。
試驗(yàn)地位于銀川市北郊20 km 的南梁農(nóng)場,該地區(qū)主要?dú)夂蝾愋蜑橹袦貛Ц珊禋夂?,年均氣溫?.9~8.6℃之間,晝夜溫差大,年均降雨量約為145.5 mm,無霜期在129~177 d 之間。樣地土壤類型為砂壤土。
在實(shí)驗(yàn)區(qū)域內(nèi)選取連續(xù)種植10年的枸杞田和無枸杞種植史的設(shè)施農(nóng)業(yè)用地,分別將其設(shè)置為連作及對照樣地。連作及對照枸杞園土壤基礎(chǔ)理化性質(zhì)如下:pH 分別為5.2 和6.4;電導(dǎo)率199.7 和1006 μS·cm-1;有機(jī)碳9.7 和12.9 g·kg-1;全氮0.9和0.8 g·kg-1;堿解氮46.7 和30.1 mg·kg-1;銨態(tài)氮13.9 和13.7 mg·kg-1;硝態(tài)氮30.3 和20.6 mg·kg-1;全磷0.9 和0.8 g·kg-1;有效磷48.9 和12.3 mg·kg-1;速效鉀260.7 和309.1 mg·kg-1。
挖去連作樣地枸杞,連作及對照樣地經(jīng)整理后于2012年4月種植一年生枸杞苗(寧杞1 號),行距3 m、株距1 m,樣方面積666 m2,枸杞田間管理與前期研究一致[6]。2016年8月采集寧夏枸杞根際、非根際土壤樣本(最后一茬枸杞果實(shí)采摘后)。實(shí)驗(yàn)田內(nèi)隨機(jī)設(shè)置3 個(gè)樣方(10 m×10 m),每個(gè)樣方內(nèi)隨機(jī)選取10 株枸杞,利用土鉆在距離枸杞莖20 cm 處采集0~20 cm 深度土壤;同時(shí)挖取枸杞根(直徑小于2 mm),抖去浮土后放入無菌離心管。合并同一樣方內(nèi)的根際土及非根際土,置于冰盒保存。非根際土過2 mm 篩去雜物并充分混勻,留取1 g鮮土于-20℃保存用于提取DNA;剩余土樣經(jīng)風(fēng)干后用于土壤理化性質(zhì)分析[6]。根際土壤采用納小凡等[6]的方法獲得。土壤DNA 按照Power SoilTMDNA提取試劑盒說明提取。所提DNA 經(jīng)瓊脂糖凝膠電泳檢測合格后分為三份,一份用于前序研究對 16S rDNA V4 區(qū)段的測序[6],一份用于本研究中真菌ITS2 區(qū)段的測序,一份保留至-20℃用于qRT-PCR測定16S rDNA 和18S rDNA 基因豐度。
采用兩步擴(kuò)增法對真菌的ITS2 區(qū)段進(jìn)行擴(kuò)增,第一步擴(kuò)增法引物序列分別為 5′-TCCGTAGGTGA ACCTGCGC-3′和 5′-TCCTCCGCTTATTGATATGC- 3′。擴(kuò)增產(chǎn)物經(jīng)稀釋后用于第二步擴(kuò)增反應(yīng)的模板,引物序列分別為 5′-GCATCGATGAAGAACGCAG C-3′和5′-TCCTCCGCTTATTGATATGC-3′。擴(kuò)增條件為:98℃,預(yù)變性1 min;98℃,變性10 s;50℃,退火30 s;72℃,延伸1 min;30 個(gè)循環(huán)。72℃,延伸5 min。PCR 擴(kuò)增完成后,應(yīng)用 Illumina HiSeq 2500 基因組分析平臺對真菌基因ITS2 區(qū)段進(jìn)行測序。
所提土壤基因組DNA 中16S rDNA(細(xì)菌)和18S rDNA(真菌)基因豐度利用QuantStudio?5實(shí)時(shí)定量熒光PCR 系統(tǒng)(Thermo Fisher Scientific,USA)進(jìn)行。所用細(xì)菌引物為Eub338(ACTCCTAC GGGAGGCAG CAG)和Eub518(ATTACCGCGGCT GG),真菌引物為Fu18s1(GGAAACTCAGGTCCA GA)和Nu-SSU-1536(ATTGCAATGCYCTATCCC CA)。PCR 反應(yīng)體系體積為20 μL,其中包括SYBRPremix Ex Taq(Takara,大連,中國)10 μL、上下游引物各1 μL(10μmol·L-1)、DNA 模板2 μL、ROX 0.4 μL、無菌去離子水5.6 μL。細(xì)菌用兩步法進(jìn)行擴(kuò)增,擴(kuò)增程序?yàn)椋?5℃預(yù)變性30 s,在95℃變性5 s, 60℃退火和延伸34 s。真菌用三步法進(jìn)行擴(kuò)增,擴(kuò)增程序如下:95℃預(yù)變性30 s,95℃變性5 s,50℃退火30 s,72℃延伸30 s,循環(huán)40 次。細(xì)菌和真菌溶解曲線按照QuantStudio 5 實(shí)時(shí)熒光定量檢測系統(tǒng)默認(rèn)程序分析:95℃變性15 s,60℃復(fù)性60 s,95℃復(fù)性15 s。將含16S rDNA 和18S rDNA 基因片段的重組質(zhì)粒DNA 梯度稀釋,并分別以此為模板制備16S rDNA 和18S rDNA 基因擴(kuò)增標(biāo)準(zhǔn)曲線,計(jì)算土壤DNA 樣本中相應(yīng)基因的絕對含量。
不同樣品間數(shù)據(jù)的顯著性分析由SPSS(22.0)中的單因素方差分析完成(Turkey test,n=3)。真菌營養(yǎng)型預(yù)測利用FUNGuild數(shù)據(jù)庫進(jìn)行對比檢索[16]。土壤理化因子與土壤真菌群落結(jié)構(gòu)變化的關(guān)系利用Canoco(5.0)中基于距離的冗余分析(db-RDA) 完成。作圖使用Origin 9.0 軟件。
實(shí)時(shí)熒光定量PCR 結(jié)果表明,連作樣地根際及非根際土壤細(xì)菌和真菌豐度分別顯著高于對照樣地根際及非根際土壤(P<0.05,表1)。其中,連作樣地再植枸杞根際細(xì)菌和真菌豐度分別為對照樣地再植枸杞根際的4.8 倍和7.3 倍,而非根際土壤則分別為對照的2.1 倍和5.9 倍(表1)。進(jìn)一步分析發(fā)現(xiàn),對照樣地根際及非根際土壤的細(xì)菌/真菌比均高于連作樣地,其比例分別為對照樣地的1.5 倍(P>0.05)和2.7 倍(P<0.05,表1)。
表1 實(shí)時(shí)熒光定量PCR 法分析枸杞根際/非根際細(xì)菌和真菌豐度變化 Table1 Abundances of soil bacteria and fungi in the bulk soil and rhizosphere of L.barbarum detected by real-time PCR
對真菌ITS2 區(qū)段的高通量測序結(jié)果表明,雖然非根際土壤真菌群落的豐度和均一度指數(shù)均高于枸杞根際土壤,但與對照樣地相比,連作并不影響再植枸杞根際及非根際土壤真菌群落的α多樣性指數(shù)(P>0.05,圖1a、圖1b)。多維尺度分析(NMDS)表明,連作樣地與對照樣地再植寧夏枸杞根際及非根際土壤的真菌群落均存在明顯差異。但與對照樣地相比,連作樣地再植寧夏枸杞根際真菌群落結(jié)構(gòu)與非根際土壤間的差異明顯降低(圖1c)。
測序結(jié)果表明,所測定土壤樣品中的優(yōu)勢真菌門分別為子囊菌門(Ascomycota)、擔(dān)子菌門(Basidiomycota)、接合菌門(Zygomycota)、壺菌門(Chytridiomycota)、球囊菌門(Glomeromycota),這些真菌占各樣品總序列數(shù)的99.1%~99.9%(表2)。統(tǒng)計(jì)分析發(fā)現(xiàn),接合菌門在連作地枸杞根際土壤中的相對豐度較對照樣地顯著增加(P<0.05,表2),但球囊菌門的相對豐度則顯著降低(P<0.05,表2)。
圖1 再植枸杞根際/非根際土壤真菌群落多樣性分析 Fig.1 Diversity of the soil fungal communities in the bulk soil and rhizosphere of replanted L.barbarum
表2 再植枸杞根際/非根際優(yōu)勢真菌群落相對豐度變化 Table2 Relative abundances of the dominant fungal phyla in bulk soils and rhizosphere of replanted L.barbarum
進(jìn)一步對相對豐度排名前30 的真菌屬分析發(fā)現(xiàn),連作顯著促進(jìn)了再植枸杞根際土壤中毛霉屬和Gibellulopsis的相對豐度(P<0.05,表3)。此外,連作樣地再植枸杞根際和/或非根際土壤中鐮刀菌屬、葡柄霉屬和粉紅螺旋聚孢霉屬的相對豐度均高于對照樣地,但結(jié)果不顯著(P>0.05,表3)。
表3 再植枸杞根際/非根際優(yōu)勢真菌屬豐度變化 Table3 Relative abundance of dominant fungal genera in the bulk and rhizosphere of replanted L.barbarum
續(xù)表
將測序結(jié)果與FUNGuild 的真菌數(shù)據(jù)庫進(jìn)行比對,在測定的所有土壤樣本中主要發(fā)現(xiàn)六種功能真菌(表4)。統(tǒng)計(jì)分析表明,連作地再植枸杞根際土壤中叢枝菌根真菌的豐度較對照樣地再植枸杞根際土壤顯著降低(P<0.05),其他五種功能真菌的相對豐度則無顯著變化(表4)。
表4 連作及對照條件下再植枸杞根際/非根際不同功能真菌的相對豐度變化 Table4 Relative abundances of soil fungi in the bulk soils and rhizosphere of replanted L.barbarum relative to functional group and treatment
相關(guān)性分析表明,枸杞根際土壤中接合菌門的相對豐度與土壤電導(dǎo)率和硝態(tài)氮顯著負(fù)相關(guān),而與有效磷含量顯著正相關(guān);壺菌門的相對豐度與土壤全磷含量顯著負(fù)相關(guān);球囊菌門的相對豐度與硝態(tài)氮含量顯著正相關(guān),而與有效磷含量顯著負(fù)相關(guān)(P<0.05,表5)。非根際土壤中優(yōu)勢真菌門的相對豐度與土壤因子間均無顯著相關(guān)關(guān)系(P>0.05,表5)。從真菌功能分類看,枸杞根際土壤中植物病原菌的相對豐度與土壤全磷含量顯著負(fù)相關(guān);腐生真菌的相對豐度與土壤電導(dǎo)率及硝態(tài)氮含量顯著負(fù)相關(guān);叢枝菌根真菌的相對豐度與硝態(tài)氮含量顯著正相關(guān),而與有效磷含量顯著負(fù)相關(guān)。而在非根際土壤中,植物病原菌的相對豐度與土壤有機(jī)碳含量顯著正相關(guān)(P<0.05,表5)。
進(jìn)一步利用基于距離的冗余分析(db-RDA)發(fā)現(xiàn),土壤有效磷含量、電導(dǎo)率、速效氮含量及pH分別解釋了寧夏枸杞根際土壤真菌群落變化的33.7%、30.0%、27.1%和26.7%(P<0.05,表6)。土壤有效磷和硝態(tài)氮含量分別解釋了寧夏枸杞非根際土壤真菌群落變化的34.1%和29.3%(P<0.05),其他土壤因子均無顯著影響(表6)。
實(shí)時(shí)熒光定量PCR 的結(jié)果證實(shí),連作地再植枸杞根際和非根際土壤中細(xì)菌和真菌的豐度均顯著高于對照樣地(表1),這一結(jié)果說明傳統(tǒng)模式下長期種植枸杞能夠顯著促進(jìn)土壤微生物種群的繁殖。導(dǎo)致這一結(jié)果的原因可能有以下三點(diǎn):首先,長期施用肥料及枸杞落葉返還土壤能夠通過提供營養(yǎng)和碳源、調(diào)整土壤C/N 比及改變土壤微環(huán)境等方式促進(jìn)土壤微生物群落的生長和繁殖[17]。其次,Bouskill等[18]研究發(fā)現(xiàn),前期預(yù)處理能夠增強(qiáng)土壤微生物群落對后期相同類型環(huán)境脅迫的耐受性。經(jīng)過長期的定向選擇,枸杞連作地土壤微生物已適應(yīng)枸杞根際環(huán)境,進(jìn)而在再植枸杞根際土壤大量定植。此外,連作地土壤可能通過影響再植枸杞生長及其生理狀態(tài)[6],調(diào)整枸杞根系與土壤微生物間的互作關(guān)系,導(dǎo)致枸杞根系對根際土壤微生物的控制力降低而使根際微生物過度繁殖。對再植枸杞根際/非根際真菌群落結(jié)構(gòu)的分析發(fā)現(xiàn),連作地再植枸杞根際真菌群落與非根際真菌群落結(jié)構(gòu)較為相似(圖1C),證實(shí)連作可能弱化了再植枸杞根系對土壤真菌的選擇作用。然而就連作如何影響枸杞與土壤微生物間的相互作用仍需深入研究。
由于土壤細(xì)菌和真菌群落在底物利用、功能等方面的差異,決定了這兩種微生物群落在應(yīng)對環(huán)境變化時(shí)表現(xiàn)出不同的響應(yīng)機(jī)制[19]。本研究結(jié)果也表明,連作對再植枸杞根際/非根際土壤真菌群落的促進(jìn)作用高于細(xì)菌,導(dǎo)致連作地再植枸杞根際和非根際土壤中細(xì)菌/真菌的比例均明顯低于對照樣地(表1),使再植枸杞根際土壤微生物環(huán)境偏向于真菌型。近期,Grosso 等[20]研究發(fā)現(xiàn)土壤pH 和植物凋落物的C/N 是影響土壤細(xì)菌和真菌群落平衡的主要因素。其中,低pH 和高C/N 均有利于真菌生長。因此,長期施用氮肥導(dǎo)致的土壤pH 降低[6],可能是連作影響再植枸杞根際和非根際土壤細(xì)菌/真菌比例平衡的主要原因之一。
叢枝菌根真菌(Arbuscular mycorrhizal fungi,AMF)是土壤中分布極其廣泛的一類真菌,它能夠與約75%的陸生植物形成互利共生的聯(lián)合體[21],因其具有幫助植物吸收氮、磷等元素、提高植物對土傳病害和干旱等生物及非生物脅迫的耐受性等功能而受到廣泛關(guān)注[22-23]。目前研究發(fā)現(xiàn),叢枝菌根真菌主要以球囊菌門為主,因此有研究直接將球囊菌門視為叢枝菌根真菌[21,24]。本研究的測序及真菌功能預(yù)測結(jié)果表明,球囊菌門和叢枝菌根真菌的相對豐度在不同土壤樣本中均保持一致(表2和表4)。進(jìn)一步分析發(fā)現(xiàn),連作地再植枸杞根際土壤中叢枝菌根真菌的相對豐度較對照樣地顯著降低(表4),且其豐度與硝態(tài)氮含量顯著正相關(guān),而與有效磷含量則顯著負(fù)相關(guān)(表5)。這些結(jié)果說明,長期施肥可能影響了叢枝菌根真菌與再植枸杞根系間的共生關(guān)系。由此推測,叢枝菌根真菌豐度的減少可能進(jìn)一步影響了再植枸杞對生物及非生物脅迫的抗性,從而導(dǎo)致再植枸杞生長減緩。
?
表6 土壤理化因子對寧夏枸杞根際/非根際真菌群落結(jié)構(gòu)的影響 Table6 Effects of soil properties on composition of the soil fungal communities in the bulk soil and rhizosphere of L.barbarum
前期對不同種植年限枸杞根際土壤微生物群落 的分析發(fā)現(xiàn),鐮刀菌屬、李維菌屬及枝孢屬等已知植物病原菌的相對豐度在不同樹齡枸杞間均發(fā)生顯著變化[7]。在本研究中,F(xiàn)UNGuild 數(shù)據(jù)庫檢索結(jié)果顯示植物致病菌的總體相對豐度在不同樣品中無顯著差異(表4)。然而,就單一植物病原菌而言,連作地再植枸杞根際土壤中已知植物病原菌鐮刀菌屬、擬青霉屬、葡柄霉屬和粉紅螺旋聚孢霉屬的相對豐度均高于對照樣地(表3)。雖然目前就這些植物病原菌是否導(dǎo)致枸杞產(chǎn)生病害尚不清楚,但這些病原真菌相對豐度的增加說明連作地再植枸杞根際土壤微生物環(huán)境較對照樣地更差。此外,這一結(jié)果也間接證明連作可能干擾了再植枸杞根系對植物病原菌的抑制作用。
前人研究表明長期施用化肥所導(dǎo)致的土壤酸化、土壤次生鹽堿化等效應(yīng)均能夠直接或間接影響土壤微生物群落的結(jié)構(gòu)和組成[21,25]。前序基于相同條件下對細(xì)菌群落的研究發(fā)現(xiàn),土壤pH 和有效磷含量是驅(qū)動枸杞非根際細(xì)菌群落結(jié)構(gòu)變化的主要因素[6]。本研究進(jìn)一步表明,土壤pH、電導(dǎo)率、硝態(tài)氮和有效磷含量是解釋再植枸杞根際真菌群落的主要因素,其中硝態(tài)氮和有效磷含量也是驅(qū)動非根際土壤真菌群落的主要土壤因子(表6)。這些結(jié)果與前人研究結(jié)果一致[21]。在實(shí)際生產(chǎn)中,長期施用氮肥是導(dǎo)致枸杞園土壤酸化的主要原因,并且氮肥和磷肥的長期使用又能夠?qū)е碌?、磷元素過量。因此在傳統(tǒng)種植模式下,大量、長期使用化肥所導(dǎo)致的土壤理化性質(zhì)改變是驅(qū)動枸杞園土壤細(xì)菌/真菌群落及再植枸杞根際真菌群落結(jié)構(gòu)演替的主要因子。鑒于土壤細(xì)菌和真菌群落在多數(shù)土壤生態(tài)過程中具有重要作用,由此可以推測,長期使用化肥對土壤微生物群落的影響可能是導(dǎo)致枸杞連作障礙產(chǎn)生的主要原因之一。
此外,在長期種植過程中枸杞凋落物及根系分泌物的不斷輸入也可能導(dǎo)致枸杞園土壤微生物群落結(jié)構(gòu)發(fā)生改變。本研究也觀察到TOC 含量與植物病原真菌的相對豐度存在顯著正相關(guān)關(guān)系(表5),但對枸杞凋落物如何驅(qū)動枸杞園土壤微生物群落演替以及長期施肥如何影響土壤微生物群落對枸杞凋落物的利用和降解等重要過程仍不清楚。探索這些土壤生物學(xué)過程對于理解枸杞連作障礙形成機(jī)制及改進(jìn)枸杞種植技術(shù)具有重要理論價(jià)值。
長期連作顯著促進(jìn)再植枸杞根際及非根際土壤細(xì)菌和真菌的群落豐度;且連作對真菌的促進(jìn)作用強(qiáng)于細(xì)菌,導(dǎo)致枸杞根際及非根際土壤微生物環(huán)境偏向于真菌型。測序結(jié)果表明,再植枸杞根際土壤的優(yōu)勢真菌分別為子囊菌門、擔(dān)子菌門、接合菌門、壺菌門和球囊菌門,其中連作地再植枸杞根際土壤中球囊菌門(叢枝菌根真菌)的相對豐度較對照樣地顯著降低,而接合菌門的相對豐度則顯著增加。db-RDA 分析發(fā)現(xiàn),土壤pH、電導(dǎo)率、硝態(tài)氮和有效磷是驅(qū)動再植枸杞根際土壤真菌群落結(jié)構(gòu)變化的主要因素,而硝態(tài)氮和有效磷含量是影響非根際土壤真菌群落結(jié)構(gòu)變化的主要因子。