不同樹齡桑根際土壤原生生物群落組成多樣性

王子夜 王志剛 閻愛華

（1. 河北農(nóng)業(yè)大學(xué)林學(xué)院，保定 071000；2. 河北省林木種質(zhì)資源和森林保護(hù)重點實驗室，保定 071000；3. 河北省城市森林健康技術(shù)創(chuàng)新中心，保定 071000）

原生動物種類繁多，廣泛存在于海洋、森林、農(nóng)田等各生態(tài)系統(tǒng)，通過土壤初級生產(chǎn)、菌物捕食、動植物餌料和病原物等多層次作用調(diào)節(jié)土壤生產(chǎn)力，控制微生物群落組裝，影響植物生長，其種群落結(jié)構(gòu)的變化往往會影響到土壤-植物系統(tǒng)微食物網(wǎng)及生態(tài)系統(tǒng)結(jié)構(gòu)的變動，對于維持土壤能量養(yǎng)分循環(huán)的強(qiáng)度和穩(wěn)定性尤其重要［1-3］。研究表明，溫度、地理環(huán)境、海拔、土壤性質(zhì)甚至土壤大型動物都會對原生生物多樣性產(chǎn)生影響，而地上部植物功能類群主導(dǎo)了原生生物種群多樣性建成［4-5］。但當(dāng)前土壤微生物組研究的重點集中在細(xì)菌、真菌和單類群原生生物類別的互作關(guān)系，探討完整的土壤原生生物組在土壤生態(tài)系統(tǒng)中的種群分布、關(guān)鍵功能的研究仍很缺乏［6-7］。

林木生長周期長，營林措施頻繁，為應(yīng)對不同環(huán)境挑戰(zhàn)，發(fā)展出獨特的生理適應(yīng)性與緩解策略，土壤作為生態(tài)環(huán)境的重要組成部分，了解原生生物群落在土壤中產(chǎn)生和維持的機(jī)制對于理解其生物功能，維持林木-土壤良性循環(huán)，促進(jìn)林木生產(chǎn)力提高至關(guān)重要［8］。樹齡不僅影響樹體形態(tài)建成、營養(yǎng)獲取和免疫抗逆能力，還與土壤理化結(jié)構(gòu)和和水肥營養(yǎng)狀況相關(guān)，隨著樹齡增長，林木共生體（symbiont）內(nèi)生和外生真菌、細(xì)菌的種群和功能多樣性在時間空間都上都表現(xiàn)出差異［9-11］，但目前作為林木根際土壤微生態(tài)三級營養(yǎng)層次的關(guān)鍵環(huán)節(jié)，樹齡對土壤原生生物種群多樣性和群落構(gòu)建的影響仍缺乏系統(tǒng)研究［12-13］。

桑（Morus albaL.）是中國傳統(tǒng)藥食兼用的資源植物，也被廣泛用于土壤生態(tài)修復(fù)，具有良好的生態(tài)效益和社會效益。河北省桑樹種植歷史悠久，東南部古桑樹群為研究桑土壤原生生物群落多樣性演化提供了寶貴的資源。本研究按照空間換時間法，以同時期同桑園同品種的根際土壤樣品替換傳統(tǒng)不同時間采集不同樹齡根際土壤樣品為試驗材料，通過Illumina高通量測序和生物信息學(xué)分析，探討古桑園不同樹齡桑（10a、80a和200a）根際土壤原生生物種類多樣性差異和功能類群演變，明確影響原生生物群落構(gòu)建和穩(wěn)定性的不同環(huán)境因子和生物因子，為進(jìn)一步探討原生生物種群分化機(jī)制和生態(tài)服務(wù)功能提供理論依據(jù)。

1 材料與方法

1.1 材料

1.1.1 試驗地概況 樣品采自河北農(nóng)業(yè)大學(xué)營子定位站，地理位置坐標(biāo)為116°2'5''E，37°54'13''N。屬暖溫帶濕潤大陸性季風(fēng)氣候區(qū)，又屬冀南暖濕輕干旱區(qū)。年平均氣溫12.7℃，無霜期187 d，年均日照2.78 h，年均降水543 mm。土壤由古河流沖積形成，層理清晰，土層深厚，疏松涵潤，園區(qū)統(tǒng)一常規(guī)管理。

