長期施加鉀肥對小麥產量、氮磷鉀吸收量及根際土壤微生物多樣性的影響

張翔 黃媛媛 宋聰 賈振華 賈良良 宋水山

摘要：為了研究長期施加鉀肥對小麥產量、氮磷鉀吸收量以及小麥根際土壤微生物多樣性的影響，利用高通量測序技術，對河北省農林科學院辛集市馬蘭農場鉀肥定位試驗田的8個小麥根際土壤樣品進行測序，得到2組樣品的微生物群落組成信息，結合河北省農林科學院資源環境研究所提供的小麥產量和品質數據，進行相關性分析。結果表明，在NP組（不施加鉀肥）優勢細菌菌群集中在疣微菌門（Verrucomicrobia），優勢真菌菌群集中在毛殼菌科（Chaetomiaceae）、錐蓋傘屬（Conocybe）以及被孢霉目 （Mortierellales）。而NPK組（施加鉀肥）的優勢細菌菌群集中在假諾卡氏菌科（Pseudonocardiaceae）、Terrimonas，優勢真菌菌群集中在光柄菇科（Pluteaceae）、花褶傘屬（Panaeolus）。在細菌屬水平上，2組間差異較大的細菌菌屬有15種，其中NP組顯著增多的是Terrabacter、長繩菌屬（Longilinea）、硝化螺旋菌屬（Nitrospira）、狹義梭菌屬（Clostridium sensu stricto）、Mangroviflexus、Rhizorhabdus、Tahibacter，而NPK組顯著增多的是Methyloceanibacter、Terrimonas、德克斯氏菌屬（Derxia）、倫茨氏菌屬（Lentzea）、Chryseolinea、間孢囊菌屬（Intrasporangium）、Maricaulis、纖維微桿菌屬（Cellulosimicrobium）。在真菌屬水平上，差異顯著的只有NP組中較多的毛殼菌屬（Chaetomium）和球腔菌屬（Phaeosphaeria）。結合小麥產量和養分數據可以得出，鉀肥可以提高小麥的產量和氮磷鉀吸收量，而參與小麥吸收鉀元素的菌群應該有Ohtaekwangia、Gp3、其他菌屬和一些未知細菌菌屬，以及毛球殼科的一個未知菌屬（unclassified_Lasiosphaeriaceae）、未知真菌菌屬（unclassified_Fungi）和地星科的一個未知菌屬（unclassified_Geastraceae）等真菌菌屬。

關鍵詞：鉀肥；小麥根際土壤；微生物多樣性；高通量測序；冗余分析

中圖分類號：S512.106? 文獻標志碼：A

文章編號：1002-1302（2024）03-0254-07

土壤是農業的根基，是農業生產重要的組成部分。土壤中的微生物更是參與了農作物生長發育的各個階段［1］。土壤微生物主要包括真菌、細菌、古菌、病毒、原生動物以及藻類［2］。它們在土壤的碳循環、氮循環、硫循環以及有機質的分解和轉化過程中扮演著重要的角色［3］。土壤微生物種類繁多、功能多樣，是衡量土壤健康質量的重要指標［4］。在農業研究上，耕地土壤微生物已經受到了很多的關注，其中利用高通量二代測序技術了解土壤微生物多樣性已經得到了廣泛的應用［5］。該技術可以檢測出土壤中的所有微生物的相對豐度，通過對土壤中所有微生物的相對豐度可以分析出耕地土壤中的微生物群落的組成情況，對耕地土壤的改良修復和肥料的施加提供科學的數據指導［6-7］。

小麥是重要的糧食作物，在我國有著廣泛的種植［8-9］。在小麥生長發育過程中鉀是必不可少的營養元素之一，參與了小麥體內各種生理生化過程［10］。在農業生產中施用鉀肥是解決土壤中鉀元素不足以及提高農作物產量和品質的重要途徑［11］。研究表明，在氮磷肥的基礎上施加不同量的鉀肥均可以提高小麥的產量［12］。隨著我國經濟水平的提高、人口的增長，對糧食的需求也日趨提高［13］。我國已經成為世界主要的鉀肥消費國［14］。所以，合理地利用鉀肥既可以降低農業生產成本，又可以節省鉀礦資源［15］，對我國農業可持續發展有著重要的意義［16］。本研究利用高通量二代測序技術對河北省農林科學院馬蘭農場試驗站的鉀肥長期定位試驗田的土壤微生物進行測序，分析研究土壤微生物與鉀肥施加量和產量之間的關系，以期為小麥的科學施肥和耕地土壤修復改良提供科學依據。

