李巍 黃彪 黃敏敏 劉文靜 郭嘉 林香信

摘要：【目的】研究不同土地利用方式下土壤氨氧化微生物多樣性的季節(jié)響應(yīng)，以期了解武夷山不同土壤利用類型氨氧化微生物群落特征和季節(jié)變化?！痉椒ā吭诟＝ㄊ∥湟纳降某＞G闊葉林區(qū)域，分別采集同一緯度下未被開發(fā)利用的原始林地和已被開發(fā)利用的茶園土壤，利用Illumina Hiseq高通量測序技術(shù)對土壤氨氧化微生物的相對豐度及多樣性進(jìn)行分析，同時測定土壤理化性質(zhì)，并對二者進(jìn)行關(guān)聯(lián)性分析?！窘Y(jié)果】（1）武夷山常綠闊葉林土壤含有豐富的氨氧化微生物資源，不同土地利用方式下土壤氨氧化微生物在Family水平相似，主要包括氨氧化古菌界的Unclassified CandidatusⅣitrososphaera、 Unclassified Crenarchaeota、 Unclassified Thaumarchaeota、Nitrosopumilus、Others;氨氧化細(xì)菌界的Nitrosomonas、Nitrosospira

Unclassified Nitrosomonadaceae、Others，其中茶園擁有獨(dú)特的Ⅳitrosovibrio（亞硝化葉菌屬）存在。（2）武夷山林地和茶園氨氧化微生物相對豐度有顯著變化且隨季節(jié)變化顯著。（3）土壤氨氧化微生物群落結(jié)構(gòu)與環(huán)境因子的關(guān)聯(lián)性分析表明，不同的環(huán)境因子對土壤氨氧化細(xì)菌和土壤氨氧化古菌的影響不同，且不同季節(jié)的主要影響因子也存在差異，有效鉀是影響氨氧化細(xì)菌變化的主要環(huán)境因子，而銨態(tài)氮是影響氨氧化古菌的主要環(huán)境因子?！窘Y(jié)論】隨季節(jié)變化武夷山常綠闊葉林土壤氨氧化微生物多樣性呈現(xiàn)出規(guī)律性變化，不同土地利用方式對氨氧化微生物種類影響不大，但對其種群的相對豐度及季節(jié)相對豐度影響較大。

關(guān)鍵詞：Illumina Hiseq高通量測序;不同土地利用方式;氨氧化細(xì)菌;氨氧化古菌

中圖分類號：S 154

文獻(xiàn)標(biāo)志碼：A

文章編號：1008-0384（ 2020） 10-1119-12

0 引言

【研究意義】武夷山具有世界同緯度上現(xiàn)存面積最大、保存最完整的地帶性常綠闊葉林群落，是武夷山分布面積最廣、最具代表性的森林生態(tài)系統(tǒng)，具有豐富的生物資源和優(yōu)良的自然環(huán)境，是研究生物多樣性的優(yōu)勢自然基地[1]。氮素既是樹木生長最主要的限制因子，直接影響土壤質(zhì)量以及森林植被的生產(chǎn)力，又會影響生態(tài)環(huán)境的變化（如酸雨、溫室效應(yīng)、硝態(tài)氮淋溶導(dǎo)致的水體富營養(yǎng)化等）[2-7]。硝化作用是土壤氮循環(huán)的關(guān)鍵環(huán)節(jié)之一，對土壤氮素轉(zhuǎn)化、吸收等過程有重要意義[8]。土壤微生物在土壤硝化過程中發(fā)揮著決定性作用，其中氨氧化細(xì)菌（ Ammonia-oxidizing bacteria，AOB）和氨氧化古菌（ Ammonia-oxidizing archaea，AOA）是土壤氮循環(huán)的重要驅(qū)動者，也是硝化速率限制性步驟的關(guān)鍵微生物[9-16]。因此，在林地生態(tài)系統(tǒng)中氮循環(huán)過程與土壤氨氧化微生物多樣性密切相關(guān)[17]?！厩叭搜芯窟M(jìn)展】近年來，國內(nèi)外開展了對不同森林類型土壤氨氧化微生物群落多樣性的相關(guān)研究，發(fā)現(xiàn)在森林土壤中氨氧化古菌是最豐富的氨氧化微生物，土壤中氨氧化古菌的amoA基因數(shù)量比氨氧化細(xì)菌高出2～3個數(shù)量級[18]，但也有研究證明氨氧化古菌主導(dǎo)了蘇格蘭典型農(nóng)田土壤氨氧化過程[19]。不同利用方式對草原氨氧化微生物的研究中發(fā)現(xiàn)放牧和圍封樣地的氨氧化古菌豐度均高于氨氧化細(xì)菌豐度[20]。四川南充市不同利用類型的土壤理化性質(zhì)差異顯著，氨氧化微生物群落結(jié)構(gòu)存在明顯的分異特征[21]。前人研究表明，耕作方式、施肥和種植不同作物等都會影響土壤的氨氧化微生物群落結(jié)構(gòu)[22]。隨著經(jīng)濟(jì)發(fā)展的需要，許多林地被開發(fā)利用，使土地利用方式發(fā)生了變化，造成土地退化和土壤生態(tài)功能的改變?nèi)找鎳?yán)重，土壤利用或管理方式的變化也逐漸引起了廣大學(xué)者的關(guān)注[21-25]。武夷山林地改茶園的現(xiàn)象較為普遍，且據(jù)Zulkamaen等[26]的研究發(fā)現(xiàn)在紅壤地區(qū)由于土地利用和管理方式的改變會降低土壤有機(jī)質(zhì)含量，同時會加劇土壤酸化，茶園的AOB豐度大于森林?！颈狙芯壳腥朦c(diǎn)】目前針對酸性森林土壤中氨氧化古菌對硝化過程貢獻(xiàn)的研究雖有相關(guān)報道[17，20，22]，但對于不同土壤利用類型下森林土壤氨氧化群落結(jié)構(gòu)的變化研究報道甚少。土地利用方式的改變是如何改變土壤氨氧化微生物，從而影響硝化作用和反硝化作用的生理機(jī)制還有待進(jìn)一步研究?！緮M解決的關(guān)鍵問題】本研究采用土壤分析化學(xué)和Illumina Hiseq高通量測序技術(shù)相結(jié)合，研究武夷山常綠闊葉林不同土壤利用類型對土壤氨氧化微生物群落多樣性的影響及其季節(jié)動態(tài)變化規(guī)律，探討不同環(huán)境因子對土壤氨氧化微生物群落多樣性及季節(jié)動態(tài)變化的影響，以期了解氨氧化微生物在武夷山常綠闊葉林森林生態(tài)系統(tǒng)氮循環(huán)中的作用，為維持武夷山自然生態(tài)系統(tǒng)平衡提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究樣地選自武夷山市星村鎮(zhèn)四新村附近常綠闊葉林下兩種土壤利用類型（山地和茶園）的地塊（27°35'24.23”N，1 17°45'55.43”E），淘拔450—580m.該區(qū)域?qū)俚湫偷膩啛釒Ъ撅L(fēng)濕潤氣候，年均溫17.0～18.4℃，年均相對濕度75%～84%，年降水量1 800mm，年均日照時數(shù)1 910.2 h，無霜期227～246 d[1]。土壤以黃壤和紅壤為主，武夷山常綠闊葉林未開發(fā)利用原始林地的森林植被種類組成復(fù)雜，呈復(fù)層結(jié)構(gòu)，優(yōu)勢種為米櫧（Castanopsis carlesii）、甜櫧（Castanopsis eyrei）、木荷（Schima superba）等;已被開墾的茶園區(qū)域，常年種植肉桂、水仙、105、304等品種的巖茶。

