馮帆，楊秀清，閆海冰

(山西農(nóng)業(yè)大學 林學院，山西 太谷 030801)

西伯利亞白刺 (NitrariasibiricaPall)，又名小果白刺，為蒺藜科(Zygophylaceae)白刺屬(Nitraria)落葉小灌木[1]，該樹種根系發(fā)達，多分枝、生長快、果實兼具食用和藥用價值，枝、葉可作飼料，在我國西北部與北部地區(qū)的新疆、青海、寧夏、甘肅以及內(nèi)蒙古等地區(qū)的鹽堿地上多有分布，因此被認為是鹽堿地生物改良的優(yōu)良經(jīng)濟型樹種資源[2,3]。中國目前17個省區(qū)都有鹽漬土資源分布，總面積逾5億畝。全球鹽堿地面積也在逐年擴大，鹽堿地己成為全球重要的后備土地資源[4]。利用耐鹽堿植物改良土地鹽堿化是提高鹽漬化土地利用效率，并對土壤起持久改善作用的一項合理途徑。對西伯利亞白刺進行耐鹽機理研究，可為白刺有效應用于鹽堿地改良和維持生態(tài)平衡提供重要依據(jù)，同時對白刺經(jīng)濟價值的發(fā)揮和擴大其資源利用具有重要意義。

關(guān)于白刺耐鹽機理的研究，楊升等[5]通過測定鹽脅迫后西伯利亞白刺幼苗體內(nèi)脯氨酸、可溶性糖、甜菜堿及礦質(zhì)元素(尤其是Na+、K+和 Cl-)含量的變化，從生理角度揭示了這幾種滲透調(diào)節(jié)物質(zhì)對白刺耐鹽能力的影響；張麗等[6]通過人工施鹽的脅迫處理，從生理生化的角度表明白刺在鹽脅迫下膜保護和修復能力均顯著提高，體內(nèi)有機滲透物質(zhì)增加，降低細胞內(nèi)水勢，同時還可以通過葉片胞內(nèi)的離子區(qū)隔化減少鹽分對重要細胞器的損傷；程鐵龍等[7]基于蛋白組學研究了鹽脅迫后的西伯利亞白刺葉片,共檢測到71種蛋白質(zhì)的表達發(fā)生了顯著變化，這些蛋白質(zhì)參與的過程包括光合作用、碳水化合物和能量代謝，膜運輸、信號轉(zhuǎn)導及氧化還原反應等；Chen J等[8]鑒定出鹽脅迫誘導產(chǎn)生的蛋白質(zhì)，得出鹽脅迫下光合作用的降低與參與光反應和卡爾文循環(huán)的酶、蛋白下調(diào)有關(guān)的結(jié)論。這些研究在不同水平上很好地揭示了西伯利亞白刺的耐鹽機制，尤其是鹽脅迫下體內(nèi)蛋白質(zhì)動態(tài)表達規(guī)律的分析將對西伯利亞白刺耐鹽機理的認識深入到分子生物學的水平。但是，僅基于蛋白組學的研究并不能反映鹽脅迫下白刺其他耐鹽相關(guān)代謝產(chǎn)物的動態(tài)變化，對西伯利亞白刺耐鹽相關(guān)代謝途徑的認識尚屬空白[9]?；诖x組學的研究將有助于更深入的了解適應脅迫環(huán)境的分子機制。

鑒于此，本實驗基于代謝組學的研究分析了西伯利亞白刺鹽脅迫前后差異代謝物種類及含量的變化，找到參與脅迫應答的代謝途徑，深入揭示其適應鹽脅迫的代謝生物學機制，為西伯利亞白刺在鹽堿地生物改良中的科學的資源利用提供基礎依據(jù)。

1 材料與方法

1.1 材料來源及實生苗培育

白刺種子采自內(nèi)蒙古巴彥淖爾市磴口縣境內(nèi)烏蘭布和沙漠東北部(海拔1 036 m,N40°32′,E106°28′)白刺產(chǎn)地的野生種。果實采回后揉搓濾去果汁、果皮，洗凈晾干，置于涼爽、通風、干燥條件下貯藏備用。

