楊建 湯華成,3,4 曹冬梅,3,4 崔航 婁雨豪 王冀菲 張東杰

（1 黑龍江八一農(nóng)墾大學(xué)食品學(xué)院，163319，黑龍江大慶；2 黑龍江省農(nóng)產(chǎn)品加工與質(zhì)量安全重點(diǎn)實(shí)驗(yàn)室，163319，黑龍江大慶；3 黑龍江省雜糧加工及質(zhì)量安全工程技術(shù)研究中心，163319，黑龍江大慶；4 國(guó)家雜糧工程技術(shù)研究中心，163319，黑龍江大慶；5 北大荒現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)省級(jí)培育協(xié)同創(chuàng)新中心，163319，黑龍江大慶）

紅小豆是我國(guó)重要的雜糧作物之一[1]，其苗期易受雜草影響，導(dǎo)致產(chǎn)量下降，因此，合理使用除草劑既可以控制雜草，又可以減少藥害，達(dá)到有效種植，確保紅小豆增收增產(chǎn)。氟磺胺草醚（fomesafen，F(xiàn)SA）又稱虎威，屬于二苯醚類除草劑，其主要作用于原卟啉原氧化酶[2]，因其除草效率高、選擇性好[3]，可被雜草葉片及根部吸收，抑制雜草光合速率及對(duì)有機(jī)物和能量汲取，從而達(dá)到除草效果[4]，目前已廣泛用于大豆[5]、番茄[6]、小麥[7]、草莓[8]和玉米[9]等田間闊葉類雜草的防除。研究[10]表明，用25%FSA·烯草酮乳油對(duì)綠豆和紅小豆田間闊葉雜草防除效果較好。紅小豆田間噴施25%FSA 能使紅小豆增產(chǎn)[11]。但田間噴施FSA 影響大豆幼苗的光合速率[12]，抑制幼苗蔗糖代謝和碳代謝過程，阻礙同化產(chǎn)物的合成和運(yùn)輸，影響大豆根瘤的能量供應(yīng)[13-14]。

代謝組學(xué)是系統(tǒng)生物學(xué)的一個(gè)分支，在研究環(huán)境因素作用下的特定代謝型具有較高的靈敏度，在環(huán)境變化對(duì)機(jī)體的影響研究等方面具有優(yōu)勢(shì)[15-16]。前人[17]研究表明，采用LC-MS 技術(shù)分析樂果對(duì)小鼠的影響，在尿液中篩選出12 種差異代謝物，在血漿中篩選出13 種差異代謝物，且代謝物二甲基硫代磷酸和二甲基二硫代磷酸在所有的差異代謝物中最為顯著。利用LC-MS 方法檢測(cè)吡蟲啉對(duì)蜜蜂羽化出房時(shí)代謝物的影響，共鑒定出18 種差異代謝物，主要調(diào)控鳥苷酸代謝和能量代謝等途徑，對(duì)蜜蜂的生長(zhǎng)發(fā)育產(chǎn)生顯著影響[18]。利用代謝組學(xué)技術(shù)研究了惡唑烷酮對(duì)玉米幼苗代謝產(chǎn)物的影響，檢測(cè)到90 多種代謝產(chǎn)物，其中24 種差異顯著，這些差異代謝物參與玉米幼苗的能量和蛋白質(zhì)代謝[19]。

目前，已有文獻(xiàn)[10-11]報(bào)道，F(xiàn)SA 在紅小豆田間具有良好的除草效果，但其在田間施用后對(duì)紅小豆幼苗生長(zhǎng)代謝影響的相關(guān)研究未見報(bào)道。因此，本研究以紅小豆幼苗為研究對(duì)象，利用LC-MS 技術(shù)探究FSA 脅迫紅小豆幼苗生長(zhǎng)過程中代謝物及代謝通路變化情況，為紅小豆田間使用FSA 的安全性評(píng)價(jià)提供理論參考。

1 材料與方法

1.1 試驗(yàn)材料

試驗(yàn)紅小豆品種為珍珠紅。試驗(yàn)試劑為25%氟磺胺草醚水劑；甲醇、乙腈、甲酸和異丙醇均為色譜純，2-氯-L-苯丙氨酸（純度≥98%）。