1.1.2 樣品的采集及處理 本研究以同桑園的幼齡樹（10a）、中齡樹（80a）和老齡樹（200a）根際土壤為材料，每株4個方位分別取樣10-30 cm土層根際土壤采集，先抖落不含根系的土壤，然后選取松散黏附在根系表面1-4 mm范圍內(nèi)的土壤為根際土的土壤裝入無菌袋內(nèi)，每株土樣充分混合為一個樣本［14］，取樣后樣品于-80℃的冷凍室中冷凍保存。10a和80a隨機(jī)取樣5 株，由于200a樹齡為珍稀古樹故僅取1 株，用于土壤原生生物高通量測序及土壤理化性質(zhì)測定。

1.2 方法

1.2.1 土壤理化性質(zhì)測定 有機(jī)質(zhì)：以硫酸亞鐵為標(biāo)準(zhǔn)溶液，將容量法的滴定手段改為光度測定手段，進(jìn)行土壤有機(jī)質(zhì)585 nm處比色的測定［15］。

堿解氮：在擴(kuò)散皿中，用堿性溶液水解土壤，使易水解態(tài)氮（潛在有效氮）堿解轉(zhuǎn)化為NH3，NH3擴(kuò)散后為硼酸吸收。硼酸吸收液中的NH3再用標(biāo)準(zhǔn)酸滴定，然后計算土壤中速效氮的含量［14］。

速效磷、速效鉀、氯離子和鈉離子分別用鉬銻抗比色法測定、火焰光度法測定、硝酸銀滴定法和ICP-AES［16］。

全鹽含量采用電導(dǎo)法（LY/T 1251-1999）［17］。

pH值：稱土10 g，放入50 mL燒杯中，加入蒸餾水25 mL用攪拌器攪拌1 min，使土體充分散開，放置半小時然后用酸度計測定［15］。

1.2.2 土壤原生生物基因組DNA提取 基因組DNA提取方法參考Edwards等［18］。利用1%瓊脂糖凝膠檢測DNA的長度和完整性，利用NanoDrop檢測DNA的濃度和純度。

1.2.3 DNA擴(kuò)增及測序 PCR擴(kuò)增和高通量測序：通過Illumina MiSeq平臺，利用原生動物18S通用引物V9：上游引物序列1389F（5'-TTGTACACACCGCCC-3'）和下游引物序列1510R（5'-CCTTCYGCAGGTTCACCTAC-3'）進(jìn)行序列測定，目的片段長度為120 bp［19］。PCR擴(kuò) 增體系（20 μL）：5×FastPfu Buffer 2 μL，2.5 mmol/LdNTPs 2 μL，正 反 引 物（5 μmol/L）各0.8 μL，F(xiàn)astPfu Polymerase（2.5 U/μL）0.4 μL，BSA 0.2 μL，模板10 ng，加ddH2O補(bǔ)充至20 μL。反應(yīng)程序為95℃預(yù)變性3 min后，第一輪擴(kuò)增27個循環(huán)（95℃變性30 s，55℃復(fù)性30 s，72℃延伸45 s），最后72℃延伸10 min，第二輪擴(kuò)展13個循環(huán)，其余條件均不變，10℃保溫。擴(kuò)增后PCR產(chǎn)物使用濃度1.8%的瓊脂糖進(jìn)行電泳檢測（北京博美富鑫科技有限公司）對完整性進(jìn)行檢測。質(zhì)檢合格樣本構(gòu)建目標(biāo)區(qū)域PCR體系，之后進(jìn)行產(chǎn)物純化和PCR擴(kuò)增，再進(jìn)行切膠回收。對于構(gòu)建好的文庫使用采用Illumina公司的TruSeq Nano DNA LT Library Prep Kit制備測序文庫（上海派森諾生物科技有限公司）。