1 材料與方法

1.1 試驗概況

試驗于2018年在河北省辛集市馬蘭村馬蘭農場進行，地處115.201°～115.206°E、37.995°～37.997°N，屬暖溫帶半濕潤季風大陸氣候，海拔在25.0～37.8 m之間，年平均氣溫10～20 ℃，年平均降水量525.7 mm，年平均日照時數2 554.9 h。

馬蘭農場試驗站鉀肥定位試驗由河北省農林科學院資源環境研究所設計，試驗田從1992年開始小麥—玉米輪作，共5個處理，本次試驗取其中2個處理作為試驗對象，每個處理4次重復，小區面積為50 m2（0.075畝），秸稈全部清除。各處理同量施用氮磷肥，每季施N 225 kg/hm2、P2O 5 90 kg/hm2。施加鉀肥的處理底肥施150 kghm2 K2O 。田間施肥情況見表1。

1.2 樣品的采集

本次試驗所用土壤樣品及小麥樣品均在2018年6月采集完成。

本研究所用土壤樣品采用“5點采樣法”采集，在每塊地的4個角以及中間位置進行土壤樣品的采集，然后進行混樣。本研究所用土壤樣品為小麥的根際土壤，根際土壤距離小麥根部最近，土壤菌群與小麥根部交流更頻繁。采集時先將小麥根部的大塊土壤抖掉，然后用軟毛刷將附著在小麥根系表面的土壤刷到取樣袋中。將所取土壤樣品放入冰盒中帶回實驗室，存入-80 ℃冰箱中待測。

1.3 土壤微生物總DNA提取及PCR擴增

取500 mg土壤樣品放入滅過菌的2 mL離心管中，然后用OMEGA公司的E.Z.N.ATM Mag-Bind Soil DNA Kit提取試劑盒對土壤樣品進行DNA提取。土壤DNA提取后，利用引物341F（5′-CCTACGGGNGGCWGCAG-3′）和805R（5′-GACTACHVGGGTATCTAATCC-3′）對土壤中細菌16S rDNA的V3～V4區域進行擴增。同時，利用引物ITS1F（5′-CTTGGTCATTTAGAGGAAGTAA-3′）和ITS2R（5′-GCTGCGTTCTTCATCGATGC-3′）對土壤中真菌ITS區域的ITS1～ITS2區域進行擴增。

1.4 上機測序

將PCR產物送到生工生物工程（上海）股份有限公司高通量測序部，利用Illumina MiSeq高通量測序平臺對樣品進行測序。

1.5 數據處理

MiSeq測序序列中含有樣本標簽序列，以及測序時加入的引物和接頭序列。首先需要去除引物接頭序列，再根據測序序列之間的重疊關系，將成對的序列拼接成一條序列，然后按照樣品標簽序列識別并區分樣品得到各樣本數據，最后對各樣本數據的質量進行質控過濾，得到各樣本的有效數據［17］。

使用USEARCH去除預處理后序列中非擴增區域序列，而后對序列進行測序錯誤校正，并調用UCHIME進行鑒定嵌合體。隨后再將去除嵌合體的序列與數據庫代表性序列進行BLASTn比對，低于閾值的比對結果通常被認為是靶標區域外的序列，剔除掉該部分序列。

接下來，利用QIIME、UCLUST、USEARCH軟件按照97%相似性將全部優質序列進行聚類，得到代表序列和操作分類單元（OTU）表。利用Mothur等軟件進行稀釋性曲線的繪制，利用QIIME等軟件計算Chao1指數、observed species指數、goods coverage指數和Shannon指數。為了得到每個OTU的物種信息，采用RDP Classifier算法對每個OTU代表序列進行比對分析，并在各個水平注釋群落物種信息。利用R語言繪制物種組成分析和主成分分析（PCA）的圖表［18］。并根據各個樣品的物種組成信息進行LEfSe分析。

1.6 關聯性分析

利用冗余分析（RDA）法分析河北省農林科學院資源環境研究所提供的小麥相關數據與高通量數據的關聯性，篩選出與鉀肥關系密切的菌株。

2 結果與分析

2.1 測序結果分析

將8個樣品的下機數據經過預處理、去除嵌合體后，剩下的優質序列為分析所用的最終序列［19］。如表2所示，8個樣品中細菌的16S rDNA V3～V4區域的最終序列一共有557 581條，平均每個樣品有7萬條左右。OTU總數為51 801個，平均每個樣品有6 475個。如表3所示，8個樣品中真菌的 ITS1～ITS2區域的最終序列一共有514 125條，平均每個樣品有64 265條。OTU總數為6 387個，平均每個樣品有798個。從這2個表中可以看出，在測序深度相當的情況下，土壤中的細菌數量和種類遠遠大于真菌的數量和種類。