1.2 儀器與設(shè)備

試驗(yàn)所用設(shè)備如下，冷凍離心機(jī)（Centrifuge 5430 R）、多樣品冷凍研磨儀（Wonbio-96 c）、臺(tái)式快速離心濃縮干燥器（LNG-T 88）、氮?dú)獯祾邇x（JXDC-20）、UHPLC 液相色譜系統(tǒng)（Vanquish Horizon system）、質(zhì)譜儀（Q-ExactiveHF-X）。

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

1.3.1 FSA 脅迫紅小豆種植試驗(yàn)及樣品采集 挑選出籽粒大小均勻、飽滿的紅小豆種子，用去離子水洗去表面污漬后，浸泡于30%H2O2溶液中消毒5min，再用去離子水沖洗，直到泡沫完全消失，晾干后備用。根據(jù)除草劑FSA 的推薦使用量進(jìn)行紅小豆萌發(fā)試驗(yàn)。每盆種50 粒，每3 盆為1 個(gè)處理，每個(gè)處理3 個(gè)重復(fù)，以未噴藥為空白對(duì)照，待紅小豆苗長(zhǎng)到2～3 葉時(shí)用除草劑FSA（0.3mL/m2）進(jìn)行噴施，噴藥次數(shù)為2 次，間隔時(shí)間為14d，第2 次噴藥結(jié)束7d 后，對(duì)紅小豆幼苗進(jìn)行采樣（即紅小豆幼苗為5 葉時(shí)），采集未噴施FSA 的紅小豆幼苗為Z-2-ZZ 組，噴施FSA 的紅小豆幼苗為Z-2-ZZ-2組，所有采集樣品置于-20℃冰箱中儲(chǔ)存?zhèn)溆谩?/p>

1.3.2 樣品處理 準(zhǔn)確稱取50mg 樣品至2mL 離心管中，加入一顆直徑為6mm 的研磨珠和400μL提取液[甲醇:水=4:1（v:v）]，在冷凍組織研磨儀研磨6min（-10℃，50Hz），再用低溫超聲提取30min（5℃，40kHz），將樣品于-20℃條件下靜置30min，離心15min（13 000g，4℃），移取上清液至帶內(nèi)插管的進(jìn)樣小瓶中上機(jī)分析。最后每個(gè)樣本分別移取20μL 上清液混合后作為質(zhì)控樣本。

1.3.3 LC-MS 檢測(cè) 色譜條件參考林立銘等[20]的檢測(cè)條件。質(zhì)譜條件參照王琪琪[21]的檢測(cè)條件并稍作修改，毛細(xì)管溫度325℃，加熱溫度425℃，分辨率7500MS2。

1.4 數(shù)據(jù)處理

使用ProgenesisQI（美國(guó)）軟件進(jìn)行峰提取等預(yù)處理，并將提取的特征峰在公共的代謝組數(shù)據(jù)庫(kù)KEGG 和HMDB 進(jìn)行搜庫(kù)鑒定，匹配出的代謝集用ropls（R）軟件進(jìn)行多元統(tǒng)計(jì)分析，結(jié)合Student’s t 檢驗(yàn)的P＜0.05 和正交偏最小二乘―判別分析（OPLS-DA）中VIP＞1 篩選出差異代謝物，最后對(duì)代謝物進(jìn)行KEGG 通路富集分析。

2 結(jié)果與分析

2.1 FSA 脅迫下紅小豆幼苗主成分分析（PCA）

PCA 圖直觀地反映樣本中代謝物的豐富度情況，樣本間位置越近則越相似，越遠(yuǎn)則反之[22]。對(duì)Z-2-ZZ-2 組與Z-2-ZZ 組進(jìn)行差異代謝物主成分分析，并繪制得分圖。由圖1 可知，2 組樣本除少數(shù)點(diǎn)外，其余都在95%置信圈內(nèi)，Z-2-ZZ-2 組與Z-2-ZZ 組的主成分明顯不同，分別位于置信圈的左右兩側(cè)。從分布模式來(lái)看，正離子模式下，置信圈內(nèi)右側(cè)的Z-2-ZZ-2 組樣本存在部分重疊現(xiàn)象，可能是因?yàn)橄嗤幚順颖镜拇x物具有相似性造成的；負(fù)離子模式下，置信圈左側(cè)6 個(gè)平行樣品間的距離更近，說明6 個(gè)樣本代謝物鑒定結(jié)果十分相近。