挑選飽滿種粒于0.5‰KMnO4溶液中避光浸泡2 h，后于200 mg.L-1ABT1號溶液中預處理30 min，洗凈后播于16 cm×14 cm規(guī)格裝有純凈濕沙的營養(yǎng)缽中，每缽播種70粒，種植深度4 cm，共計6缽，置于光照培養(yǎng)箱中培養(yǎng)。培養(yǎng)箱設置光照時間14 h，光照強度1 200 Lux，溫度27 ℃。試驗過程中澆灌蒸餾水，澆水量為100 mL，每隔3天澆灌一次。

1.2 鹽脅迫處理及取樣

播種12天后種子發(fā)芽率為68.1%，待幼苗生長50天進行脅迫前取樣，此時6個營養(yǎng)缽中正常生長的西伯利亞白刺幼苗共255株。每一營養(yǎng)缽中隨機剪取白刺幼苗葉片100 mg，6缽共計6次重復，取下樣品用錫箔紙包裹、編號，浸泡于液氮7 min后迅速取出置于-80 ℃冰箱備測。取樣后立即進行脅迫處理，用300 mM NaCl溶液澆灌幼苗，每盆澆灌100 mL，每3天澆灌一次。為防止鹽分流失，營養(yǎng)缽底部墊有塑料托盤，澆灌后滲出的溶液重新倒回營養(yǎng)缽內(nèi)，同時定時記錄白刺的生長狀況。脅迫后第6天幼苗頂端的嫩葉邊緣出現(xiàn)微黃甚至萎蔫，此時取脅迫后的樣品，取樣方法與脅迫前一致，對所取樣品編號備測。

1.3 代謝物的提取及檢測

取樣本25±1 mg于2 mL EP管中，加入300 μL提取液(甲醇∶水=3∶1)，再加入5 μL核糖醇，渦旋30 s后低溫離心15 min(12000 rpm轉(zhuǎn)速)，移取200 μL上清液于1.5 mL EP管中，真空濃縮干燥提取物，向干燥后的代謝物中加入30 μL甲氧胺鹽試劑(甲氧胺鹽酸鹽，溶于吡啶20 mg·mL-1)，混勻后放入烘箱中80 ℃孵育30 min，再向每個樣品中加入40 μL BSTFA(含有1% TMCS,v/v)，70 ℃孵育1.5 h，冷卻至室溫，加入5 μL FAMEs。使用配有Agilent DB-5MS毛細管柱的Agilent 7890氣相色譜-飛行時間質(zhì)譜聯(lián)用儀(30 m×250 μm×0.25 μm,J&W Scientific,Folsom,CA,USA)對樣品進行GC-TOFMS檢測。

1.4 數(shù)據(jù)處理及統(tǒng)計分析

使用ChromaTOF軟件(V 4.3x，LECO)對質(zhì)譜數(shù)據(jù)進行峰提取、基線矯正、解卷積、峰積分、峰對齊等分析[10]。在物質(zhì)定性分析時，使用LECO-Fiehn Rtx5數(shù)據(jù)庫，包括質(zhì)譜匹配及保留時間指數(shù)匹配[11]，采用多元變量統(tǒng)計分析方法對代謝物進行篩選，并找出差異代謝物，利用富集分析和拓撲分析對差異代謝物所在通路進行綜合分析，進一步確定西伯利亞白刺耐鹽性相關(guān)差異代謝物及耐鹽代謝通路。

2 結(jié)果與分析

2.1 鹽脅迫下白刺葉片代謝物改變的應答及層次聚類分析

鹽脅迫后西伯利亞白刺葉片中鑒定出差異代謝物共108種。除己內(nèi)酰胺和黃烷腺嘌呤衍生物2種代謝物濃度顯著下調(diào)外(下調(diào)濃度分別為3.18倍和4.22倍)，其余106種差異代謝物包括51種氨基酸類、22種糖類、11種脂肪酸類、8種有機酸類、7種醇類、5種生物堿以及2種維生素濃度均顯著上調(diào)(上調(diào)范圍在0.54倍和13.80倍之間)。其中上調(diào)濃度最大的代謝物為苯丙氨酸，上調(diào)濃度最小的代謝物為6-羥基己酸二聚物。顯著上調(diào)的前10個差異代謝物見表(表1)。將已定性的108種差異代謝物進行層次聚類分析(圖1)可知：層次聚類分析將這些代謝物主要分成了A和B兩大類，其中A類包括生物堿和堿基衍生物，B類主要包括大部分氨基酸類、糖類、脂肪酸類、醇類、有機酸類、生物堿等。

表1 鹽脅迫下西伯利亞白刺葉片中差異代謝物Table 1 Differentially expressed metabolites of Nitraria sibirica Pall under salt stress