1.2 土壤樣品采集

采集同一生境下兩種不同土地利用方式的土壤為研究對象：未開發(fā)利用的林地（SD）和相鄰開發(fā)利用的茶園（CY），茶園每年施用復(fù)合肥一次。根據(jù)不同季節(jié)設(shè)定取樣時間為秋季2015年9月18日、冬季2015年12月22日，春季2016年3月29日、夏季2016年6月21日，共采樣4次，分別采集不同土壤利用類型的表層土壤（0～10 cm），每種利用方式土壤選取3個樣點(diǎn)，每個樣點(diǎn)采用隨機(jī)取點(diǎn)法取樣，混合均勻后取1kg裝入滅菌袋中，在4℃冷藏箱保存并盡快帶回實驗室，挑根去石礫，將土樣分成2份，一份置于超低溫冰箱（一70℃）中保存，用于土壤微生物多樣性分析;另一份自然風(fēng)干后過篩，用于土壤理化性質(zhì)分析。

1.3 土壤理化性質(zhì)測定

采用凱氏定氮法測定土壤全氮;氯化鉀浸提.靛酚藍(lán)比色法測定土壤銨態(tài)氮;重鉻酸鉀氧化.外加熱法測定土壤總有機(jī)碳;冷硝酸浸提一火焰光度法測定土壤有效鉀;電位法測定土壤pH;烘干法測定土壤濕度;土壤溫濕度計測定土壤溫度。具體測定方法參考《土壤農(nóng)化分析》[27]。

土壤硝化勢的測定參照宋亞娜等[13-14]在稻田測定土壤硝化勢的方法，具體方法：稱取5 g鮮土置于50 mL離心管中，加20 mL含（NH4） 2S04濃度為1 mmol.L^-1的磷酸緩沖液，黑暗的條件下保持在室溫培養(yǎng)24 h，加入質(zhì)量濃度為10 g.L^-1KC103來抑制亞硝酸鹽氧化;培養(yǎng)完成后加入5 mL 2 mol.L^-1的KCI，提取NO2-N后離心，再用格里斯試劑比色法測定提取液中NO2_N濃度（單位：μg.g^-1.h^-1干土）。

1.4 土壤氨氧化微生物群落多樣性分析

1.4.1 土壤DNA提取和測序 將上述-70℃保存的土壤樣品進(jìn)行宏基因組DNA提取，對檢測合格的DNA樣品進(jìn)行高通量測序。土壤氨氧化細(xì)菌的測序以Illumina Hiseq 2 500 PE250為測序平臺，以土壤總DNA作為模板，土壤氨氧化細(xì)菌基因利用特異引物對AOB+ad-1（5 '-TCCCTACACG ACGCTCTTCCGATCTGGGGTTTCTACTGGTGG T-3）和AOB+ad-2（5-AGACGTGTGCT CTTCCGAT CTCCCCTCKGSAAAG CCTTCTTC-3）進(jìn)行擴(kuò)增;土壤氨氧化古菌基因利用特異引物對AOA+ad-I（5，_TCCCTACACGA CGCTCTTCCATCTSTAATGGTCTGGC TTAGACG-3）和AOA+ad-2（5-AGACGTGTGCTCT TCCGATCTGCGGCCATCC ATCTGTATGT-3'）進(jìn)行擴(kuò)增，引物由北京百邁客生物科技有限公司合成。

