向芬 李維 劉紅艷 銀霞 曾澤萱 周凌云

（湖南省農(nóng)業(yè)科學院茶葉研究所，長沙 410125）

氮肥投入是提高茶樹產(chǎn)量的有利保障，但過量長期施用氮肥反而會造成茶葉減產(chǎn)，土壤板結等問題，亦會對土壤微生態(tài)造成較大的影響［1］，進而影響土壤微生物的活性與群落結構［2］。土壤微生物參與土壤有機質(zhì)分解和養(yǎng)分轉(zhuǎn)化［3-4］，間接促進了土壤物質(zhì)養(yǎng)分循環(huán)［5-6］，是土壤環(huán)境質(zhì)量早期評價的指標［7］。

不同品種茶園的根際土壤細菌種類數(shù)、個體數(shù)各不相同［8］。茶園不同施肥模式比較發(fā)現(xiàn)，與不施肥對照比較，施肥模式的根際土壤細菌、真菌、放線菌等微生物數(shù)目增加，根際土壤微生物量碳、氮和土壤酶活性亦增加［9］，雖然N肥的施用能顯著提高非根際土壤微生物量碳氮含量［10-11］，有機肥施用能提高茶園土壤微生物區(qū)系的多樣性［12］，但對于不同減氮模式處理下的茶園土壤細菌的活性及群落結構鮮見報道，而土壤中90%微生物為細菌［13］，可能與土壤中細菌種類繁雜，分離培養(yǎng)周期長，難以比較全面的檢測土壤細菌的類別有關［14-15］。研究者前期對不同減氮模式下茶園土壤理化性質(zhì)、茶園產(chǎn)量、氮肥農(nóng)學效率進行了研究［16］，發(fā)現(xiàn)減氮處理B能保證產(chǎn)量，氮肥農(nóng)學效率最高，土壤養(yǎng)分利用率高?；诤昊蚪M的高通量測序方法無需單獨對細菌分離培養(yǎng)，可對樣本中所有細菌種類進行分析，現(xiàn)已被廣泛應用于土壤微生物群落研究［1，14，17］。為了探究不同減氮模式對茶樹地下部分細菌菌群組成的影響，本文通過高通量測序方法測定不同施氮量茶樹地下部分的細菌結構及組成，為進一步明確氮代謝相關細菌活動規(guī)律及田間氮肥用量標準提供參考。

1 材料與方法

1.1 材料

試驗地位于湖南省茶葉研究所馬坡嶺基地茶園，地理位置113°04.512′ E，28°12.397′ N，地勢平坦，試驗地土壤為第四紀紅壤，土層深厚，茶園土壤0-20 cm土層基本情況：全碳15.9 g/kg，全氮1.7 g/kg，堿解氮162 mg/kg，有效磷44.8 mg/kg，有效鉀264 mg/kg，碳氮比為9.35。

1.2 方法

1.2.1 試驗設計 茶園連續(xù)4年管理水平保持一致，于2013年10月18日-2016年10月18日對茶園連續(xù)4年進行不同施氮量處理，試驗設16 kg/667 m2（減氮55.6%處理A），26 kg/667 m2（減氮27.8%處理B），36 kg/667 m2（常規(guī)施肥處理C）、不施氮（CK）4個處理水平，每個處理小區(qū)200 m2，每個處理3次重復。具體施肥時以菜籽餅、硫酸鉀型復合肥為基肥，尿素為追肥的方式進行施用，施用方式見表1。

表1 試驗茶園施肥方式Table 1 Fertilizer application methods in tea garden

1.2.2 樣品采集 從茶園施肥處理第4年開始進行樣品采集，非根際土壤樣品、根際樣品采集于2017年10月8日進行，用土鉆（直徑2.5 cm）在每個小區(qū)按Z型5個點（非根際0-20 cm土壤，非根際20-40 cm土壤）進行取樣，每個處理3個重復，非根際土壤樣品在24 h內(nèi)去除茶樹根系，置于-4℃保存，用于非根際土壤細菌菌落結構的測定。離茶樹根頸部25-45 cm處，深0-40 cm，寬50 cm挖取茶樹吸收根，盡量減少對根系的損傷，用軟毛刷清掃根系上的土壤，即根際土。

1.2.3 樣品測定 葉面積指數(shù)于2017年4月采用LAI-2200C冠層分析系統(tǒng)對各減氮處理葉面積指數(shù)進行測定。土壤碳通量：于2017年5月采用Soilbox-343 便攜式土壤呼吸系統(tǒng)（VAISALA，芬蘭）進行測定。將待測區(qū)域內(nèi)植被地上部全部清理干凈，并揀去較多的枯枝落葉。將非透明呼吸室罩在待測區(qū)域，呼吸室與土壤的接觸邊緣用干土密封，每小區(qū)隨機選取3 個點進行測定。根據(jù)單位時間箱內(nèi)CO2濃度的變化計算土壤呼吸速率［18］。

