王曉丹 馬 花 官 梅 官春云

(湖南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院/國(guó)家油料改良中心湖南分中心，湖南 長(zhǎng)沙 410000)

20世紀(jì)中后期以來(lái)，我國(guó)化肥工業(yè)快速發(fā)展，促使我國(guó)糧食產(chǎn)量不斷增加[1]。然而化肥的大量施用造成了江河湖泊水域富營(yíng)養(yǎng)化[2]，土壤板結(jié)、有機(jī)碳含量下降，農(nóng)田氮素遷移[3]，臭氧層破壞，自然災(zāi)害頻發(fā)[4]等一系列問(wèn)題，使可持續(xù)發(fā)展面臨巨大的壓力和挑戰(zhàn)[5]。作為一種綠色防控技術(shù)，種植綠肥可減少我國(guó)化肥農(nóng)藥使用量[6]。因此培育優(yōu)質(zhì)綠肥品種，發(fā)展綠肥產(chǎn)業(yè)對(duì)保護(hù)生態(tài)環(huán)境及農(nóng)業(yè)綠色可持續(xù)發(fā)展具有重要意義。

綠肥油菜作為近年來(lái)迅速興起的有機(jī)肥，其種植區(qū)域分布廣泛，適應(yīng)性強(qiáng)[7]，且油菜干物質(zhì)積累量大，翻壓可提高土壤養(yǎng)分含量、增強(qiáng)土壤肥力、改善土壤物理性狀[8-9]。研究表明，油菜種植可改善細(xì)菌群落多樣性，增加土壤有益菌群落數(shù)量[10]。油菜的碳氮比適中，較紫云英更有利于被土壤微生物分解，養(yǎng)分易被下茬作物吸收[6,11]；油菜中硫苷水解物可抑制雜草[12]、土壤病原菌[13]及土壤線(xiàn)蟲(chóng)[14]等，是區(qū)別于其他綠肥作物的重要特征之一。同時(shí)，綠肥油菜可結(jié)合不同的農(nóng)作制度需求，如長(zhǎng)江流域游的“稻-稻-油”[15]、“蓮-油”輪作[16]等，具有獨(dú)特的應(yīng)用價(jià)值和廣闊發(fā)展前景。

養(yǎng)分特性是衡量綠肥作物優(yōu)劣的重要指標(biāo)，也是決定綠肥對(duì)土壤和后茬作物影響大小的重要基礎(chǔ)[17]。研究表明，同一生態(tài)區(qū)不同綠肥作物的養(yǎng)分特性及綠肥作物對(duì)土壤理化性質(zhì)與微生物的影響差異明顯[18]。本研究以?xún)蓚€(gè)綠肥油菜品種油肥1號(hào)[19]、油肥2號(hào)[20]為對(duì)照，檢測(cè)12個(gè)油菜品系植株養(yǎng)分含量、土壤養(yǎng)分含量及土壤微生物類(lèi)群，比較品種間差異，以期篩選出優(yōu)質(zhì)綠肥種質(zhì)資源，促進(jìn)綠肥油菜育種。

1 材料與方法

1.1 試驗(yàn)地概況

試驗(yàn)地點(diǎn)位于湖南省長(zhǎng)沙市湖南農(nóng)業(yè)大學(xué)耘園基地(113°07′E、28°18′N(xiāo))，屬亞熱帶季風(fēng)濕潤(rùn)氣候，年平均氣溫17.2℃，年積溫為5 457℃，年均降水量 1 361.6 mm。 試驗(yàn)田施行“稻-油”輪作，土壤pH值為5.6，有機(jī)質(zhì)、全氮、全磷、速效鉀、速效磷的含量分別23.80、1.40、0.70、0.28、0.01 g·kg-1。

1.2 試驗(yàn)設(shè)計(jì)

試驗(yàn)采用田間隨機(jī)區(qū)組設(shè)計(jì)，油肥1號(hào)、油肥2號(hào)和12個(gè)油菜品種(系)見(jiàn)表1。小區(qū)面積為2×4 m2，種植密度為15萬(wàn)株·hm-2，每個(gè)品系設(shè)3次重復(fù)，共42個(gè)小區(qū)。于2018年9月27日播種，條播，田間管理按大田油菜栽培方法進(jìn)行。

表1 供試油菜品種(系)Table 1 List of tested oilseed rape varieties (lines)