2.2 測序深度分析

Shannon指數是反映土壤微生物多樣性的指數。由圖1和圖2可以看出，隨著測序深度的增加各個樣品的曲線趨于平緩，說明隨著測序深度的增加，新增微生物種類明顯減少，所以當前的測序深度已經可以足夠反映土壤中所有細菌和真菌的真實情況。

2.3 2組間PCA結果

PCA是利用降維方法分析各個土壤樣品中菌群的相似性。PCA中距離越近說明樣品間的相似性越高。圖3和圖4分別展示了這8個樣品中細菌和真菌的PCA三維圖。從圖中可以看出，NP組和NPK組2組樣品間的距離要遠遠大于組內的距離，說明它們的組間的菌群差異要大于組內差異，也就說明了這樣的分組方式有意義，可以進一步進行組間菌群差異的分析。

2.4 2組間LEfSe結果

LEfSe是對2組間差異菌群進行分析的方法。如5和圖6所示，分別是2組間細菌和真菌差異菌群的LEfSe分支圖，其中紅色是NP組的優勢菌群，綠色是NPK組的優勢菌群?？梢钥闯?，NP組的優勢細菌菌群集中在了疣微菌門（Verrucomicrobia），優勢真菌菌群集中在了毛殼菌科（Chaetomiaceae）、錐蓋傘屬（Conocybe）以及被孢霉目（Mortierellales）。

而NPK組的優勢細菌菌群集中在了假諾卡氏菌科（Pseudonocardiaceae）、Terrimonas，優勢真菌菌群集中在了光柄菇科（Pluteaceae）、花褶傘屬（Panaeolus）。

2.5 2組間屬水平差異微生物菌群分析

屬水平是高通量測序技術能達到的測量較準確的一個分類水平。圖7和圖8展示的是2組樣品中土壤微生物在屬水平上細菌和真菌相對豐度差異較大的菌屬，其中紅色部分是相對豐度有顯著差異的菌屬。從圖7中可以看出，2組間差異較大的細菌菌屬有15種，其中NP組顯著增多的是地桿菌屬（Terrabacter）、長繩菌屬（Longilinea）、硝化螺旋菌屬（Nitrospira）、狹義梭菌屬（Clostridium sensu stricto）、Mangroviflexus、Rhizorhabdus、Tahibacter，而NPK組顯著增多的是Methyloceanibacter、Terrimonas、德克斯氏菌屬（Derxia）、倫茨氏菌屬（Lentzea）、Chryseolinea、間孢囊菌屬（Intrasporangium）、Maricaulis、纖維微桿菌屬（Cellulosimicrobium）。圖8展示了2組間屬水平上差異較大的真菌菌群，其中差異顯著的只有NP組中較多的毛殼菌屬（Chaetomium）和球腔菌屬（Phaeosphaeria）。

產生這些差異菌屬的原因，可能是由于長期的鉀肥施加使土壤中微生物群落的結構發生了改變。因為微生物群落之間存在著聯系，所以鉀肥施加后，一些能夠利用、分解鉀肥的微生物以及其他相關微生物就會聚集。

2.6 小麥產量、氮磷鉀吸收量與土壤微生物關聯性分析

表4中展示的是2組處理小麥籽粒產量和小麥氮磷鉀吸收量，該數據由河北省農林科學院資源環境研究所提供。從表4中可以看出，施加鉀肥后，小麥的籽粒產量和小麥中氮磷鉀的成分都有了明顯的提高，說明鉀肥對小麥的生長發育和小麥的產量和品質都有一定的影響。

圖9和圖10展示了屬水平上細菌和真菌相對豐度與小麥產量和氮磷鉀吸收量的相關性。從圖9可以看出，鉀肥的施加與小麥籽粒產量和氮磷鉀的吸收量成呈正相關趨勢，同時也與Ohtaekwangia、Gp3、其他菌屬以及一些未知細菌菌屬相對豐度呈正相關的趨勢。未知菌屬中是一些目前沒有命名的菌屬，這些菌屬中應該有與鉀肥有密切聯系的菌屬。由圖10可知，鉀肥的施加與毛球殼科的一個未知菌屬（unclassified_Lasiosphaeriaceae）、未分類真菌菌屬（unclassified_Fungi）以及地星科的一個未知菌屬（unclassified_Geastraceae）等真菌菌屬相對豐度呈正相關的趨勢。與鉀肥施加呈正相關的這一系列細菌和真菌有可能參與了鉀肥的分解以及小麥對鉀的吸收過程。