將茶園茶樹根際土壤、非根際0-20 cm土壤、非根際20-40 cm土壤使用TIANamp Bacteria DNA kit核酸提取試劑盒提取樣本中總DNA。對細菌16S rRNA基因V3-V4高變區(qū)進行PCR擴增（Bio-rad T100 梯度PCR儀），引物序列為：341F：CCTAYGGGRBGCASCAG，806R：GGACTACNNGGGTATCTAAT。酶和緩沖液：使用 New England Biolabs 公 司 的 Phusion? High-Fidelity PCR Master。PCR 擴 增 體 系（30 μL）：2×taq PCR mix 15 μL，Primer F（10 μmol/L）3 μL，Primer R（10 μmol/L）：3 μL，gDNA 2.5 μL，H2O 6.5 μL。PCR反應程序：95℃預變性 5 min；34個循環(huán)包括（94℃，1 min；57℃，45 s；72℃，1 min）；72℃最后延伸10 min。將產(chǎn)物送交北京諾禾致源生物信息科技有限公司采用IonS5TMXL高通量測序平臺進行測序。先對 PCR產(chǎn)物進行純化與定量，再使用 Thermofisher 公司的 Ion Plus Fragment Library Kit 48 rxns 建庫試劑盒進行文庫的構建，構建好的文庫檢測合格后，使用Thermofisher 的Ion S5TMXL 進行上機測序。

1.2.4 序列數(shù)據(jù)處理與數(shù)據(jù)分析 基于IonS5TMXL測序平臺，利用單端測序（Single-End）的方法，構建小片段文庫進行單端測序。通過對Reads剪切過濾，操作單元分類數(shù)（Operational taxonomic units，OTUs）聚類，并進行物種注釋及豐度分析，可以揭示樣品物種構成，應用VENN圖可以用來統(tǒng)計樣本中共有的和獨有的OTUs數(shù)目，采取BLAST將代表序列與GreenGene數(shù)據(jù)庫比對，對每個OTU進行物種分類注釋，分別在門、綱、目、科、屬的分類水平上統(tǒng)計樣本的菌群組成，進一步的α多樣性分析（Alpha Diversity）及β多樣性分析（Beta Diversity）可以挖掘樣品之間的差異。為進一步挖掘分組樣品間的群落結構差異，選用T-test統(tǒng)計分析方法對分組樣品的物種組成和群落結構進行差異顯著性檢驗。采用DPS7.05軟件進行微生物豐富度指數(shù)與多樣性指數(shù)統(tǒng)計分析。

2 結果

2.1 減氮處理對茶樹冠層的影響

茶樹為葉用經(jīng)濟作物，葉面積指數(shù)能較好的反應作物產(chǎn)量。從圖1可以看出隨著施氮量的增加茶樹葉面積指數(shù)增加，與對照不施肥比較，減氮處理B、常規(guī)施肥C的葉面積指數(shù)顯著、極顯著增加，但減氮處理B與常規(guī)施肥C的葉面積指數(shù)無顯著性差異。

圖1 減氮處理對茶樹葉面積指數(shù)的影響Fig.1 Effect of nitrogen reduction on tea leaf area index

2.2 測序結果質(zhì)量分析

根據(jù)圖2稀釋性曲線看出，各樣品的序列數(shù)據(jù)量達到30 000條時，OTUs基本趨于飽和，表明測序獲得了足夠多的數(shù)據(jù)量，樣本的測序數(shù)據(jù)量基本覆蓋所有類群的信息。圖3根際樣品曲線在橫軸上的跨度最小，表示物種的豐富度最小，而非根際土壤細菌的代表曲線相對更平緩，表示其中分布的物種更均勻。

圖2 不同土壤樣品測序稀釋曲線圖Fig.2 Rarefaction analysis of 16S libraries from different soil samples

圖3 不同土壤樣品的Rank Abundance曲線Fig.3 Rank Abundance curves of different soil samples

2.3 減氮處理對各土層中細菌多樣性指數(shù)的影響

Chao1指數(shù)和Ace指數(shù)是衡量樣品中物種豐度的重要指標，而香農(nóng)（Shannon）指數(shù)和辛普（Simpson）指數(shù)常用來衡量樣品物種多樣性，Shannon指數(shù)值越大，Simpson指數(shù)值越小，則說明樣品中物種群落多樣性越高。由表2可知，與對照比較，施肥處理的茶樹根際土壤中的細菌物種豐度增加；而非根際土壤中細菌菌落的多樣性和豐度降低，細菌種群數(shù)在非根際土壤中比根際土壤高出44.50%-131.08%，非根際土壤中施肥處理以常規(guī)施肥處理的土壤細菌豐度較高，并以減氮較多的處理A最低，而根際土壤細菌豐度則以常規(guī)施肥處理（36 kg/667 m2）最低，減氮處理B（26 kg/667 m2）最高，但差異不顯著。其中細菌多樣性指標Shannon指數(shù)顯示減氮處理B的土壤細菌多樣性在根際土壤與非根際0-20 cm土壤中最高，而在非根際20-40 cm土壤則是常規(guī)施肥處理C的細菌多樣性顯著增加（P<0.05）。