注：Peak：定性得到的物質(zhì)名稱；Similarity：定性分析中該物質(zhì)與標準庫中物質(zhì)的匹配程度，取值[0,1000],數(shù)值越高說明定性出的物質(zhì)越準確；Rt：該物質(zhì)的色譜保留時間；LOG_FOLDCHANGE：該物質(zhì)對比兩組實驗間倍數(shù)關(guān)系以2為底的對數(shù)值，正號表示脅迫后較脅迫前的差異代謝物含量是顯著上調(diào)關(guān)系，數(shù)值越大濃度上調(diào)越多。
Note：Peak：A qualitative description of the substance；Similarity：The matching degree between the substance and the standard substance in the qualitative analysis，The range is 0 to 1000，The higher the value,the more accurate the substance；Rt：The chromatographic retention time of the substance；LOG_FOLDCHANGE：The substance was compared between the two groups of experiments with the logarithm of base 2，The positive sign indicates that the difference of metabolite content after stress is significantly higher than that before stress，The higher the value,the higher the concentration.

圖1 鹽脅迫下西伯利亞白刺差異代謝物聚類分析圖Fig.1 Hierarchical clustering analysis of differentially expressed metabolites of Nitraria sibirica Pall under salt stress

2.2 鹽脅迫后西伯利亞白刺差異代謝物的代謝通路分析

鹽脅迫后西伯利亞白刺葉片中濃度上調(diào)的106種差異代謝物經(jīng)過富集分析和拓撲分析被注釋到46條代謝通路中(圖2)，其中參與到氨基酸生物合成代謝途徑(圖3)和苯丙氨酸代謝途徑(圖4)中的差異代謝物種類最多。被注釋到氨基酸生物合成代謝途徑中的與鹽脅迫相關(guān)的差異代謝物有6種，分別為檸康酸、L-亮氨酸、3-異丙基蘋果酸、L-纈氨酸、L-異亮氨酸、丙酮酸；被注釋到苯丙氨酸代謝途徑中的與鹽脅迫相關(guān)的差異代謝物也有6種，分別為苯丙氨酸、L-半胱氨酸、纈氨酸、賴氨酸、異亮氨酸、天冬酰胺。乳糖、葡糖糖、木糖、阿洛糖主要參與糖降解、糖異生及半乳糖代謝過程。參與脂肪酸合成過程的代謝物主要包括棕櫚酸、硬脂酸和肉豆蔻酸等飽和脂肪酸以及亞麻酸、亞油酸等不飽和脂肪酸。檸檬酸、α-酮戊二酸、蘋果酸主要參與TCA循環(huán)過程。此外，其余差異代謝物還參與苯丙類生物合成途徑、嘌呤代謝、二羧酸鹽代謝、戊糖和葡聚糖酸酯間轉(zhuǎn)化代謝及氧化磷酸化過程等。同時，谷甾醇作為植物固醇的一種重要物質(zhì)，其含量也顯著上調(diào)，有利于維持膜的完整性。

3 討論與結(jié)論

本實驗鑒定到氨基酸類差異代謝物共51種。氨基酸是合成結(jié)構(gòu)蛋白與應激蛋白的重要小分子物質(zhì)，在植物的生命活動中承擔著極為重要的生理作用[13]。在正常生長狀態(tài)下，氨基酸、有機酸的含量往往很低，只有在逆境條件下這些有機物的含量才會上升，但他們本身并不會對植物的生長造成傷害[14]。Hare等[15]研究了鹽脅迫條件下脯氨酸含量的變化，結(jié)果表明脯氨酸的積累有助于植物抵抗逆境環(huán)境。本實驗結(jié)果顯示，鹽脅迫后西伯利亞白刺葉片中51種氨基酸差異代謝物的含量較脅迫前均顯著上調(diào)，其中，脯氨酸、苯丙氨酸、纈氨酸、丙氨酸、絲氨酸、谷氨酸、異亮氨酸等含量的上調(diào)范圍均在2.07倍以上。脯氨酸的合成與降解會促進細胞質(zhì)基質(zhì)、線粒體的氧化還原反應及腺苷酸的合成，同時還可以作為碳、氮及其他細胞抵抗逆境所需能源物質(zhì)的還原劑，這會提供給西伯利亞白刺葉片一個較高的逆境抵抗力，表明氨基酸在西伯利亞白刺體內(nèi)抵抗鹽脅迫、維持穩(wěn)態(tài)方面發(fā)揮了重要的作用。