1.4.2土壤微生物數(shù)據(jù)分析 土壤氨氧化微生物序列分析利用QIIME software package和UPARSE pipeline完成。利用QIIME的默認(rèn)參數(shù)過濾掉低質(zhì)量、具有嵌合體結(jié)構(gòu)的序列;然后利用UPARSE軟件生成OTU（ Operational taxonomic units），利用RDP從每個OTU中選擇一條序列作為OTU代表;根據(jù)獲得的OTU數(shù)據(jù)，以測序樣本的序列數(shù)為橫坐標(biāo)，OTU值作為縱坐標(biāo)繪出每個樣品的稀釋曲線，同時根據(jù)樣本中所共有和獨(dú)有的OTU數(shù)目作OTU分布Venn圖。然后對數(shù)據(jù)進(jìn)行冗余分析（ Redundancy analysis，RDA）、直接梯度分析，從統(tǒng)計學(xué)角度評價一個或一組變量與另一組多變量數(shù)據(jù)之間的關(guān)系。同時，通過顏色梯度來代表數(shù)據(jù)矩陣中數(shù)值的大小并根據(jù)物種或樣品豐度相似性進(jìn)行聚類形成熱圖（ Heat map）展示。

1.5 土壤理化性質(zhì)與微生物多樣性的相關(guān)性分析

利用Excel 2010和R語言進(jìn)行數(shù)據(jù)分析，通過主成分分析（ Principal Component Analysis，PCA）將多組數(shù)據(jù)的差異反映在二維坐標(biāo)圖上，坐標(biāo)軸能夠最大反映方差的兩個特征值，兩個樣品距離越近表示這兩個樣品的組成越相似。

2 結(jié)果與分析

2.1 不同土地利用方式對土壤氨氧化微生物豐度的影響

武夷山林地土壤中的氨氧化細(xì)菌變化具有一定的規(guī)律性（圖1-A），在春季，亞硝化單胞菌屬Nitrosomonas是氨氧化細(xì)菌中的優(yōu)勢屬，占50.08%，但至夏季時驟然減少，僅占0.98%，秋季低至0.000 695%，冬季略微上升至0.001 16%。在春季，相對豐度第二位為Nitrosomonadaceae的其他屬，占28.7g%，但在夏季的土壤樣品中達(dá)到90.92%，在秋冬季節(jié)分別為84.03%和77.94%。在春季相對豐度第三位的細(xì)菌是亞硝化螺菌屬Nitrosospira，占2.50%，該屬在夏、秋、冬季的豐度也相對較低，分別為4.02%、9.33%、9.20%。而亞硝化葉菌屬Nitrosovibrio在四季的土壤中均未被檢出。Others部分包括了大部分未歸屬的Proteobacteria及在界分類中未能歸屬的OTU，其在四季的相對豐度分別為18.62%、4.09%、6.63%、12.86%。

古菌界相對豐度變化比較穩(wěn)定（圖1-B），Unclassified Candidatus Nitrososphaera四季相對豐度分別為0.061%、0.021%、0.051%、0.042%;泉古菌屬Unclassified Crenarchaeota的四季相對豐度分別為1.24%、1.72%、1.38%、1.51%; UnclassifiedThaumarchaeota的四季相對豐度為4.68%、2.77%、2.8g%、3.43%，Nitrosopumilus屬相對豐度極低，僅在春季被檢出，相對豐度為0.000 167%。

武夷山茶園土壤中的氨氧化微生物在屬水平上的相對豐度較為穩(wěn)定（圖1）。細(xì)菌界中Nitrosomonadaceae下的其他屬是氨氧化細(xì)菌中的優(yōu)勢屬，在4個季節(jié)的樣本中分別占73.75%、65.62%、60.95%、68.61%，即該屬在春、夏、秋三季不斷減少，在冬季回升。與林地土壤不同的是，Nitrosospira是茶園土壤中的第二大屬，其變化規(guī)律與Nitrosomonadaceae下的其他屬相反，在春、夏、秋三季不斷增加，但在冬季降低。Nitrosomonas在四季中的相對豐度逐漸增加，分別為0.000 2%，0.01%，0.13%，4.74%.Nitrosovibrio的相對豐度是最低的，其規(guī)律為冬、春、夏季逐漸降低，秋季豐度最高，分別為0.005 6%、0.007 7%、0.0002 2%、0.018%。Others在四季的相對豐度分別為0.25%、0.055%、0.008 2%、0.022 8%，即在冬春季節(jié)逐漸增加，夏秋季減少。

武夷山茶園土壤中的氨氧化古菌的UnclassifiedCrenarchaeota與Unclassified Candidatus Nitrososphaera在秋、冬、春三季逐漸下降，但在夏季快速下降，變化規(guī)律與Others相反，僅在夏季表現(xiàn)出了較高的相對豐度，其中，Unclassified Crenarchaeota的四季相對豐度為1.42%、13.48%、3.32%、1.44%， UnclassifiedCandidatusNitrososphaera的相對豐度分別為0.028%、7.48%、0.30%、0.24%。

2.2不同土地利用方式對土壤氨氧化微生物組成的影響

根據(jù)武夷山林地土壤樣本不同分組繪制成OTU-Venn（韋恩）圖，通過對土壤樣品測序所得OTU（圖2）進(jìn)行比對，結(jié)果表明：林地土壤樣品中，春季細(xì)菌界特有OTU 16個，夏季23個，秋季35個，冬季78個;春季古菌界特有OTU數(shù)為8個，夏季為2個，秋季1個，冬季2個。韋恩圖能夠在一定程度上從OTU層面反映出土壤菌群的季節(jié)多樣性，利于發(fā)現(xiàn)特殊物種。從研究結(jié)果中可以看出，春季土壤古菌界中發(fā)現(xiàn)1條特有OTU，歸屬于Nitrosopumilus。由于所得OTU在RDP分類學(xué)數(shù)據(jù)庫分類時大部分物種無法歸類到屬或種，導(dǎo)致特殊OTU所代表的物種無法被精確挖掘，但特殊OTU的數(shù)量可證明不同季節(jié)的樣本中存在不同的微生物。其中，春季細(xì)菌界中特有OTU中有10個被歸屬為Nitrosomonas，但也有5個OTU屬于Others，1條屬于Nitrosomonadaceae;古菌界中則有7條屬于Others，1條屬于Nitrosopumilus。