1.2.4 數(shù)據(jù)分析 （1）基于有效數(shù)據(jù)進(jìn)行OTUs（operational taxonomic units）聚類和物種分類分析，使用 Trimmomatic v0.33軟件，對測序得到的 Raw Reads進(jìn)行過濾；利用cutadapt 1.9.1軟件進(jìn)行引物序列的識別與去除，得到不包含引物序列的Clean Reads；使用Usearch v10軟件雙端序列拼接，通過overlap對每個樣品的Clean Reads 進(jìn)行拼接，根據(jù)不同區(qū)域的長度范圍對拼接后數(shù)據(jù)進(jìn)行長度過濾；利用 UCHIME v4.2軟件，鑒定并去除嵌合體序列，得到最終有效數(shù)據(jù)（effective reads）。并將OTU和物種注釋結(jié)合，人工剔除測序結(jié)果中植物、真菌和動物等序列，得到每個樣品中原生生物OTUs和物種注釋等基本分析結(jié)果。

（2）應(yīng)用Mothur軟件對根際土壤中原生生物群落在各個分類水平上的豐度劃分alpha多樣性指數(shù)（Chao1、Shannon、Simpson 指數(shù)以及覆蓋度），beta多樣性通過主成分分析（principal component analysis，PCA）等統(tǒng)計分析，獲得群落樣本的分布特征［20］。

（3）采用Excel 2003進(jìn)行數(shù)據(jù)處理，采用SPSS 17.0進(jìn)行Duncan氏差異顯著性檢驗和Pearson相關(guān)性分析。計算原生生物與真菌、細(xì)菌（為本實驗室前期獲得數(shù)據(jù)，尚未公開發(fā)表）之間的Spearman等級相關(guān)系數(shù)，對其中rho>0.6且P<0.01的相關(guān)優(yōu)勢屬構(gòu)建關(guān)聯(lián)網(wǎng)絡(luò)，并導(dǎo)入Cytoscape（http://www.cytoscape.org/）軟件進(jìn)行可視化［21］。

2 結(jié)果

2.1 多樣性分析

2.1.1 α多樣性分析 200 a原生動物獲得最多的有效序列數(shù)，但產(chǎn)生的OTUs最少?；贠TU數(shù)及相對豐度，結(jié)果表明果桑不同樹齡多樣性隨樹齡增大而降低，不同樹齡Shannon、和Chao1指數(shù)表明原生生物隨樹齡增大其多樣性和豐富度降低（表1）。

表1 土壤原生生物數(shù)據(jù)統(tǒng)計及alpha多樣性Table 1 Soil protist data statistics and alpha diversity

2.1.2 β多樣性分析 基于物種豐度進(jìn)行主成分分析（principal component analysis，PCA），在OTU水平上對β多樣性進(jìn)行評估，分析不同樣本組成結(jié)構(gòu)的相似性和差異關(guān)系。基于Bray-Curtis距離算法，構(gòu)建不同樹齡原生生物的層級聚類樹（圖1），各樣點聚集表明樹齡對原生生物群落β多樣性影響較小。

圖1 根際土壤原生生物主成分分析Fig. 1 Principal component analysis（PCA）of rhizosphere soil protist

2.2 門和屬水平土壤原生生物群落豐富度分析

不同樹齡根部土壤原生生物在門水平上Chlorophyta、Kinetoplastida、Bacillariophyta、Rotaliida、Oomycetes、Petalomonadida、Blastocystidae被 定 義為優(yōu)勢門（平均相對豐度>1）。Bacillariophyta相對豐度隨樹齡的增大降低，綠藻門Chlorophyta、Petalomonadida和Blastocystidae相對豐度升高。由圖2可知原生生物群落優(yōu)勢門種類基本一致，但它們相對豐度存在差異。

圖2 門水平物種統(tǒng)計圖Fig. 2 Statistics at phylum level

在屬水平上，各樹齡中Poteriospumella、Vermamoeba、Paracercomonas、吻滴蟲屬Rhynchomonas、Synura、Neobodo、Spumella、吻波豆蟲Rhynchobodo、卵囊藻Oocystis、Petalomonas、棘阿米巴屬Acanthamoeba和Blastocystis為優(yōu)勢屬（平均相對豐度>1）（圖3）。10a中獨有屬18 個，80a獨有12 個，200a獨有2 個。利用單因素ANOVA檢驗發(fā)現(xiàn)在10a-80a、10a-200a和80a-200a中 相 對 豐 度 差 異顯著的屬分別為Paracercomonas、Petalomonas和Acanthamoeba、Rhynchomonas，結(jié)果表明土壤原生生物群落組成和相對豐度具有明顯樹齡變異（圖4）。