圖2 鹽脅迫下西伯利亞白刺代謝通路分析氣泡圖Fig.2 Metabolic pathway analysis of Nitraria sibirica Pall under salt stress注：氣泡圖中每一個氣泡代表一個代謝通路，氣泡所在橫坐標和氣泡大小代表該通路在拓撲分析中的影響因子，氣泡越大影響因子越大；氣泡所在縱坐標和顏色表示富集分析的P值，顏色越深P值越小，富集程度越顯著。Note：Each bubble represents a metabolic pathway in the bubble map,The abscissa and bubble size represent the influencing factors of the path in topological analysis,The bigger the bubble,the bigger the impact factor;The ordinate and color of the bubble represent the P values for enrichment analysis,The darker the color,the smaller the P value,The higher the enrichment degree is.

圖3 氨基酸生物合成途徑Fig.3 Biosynthesis of amino acids pathway

圖4 苯丙氨酸代謝途徑Fig.4 Phenylalanine metabolism pathway

植物體內(nèi)糖類的形成、轉(zhuǎn)運、貯藏和消耗與植物生理代謝對外界的響應有緊密地聯(lián)系，逆境條件會影響植物體內(nèi)糖代謝途徑，因此，可以用糖代謝水平來衡量環(huán)境對植物的脅迫程度及植物對逆境環(huán)境適應性的強弱[16]。黃相玲等[17]研究了小葉欖仁幼苗對鹽堿環(huán)境的適應機理，結(jié)果表明，在0.6%～1.0%鹽脅迫下，可溶性糖含量隨鹽濃度的升高呈遞增趨勢。本實驗中，脅迫后西伯利亞白刺葉片中包括乳糖、阿洛糖、赤蘚糖、海藻糖、脫水半乳糖、景天庚酮糖等在內(nèi)的22種糖類的含量較脅迫前均顯著上調(diào)，這些糖類代謝物主要參與糖降解、糖異生及半乳糖代謝途徑，在參與西伯利亞白刺的鹽脅迫應答過程中，除了可以作為滲透調(diào)節(jié)物質(zhì)外，所參與的這些代謝過程中還會釋放大量ATP可為白刺抵抗鹽脅迫環(huán)境提供重要能量。

脂肪酸是植物體內(nèi)一類重要的生物活性物質(zhì)，其不僅是細胞膜脂的重要組成成分，還是生物體貯藏能量的主要來源[18]。張婧[19]通過比較中性鹽脅迫下兩種大豆的生長參數(shù)，表明大豆的耐鹽能力主要是依靠脂肪酸等小分子滲透調(diào)節(jié)物質(zhì)的積累。在脅迫后西伯利亞白刺葉片中檢測到11種脂肪酸，其中棕櫚酸和亞麻酸的濃度上調(diào)最為顯著，參與的代謝途徑主要是脂肪酸合成過程，通過合成脂肪酸來維持膜結(jié)構(gòu)的完整性，從而提高西伯利亞白刺的耐鹽性。

TCA循環(huán)不僅可為植物生命活動提供能量，而且還將糖類、脂類及蛋白質(zhì)代謝緊密聯(lián)系在一起，同時，TCA循環(huán)產(chǎn)生的許多中間產(chǎn)物也是植物體內(nèi)許多重要物質(zhì)的合成原料[20]。代謝通路分析可知，鹽脅迫條件下西伯利亞白刺葉片中檸檬酸、異檸檬酸、α-酮戊二酸、琥珀酸及蘋果酸等代謝物含量均顯著上調(diào)，表明鹽脅迫條件下白刺葉片中TCA循環(huán)顯著增強，這為其抵抗鹽脅迫提供了更多的能量，同時又將代謝組的多個耐鹽相關(guān)代謝物聯(lián)系起來，成為西伯利亞白刺耐鹽機制的重要組成部分。

鹽脅迫條件下西伯利亞白刺葉片中氨基酸代謝、糖代謝、脂肪酸代謝均參與滲透調(diào)節(jié)，TCA循環(huán)過程也顯著增強，說明西伯利亞白刺不僅依靠糖代謝及TCA循環(huán)過程產(chǎn)生能量，還通過積累滲透調(diào)節(jié)物質(zhì)以適應不良環(huán)境，西伯利亞白刺耐鹽機制較為完善。