1.3 試驗(yàn)方法

在作物苗期、盛花期分別進(jìn)行植株和土壤取樣，參照王營(yíng)營(yíng)[21]的方法測(cè)定植株指標(biāo)(全氮、全磷、全鉀含量)與土壤養(yǎng)分指標(biāo)(全氮、速效鉀含量)。土壤樣品委托北京諾禾致源生物信息科技有限公司，基于IonS5TMXL測(cè)序平臺(tái)，利用單端測(cè)序(Single-End)的方法進(jìn)行16S(V3+V4)區(qū)域擴(kuò)增子測(cè)序。

1.4 數(shù)據(jù)分析

采用SPSS 24.0軟件進(jìn)行相關(guān)性與主成分分析，并用Origin 2019軟件作圖，差異顯著性檢驗(yàn)采用LSD法，P<0.05時(shí)認(rèn)為差異達(dá)到顯著水平。高通量數(shù)據(jù)的生物信息學(xué)分析通過(guò)諾和云平臺(tái)進(jìn)行。

2 結(jié)果與分析

2.1 植株養(yǎng)分特性差異

14個(gè)油菜品種(系)的氮(N)、磷(P)、鉀(K)含量測(cè)定結(jié)果見(jiàn)表2，相比對(duì)照，多數(shù)品系苗期養(yǎng)分含量差異顯著。在氮含量方面，品系編號(hào)為4、10、12、11、9、7的植株苗期全氮含量均高于CK1和CK2；在磷含量方面，除品系5外，其他品系的苗期全磷含量均高于CK1和CK2，其中品系6全磷含量顯著高于其他品系(P<0.05)；在鉀含量方面，品系10、9、4、12苗期全鉀含量顯著高于CK1和CK2。相比對(duì)照，多數(shù)品系盛花期養(yǎng)分含量存在顯著差異。在氮含量方面，品系11、13、4、8、5植株盛花期全氮含量均高于CK2；品系8植株盛花期全磷含量最高，除品系7和1外，其他參試品系均高于CK1；品系1和9盛花期全鉀含量高于CK1和CK2，其中品系1顯著高于其他品種(系)(P<0.05)?？傮w而言，品系4和11植株養(yǎng)分含量較高。

14個(gè)油菜品種(系)氮(N)、磷(P)、鉀(K)單株積累量測(cè)定結(jié)果表明(表2)，苗期品系5植株的氮、磷、鉀積累量均顯著高于其他品種(系)(P<0.05)；盛花期時(shí)，品系11、4、10、8、2、3、12、5的全氮單株積累量高于CK1和CK2，品系11號(hào)顯著高于其他品種(系)(P<0.05)，品系8、10、11、3、2、4、5的全磷單株積累量較高，品系8高于其他品種(系)(P<0.05)，品系11、4、2、1、3、9的全鉀單株積累量高于CK1和CK2，品系11號(hào)顯著高于其他品種(系)(P<0.05)。結(jié)合植株養(yǎng)分含量分析表明，油菜養(yǎng)分含量苗期普遍比盛花期高，平均養(yǎng)分積累量則是苗期普遍低于盛花期。以盛花期養(yǎng)分積累量作為篩選指標(biāo)，品系4、11、3、5表現(xiàn)較優(yōu)。

表2 不同品種(系)植株養(yǎng)分含量與養(yǎng)分積累量差異Table 2 Nutrient contents and nutrient accumulation of per plant in different varieties

表3 不同品種(系)土壤養(yǎng)分含量差異Table 3 difference of soil nutrient content among different varieties (lines)

2.2 不同生育期/不同品種(系)間土壤養(yǎng)分及微生物特性

2.2.1 不同生育期/不同品種(系)間土壤養(yǎng)分含量差異 14個(gè)油菜品種(系)的土壤養(yǎng)分含量測(cè)定結(jié)果表明(表3)，在苗期，品系11土壤氮含量，品系3、10、4、1、2、9土壤速效磷含量，品系4、11、6、1土壤速效鉀含量分別高于CK1和CK2中土壤氮、速效磷、速效鉀含量，其中品系4土壤中速效鉀含量顯著高于兩對(duì)照(P<0.05)。盛花期，品系10、12、5土壤全氮含量較高，與CK1差異顯著(P<0.05)；品系12、2、1土壤速效磷含量顯著高于CK1和CK2(P<0.05)；品系11、6、9、10、4、5土壤速效鉀含量較高，且品系11顯著高于CK1和CK2(P<0.05)。多數(shù)品系土壤氮含量在苗期和盛花期以及品種(系)間無(wú)顯著差異，不同品種(系)土壤全氮、速效磷和速效鉀含量在苗期和盛花期的變化情況不一致，整體而言，品系10、11土壤養(yǎng)分含量較高。