表2 各土壤樣品微生物豐富度指數(shù)與多樣性指數(shù)Table 2 Diversity and richness index of bacteria from twelve samples

2.4 各減氮處理土壤樣品細菌的Venn圖比較

根際土壤、非根際0-20 cm土壤及非根際20-40 cm土壤共有OTUs分別為505、854及835個。其中根際土壤、非根際0-20 cm土壤、非根際20-40 cm土壤的OUTs分布隨著施氮量增加逐漸降低，其中非根際20-40 cm土壤各處理分別占31.58%、27.86%、26.23%（圖4）。由圖5可知，茶樹地下部分各樣品在門水平上的優(yōu)勢菌基本為氮代謝相關菌類。其中共有3個優(yōu)勢類群，分別為變形菌門（Proteobacteria）、放線菌門（Actinobacteria）、綠彎菌門（Chloroflexi），其中相對富集在根際的優(yōu)勢菌為藍藻菌門（Cyanobacteria）；非根際土壤中的則為酸桿菌門（Acidobacteria），另外3種分別為擬桿菌門（Bacteroidetes）、芽單胞菌門（Gemmatimonadetes）、糖化菌門（Saccharibacteria）。在屬的水平上（圖6），根際與非根際土壤細菌優(yōu)勢菌差異明顯，其中根際為旁氏伯克霍爾德菌（Burkholderia Paraburkholderia）、放線菌（Actinospica）、熱孢子絲菌（Thermosporothrix）3種，而非根際土壤則為水蛭桿菌（Mizugakiibacter）、嗜酸棲熱菌屬（Acidothermus）2種。并發(fā)現(xiàn)常規(guī)施肥處理C在非根際0-20 cm土壤與非根際20-40 cm土壤有兩個特異菌，分別為克羅斯耶拉菌（Crossiella）與羅爾斯頓菌（Ralstonia）。并從圖7的組間硝化螺旋菌門顯著性差異比較發(fā)現(xiàn)，非根際0-20 cm土壤的硝化螺旋菌門豐度較高，與對照比較，根際土壤各處理之間無顯著性差異；根際土壤的各施肥處理與非根際0-20 cm土壤的CK處理差異極顯著；非根際0-20 cm土壤的減氮處理A、B與CK差異極顯著，非根際20-40 cm土壤的減氮A、B處理與非根際0-20 cm土壤CK處理差異極顯著。

圖4 不同土壤樣品OTUs分布的Venn圖Fig.4 Venn diagrams of OTUs distribution of different soil samples

圖5 門水平的細菌組成聚類堆積圖Fig.5 Cluster stacking graph of bacteria composition at phylum level

圖6 屬分類水平的細菌組成聚類堆積圖Fig.6 Cluster stacking graph of bacteria composition at genus level

圖7 組間硝化螺旋菌門顯著性差異統(tǒng)計圖Fig.7 Statistical chart of significant difference with Nitrospirae between groups

2.5 減氮處理土壤樣品細菌群落結構的相似性 分析

采用主坐標分析方法對各處理組土壤樣品的細菌在屬分類水平上進行分析，評估細菌群落間的差異。從圖8可以看出，主成分1（PC1）和主成分2（PC2）對固氮細菌群落結構變異的解釋量分別為35%和19.34%，A20、CK20與C40、CK40樣本相較與其他樣本而言分布較為分散，樣品間距離較遠。

圖8 各土壤樣品細菌群落主坐標分析Fig.8 PCoA analysis of bacteria in different soil samples

2.6 減氮處理對茶園土壤碳通量的影響

土壤碳通量變化與土壤微生物變化密切相關，對不同減氮處理后茶園的土壤碳通量進行測定，結果如圖9所示。以減氮處理B的土壤碳通量較高，為2.70 μmol/m2/s；常規(guī)施肥 C 處理次之，為2.63 μmol/m2/s；對照處理最低，為2.02 μmol/m2/s；各處理間差異不顯著。

圖9 減氮處理對茶園土壤碳通量的影響Fig.9 Effect of nitrogen reduction on soil carbon flux in tea garden

3 討論