根據(jù)武夷山茶園土壤微生物文庫中的OTU繪制成韋恩圖（圖3），比較4個季節(jié)特有的OTU可知，細(xì)菌界中，春季擁有4個特殊的OTU，夏季10個，秋季1個，冬季11個。其中，春季特殊OTU中3個屬于Others，1個屬于Nitrosomonas;夏季的10個OTU中有6個被分類為Others，1個為Nitrosomonas，3個屬于Nitrosomonadaceae的其他屬;秋季僅含的1個特殊的OTU歸屬于Others;冬季細(xì)菌界包含11個特殊OTU。古菌界中所含的特有OTU數(shù)量較少，其中，春季1個，屬于Others;夏季3個，2個OTU屬于Others，1個屬于Unclassified Crenarchaeota;秋季1個，屬于Others;冬季3個，屬于Others。雖然從RDP分類學(xué)數(shù)據(jù)庫得到的分類結(jié)果無法識別出不同樣本中氨氧化微生物的特殊核心微生物，但是從所得Venn圖可知，不同樣本中存在不同的微生物。

2.3不同土地利用方式對土壤氨氧化微生物群落結(jié)構(gòu)的影響

根據(jù)組間PCA分析結(jié)果，在第一、第二主成分上的累積貢獻(xiàn)率大于90%，代表性強(qiáng)，在細(xì)菌界PCA圖4-A中，茶園土壤與林地土壤組內(nèi)、組間的樣本離散性都很高，說明樣本差異性受季節(jié)及采集地影響很大。在古菌界PCA圖.B中，茶園土壤組內(nèi)樣本離散性很高，說明樣本差異性受季節(jié)影響很大，但林地土壤樣本的離散性較小，即具有重復(fù)性，說明樣本內(nèi)的古菌變化受季節(jié)影響較小。林地與茶園古菌組間離散型高，說明不同采集地獲得的樣本古菌組成差異性大。

從林地土壤及茶園土壤樣品屬分類水平上的Heat map聚類分析（圖5），可以反映出不同季節(jié)、不同采集地樣本群落組成的相似性和差異性。由樣本細(xì)菌界熱圖聚類結(jié)果可知，夏季林地土壤樣本與秋季林地土壤樣本具有最高的相似性，春季林地土壤與其他土壤樣本具有最大的差異性。在所有樣本中，春季林地土壤含有最高Nitrosomonas及Others豐度;夏季林地土壤中Unclasified Nitrosomonadaceae相對豐度最高;Nitrosospira在秋季茶園土壤中相對豐度最高;Nirosovibrio在冬季茶園土壤中相對豐度最高。由樣本古菌界的熱圖可知，夏季林地、秋季林地、春季茶園3個土壤樣本古菌相對豐度具有很高的相似性;夏季茶園與其他土壤樣品古菌相對豐度差異性最高，Nitrosopumilus，Unclassified CandidatusNitrososphaera，Unclassified Crenarchaeota在夏季茶園土壤樣本中相對豐度最高，Unclassified Thaumarchaeota在春季林地土壤中的豐度最高。

在氨氧化細(xì)菌界中，由于土壤耕作的利用方式不同，在相同季節(jié)，茶園土壤中Nitrosospira相對豐度總是高于林地土壤，但四季豐度變化趨勢相同，即氣候變化對Nitrosospira的影響是相同的。此外，Nitrosovibrio僅存在茶園土壤中，而不存在林地土壤中。特別的是，林地土壤中Nitrosomonas的相對豐度在春、夏季比茶園土壤多，但在秋冬季節(jié)比茶園土壤少，而茶園土壤中Nitrosomonas的相對豐度較穩(wěn)定;Unclassified Nitrosomonadaceae在林地土壤中的相對豐度波動較大，在茶園土壤中處于較穩(wěn)定狀態(tài)。因此，在對比相同季節(jié)不同采集地樣品后可知，土壤中Nitrosomonas與Unclassified Nitrosomonadaceae的相對豐度與土地利用方式具有很大關(guān)系。

在氨氧化古菌界中，在相同季節(jié)，茶園土壤中的Nitrosopumilus的相對豐度始終比林地土壤高;林地土壤中Unclassified Thaumarchaeota相對豐度始終比茶園土壤高;在夏季時，茶園土壤中UnclassifiedCrenarchaeota及Unclassified Candidatus Nitrososphaera的相對豐度比林地土壤高出很多，而在其他季節(jié)則與林地土壤的相對豐度相似，因此，這些微生物在夏季相對豐度上升與茶園土壤的利用方式有很大關(guān)聯(lián)。

2.4 土壤氨氧化微生物群落結(jié)構(gòu)與環(huán)境因子的關(guān)聯(lián)性分析

氨氧化微生物受多種環(huán)境因素影響，這些因素包括土壤中銨態(tài)氮濃度、pH值、溫度等[28-30]。因此，本研究測定了武夷山常綠闊葉林下林地和茶園的有機(jī)質(zhì)、銨態(tài)氮、速效磷、有效鉀和pH值5項基本土壤理化性質(zhì)（表1），對5個環(huán)境因子分別進(jìn)行統(tǒng)計分析，得到其重要性排序，以評估每個備選環(huán)境變量對相應(yīng)變量的解釋貢獻(xiàn)，其中，環(huán)境因子對細(xì)菌群落在約束性排序的前兩軸的累積解釋量為61.68%，對古菌的累積解釋量為81.60%，因此，環(huán)境因子的變化能夠解釋林地土壤中氨氧化微生物的變化（表2）。