圖1 紅小豆幼苗Z-2-ZZ-2 組和Z-2-ZZ 組主成分分析得分圖Fig.1 Scores of principal component analysis of adzuki bean seedlings in Z-2-ZZ-2 and Z-2-ZZ groups

2.2 FSA 脅迫紅小豆幼苗正交偏最小二乘判別分析（OPLS-DA）

OPLS-DA 可剔除PCA 中不相關(guān)的變量，計(jì)算出VIP 值，用于差異代謝物的篩選[23]。通過對(duì)Z-2-ZZ-2 組和Z-2-ZZ 組進(jìn)行OPLS-DA 置換檢驗(yàn)分析，由圖2 可知，正離子模式下，R2X累積=0.75，R2Y累積=1，Q2=0.997，負(fù)離子模式下，R2X累積=0.79，R2Y累積=1，Q2=0.998，都比較接近于1，說明模型的預(yù)測(cè)能力和穩(wěn)定性較好。經(jīng)過200 次置換檢驗(yàn)后，正負(fù)離子模式下的R2和Q2的回歸直線均與縱坐標(biāo)相交且Q2分別為-0.201 和-0.1971，說明模型具有有效性和可靠性，不存在過擬合現(xiàn)象。

圖2 紅小豆幼苗Z-2-ZZ-2 組和Z-2-ZZ 組OPLS-DA 置換檢驗(yàn)結(jié)果Fig.2 OPLS-DA replacement test results of adzuki bean seedlings in Z-2-ZZ-2 and Z-2-ZZ groups

2.3 FSA 脅迫下紅小豆幼苗差異代謝物數(shù)量和類別鑒定結(jié)果

由圖3 可知，紅色點(diǎn)表示差異代謝物的表達(dá)量呈上調(diào)趨勢(shì)，藍(lán)色點(diǎn)表示差異代謝物的表達(dá)量呈下調(diào)趨勢(shì)，越往左右兩邊，差異越顯著。2 種離子模式下分析出的代謝物數(shù)量極多，但多數(shù)僅有變化趨勢(shì)，差異不顯著。

圖3 FSA 脅迫下紅小豆幼苗差異代謝物火山圖Fig.3 Volcanic figure of different metabolites of adzuki bean seedlings under FSA stress

Z-2-ZZ-2 與Z-2-ZZ 相比，從數(shù)量上看，在P＜0.05、VIP＞1 的條件下共篩選出236 種差異代謝物，如表1 和表2 可知，正離子模式下，分析出有名稱的代謝物106 種，其中上調(diào)50 種，下調(diào)56種；負(fù)離子模式下，分析出代謝物130 種，其中上調(diào)42 種，下調(diào)88 種。從分類上看，正離子模式下，上調(diào)的差異代謝物主要是苯丙烷和聚酮化合物、有機(jī)化合物、脂質(zhì)和類脂分子等，下調(diào)的差異代謝物主要是有機(jī)酸及其衍生物、脂質(zhì)和類脂分子、有機(jī)雜環(huán)化合物、核苷、核苷酸及其類似物、苯丙烷和聚酮化合物等；負(fù)離子模式下，上調(diào)的差異代謝物主要是苯丙烷和聚酮化合物、脂質(zhì)和類脂分子、有機(jī)氧化合物以及有機(jī)雜環(huán)化合物等，下調(diào)的差異代謝物主要是有機(jī)酸及其衍生物、脂質(zhì)和類脂分子、有機(jī)氧化合物、核苷、核苷酸及其類似物、苯丙烷和聚酮化合物等。

表1 FSA 脅迫下紅小豆幼苗差異代謝物定性結(jié)果（正離子)Table 1 Qualitative results of different metabolites in adzuki bean seedlings under FSA stress(cation)種species

表2 FSA 脅迫下紅小豆幼苗差異代謝物定性結(jié)果（負(fù)離子)Table 2 Qualitative results of different metabolites in adzuki bean seedlings under FSA stress(anion)種species