Illumina Miseq測序技術能夠克服傳統(tǒng)檢測方法的缺陷，具有通量高、錯誤率低、成本低、流程自動化、速度快等優(yōu)勢。已被廣泛用于微生物群落結構研究。本研究采用高通量Illumina Miseq測序技術對茶樹地下部分微生物進行系統(tǒng)分析，能真實全面地了解土壤及根際微生物的組成情況。施氮能促進茶樹生長，茶葉葉面積指數(shù)增加，提高茶葉產(chǎn)量，本研究結果表明施氮能提高茶樹的茶葉葉面積指數(shù)。同時發(fā)現(xiàn)減氮處理B的茶園土壤固氮細菌組成增加，茶樹地下部分的根際微生物多樣性增加。從而推測施氮量為26 kg/667 m2（減氮處理B）的處理有利于土壤中氮素循環(huán)平衡與氮素固定［19］，也有利于維持土壤微生物平衡。本研究結果表明施肥處理后茶園非根際土壤細菌多樣性與豐度均降低，推測氮素促使茶園土壤的變形菌門數(shù)量變多，從而改變了土壤的細菌群落構成，這與李鳳霞等［20］的研究結果一致。各施肥處理茶樹地下部分在門水平上的優(yōu)勢菌基本為氮代謝相關菌類，如優(yōu)勢類群之一的變形菌門以及在根際的優(yōu)勢菌藍藻菌門（Cyanobacteria）［1］。其中施氮量為26 kg/667 m2（減氮處理B）的處理的Chao1指數(shù)、Ace指數(shù)和Shannon 指數(shù)結果表明0-20 cm土壤中細菌豐度較高，尤其是多樣性最高，而在非根際20-40 cm土壤中則是施氮量36 kg/667 m2（常規(guī)施肥C處理）的細菌多樣性顯著增高，這說明適量施氮的茶園土壤細菌種類較多，并主要集中在養(yǎng)分利用的根際及周邊的土壤中，與吳海寧等［1］的結果一致，亦能和研究者課題組前期研究結果減氮處理B能保證常規(guī)施肥的產(chǎn)量，各施肥處理中以減氮處理B氮肥農(nóng)學效率最高，土壤養(yǎng)分剩余較少，氮肥的利用效率較高［16］相互印證。進一步印證合理施氮條件下茶園土壤微生物的種群構成、分布與數(shù)量及土壤養(yǎng)分具有良好的調(diào)節(jié)作用［21］。而施氮量36 kg/667 m2（常規(guī)施肥C處理）在非根際20-40 cm土壤中的多樣性較高，可能是由于在其非根際土壤中多了兩種特異菌：克羅斯耶拉菌與羅爾斯頓菌。主成分分析結果顯示CK40、C40樣本相較與其它樣本而言分布較為分散，樣品間距離較遠，表明該土層攜帶更為豐富的物種多樣性，與Shannon指數(shù)的結果一致，其余樣品間的點距離比較近，說明它們的固氮細菌群落結構比較相似。

組間硝化螺旋菌門顯著性差異統(tǒng)計情況表明，非根際土壤的硝化螺旋菌門較根際土壤豐度高，非根際土壤中以0-20 cm土壤的硝化螺旋菌門豐度較高，可能是由于茶園常年耕作的深度在20 cm左右，茶園土壤0-20 cm的土層較疏松，同時茶園的養(yǎng)分施用也集中在此土層，所以該土層的土壤微生物特別是固氮微生物較多。與對照比較，根際土壤各處理之間無顯著型差異；根際土壤的各施肥處理與非根際0-20 cm土壤的CK處理差異極顯著；非根際0-20 cm土壤的減氮處理A、B與CK差異極顯著，非根際20-40 cm土壤的減氮A、B處理與非根際0-20 cm土壤CK處理差異極顯著，而非根際土壤的常規(guī)施肥處理C與對照差異不顯著，亦說明常規(guī)施肥處理C的土壤中豐度較高，可能是由于硝化螺旋菌門較多，說明可能常規(guī)施肥處理C硝化作用較多，而茶樹偏好銨態(tài)氮［22］，不利于其養(yǎng)分的有效利用。土壤碳通量的結果顯示以減氮處理B的土壤二氧化碳通量指數(shù)最高，與高通量測序結果土壤細菌多樣性指數(shù)分析的結果基本一致。

4 結論

本研究利用高通量測序技術明確減氮處理B下的細菌種類較多，細菌豐度較高，并主要集中在養(yǎng)分吸收利用的根際及周邊土壤中，有利于茶園養(yǎng)分的吸收，土壤二氧化碳通量結果基本一致。雖然常規(guī)施肥處理C下的非根際20-40 cm土壤中細菌多樣性較多，但其中一部分為克羅斯耶拉菌與羅爾斯頓菌以及促進硝化作用的硝化螺旋菌門，不利于養(yǎng)分的吸收與轉(zhuǎn)化，長期過量施用氮肥對茶園土壤生態(tài)系統(tǒng)不利。這為進一步探討氮肥施用效率機制及開發(fā)微生物肥料奠定了基礎。