對環(huán)境因子和微生物物種之間進(jìn)行相關(guān)系數(shù)和重要性分析（表3），結(jié)果表明，環(huán)境因子對氨氧化細(xì)菌的重要性由大到小依次為有效鉀>pH值>總有機(jī)碳>銨態(tài)氮>速效磷;對氨氧化古菌結(jié)構(gòu)影響重要性大小依次為銨態(tài)氮>有效鉀>pH值>總有機(jī)碳>速效磷。由此可知，環(huán)境中的有效鉀、pH值、總有機(jī)碳與氨氧化微生物的相關(guān)性較大，相關(guān)系數(shù)均大于0.2;銨態(tài)氮與古菌的相關(guān)性最大，相關(guān)系數(shù)為0.457 8，但其與細(xì)菌的相關(guān)性最小僅為0.066 9;速效磷與細(xì)菌豐度變化的相關(guān)系數(shù)0.056 4，與古菌的相關(guān)系數(shù)僅0.047 3，所以氨氧化微生物受土壤速效磷影響較小。

對武夷山常綠闊葉林下不同土地利用類型土壤中氨氧化細(xì)菌與古菌群落結(jié)構(gòu)與季節(jié)變化中環(huán)境因子的關(guān)聯(lián)性分別進(jìn)行RDA分析（圖6）的結(jié)果亦表明，不同的環(huán)境因子對氨氧化細(xì)菌和古菌的影響不同，而且不同季節(jié)的主要影響因子也存在差異。對于氨氧化細(xì)菌（圖6-A）來說，相關(guān)性最大的環(huán)境因子是有效鉀;而對于氨氧化古菌（圖6-B），相關(guān)性最大的環(huán)境因子是銨態(tài)氮。對于林地樣本來說，氨氧化細(xì)菌和古菌的物種結(jié)構(gòu)相似性季節(jié)變化一致，均是夏、秋、冬三季相似性高，而春季則與其他3個季節(jié)相似性較低;而對于茶園樣本來說，氨氧化細(xì)菌的物種結(jié)構(gòu)為春夏兩季相似性高，而與秋冬兩季節(jié)相似性較低，氨氧化古菌的物種結(jié)構(gòu)為春、秋、冬三季相似性高，而夏季則與其他3個季節(jié)相似性較低。

從RDA（圖6-A）中可以看出Nitrosospira菌屬與Apo、Aph呈正相關(guān)，而與pH、An、Oc呈負(fù)相關(guān);Nitrosomonas和Others菌屬與pH、An、Oc呈正相關(guān)，而與Apo、Aph負(fù)相關(guān)，分析物種垂直投影于環(huán)境變量延長線上距離可以看出，pH值上升對Nitrosomonas的生長最有利，UnclassifiedNitrosomon-adaceae適合相對低的pH值;有效鉀對秋季茶園土壤氨氧化細(xì)菌群落結(jié)構(gòu)的差異性具有較大貢獻(xiàn)。從物種之間的卡方距離來看，Others與Nitrosomonas親疏關(guān)系較近，可能存在互惠的生存關(guān)系，而與其他菌屬親疏關(guān)系較遠(yuǎn)，可能存在競爭關(guān)系。

從RDA（圖6-B）中可以看出，UnclassifiedCandidatusNitrososphaera、Unclassified Crenarchaeota，Nitrosopumilus均與An呈正相關(guān)，而與Apo、pH呈負(fù)相關(guān);Unclassified Thaumarchaeota與pH、Oc呈正相關(guān)，而與Apo、Aph呈負(fù)相關(guān)。分析物種垂直投影于環(huán)境變量延長線上距離可以看出，UnclassifiedCandidatus Nitrososphaera適合相對較高銨態(tài)氮的環(huán)境，而Others適合相對低的銨態(tài)氮的環(huán)境。銨態(tài)氮對春季林地土壤氨氧化古菌群落結(jié)構(gòu)的差異性具有較大貢獻(xiàn)。從物種之間的卡方距離來看，UnclassifiedCandidatus Nitrososphaera、Unclassified Crenarchaeota與Nitrosopumilus親疏關(guān)系較近，可能存在互惠的生存關(guān)系，而與其他菌屬親疏關(guān)系較遠(yuǎn)，可能存在競爭關(guān)系。

硝化勢是研究氨氧化微生物活動的重要指標(biāo)，從圖6-C中可以看出，Nitrosospira及Nirosovibrio適合在硝化勢較低的環(huán)境中生長，即在硝化勢較高的環(huán)境中可能因為無法適應(yīng)環(huán)境而失去對硝化作用的貢獻(xiàn)，而Unclassified Nitrosomonadaceae則在較高的硝化勢環(huán)境中依然能夠完成硝化作用。夏季林地的硝化勢最高，因此夏季林地土壤中細(xì)菌的結(jié)構(gòu)容易形成較高的硝化勢。從圖6-D中可以看出，Nitrosopumilus適合在較低的硝化勢環(huán)境中生長，而Other適合在較高的硝化勢環(huán)境中生長。秋季林地樣本的硝化勢是最高的，說明秋季林地土壤中古菌結(jié)構(gòu)適合造成較高的硝化勢。同時，針對夏季林地土壤中的細(xì)菌結(jié)構(gòu)和秋季林地土壤中的古菌結(jié)構(gòu)進(jìn)行分析，細(xì)菌中的Unclassified Nitrosomonadaceae和古菌中的Other可能是形成高硝化勢的主要微生物，抑或可能是這些時期的群落結(jié)構(gòu)容易使菌群對土壤產(chǎn)生較高的硝化勢，而不是單一微生物引起的。然而硝化勢的形成是由微生物之間的緩慢的更替形成，若無低硝化勢時期的微生物群落結(jié)構(gòu)，可能也無法形成高硝化勢的微生物群落結(jié)構(gòu)。

3 討論