圖3 屬水平物種統(tǒng)計圖Fig. 3 Statistics at genus level

圖4 獨有和共有原生生物種類數(shù)Fig. 4 Number of unique and common protist species

2.3 基于影響群落組成土壤環(huán)境因素的分析

根據(jù)差異顯著分析結(jié)果表明：各樹齡速效鉀含量差異顯著，80a vs 200a間速效磷、氯離子含量差異顯著，pH值10a vs 80a差異極顯著，堿解氮含量、有機(jī)質(zhì)含量、全鹽量和鈉離子隨樹齡變化不顯著（圖5）

圖5 土壤理化指標(biāo)差異顯著性分析Fig. 5 Significance analysis of differences in soil physical and chemical indicators

根際土壤原生生物群落多樣性、豐富度（Chao指數(shù)、Simpson指數(shù)、Shannon指數(shù)）、均勻度與樹齡和氯離子含量極顯著負(fù)相關(guān)，Simpson指數(shù)與全鹽量顯著負(fù)相關(guān)，覆蓋度與樹齡、速效鉀和堿解氮含量顯著正相關(guān)；同時各多樣性指數(shù)與速效磷、有機(jī)質(zhì)、鈉離子含量和pH值相關(guān)但不顯著（表2）。

表2 土壤理化性狀與土壤原生生物多樣性指數(shù)的相關(guān)性分析Table 2 Correlation analysis between soil physical and chemical properties and biodiversity indexes of soil protist

2.4 基于細(xì)菌、真菌的關(guān)聯(lián)分析

使用Mothur軟件對各樣本中原生生物與細(xì)菌、真菌在rho > 0.6且P值< 0.01相關(guān)的屬構(gòu)建關(guān)聯(lián)網(wǎng)絡(luò)，結(jié)果表明關(guān)聯(lián)最多且互為共生關(guān)系的原生生物為破囊壺菌Aplanochytrium和微藻Chrysoxys（分別有22 個關(guān)聯(lián)屬），其次是Chlamydaster（19 個關(guān)聯(lián)屬）、鞭毛蟲Dimastigella（18 個關(guān)聯(lián)屬）和Goniomonas（17 個關(guān)聯(lián)屬）；與原生生物關(guān)聯(lián)最多的是Aplanochytrium（7 個關(guān)聯(lián)屬）；Chlamydaster、Chrysoxys、Dimastigella、Goniomonas、Rhogostoma、Rosculus和Uroleptus互為共生關(guān)系。絲壺菌屬Hyphochytrium與細(xì)菌關(guān)聯(lián)最多（4 個關(guān)聯(lián)屬），分別是動球菌屬Kineococcus、水生角質(zhì)菌屬Pilimelia、Pseudorhodoferax和疣孢菌屬Verrucosispora。與真菌關(guān)聯(lián)最多（15 個關(guān)聯(lián)屬），分別是Flammulina、腹孢 菌 屬Gastrosporium、Genolevuria、Kalaharituber、Kondoa、Mycena、Perenniporia、 紅 酵 母 屬Rhodotorula、莢 孢 腔 菌 屬Sporormia、匍 柄 霉屬Stemphylium、Termitomyces、Trichophaeopsis、Tuber、Xerochrysium。在 各 原 生 生 物 集 團(tuán) 中 與大量子囊菌門、擔(dān)子菌門、球囊菌門等真菌共生。Aplanochytrium、Chlamydaster、Dimastigella、Goniomonas、Eocercomonas和絲壺菌屬僅存在于10齡根部土壤中，Chrysoxys僅存在于80 齡根部土壤中。Pervetustus與細(xì)菌Ramicandelaber為競爭關(guān)系，Pervetustus在10 齡到80 齡相對豐度上升到200 齡消失，而Ramicandelaber在10 齡到80 齡相對豐度降低，200 齡豐度上升（圖6）。

圖6 原生生物與真菌、細(xì)菌網(wǎng)絡(luò)關(guān)聯(lián)圖Fig. 6 Correlation diagram of protis，fungi and bacteria network