2.2.2 不同生育期/不同品種(系)間土壤微生物差異 對(duì)14個(gè)品種(系)苗期與盛花期土壤樣本進(jìn)行16S測(cè)序，質(zhì)控后按照97%相似性進(jìn)行操作分類(lèi)單元(operational taxonomic unit，OTU)聚類(lèi)分析(圖1-A、B)，結(jié)果顯示苗期共享OTU數(shù)為1 658，占OTU總數(shù)的43.96%，花期共享OTU數(shù)為1 480，占OTU總數(shù)的40.29%，表明各品種(系)的種群具有一定差異。根據(jù)SILVA數(shù)據(jù)庫(kù)(https://www.arb-silva.de/)注釋結(jié)果，苗期各品種(系)土壤細(xì)菌在門(mén)(Phylum)水平上豐度排名前10的物種(圖1-C)分別為變形桿菌(Proteobacteria)、酸桿菌(Acidobacteria)、放線(xiàn)菌(Actinobacteria)、綠彎菌(Chloroflexi)、擬桿菌門(mén)(Bacteroidetes)、厚壁菌(Firmicutes)、芽單胞菌(Gemmatimonadetes)、疣微菌(Verrucomicrobia)、硝化螺旋菌(Nitrospirae)、藍(lán)細(xì)菌(Cyanobacteria)，盛花期其中八類(lèi)與苗期一致，另外兩類(lèi)為奇古菌(Thaumarchaeota)與未知門(mén)類(lèi)微生物(圖1-D)?？傮w而言，同一時(shí)期微生物門(mén)水平的組成無(wú)明顯差異，不同時(shí)期微生物門(mén)水平的組成差異不大，但綠肥的不同品種(系)和種植的不同時(shí)期會(huì)影響不同門(mén)類(lèi)微生物的相對(duì)含量。

2.3 油菜植株與土壤養(yǎng)分性質(zhì)以及土壤微生物相關(guān)性分析與主成分分析

2.3.1 油菜植株與土壤養(yǎng)分性質(zhì)以及土壤微生物相關(guān)性分析 植株養(yǎng)分、土壤微生物和土壤養(yǎng)分間相關(guān)性分析結(jié)果表明(表4)，在苗期和盛花期，植株養(yǎng)分含

注：A、B分別為苗期(M)、花期(H)不同品種(系)微生物群落種類(lèi)Venn圖； C、D分別為苗期、花期不同品種(系)在門(mén)水平的差異。Note: A and B are microbial community types of varieties (lines) at seedling stage and flowering stage respectively(Venn graph). C and D are the differences of varieties (lines) at seedling stage and flowering stage respectively(at the phylum level).圖1 不同生育期/不同品種(系)間土壤微生物差異Fig.1 Differences of soil microorganisms in different growth stages/varieties (lines)

表4 油菜植株養(yǎng)分和主要微生物類(lèi)群與土壤養(yǎng)分性質(zhì)的相關(guān)系數(shù)Table 4 Correlation coefficients of plant nutrients and main microbial groups with soil nutrient properties

量與土壤養(yǎng)分含量均無(wú)顯著相關(guān)性，單株氮含量與土壤全氮含量分別存在極顯著和顯著正相關(guān)，同時(shí)，盛花期單株磷含量與土壤磷含量之間呈顯著正相關(guān)。苗期與盛花期的部分微生物類(lèi)群與土壤全氮、速效磷、速效鉀養(yǎng)分含量之間呈現(xiàn)正或負(fù)相關(guān)性，尤其是厚壁菌門(mén)與疣微菌門(mén)，在苗期與盛花期與土壤全氮含量均有顯著負(fù)相關(guān)性。由此可知，植株單株養(yǎng)分含量與土壤微生物是影響土壤養(yǎng)分含量的重要因素。

2.3.2 油菜植株與土壤養(yǎng)分性質(zhì)以及土壤微生物主成分分析 據(jù)2.3.1分析，植株單株養(yǎng)分積累量與部分微生物類(lèi)群有正或負(fù)相關(guān)性，將各品種(系)植株單株養(yǎng)分含量和土壤養(yǎng)分含量及土壤厚壁菌和疣微菌類(lèi)群標(biāo)準(zhǔn)化后進(jìn)行主成分分析。表5列出了前5個(gè)主成分，可代表14個(gè)油菜品種(系)82.87%的養(yǎng)分信息。第一主成分的貢獻(xiàn)最大，達(dá)28.85%，苗期土壤N含量及植株養(yǎng)分積累量，占較大正向負(fù)荷；第二主成分占20.12%，盛花期養(yǎng)分積累量占較大正向負(fù)荷；第三主成分中土壤速效鉀含量占較大正向負(fù)荷；第四、第五主成分占較大正向負(fù)荷的是土壤微生物含量?？傮w來(lái)看，第一、第二主成分為植株養(yǎng)分特性因子，第三、四、五因素主要為土壤養(yǎng)分與微生物因子，其中植株養(yǎng)分積累量是綠肥油菜種質(zhì)資源養(yǎng)分特性評(píng)價(jià)的主導(dǎo)因子。

將各品種(系)植株和土壤養(yǎng)分特性標(biāo)準(zhǔn)化數(shù)值與前5個(gè)主成分的特征向量相乘，并加權(quán)得到各主成分得分及排名(表6)。與兩對(duì)照品種相比，在植株養(yǎng)分特性方面，以品系11、4、10、2表現(xiàn)較優(yōu)，在土壤養(yǎng)分及微生物特性方面，以品系11、6、4、1較優(yōu)，整體分析來(lái)看，品系11、4優(yōu)于CK1和CK2。

3 討論

綠肥養(yǎng)分特性對(duì)土壤理化性質(zhì)和后茬作物生長(zhǎng)具有重要意義。本研究中苗期和盛花期油菜植株養(yǎng)分含量及單株積累量分析結(jié)果表明，盛花期N、P、K含量均低于苗期，與盛花期植株部分養(yǎng)分由莖葉向生殖器官轉(zhuǎn)移有關(guān)[22]。而盛花期單株養(yǎng)分積累量高于苗期，表明在油菜營(yíng)養(yǎng)生長(zhǎng)過(guò)程中，單株鮮草產(chǎn)量不斷增加，進(jìn)而使單株養(yǎng)分含量不斷積累，這與王佳銳[23]的研究結(jié)果一致。對(duì)于翻壓綠肥，翻壓養(yǎng)分總量即養(yǎng)分累積量比植株養(yǎng)分含量更為重要，且油菜在盛花期進(jìn)行翻壓可提供更多植株養(yǎng)分，以盛花期養(yǎng)分積累量作為篩選指標(biāo)，則品系4、11、3、5表現(xiàn)較優(yōu)。另一方面，盡管苗期植株單株積累量略低，但苗期營(yíng)養(yǎng)生長(zhǎng)旺、莖葉產(chǎn)量較高，栽培時(shí)間短，在翻壓后易腐熟，以苗期養(yǎng)分積累量作為篩選指標(biāo)，則品系5、11、8表現(xiàn)較優(yōu)。

表5 植株和土壤養(yǎng)分特性的主成分分析Table 5 Principal component analysis of plant and soil nutrient characteristics

表6 基于主成分分析的綠肥油菜植株和土壤養(yǎng)分與微生物特性的綜合比較Table 6 Comprehensive comparison of nutrient characteristics via principal component analysis

油菜在不同生育期對(duì)土壤養(yǎng)分的吸收和利用率不同[24]，不同品種生長(zhǎng)情況亦對(duì)土壤養(yǎng)分有影響。研究表明，油菜對(duì)氮素的吸收量在抽薹前、抽薹至終花、角果發(fā)育期分別約占45%、45%、10%[25]，同時(shí)氮肥利用率低在油菜生產(chǎn)中亦普遍存在[26]。本研究中，植株氮含量在苗期和盛花期均較高，而相應(yīng)的土壤氮含量在苗期和盛花期以及品種(系)間均變化差異較小，表明油菜植株需氮量較大，但氮吸收率較低。盡管如此，油菜所吸收的80%以上養(yǎng)分能夠以各種形式還田[25]。肖小軍等[27]研究表明，在油菜秸稈還田條件下，僅需少量N肥即可滿(mǎn)足后茬水稻生長(zhǎng)需求。另一方面，油菜根系可分泌多種有機(jī)酸和酶類(lèi)物質(zhì)，尤其是酸性磷酸酶的分泌，活化釋放土壤中難溶態(tài)磷，從而提高根際土中速效磷的含量[28-29]，本試驗(yàn)中土壤苗期速效磷的含量低于盛花期，可能與盛花期根系發(fā)達(dá)，釋放了更多難溶態(tài)磷有關(guān)。整體而言，不同品種(系)氮磷鉀含量在苗期和盛花期的變化情況不一致，推測(cè)與不同品種(系)的生長(zhǎng)發(fā)育差異性有關(guān)。以土壤養(yǎng)分為篩選指標(biāo)，則品系11、10、4表現(xiàn)較佳。

不同生育期土壤微生物組成類(lèi)群不同，不同品系間同一類(lèi)群占比不同。本研究中，苗期與盛花期微生物豐度及微生物類(lèi)群存在一定差異，且同一微生物類(lèi)群在不同品系間含量亦有一定差異。同時(shí)，相關(guān)性分析表明厚壁菌和疣微菌門(mén)類(lèi)群在苗期與盛花期與土壤養(yǎng)分有顯著或極顯著相關(guān)性。疣微菌門(mén)是一門(mén)被劃分出不久的細(xì)菌類(lèi)群，只有少數(shù)被識(shí)別的種類(lèi)，Aguirre-von-Wobeser等[30]研究表明，與對(duì)照土壤相比，疣微菌門(mén)是唯一一個(gè)在玉米(ZeamaysL.)根部土壤中的豐度顯著增加的門(mén)類(lèi)，盡管其功能研究尚不明確，但在水稻(OryzasativaL.)[31]、韭菜(Alliumporrum)[32]、甘蔗(Saccharumofficinarum)[33]中已經(jīng)被證明該類(lèi)菌可利用植物性碳源促進(jìn)生長(zhǎng)。且疣狀菌可以被描述為寡養(yǎng)型，當(dāng)?shù)^(guò)量時(shí)，其豐度降低[34-35]。據(jù)此可調(diào)控土壤施肥量，避免因施肥量過(guò)多而造成的浪費(fèi)與負(fù)面環(huán)境影響。厚壁菌是最常見(jiàn)的植物根際細(xì)菌之一[36]。據(jù)報(bào)道，厚壁菌門(mén)可培養(yǎng)性很高，在可培養(yǎng)的細(xì)菌種類(lèi)中占14%[37]。Polyanskaya等[38]從豌豆(Pisumsativum)根際分離到一株厚壁菌門(mén)的梭狀芽胞桿菌(Clostridiumsp.)，對(duì)大麥(Hordeumvulgare)和黃瓜(Cucumissativus)的生長(zhǎng)均有促進(jìn)作用，且能抑制番茄(Lycopersicumesculentum)根系有害菌的生長(zhǎng)[39]，并提高芒草(Miscanthussinensis)對(duì)鹽脅迫的耐受性[40]。此外，厚壁菌參與細(xì)胞壁生物合成和跨膜轉(zhuǎn)運(yùn)等[41]，且具有磷酸鹽增溶和固氮等多種促進(jìn)植物生長(zhǎng)的特性[42]，以疣微菌和厚壁菌豐度作為篩選綠肥油菜的評(píng)價(jià)因子進(jìn)一步研究，對(duì)于發(fā)展可持續(xù)農(nóng)業(yè)具有重要意義。

4 結(jié)論

本研究發(fā)現(xiàn)，種植油菜綠肥會(huì)改變土壤養(yǎng)分含量與微生物群落豐度。與對(duì)照相比，12個(gè)油菜品系中，11號(hào)(辰溪臘油菜)、4號(hào)(紫葉芥)品系綜合表現(xiàn)突出，是良好的綠肥種質(zhì)。此外，除植株養(yǎng)分積累量與土壤養(yǎng)分外，土壤微生物中疣微菌和厚壁菌豐度同樣可用于篩選綠肥種質(zhì)。