土壤氨氧化微生物在武夷山常綠闊葉林森林生態(tài)系統(tǒng)的氮循環(huán)過程中具有重要的意義。本文結(jié)合了土壤化學(xué)和現(xiàn)代分子生態(tài)學(xué)研究方法對武夷山常綠闊葉林生態(tài)系統(tǒng)下不同土地利用類型的氨氧化微生物的群落多樣性進(jìn)行了初步分析，研究結(jié)果表明：武夷山常綠闊葉林原始林地土壤及同域下茶園土壤中氨氧化微生物種類相似，在茶園土壤中僅發(fā)現(xiàn)亞硝化葉菌屬Nitrosovibrio被檢出，但豐度極低。這與廣東省國家級自然保護(hù)區(qū)的研究結(jié)果不同，在廣東南嶺亞熱帶森林生態(tài)系統(tǒng)中所有的氨氧化細(xì)菌都屬于亞硝化螺菌屬Nitrosospira[31]，說明武夷山常綠闊葉林土壤的氨氧化細(xì)菌微生物的資源更加豐富。

武夷山林地和茶園中氨氧化微生物在不同的季節(jié)其群落結(jié)構(gòu)的相對豐度大不相同，分別在以下季節(jié)中達(dá)到最高的相對豐度，土壤氨氧化細(xì)菌的林地樣本中春季Nitrosomonas占比較大，夏、秋、冬季均是Unclasified Nitrosomonadaceae相對豐度占比最大，而茶園樣本中4個季節(jié)均為Unclasified Nitroso-monadaceae相對豐度占比最大;土壤氨氧化古菌中林地樣本除Others屬外林地夏、秋季UnclassifiedThaumarchaeota占比較大，而茶園夏、秋季UnclassifiedCrenarchaeota占比較大。

對中國桃園研究站[32]3種土地利用類型中（農(nóng)田、天然林、茶園）的研究結(jié)果表明，土地利用類型是引起氨氧化微生物的群落結(jié)構(gòu)和相對豐度顯著變化的主要因素之一，前人的研究結(jié)果[16， 33-36]與本文的研究結(jié)果一致。本研究也發(fā)現(xiàn)不同采集地的樣本結(jié)構(gòu)差異較大，武夷山常綠闊葉林下除春季茶園土壤樣本與夏、秋季林地古菌的相對豐度相似外，其他樣本的均無相似性，其中夏季茶園與其他土壤樣品古菌相對豐度差異性最高，Nitrosopumilus，Unclassified Candidatus，Nitrososphaera.UnclassifiedCrenarchaeota在夏季茶園土壤樣本中相對豐度最高，Unclassified Thaumarchaeota在春季林地土壤中的豐度最高。在對比相同季節(jié)不同采集地樣品后可知，土壤中Nitrosomonas與UnclassifiedNitrosomonadaceae的相對豐度與土壤利用方式具有很大關(guān)系。

不同的環(huán)境因子對氨氧化細(xì)菌和古菌的影響不同，而且不同季節(jié)的主要影響因子也存在差異。有效鉀是影響氨氧化細(xì)菌變化的主要環(huán)境因子;而銨態(tài)氮是影響氨氧化古菌變化的主要環(huán)境因子;速效磷對氨氧化古菌和細(xì)菌的影響均最小。有研究也證明[37]AOA和AOB在中國亞熱帶森林的酸性土壤中差異分布，其受土壤理化特性的影響也不同[37]。Nitrosospira菌屬與銨態(tài)氮呈負(fù)相關(guān);Nitrosomonas菌屬與銨態(tài)氮呈正相關(guān)，這與曹彥強(qiáng)[38]等的研究結(jié)果也一致，Nitrosomonas菌屬適合在高濃度銨環(huán)境下存在，而Nitrosospira菌屬在自然環(huán)境中經(jīng)常被發(fā)現(xiàn)對亞硝酸鹽的親和力較高。對于林地樣本來說，氨氧化細(xì)菌和古菌的物種結(jié)構(gòu)相似性一致，均是夏、秋、冬三季相似性高，而春季則與其他3個季節(jié)相似性較低;而對于茶園樣本來說，氨氧化細(xì)菌的物種結(jié)構(gòu)為春、夏兩季相似性高，而與秋冬兩季節(jié)相似性較低，氨氧化古菌的物種結(jié)構(gòu)為春、秋、冬三季相似性高，而夏季則與其他3個季節(jié)相似性較低。