2.4 FSA 脅迫下紅小豆幼苗差異代謝物代謝通路分析

植物的生長(zhǎng)伴隨著復(fù)雜的代謝過程，由多種小分子代謝物共同調(diào)控，并不能從單一物質(zhì)含量變化對(duì)整體變化進(jìn)行判斷，因此需進(jìn)一步對(duì)其代謝途徑進(jìn)行具體的分析[24]。將篩選出的差異代謝物與KEGG 通路數(shù)據(jù)庫(kù)信息進(jìn)行匹配，對(duì)匹配到的代謝通路進(jìn)行拓?fù)浞治霾⒗L制氣泡圖。由圖4 可知，正離子模式下，共分析出16 條代謝通路；負(fù)離子模式下，共分析出12 條代謝通路，再結(jié)合Impact 值和P＜0.05 綜合分析，篩選出差異顯著的代謝通路以及參與相應(yīng)通路的差異代謝物。如表3 和表4 可知，正離子模式下，富集到5 條差異顯著的代謝通路，分別是嘧啶代謝、異黃酮生物合成、嘌呤代謝、半乳糖代謝和精氨酸和脯氨酸代謝，映射出13 種差異代謝物，其中差異代謝物尿嘧啶和胞嘧啶參與嘧啶代謝的代謝途徑，2,7,4′-三羥基異黃酮、甘氨酸和毛蕊異黃酮參與了異黃酮生物合成，鳥嘌呤、次黃嘌呤和β-谷甾酮參與嘌呤代謝，甘露三糖和蔗糖參與了半乳糖代謝，4-（谷氨酰胺）丁酸酯、對(duì)香豆酰腐胺和（2 S）-2-（3-羧基丙酰氨基）-5-氧代戊酸參與了精氨酸和脯氨酸代謝；負(fù)離子模式下，富集到2 條差異顯著的代謝通路，映射出5 種差異代謝物，分別是花青素生物合成和黃酮類生物合成，其中差異代謝物矢車菊素-3-葡萄糖苷和花葵素參與了花青素生物合成，根皮苷、（2 R,3 R）-3,4′,7-三羥基黃烷酮、花葵素和香樹脂參與了黃酮類生物合成。

表3 FSA 脅迫下紅小豆幼苗代謝通路富集結(jié)果（正離子）Table 3 Enrichment of metabolic pathway in adzuki bean seedlings under FSA stress(cation)

圖4 FSA 脅迫下紅小豆幼苗代謝途徑拓?fù)浞治鰵馀輬DFig.4 Bubble diagram of topological analysis of metabolic pathway of adzuki bean seedlings under FSA stress

由圖4 可知，顏色越紅，P值越小，說明通路越顯著，氣泡越大，通路越重要。正離子模式下，嘧啶代謝途徑、異黃酮生物合成途徑和嘌呤代謝途徑的氣泡相對(duì)較大，說明這3 條通路在分析FSA對(duì)紅小豆幼苗代謝途徑的影響較為重要，其次，半乳糖代謝和精氨酸和脯氨酸代謝的氣泡相對(duì)較小，說明FSA 對(duì)紅小豆幼苗代謝途徑的影響也有一定的重要性；負(fù)離子模式下，花青素生物合成途徑的氣泡較大，說明FSA 對(duì)紅小豆幼苗的花青素生物合成影響較大，通路較為重要，黃酮類生物合成途徑的氣泡稍小，說明FSA 對(duì)紅小豆幼苗中黃酮類物質(zhì)的合成有顯著影響。

3 討論