3 討論

原生生物對各種土壤-植物系統(tǒng)表現(xiàn)出積極的動態(tài)響應(yīng)和生物地理屬性［22］，對森林、草原、稻田的研究都表明，植被類型與土壤理化性質(zhì)對土壤生物區(qū)系的影響大于植被的豐富度，且與微生物群落組成具有很強(qiáng)的相關(guān)性［4，23］。Ceja等［24］對開關(guān)草（Panicum virgatum）的研究表明，隨著時間的推移，植物的生理學(xué)指標(biāo)水平發(fā)生變化，而根際土壤原生動物多樣性減少。樹木生命周期長，個體空間更大，隨著時間積累，單個樹種對土壤中的許多生化和生態(tài)過程的驅(qū)動顯著強(qiáng)于多樹種，并且表現(xiàn)出樹種特異性［25］。前人研究表明全球不同生態(tài)區(qū)土壤內(nèi)Cercozoa和Ciliophor等捕食者（consumers）原生生物類群占支配地位［26］，但時間因素是否驅(qū)動土壤原生生物種類多樣性和群落結(jié)構(gòu)分化仍然未知。本研究利用空間換時間法從較長時間尺度上探討原生生物群落的結(jié)構(gòu)和組成動態(tài)變化，發(fā)現(xiàn)土壤原生生物群落組成和相對豐度具有明顯樹齡變異，不同樹齡根際原生生物的α多樣性和群落結(jié)構(gòu)存在明顯的差異，原生生物豐富度和多樣性隨樹齡增大降低?？赡苁巧ＩL盛期樹體營養(yǎng)充足，促進(jìn)了原生生物的群落多樣性和豐富度的增長；幼樹期或老齡期競爭養(yǎng)分和自身代謝能力較低，限制原生生物的生長及繁殖從而使多樣性和豐富度較低。在門水平上，在對各樹齡種類基本一致，但相對豐度存在一定差異，隨樹齡增長，提高了綠藻門等光能營養(yǎng)型原生生物的相對豐度，降低了硅藻門Bacillariophyta的相對豐度。在屬水平上，隨樹齡的增加，差異種群減少，10a中獨有屬18 個，80a獨有12 個，200a獨有2 個，Petalomonas和Acanthamoeba等捕食者原生動物的種類多樣性和豐度變異增強(qiáng)。

研究發(fā)現(xiàn)土壤pH值是原生生物多樣性、群落結(jié)構(gòu)和功能組成變化的驅(qū)動力，自養(yǎng)類原生生物在較低酸堿度條件下多樣性下降［27-28］。同時，原生生物是對施氮肥最敏感的微生物類群［29］，在堿性潮土中施肥可以顯著降低捕食性原生生物的相對豐 度［30］。本 研究10a根際土 壤pH值高于80a和200a，這與黃晶淼和王靜等對蘋果樹齡對土壤pH值的研究趨勢一致［31-32］，同時，速效鉀含量差異顯著，80a vs 200a間速效磷、氯離子含量差異顯著；氯離子含量顯著影響了原生生物群落多樣性，速效鉀和堿解氮含量的增長也顯著促進(jìn)了群落覆蓋度的增大。速效鉀和氯離子可調(diào)節(jié)土壤酸度，也證實了非生物參數(shù)是原生生物群落重要性的決定因素［33-34］。同 時，10a vs 80a、10a vs 200a和80a vs 200a中相對豐度差異顯著的屬也分別為動物性或腐生性營養(yǎng)型（Consumer）的Paracercomonas、溝滴蟲屬Petalomonas和棘阿米巴屬Acanthamoeba、吻滴蟲屬Rhynchomonas，群落多樣性也隨樹齡增大降低。這可能是樹齡影響了土壤主要營養(yǎng)和理化指標(biāo)，不同類群的原生生物對根際土壤酸堿度和營養(yǎng)狀況適應(yīng)性不同，捕食者群體或許可以作為土壤狀況的指示性生物［5］。