據(jù)Lu等人[39]的研究結(jié)果也證明，森林土壤改為茶園土壤后，其主導(dǎo)的AOA菌群從隸屬于陸地分支的Group 1.1b變?yōu)殡`屬于海洋分支的Group 1.1a-associated;氨氧化微生物菌群豐度與土壤硝化勢有相關(guān)性，如硝化潛勢與隸屬于54d9-like cluster的AOA呈正相關(guān)。這一研究結(jié)果與本研究的結(jié)論相似。在武夷山在較高的硝化勢環(huán)境中，秋季林地樣本中氨氧化古菌豐度及夏季林地中氨氧化細(xì)菌豐度都更為豐富，其中，細(xì)菌中的Unclassified Nitrosomonadac-eae和古菌中的Other均有可能是在高硝化勢的土壤中形成的主要微生物，此外，也有可能是采樣時間影響了群落結(jié)構(gòu)（即各微生物的相對豐度結(jié)構(gòu)，形成利于硝化勢升高的代謝途徑），容易使菌群對土壤產(chǎn)生較高的硝化勢，而不是單一微生物引起的。

4 結(jié)論

武夷山常綠闊葉林土壤含有豐富的氨氧化微生物資源，不同土地利用類型的氨氧化微生物種類大致相似，但茶園有獨(dú)特的Nitrosovibrio亞硝化葉菌屬存在，且氨氧化微生物的相對豐度有顯著差異;同時，氨氧化細(xì)菌和古菌菌落相對豐度隨季節(jié)變化發(fā)生明顯變化。有效鉀是影響氨氧化細(xì)菌變化的主要環(huán)境因子，而銨態(tài)氮是影響氨氧化古菌變化的主要環(huán)境因子。本研究有助于增強(qiáng)對氨氧化微生物在武夷山常綠闊葉林生態(tài)系統(tǒng)氮循環(huán)中的作用的認(rèn)識，同時為武夷山自然保護(hù)區(qū)的土壤利用開發(fā)及經(jīng)營管理提供科學(xué)依據(jù)。

參考文獻(xiàn)：

[1] 李巍，劉洋，羅欽，等.武夷山常綠闊葉林土壤微生物多樣性的季節(jié)動態(tài)[J].熱帶亞熱帶植物學(xué)報，2017， 25 （2）：115-126.

LI W， LIU Y， LUO Q， et al_ Seasonal dynamics in soilmicroorgamsms diversity of evergreen broad-leaved forest in Wuyimountains， southeastern China[J]. Journal of Tropical andSubtropical Botany， 2017， 25（2）：115-126. （in Chinese）

[2] 沙麗清，孟盈，西雙版納不同熱帶森林土壤氮礦化和硝化作用研究[J].植物生態(tài)學(xué)報，2000. 24 （2）：152-156.

SHA L Q，MENG Y， et al. Nitrification and net N mineralization rateof soils under different tropical forests in Xishuangbanna， southwestChina EJl. Acta Phvtoecologica Sinica， 2000， 24 （2）：152-156. （inChinese）

[3]李檢舟，沙麗清，王君，等，云南哀牢山中山濕性常綠闊葉林土壤氮礦化季節(jié)變化[J].山地學(xué)報，2006， 24 （2）：186-192.

LI J Z，SHA L Q，WANG J， et aL Seasonal variation of soil nitrogenmineralization in a mountane moist evergreen broad-leaved forest inAilao Mountains， SW China EJl.Journal of Mountain Science， 2006，24 （2）： 186-192. （in Chinese）

[4] 李貴才，韓興國，黃建輝，等.森林生態(tài)系統(tǒng)土壤氮礦化影響因素研究進(jìn)展[J].生態(tài)學(xué)報，2001， 21 （7）：1187-1195.

LIGC，HANXG，HUANG JH，et al.A review of affecting factors ofsoil nitrogen mineralization in forest ecosystems [J].Acta EcologicaSinica. 2001， 21 （7）： 1187-1195. （in Chinese）

SONG Y N， LIN Z M， LrN J. Composition of ammonia-oxidizingbacteria and ammonia-oxidizing archaea communities in paddy soilsof different rice cultivars [J] Chinese Journal of Eco-Agriculture，2009. 17 （6）： 1211-1215. （in Chinese）

[15]宋亞娜，林智敏，紅壤稻田不同生育期土壤氨氧化微生物群落結(jié)構(gòu)和硝化勢的變化[J].土壤學(xué)報，2010， 47 （5）：987-994.

SONG Y N. LIN Z M. Changes in community structures of ammonia-oxidizers and potential nitrification rates in red paddy soil at differentgrowth stages of rice [J].Acta Pedologica Sinica. 2010， 47 （5）：987-994. （in Chinese）

[16]葉磊，祝貴兵，王雨，等，白洋淀湖濱濕地岸邊帶氨氧化古菌與氨氧化細(xì)菌的分布特性[J].生態(tài)學(xué)報，2011， 31 （8）：2209-2215.

YE L，ZHU G B，WANG Y， et al_Abundance and biodiversity ofammonia-oxidizing archaea and bacteria in littoral wetland ofBaiyangdian Lake. North China[J]. Acta Ecologica Sinica， 2011，31（8）：2209-2215. （in Chinese）

[17]隋心，張榮濤，鐘海秀，等.森林生態(tài)系統(tǒng)中主要功能微生物的研究進(jìn)展[J].中國農(nóng)學(xué)通報，2014， 30 （28）：1-5.

SUI X， ZHANG R T，ZHONG H X. et aL Research progress on mainfunctional microorgamsms in forest ecosystems [J]. ChineseAgricultural Science Bulletin， 2014， 30 （28）：1-5.（ in Chinese）

[18] XIA W W， ZHANG C X， ZENG X W， et al.Autotrophic growth ofnitrifying community in an agricultural soil [J].The ISME Journal，2011.5 （7）：1226.

[19] SEGAL L M， MILLER D N， MCGHEE R P，et al.Bacterial andarchaeal ammonia oxidizers respond differently to long-term tillageand fertilizer management at a continuous maize site [J]. Soil andTillage Research. 2017. 168：110-1 17.

[20]黃蓉，張金波，鐘文輝，等.土地利用方式對萬木林土壤氨氧化微生物豐度的影響[J].土壤。2012， 44 （4）：581-587.

HUANG R，ZHANG J B，ZHONG W H. et al_Abundances ofammonia-oxidizing prokaryotes and gross nitrification activities inforest soils under different vegetations in a natural reserve [J]. Soils，2012， 44 （4）： 581-587. （in Chinese）

[21]朱蕊，陳清，馬成倉，等，不同利用方式對內(nèi)蒙古羊草草原氨氧化微生物豐度的影響[J].草地學(xué)報，2019. 27 （2）：437-442.

ZHU R， CHEN Q，MA C C，et al. Effects of land use pattem changeson abundance of ammonia-oxidizing nucroorganisms in LeymusChinensis grassland in Inner Mongolia [J]. Acta Agrestia Sinica，2019， 27 （2）： 437-442. （in Chinese）

[22]路璐，何燕.不同林分土壤中氨氧化微生物的群落結(jié)構(gòu)和硝化潛勢差異及其驅(qū)動因子[J].南方農(nóng)業(yè)學(xué)報，2018. 49 （11）：2169- 2176.

LU L，HE Y.The difference of ammonia-oxidizing microorganismcommunities structure and nitrification potential in soils of differentforest stands and their driving factors [J]. Journal of SouthernAgriculture， 2018， 49（ 11）： 2169-2176. （in Chinese）

[23] LAMBIN E F，MEYFROIDT P.Global land use change， economicglobalization， and the looming land scarcity [J].Proceedings of theNational Academy of Sciences of the United States of America， 20u，108 （9）：3465-3472.

[24] MOORE N. ALAGARSWAMY G， PIJANOWSKI B，et aL EastAfrican food security as influenced by future climate change and landuse change at local to regional scales [J]. Climatic Change， 2012，110（ 3/4）：823-844.

[25]劉靈芝，李景云，秦嗣軍，等，“寒富”蘋果園土壤氨氧化細(xì)菌的篩選與鑒定[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報，2015. 46 （4）：486-491.

LIU L Z，LI J Y， QIN S J，et al.Screening and identification of soilammonia-oxidizing bacteria from "Hanfu" apple orchard [J]. Journalof Shenyang Agricultural UniversitV， 2015，46 （4）：486-491. （inChinese）

[26] ZULKARNAEN N. ZHANG Y， ZHANG P，et al.Abundance ofAOA， AOB. nirS. nirK. and nosZ in red soil of China under differentland use[J]. IOP Conference Series： Earth and EnvironmentalScience. 2019， 393： 12007.

[27]鮑士旦.土壤農(nóng)化分析[M].3版.北京：中國農(nóng)業(yè)出版社，2000： 40， 81.

[28]辛亮，武傳東，曲東.長期施肥對旱地土壤中氨氧化微生物豐度和分布的影響[J].西北農(nóng)業(yè)學(xué)報，2012， 21 （6）：41-46.

XIN L，WU C D，QU D.Long-term fertilization determiningammonia-oxidizing organism abundance and distribution in dryhighland soil of loess plateau [J]. Acta Agriculturae Boreali-Occidentalis Sinica， 2012，2l（6）：41-46.（ in Chinese）

[29]宋三多，陳強(qiáng)，熊璐，等.不同沼肥處理對小麥分蘗期土壤硝化作用強(qiáng)度及氨氧化細(xì)菌和古菌群落的影響[J].麥類作物學(xué)報，2016.36 （1）：111-119.

SONG S D， CHEN Q， XIONG L，et al.Effect of biogas manure onsoil nitrification intensity and soil ammonia-oxidizing bacteria andammonia-oxidizing Archaea communities at wheat tillering stage [J].Journal of Triticeae Crops， 2016， 36（1）：111-1 19.（ in Chinese）

[30] NIU J. KASUGA L KURISU F， et al. Abundance and diversity ofammonia-oxidizing archaea and bacteria on gianular activated carbon and their fates during drinking water purification process [J] . AppliedMicrobiology and Biotechnology， 20 1 6， 100 （2 ）： 729-742.

[31] LIU B， LI Y M. ZHANG J P， et al. Abundance and diversity ofammonia-oxidizing microorganisms in the sediments of Jinshanlake EJl . Current Microbiology， 2014， 69 （ 5 ）： 75 1-757.

[32] LI M. CAO H L， HONG Y G. et aL Spatial distribution andabundances of ammonia-oxidizing archaea （AOA） and ammonia-oxidizing bacteria （AOB） in mangrove sediments EJl. AppliedMicrobiology and Biotechnology， 201 1， 89 （4 ）： 1243-1254.

[33] ADAIR K L， SCHWARTZ E. Evidence that ammonia-oxidizingarchaea are more abundant than ammonia-oxidizing bacteria insemiarid soils of northern Arizona， USA [J]. Microbial Ecology，2008. 56 （3） ： 420-426.

[34]CHEN Z M. Ecosystem functions in Nanling National Nature Reserveand protective measures concemed [J]. Shaanxi Forest Science &Technology， 2012， 61-63： 19.

[35] QIN H L. YUAN H Z. ZHANG H， et al. Ammonia-oxidizing archaea are more important than ammonia-oxidizing bacteria in nitrificationand NO3～N loss in acidic soil of sloped land[J]Biology and Fertilityof Soils，2013，49（6）：767—776.

[36]JORDAN F L，CANTERA J J L，F(xiàn)ENN M E，et al.Autotmphicammonia-0xidizing bacteria COntribute minimally to nitriflCatiOn in anitrogen—impacted forested ecosystem[J].Applied and Eenvironmental，Microbiology，2005，71（1）：197—206.

[37]HAYDEN H L，DRAKE J，IMHOF M，et al.The abundance ofnitr09en CyCle genes AmoA and nim dependS on land-uSeS and S0iltypes in South-Eastem Australia[J].Soil Biology and Biochemistry，2010.42（lO）：1774—1783.

[38]曹彥強(qiáng)，閆小娟，羅紅燕，等.不同酸堿性紫色土的硝化活性及微生物群落組成[J].土壤學(xué)報，2018，55（1）：194—202.

CAO Y Q，YAN X J，LUO H Y，et al.Nitrification activity andrnicrobial community structure in purple soils with differem pH[J]Acta Pedologica Sinica，2018，55（1）：194—202.（in Chinese）

[39]LU L，HAN W Y，ZHANG J B，et al.Nitrification of archaealammonia oxidizers in acid soils is supported by hydr01ysis ofurea[J].The ISME Journal 2012.6（10）：1978—1984.

（責(zé)任編輯：張梅）

作者簡介：李巍（1984-），女，碩士，助理研究員，主要從事土壤化學(xué)與植物營養(yǎng)研究（E-mail： liwei6055@126.com）

通信作者：林香信（1972-），男，副研究員，主要從事農(nóng)產(chǎn)品質(zhì)量安全研究（E-mail： 1353645173@qq.com）