3 討論與結論

本試驗對河北省農林科學院辛集市馬蘭農場試驗站的鉀肥定位試驗田的土壤進行采集和高通量測序分析。本次試驗的土壤樣品采集了來自鉀肥定位試驗的2個處理，分別是NP組和NPK組。NP組在整個小麥生長發育過程中不施加鉀肥，NPK組小麥生長發育過程中施加適量的鉀肥。本試驗一共采集了8個小麥根際土壤樣品，并利用Illumina MiSeq高通量測序平臺對這8個樣品中的細菌16S rDNA的V3～V4區域以及真菌ITS1～ITS2區域進行測序。其中細菌得到557 581條最終序列，真菌得到514 125條最終序列。通過聚類細菌得到 51 801 個OTU，真菌得到6 387個OTU。通過稀釋曲線分析得出所有土壤樣品的測序深度已經能夠反映土壤微生物的真實情況。通過PCA，說明NP與NPK組間差異大于組內差異，具有分組意義，可以進行組間菌群的差異性分析。

LEfSe是對組間差異菌群的分析，通過分析可以得出，NP組的優勢細菌菌群集中在了疣微菌門（Verrucomicrobia），優勢真菌菌群集中在了毛殼菌科（Chaetomiaceae）、錐蓋傘屬（Conocybe）以及被孢霉目（Mortierellales）。而NPK組的優勢細菌菌群集中在了假諾卡氏菌科（Pseudonocardiaceae）、Terrimonas，優勢真菌菌群集中在了光柄菇科（Pluteaceae）、花褶傘屬（Panaeolus）。進一步分析屬水平的差異菌群可以得出，NP組顯著增多的細菌菌屬是Terrabacter、Longilinea、Nitrospira、Clostridium sensu stricto、Mangroviflexus、Rhizorhabdus、Tahibacter，而NPK組顯著增多的細菌菌屬是Methyloceanibacter、Terrimonas、Derxia、Lentzea、Chryseolinea、Intrasporangium、Maricaulis、Cellulosimicrobium。在真菌屬水平上的差異菌群只有NP組中較多的毛殼菌屬（Chaetomium）和球腔菌屬（Phaeosphaeria）。這些差異菌群產生的原因可能是小麥田長期施加鉀肥造成的，這些差異菌群應該與土壤中鉀元素的轉化分解有關。

結合河北省農林科學院資源環境研究所提供的小麥產量和氮磷鉀吸收量的相關數據，可以看出鉀肥的施加可以提高小麥的產量和氮磷鉀吸收量。鉀肥對小麥的一系列影響，應該也與土壤微生物有著一定的聯系。通過RDA圖可以看出，小麥的產量和品質與Ohtaekwangia、Gp3、其他菌屬以及一些未知細菌菌屬的相關豐度呈正相關的趨勢，也與unclassified_Lasiosphaeriaceae、unclassified_Fungi、unclassified_Geastraceae等真菌菌屬的相關豐度呈正相關的趨勢。這些菌屬應該是直接或間接地參與了小麥在整個生長發育過程中對鉀元素的吸收。

通過以上分析可以看出，鉀肥的施加量與小麥的生長發育有著重要的聯系，這與路欣春等的研究結果［20-21］相符。結合高通量測序的分析結果可以看出，土壤微生物參與了小麥吸收鉀肥的過程，并可以有效促進小麥吸收鉀元素，這與Yadav等的研究結果［22-24］相符。

綜上所述，通過對河北省農林科學院辛集市馬蘭農場鉀肥定位試驗田小麥根際土壤的采集與土壤微生物進行分析，得出的結論是鉀肥的合理施加可以提高小麥的產量和氮磷鉀吸收量，土壤中的一系列細菌和真菌直接或間接地參與了小麥對土壤中鉀元素的吸收過程。另外，還有哪些微生物參與了鉀元素的吸收過程，以及參與的具體方式還需要進一步的試驗來進行探討。

參考文獻：

［1］查美琴，徐海東，成向榮，等. 不同林齡杉木+閩楠復層林土壤磷形態及微生物功能多樣性變化［J］. 生態學報，2020，40（19）：6938-6947.

［2］王宏杰，劉紹東，劉瑞華，等. 輪作對棉花根際土壤細菌群落的影響［J］. 生物技術通報，2020，36（9）：117-124.