植物中產(chǎn)生的次級(jí)代謝產(chǎn)物主要用于應(yīng)對(duì)不良環(huán)境條件，如黃酮類物質(zhì)、萜類化合物和生物堿等，它們對(duì)非生物的脅迫具有耐受性[25]。紅小豆田間噴施FSA，主要對(duì)紅小豆幼苗中苯丙烷和聚酮類化合物、脂質(zhì)和類脂分子類化合物以及有機(jī)酸及其衍生物類化合物產(chǎn)生顯著影響。苯丙烷和聚酮類化合物中的差異代謝物主要是黃酮類化合物，研究[26]表明，黃酮類化合物是豆類中的一種酚類化合物，具有調(diào)理細(xì)胞分化、凋亡、抗自由基等作用[27-28]，紅小豆田間噴施FSA 能刺激紅小豆幼苗產(chǎn)生次級(jí)代謝產(chǎn)物中的黃酮類物質(zhì)，促進(jìn)紅小豆幼苗的生長(zhǎng)代謝。植物根系分泌有機(jī)酸是為了適應(yīng)環(huán)境脅迫，紅小豆幼苗生長(zhǎng)過程中受除草劑FSA 的脅迫，產(chǎn)生的有機(jī)酸釋放到根際土壤中，從而改變土壤環(huán)境，例如土壤的理化性質(zhì)和微生物活性等[29-31]，為紅小豆幼苗的生長(zhǎng)提供養(yǎng)分和能量，并且有機(jī)酸及其衍生物化合物中的羧酸及其衍生物屬于萜類化合物，具有抵御FSA 脅迫的作用。

糖類物質(zhì)能為植物的細(xì)胞代謝提供碳源和能量[32]。本研究發(fā)現(xiàn)，參與半乳糖代謝途徑的差異代謝物——蔗糖的表達(dá)量顯著下調(diào)，與程茁等[14]的研究結(jié)果一致，糖代謝為氨基酸代謝提供能量[33]，因此，氨基酸代謝中的精氨酸和脯氨酸代謝也受一定的影響。嘌呤和嘧啶主要以核苷酸的形式存在，而參與精氨酸和脯氨酸代謝的差異代謝物4-（谷氨酰胺）丁酸酯為嘌呤和嘧啶核苷酸類化合物的合成提供一部分氮源[34]。嘌呤和嘧啶屬于生物堿，生物堿的合成和積累受生物和非生物的調(diào)控[35]，本研究中嘌呤代謝和嘧啶代謝途徑富集到的差異代謝物尿嘧啶、胞嘧啶、鳥嘌呤、次黃嘌呤等的表達(dá)量在FSA 脅迫下發(fā)生顯著下調(diào)，噴施FSA 抑制了紅小豆幼苗的嘌呤代謝和嘧啶代謝。此外，植物類黃酮能使植物抵御外界環(huán)境的壓迫[25,36]，噴施FSA 能促進(jìn)紅小豆幼苗中黃酮類物質(zhì)的合成?；ㄇ嗨貙儆邳S酮類化合物，是一種天然色素，幼苗的花色由各種花青素的比例決定，矢車菊素-3-葡萄糖苷和花葵素分別產(chǎn)生紅色和橙色的花[37]，但花青素生物合成容易受外界因素的影響[38]，紅小豆田間噴施FSA 顯著影響紅小豆幼苗中花青素的生物合成。

4 結(jié)論

采用LC-MS 代謝組學(xué)技術(shù)分析紅小豆幼苗受FSA 脅迫后代謝物的變化具有可行性。結(jié)果表明，F(xiàn)SA 對(duì)紅小豆幼苗代謝產(chǎn)物的種類和數(shù)量都有顯著影響，主成分分析顯示，Z-2-ZZ-2 組和Z-2-ZZ組之間存在代謝組學(xué)差異，經(jīng)多元統(tǒng)計(jì)分析后，共篩選出236 種差異代謝物，且多數(shù)呈現(xiàn)下調(diào)趨勢(shì)，說明紅小豆田間噴施FSA 對(duì)紅小豆幼苗的生長(zhǎng)代謝造成顯著影響。

KEGG 通路分析發(fā)現(xiàn)，除草劑FSA 對(duì)紅小豆幼苗的生長(zhǎng)代謝途徑有顯著影響，正離子模式下，差異顯著且最為重要的通路是嘧啶代謝和異黃酮生物合成，負(fù)離子模式下，差異顯著且最為重要的通路是花青素生物合成。本研究中篩選出7 條差異顯著的代謝通路中有3 條（異黃酮生物合成、花青素生物合成和黃酮類生物合成）與紅小豆幼苗中黃酮類化合物的合成有關(guān)，噴施FSA 對(duì)紅小豆幼苗中黃酮類物質(zhì)的影響較大，黃酮類化合物是通路分析的關(guān)鍵代謝物。