在土壤微生物三級營養(yǎng)結(jié)構(gòu)中，層次群落由自上而下（捕食者）和自下而上（營養(yǎng)）綜合效應(yīng)決定［35］。原生動物作為土壤微生物組的中心樞紐，連接著不同的細(xì)菌和真菌種群，處于土壤微生態(tài)網(wǎng)絡(luò)對樹齡動態(tài)反應(yīng)的關(guān)鍵位置，在土壤管理和調(diào)控細(xì)菌和真菌數(shù)量過程中發(fā)揮了杠桿作用［36-38］。例如異養(yǎng)型原生生物破囊壺菌Aplanochytrium和光能自養(yǎng)型卵囊菌屬Oocystis是各樹齡根際土壤中的優(yōu)勢菌屬，他們與根際真菌、細(xì)菌和原生生物之間呈現(xiàn)出復(fù)雜的互作關(guān)系：破囊壺菌在和微藻Chrysoxys之間不僅存在共生關(guān)系，還共同與棘尾蟲屬Uroleptus等捕食性原生生物存在協(xié)同共生關(guān)系，與病原細(xì)菌氣單胞菌屬Aeromonas和大量植物病原真菌如假尾孢屬Pseudocercospora、新殼二孢Neoascochyta和絲核菌屬Rhizoctonia存在正向調(diào)控；而卵囊菌Oocystis則與老傘屬Gerronema、空團(tuán)菌屬Cenococcum等植物菌根共生。原生生物通過直接取食病原體和通過捕食誘導(dǎo)細(xì)菌和真菌的分類和功能組成的變化來保護(hù)植物，控制微生物群落的多樣性和結(jié)構(gòu)，從而促進(jìn)植物生長［2，39-40］。例如近爪屬Parachela集團(tuán)中與異常球菌屬Deinococcus、鐵銹菌Ferrovibrio和菌根菌膠膜菌屬Tulasnella、病原菌暗球腔菌屬Phaeosphaeria等11個真菌共生。異常球菌屬可提高鹽敏植物的抗旱能力和對鹽的耐受性，在一定程度上可提高樹體耐鹽堿能力，至于是否抑制或促進(jìn)相關(guān)真菌的生長還未可知［41］。鐵銹菌可水解纖維素促進(jìn)土壤中的能量流動以及氮、磷和硫的循環(huán)［42-44］。Pseudofulvimonas、四鞭毛蟲屬Tetramitus和Telaepolella集團(tuán)與細(xì)菌Angustibacter和病原菌鏈蠕孢屬Dendryphion、菌屬Entyloma和叉絲殼屬Microsphaera等8屬真菌共生；絲壺菌屬和水生角質(zhì)菌屬等4個細(xì)菌屬和犁頭霉屬Absidia等6個屬真菌共生，水生角質(zhì)菌屬具有抗真菌活性［45］；其余各集團(tuán)Globigerinita、尾滴蟲屬Cercomonas、吻滴蟲屬Rhynchomonas分別與鞘藻屬Microcoleus、甲基桿菌屬Methylibium和Methylacidiphilum共生，甲基桿菌屬細(xì)菌具有較高耐鹽、促生和抗菌性能，可促進(jìn)桑的生長，抑制病蟲害發(fā)生［46］。原生生物與大量子囊菌門、擔(dān)子菌門、球囊菌門真菌，有益或潛在致病細(xì)菌共生，表明它們之間存在共同的環(huán)境偏好或潛在的生物相互作用，以互利共生的原則維持根際土壤微生態(tài)穩(wěn)定。

4 結(jié)論

不同樹齡間優(yōu)勢菌門一致但相對豐度存在差異，屬水平上種類組成及相對豐度均受到樹齡的影響。隨著樹齡的變化，樹體長勢和樹體吸收養(yǎng)分能力不同其生命活動也受到土壤理化指標(biāo)的影響，其中氯離子含量和全鹽量可顯著影響原生生物群落多樣性、豐富度和均勻性，速效鉀和氯離子可顯著影響群落覆蓋度。在土壤微生物三級營養(yǎng)結(jié)構(gòu)中，原生生物與真菌細(xì)菌構(gòu)成不同的集團(tuán)共生，表現(xiàn)出相同的環(huán)境偏好共同維持土壤微生物環(huán)境的平衡，為進(jìn)一步探討原生生物種群分化機(jī)制和生態(tài)服務(wù)功能提供理論依據(jù)。