［3］張 翔，宋水山，黃亞麗，等. 基于高通量測序分析河北省中南部地區耕地土壤細菌多樣性［J］. 華北農學報，2018，33（4）：196-203.

［4］陳建國，張楊珠，曾希柏，等. 不同配方施肥對長期缺施鉀肥的紅壤性水稻土微生物特性的影響［J］. 植物營養與肥料學報，2008，14（6）：1200-1205.

［5］薛宇昂，郭 壯，趙慧君，等. 基于MiSeq測序技術石花酒大曲中微生物多樣性解析［J］. 中國釀造，2020，39（7）：163-168.

［6］秦 濤，劉新社. 氮鉀肥配施對土壤微生物與西瓜形態建成、品質、產量的影響［J］. 江蘇農業科學，2022，50（16）：154-161.

［7］湯明堯，沈重陽，陳署晃，等. 新疆小麥、玉米的產量和氮磷鉀肥利用效率［J］. 中國農業科學，2022，55（14）：2762-2774.

［8］李 博，常旭虹，王德梅，等. 不同土壤條件下追鉀肥對小麥產量和品質的影響［J］. 農業科技通訊，2019（11）：75-79.

［9］秦永林，石曉華，賈立國，等. 有機肥替代化肥對馬鈴薯—小麥輪作體系產量及鉀肥利用率的影響［J］. 北方農業學報，2019，47（4）：47-51.

［10］李淑靖，吳金芝，黃 明，等. 旱地條件下不同地理來源小麥品種的產量與氮磷鉀積累利用差異［J］. 華北農學報，2023，38（4）：129-140.

［11］黃雪嬌，王 菲，谷守寬，等. 鉀肥及與秸稈配施對紫色土作物產量和微生物群落結構的影響［J］. 生態學報，2018，38（16）：5792-5799.

［12］錢晨誠，陳 立，馬 泉，等. 磷鉀肥施用量和方法對弱筋小麥籽粒產量和蛋白質含量及養分吸收利用的影響［J］. 植物營養與肥料學報，2023，29（2）： 287-299.

［13］唐 旭，計小江，李超英，等. 水稻—大麥長期輪作體系鉀肥效率及土壤鉀素平衡［J］. 中國農業科學，2014，47（1）：90-99.

［14］陳 娜. 鉀肥對煙草青枯病的防控效果及其根際微生物群落的影響［J］. 安徽農業科學，2018，46（22）：125-127，134.

［15］王春寧，余俊清，陳 良，等. 鉀鹽資源全球分布和我國找鉀實踐及方法探究［J］. 鹽湖研究，2007，15（3）：56-72.

［16］李丙智，王桂芳，秦曉飛，等. 沼液配施鉀肥對果園土壤理化特性和微生物及果實品質影響［J］. 中國農業科學，2010，43（22）：4671-4677.

［17］楊 敏，殷 絨，張國濤，等. 基于高通量測序技術的香格里拉葡萄酒產區根際微生物多樣性研究［J］. 云南農業大學學報（自然科學），2020，35（3）：392-400.

［18］彭玉嬌，崔學宇，譚夢超，等. 不同立地條件沙田柚果園土壤微生物多樣性分析［J］. 南方農業學報，2020，51（5）：1136-1144.

［19］楊 慧，路風銀，韓俊豪，等. 高通量測序技術分析鮮食花生貯藏過程中真菌群落的的組成和變化［J］. 河南農業科學，2023，52（2）：171-180.

［20］陸欣春，韓曉增，鄒文秀. 作物高效利用土壤磷素的研究進展［J］. 土壤與作物，2013，2（4）：164-172.

［21］張 磊，張維樂，魯劍巍，等. 秸稈還田條件下不同供鉀能力土壤水稻、油菜、小麥鉀肥減量研究［J］. 中國農業科學，2017，50（19）：3745-3756.

［22］Yadav R L，Dwivedi B S，Prasad K，et al. Yield trends，and changes in soil organic-C and available NPK in a long-term rice-wheat system under integrated use of manures and fertilisers［J］. Field Crops Research，2000，68（3）：219-246.

［23］劉 雪. 腐殖酸復合微生物肥料對土壤微生物及蔬菜生長的影響［J］. 新農業，2020（13）：41-42.

［24］王素萍，杜 雷，陳 鋼，等. 有機無機鉀肥配施對甜玉米產量、品質及土壤微生物的影響［J］. 湖北農業科學，2018，57（23）：55-60.