基于多元統(tǒng)計(jì)分析的曹妃甸區(qū)24份鹽地堿蓬成分指標(biāo)綜合評(píng)價(jià)

孟 然，張曉東，楊雅華，王秀萍，薛志忠

（河北省農(nóng)林科學(xué)院濱海農(nóng)業(yè)研究所，唐山市植物耐鹽研究重點(diǎn)實(shí)驗(yàn)室，河北唐山 063200）

土壤鹽漬化是農(nóng)業(yè)發(fā)展面臨的重要問題之一，會(huì)導(dǎo)致土壤質(zhì)地粘重、透氣性不良，表層積鹽重、上層高礦化度地下水積累多、心底土鹽分含量高，全鹽含量常在0.4%～1.8%，pH常在7～8之間[1?2]。唐山市地處渤海之濱，南部沿海擁有近1100 km2的灘涂及鹽堿地，而曹妃甸區(qū)是唐山市境內(nèi)鹽堿地程度較高、鹽漬化區(qū)域較集中的濱海地區(qū)，全區(qū)境內(nèi)80%的土地為鹽堿土壤，在我國(guó)鹽堿地中具有代表性。

鹽地堿蓬（Suaeda salsa(L.)Pall.）俗名黃須菜，為藜科堿蓬屬一年生草本鹽生植物，對(duì)鹽堿地區(qū)的土壤、氣候等環(huán)境條件具有很強(qiáng)的適應(yīng)性，是鹽堿地的典型指示植物[3]，在曹妃甸區(qū)野生資源豐富。鹽地堿蓬幼芽莖葉鮮嫩可食用，富含維生素、蛋白質(zhì)、氨基酸等營(yíng)養(yǎng)物質(zhì)，其中維生素尤以抗壞血酸含量最為突出，最高可達(dá)到86.00 mg/100 g[4]，與蔬菜中含量較高的辣椒幾乎相當(dāng)；蛋白質(zhì)與膳食纖維含量共同占干物質(zhì)總量的30%以上[5]；氨基酸含量豐富，種類齊全，必需氨基酸構(gòu)成合理[6]，是一種極具營(yíng)養(yǎng)價(jià)值的野菜，被譽(yù)為“菜中靈芝”[7]。此外，還含黃酮類、萜類、多酚類、多糖類物質(zhì)等[8]，具有降血糖、降血脂、抗腫瘤、抗氧化、免疫調(diào)節(jié)的作用[9?11]。

亞硝酸鹽超標(biāo)已經(jīng)成為影響食品安全的重要原因之一，由于生長(zhǎng)環(huán)境的特殊性，鹽地堿蓬等野菜中亞硝酸鹽含量為一項(xiàng)重要監(jiān)測(cè)指標(biāo)?，F(xiàn)有文獻(xiàn)多記載山東鹽城、東營(yíng)及山西汾河灘地區(qū)堿蓬亞硝酸鹽情況，不同地區(qū)堿蓬亞硝酸鹽含量存在較大差異，部分亞硝酸鹽含量超過了聯(lián)合國(guó)糧食農(nóng)業(yè)組織/世界衛(wèi)生組織規(guī)定和中國(guó)食品安全國(guó)家標(biāo)準(zhǔn)[12]，而關(guān)于河北地區(qū)唐山一帶堿蓬鮮有報(bào)道。隨著對(duì)鹽地堿蓬作為優(yōu)質(zhì)的蔬菜、經(jīng)濟(jì)作物和具有醫(yī)療保健價(jià)值屬性的充分認(rèn)識(shí)，堿蓬資源的開發(fā)引起人們的重視，但缺少對(duì)于唐山曹妃甸地區(qū)鹽地堿蓬成分指標(biāo)評(píng)價(jià)較為系統(tǒng)的研究。

據(jù)此，本研究以曹妃甸內(nèi)不同區(qū)域的鹽地堿蓬為研究對(duì)象，對(duì)其水分、抗壞血酸、蛋白質(zhì)、氨基酸、總黃酮、亞硝酸鹽含量進(jìn)行檢測(cè)，利用相關(guān)性分析、主成分分析以及聚類分析方法，對(duì)鹽地堿蓬成分指標(biāo)進(jìn)行綜合評(píng)價(jià)，擬篩選出綜合評(píng)價(jià)不同地區(qū)鹽地堿蓬品質(zhì)的指標(biāo)，并探究了土壤鹽分、速效養(yǎng)分與品質(zhì)之間的關(guān)系，以期為曹妃甸區(qū)鹽地堿蓬質(zhì)量評(píng)價(jià)提供可靠的分析方法與參考依據(jù)，也為進(jìn)一步加工發(fā)展與綜合利用提供數(shù)據(jù)支持。

1 材料與方法

1.1 材料與儀器

樣本 采集于唐山市曹妃甸區(qū)不同采樣點(diǎn)，共24份材料于2020年6月統(tǒng)一采集，經(jīng)河北省農(nóng)林科學(xué)院王秀萍研究員鑒定為藜科堿蓬屬鹽地堿蓬（Suaeda salsa）；蘆丁標(biāo)準(zhǔn)品（貨號(hào)：R189033，純度95%） 阿拉丁公司；磷酸、甲醇 色譜純，F(xiàn)isher Chemical公司；無水乙醇、草酸、抗壞血酸、碳酸氫鈉、2,6-二氯酚靛酚等 分析純，天津市恒興化學(xué)試劑制造有限公司；氨基酸含量檢測(cè)試劑盒（貨號(hào)：BC1570）、亞硝酸鹽含量測(cè)定試劑盒（貨號(hào)：BC1490）、雙縮脲法蛋白質(zhì)含量檢測(cè)試劑盒（貨號(hào)：BC3180）北京Solaibio科技有限公司；純凈水 屈臣氏。

AP124W電子天平 島津公司；JT-1027HTD機(jī)械超聲波清洗機(jī) 深圳市潔拓超聲波清洗設(shè)備有限公司；D-37520 osterode臺(tái)式高速冷凍離心機(jī) 德國(guó)Biofagestratos公司；TU-1810紫外可見分光光度計(jì)北京普析通用儀器有限責(zé)任公司；HH-6數(shù)顯恒溫水浴鍋 上海力辰邦儀器科技有限公司；SCIENTZ-48組織研磨機(jī) 寧波新芝生物科技股份有限公司；1200型高效液相色譜儀 Agilent公司。

1.2 實(shí)驗(yàn)方法

1.2.1 樣品的采集方法 為避免工廠污染的影響，選取曹妃甸區(qū)內(nèi)鹽地堿蓬生長(zhǎng)密集區(qū)，從不同區(qū)域的野外鹽堿荒地隨機(jī)進(jìn)行鹽地堿蓬采樣，共獲得24份樣本，分別來自于唐山市曹妃甸區(qū)一場(chǎng)、二場(chǎng)、七場(chǎng)、十一場(chǎng)、柳贊鎮(zhèn)、南堡區(qū)、大學(xué)城、工業(yè)區(qū)，具體采樣地的地理坐標(biāo)及采收日期見表1。

表1 鹽地堿蓬采樣地點(diǎn)及時(shí)間Table 1 Sampling location and time of Suaeda salsa

選擇生長(zhǎng)狀態(tài)良好、無機(jī)械損傷綠色表型的鹽地堿蓬植株，將樣本裝入牛皮紙袋，編號(hào)，每個(gè)樣點(diǎn)采集樣本3份。將每一份鹽地堿蓬樣本分成兩部分，一部分留用鮮樣待測(cè)，一部分用流動(dòng)水洗去表面泥沙，再用蒸餾水沖洗3～4次，于60 ℃烘干12 h后粉碎過0.15 mm孔徑篩，粉末備用。

在所取植株近根處采用取土器垂直取20 cm深度的耕層土壤，每個(gè)樣點(diǎn)采集樣本3份，混勻后用四分法留取1 kg土樣，裝袋帶回，經(jīng)60 ℃烘干12 h后過1 mm和0.25 mm孔徑篩備用。

1.2.2 不同區(qū)域鹽地堿蓬品質(zhì)指標(biāo)測(cè)定

1.2.2.1 基本營(yíng)養(yǎng)成分的測(cè)定 水分、抗環(huán)血酸、蛋白質(zhì)、氨基酸含量采用鮮樣測(cè)定，水分含量的測(cè)定參照GB5009.3—2010食品中水分的測(cè)定-直接干燥法[13]?？弓h(huán)血酸含量的測(cè)定參照GB 5009.86-2016《食品安全國(guó)家標(biāo)準(zhǔn)-抗壞血酸的測(cè)定》[14]。

蛋白質(zhì)含量測(cè)定采用試劑盒法，蛋白質(zhì)（mg/g 鮮重）=C標(biāo)準(zhǔn)管/（A標(biāo)準(zhǔn)管?A空白管）×（A測(cè)定管?A空白管）×V樣總/W。其中，C標(biāo)準(zhǔn)管：5 mg/mL；A標(biāo)準(zhǔn)管：標(biāo)準(zhǔn)管的吸光值；A空白管：空白管的吸光值；V樣總：樣本總體積，1 mL；W：樣本質(zhì)量，g。

氨基酸含量測(cè)定采用試劑盒法，氨基酸含量（μmol/g 質(zhì)量）=[C標(biāo)準(zhǔn)管×V標(biāo)準(zhǔn)管×△A測(cè)定管/△A標(biāo)準(zhǔn)管]×（V樣總/V樣）/W。其中，C標(biāo)準(zhǔn)管：標(biāo)準(zhǔn)品濃度，10 μmol/mL；V標(biāo)準(zhǔn)管：反應(yīng)體系中加入標(biāo)準(zhǔn)品體積，0.05 mL；V樣總：樣本總體積，1 mL；V樣：反應(yīng)體系中加入樣本體積，0.05 mL；W：樣本質(zhì)量，g。

1.2.2.2 總黃酮含量的測(cè)定 采用超聲波輔助甲醇提取法提取鹽地堿蓬總黃酮[15]，取一定質(zhì)量堿蓬莖葉粉末，采用纖維素酶與果膠酶于60 ℃水浴復(fù)合酶解30 min，用70%甲醇溶液調(diào)整至料液比1:100，于60 ℃超聲120 min，轉(zhuǎn)速8000 r/min離心15 min后取上清液。利用一次性注射器和0.45 μm針頭式過濾器將最終提取液注入至進(jìn)樣瓶中備用，待進(jìn)樣檢測(cè)。參考郭麗霞等[16]的高效液相色譜法，以蘆丁標(biāo)準(zhǔn)品制作標(biāo)準(zhǔn)曲線，以峰面積定量，通過回歸方程對(duì)樣品中的總黃酮進(jìn)行定量分析?；貧w方程為y=15467x+1.9（R2=0.997），其中y是總峰面積，x是樣品中總黃酮濃度（以蘆丁計(jì)），mg/mL。

HPLC條件：色譜柱：Dikma C18柱（4.6 mm×250 nm，5 μm）；流動(dòng)相：A（體積分?jǐn)?shù)0.2%的磷酸水溶液）：B（甲醇）＝50:50；檢測(cè)波長(zhǎng)：350 nm；柱溫：30 ℃；流速：1 mL/min；進(jìn)樣量：10 μL。流動(dòng)相均經(jīng)真空泵過濾除雜和超聲清洗機(jī)脫氣。

1.2.2.3 亞硝酸鹽含量的測(cè)定 亞硝酸鹽含量使用鮮樣測(cè)定，測(cè)定方法采用試劑盒法，亞硝酸鹽含量（μmol/g 鮮重）=0.12×[(A樣品?A空白）/（A標(biāo)準(zhǔn)?A空白）]/W。其中，A樣品：測(cè)定管的吸光值；A空白：空白管的吸光值；A標(biāo)準(zhǔn)：標(biāo)準(zhǔn)管的吸光值；W：樣本質(zhì)量，g。

1.2.3 不同區(qū)域鹽地堿蓬土壤成分測(cè)定 參照《土壤農(nóng)化分析》[17]及《土壤學(xué)實(shí)驗(yàn)指導(dǎo)》[18]的方法，采用堿解擴(kuò)散法、碳酸氫鈉法和醋酸銨法分別測(cè)定土壤速效氮、速效磷和速效鉀含量。

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

遵循平行重復(fù)原則，所有試驗(yàn)均重復(fù)3次，數(shù)據(jù)以平均值±標(biāo)準(zhǔn)差表示。采用SPSS 23.0統(tǒng)計(jì)軟件進(jìn)行描述性分析，分析各指標(biāo)的均值、極值、標(biāo)準(zhǔn)差和變異系數(shù)，并進(jìn)行單因素ANOVA方差分析，運(yùn)用Duncan檢驗(yàn)法對(duì)成分指標(biāo)、土壤鹽分及速效養(yǎng)分含量分別進(jìn)行顯著性評(píng)價(jià)；采用雙變量Pearson相關(guān)分析法對(duì)成分指標(biāo)之間、與土壤鹽分及速效養(yǎng)分之間的相關(guān)性進(jìn)行分析；同樣采用SPSS對(duì)品質(zhì)指標(biāo)進(jìn)行主成分分析，以特征值大于1為標(biāo)準(zhǔn)，以每個(gè)主成分分別對(duì)應(yīng)的方差貢獻(xiàn)率占累積方差貢獻(xiàn)率的比例作為權(quán)重[19]，構(gòu)建堿蓬品質(zhì)的評(píng)價(jià)函數(shù)；采用組間平均數(shù)聯(lián)結(jié)法，根據(jù)平方歐式距離，對(duì)不同區(qū)域鹽地堿蓬進(jìn)行聚類分析。

2 結(jié)果與分析

2.1 不同區(qū)域鹽地堿蓬成分指標(biāo)分析

為探究鹽地堿蓬成分指標(biāo)的產(chǎn)地差異，測(cè)定了不同區(qū)域樣本中的水分、氨基酸、蛋白質(zhì)、抗壞血酸、亞硝酸鹽以及黃酮含量（表2），不同區(qū)域鹽地堿蓬的成分指標(biāo)之間變異系數(shù)不同，說明成分指標(biāo)之間呈現(xiàn)不同的差異性，亞硝酸鹽含量差異最大，水分含量差異最小。南堡區(qū)水分、蛋白質(zhì)及亞硝酸鹽含量均為最高值，其中南堡區(qū)、工業(yè)區(qū)鹽地堿蓬中蛋白質(zhì)含量分別為16.315和14.825 mg/g，顯著高于其他區(qū)域（P＜0.05），二場(chǎng)鹽地堿蓬蛋白質(zhì)含量最低，為4.999 mg/g，較其他組差異顯著（P＜0.05）；二場(chǎng)、工業(yè)區(qū)鹽地堿蓬抗壞血酸含量分別為4.565和4.034 mg/100 g，顯著高于其他區(qū)域（P＜0.05），且二者具顯著差異（P＜0.05），從報(bào)道的文獻(xiàn)來看，鹽地堿蓬抗壞血酸含量從0.55到86.00 mg/100g不等[4,20]，本研究中抗壞血酸含量與文獻(xiàn)中試驗(yàn)結(jié)果有所差異，其原因可能是采樣時(shí)間、采樣部位或采樣地點(diǎn)不同從而展現(xiàn)出不同的規(guī)律；南堡區(qū)鹽地堿蓬中亞硝酸鹽含量為5.885 μg/g，顯著高于其他區(qū)域（P＜0.05），根據(jù)國(guó)家食品安全標(biāo)準(zhǔn)（GB 2762—2017），蔬菜及其加工制品中亞硝酸鹽含量應(yīng)該在20 mg/kg[12]范圍以內(nèi)，所有采集區(qū)域樣本均未超標(biāo)。不同區(qū)域鹽地堿蓬中水分、氨基酸、黃酮含量無顯著差異（P＞0.05）。

表2 不同區(qū)域鹽地堿蓬成分指標(biāo)分析（n=3）Table 2 Analysis of component indexes of Suaeda salsa from different regions (n=3)

2.2 不同區(qū)域土壤鹽分及速效養(yǎng)分含量分析

如表3所示，土壤速效氮含量變異系數(shù)最大，為65.328%，說明不同區(qū)域土壤速效氮含量差異大。其中柳贊鎮(zhèn)土壤的速效氮含量達(dá)到54.250 mg/kg，顯著高于其他區(qū)域（P＜0.05）；南堡區(qū)土壤鹽分量最高，達(dá)2.360%，除柳贊鎮(zhèn)以外，較其他區(qū)域差異顯著（P＜0.05）；此外，南堡區(qū)速效鉀含量也具有最高值，為408.051 mg/kg，顯著高于其他區(qū)域（P＜0.05），七場(chǎng)土壤速效鉀含量顯著低于其他區(qū)域（P＜0.05）；二場(chǎng)土壤中速效磷含量最高，除一場(chǎng)、柳贊鎮(zhèn)以外，較其他區(qū)域差異顯著（P＜0.05）。

表3 不同區(qū)域土壤鹽分及速效養(yǎng)分含量分析(n=3)Table 3 Analysis of soil salt content and available nutrient content in different regions (n=3)

2.3 鹽地堿蓬成分指標(biāo)與土壤鹽分及速效養(yǎng)分相關(guān)性分析

相關(guān)性分析結(jié)果如表4所示，蛋白質(zhì)含量與土壤鹽分呈極顯著正相關(guān)（P＜0.01），與土壤速效鉀含量呈極顯著正相關(guān)（P＜0.01），表明土壤鹽分及速效鉀的變化顯著影響蛋白質(zhì)的變化，進(jìn)而影響鹽地堿蓬的營(yíng)養(yǎng)品質(zhì)。有研究表明，在一定范圍內(nèi)，隨土壤鹽分的增加，翅堿蓬葉片中可溶性蛋白的含量整體呈上升趨勢(shì)[21]，與本試驗(yàn)結(jié)果一致。有學(xué)者認(rèn)為，在一定濃度鹽脅迫處理下，堿蓬植物為了自我保護(hù)和抗鹽，細(xì)胞中蛋白質(zhì)合成代謝增強(qiáng)，參與滲透調(diào)節(jié)[22]。而鉀能提高植物光合作用效率，促進(jìn)碳水化合物的形成、轉(zhuǎn)運(yùn)和貯藏，促進(jìn)氮的吸收，加速蛋白質(zhì)的合成[23]。

表4 鹽地堿蓬成分指標(biāo)與土壤鹽分及速效養(yǎng)分的相關(guān)性分析Table 4 Correlation analysis of component indexes with content of soil salt and quick-acting nutrientsof Suaeda salsa

亞硝酸鹽與土壤鹽分呈顯著正相關(guān)（P＜0.05），與土壤速效鉀含量呈顯著正相關(guān)（P＜0.05），由此可知，一定程度的鹽脅迫促進(jìn)了鹽地堿蓬中蛋白質(zhì)和亞硝酸鹽含量的增加，與表2、表3結(jié)果相對(duì)應(yīng)，表明土壤鹽分及速效鉀的變化同樣顯著影響亞硝酸鹽含量，進(jìn)而影響鹽地堿蓬的安全品質(zhì)。

總黃酮含量與土壤鹽分呈極顯著負(fù)相關(guān)（P＜0.01），表明土壤鹽分的變化顯著影響總黃酮含量的變化。黃酮類物質(zhì)屬于次生代謝產(chǎn)物，在中等鹽度脅迫下，有助于植物體內(nèi)次生代謝產(chǎn)物積累[24]。但本研究結(jié)果顯示土壤鹽分阻礙了黃酮含量的積累，盡管不同植物中，報(bào)道結(jié)果并不一致，這可能與不同植物黃酮含量差異及對(duì)環(huán)境響應(yīng)不同造成的，尚需更多研究數(shù)據(jù)支持。

2.4 鹽地堿蓬成分指標(biāo)之間的相關(guān)性分析

鹽地堿蓬成分指標(biāo)間存在相互依存和制約的關(guān)系，多樣本兩指標(biāo)之間系數(shù)絕對(duì)值越大，則這兩指標(biāo)之間的聯(lián)系越緊密[25]。由表5可知，抗壞血酸與水分含量呈極顯著負(fù)相關(guān)（P＜0.01），與蛋白質(zhì)呈顯著負(fù)相關(guān)（P＜0.05），表明抗壞血酸的變化顯著影響水分及蛋白質(zhì)的變化，進(jìn)而影響鹽地堿蓬的營(yíng)養(yǎng)品質(zhì)；黃酮含量與亞硝酸鹽含量則呈顯著負(fù)相關(guān)（P＜0.05），表明亞硝酸鹽含量的變化影響了黃酮含量，進(jìn)而影響鹽地堿蓬的功效品質(zhì)；其他指標(biāo)間相關(guān)性不顯著（P＞0.05）。綜上所述，鹽地堿蓬6個(gè)品質(zhì)指標(biāo)相互間具有一定的相關(guān)性，但是所統(tǒng)計(jì)的原始數(shù)據(jù)反映的信息存在重疊。由此進(jìn)一步采取主成分分析法對(duì)6項(xiàng)品質(zhì)指標(biāo)簡(jiǎn)化，提高鹽地堿蓬品質(zhì)評(píng)價(jià)的分析效率，同時(shí)有利于提高統(tǒng)計(jì)信息的準(zhǔn)確性。

表5 鹽地堿蓬成分指標(biāo)之間的相關(guān)性分析Table 5 Correlation analysis of component indexes of Suaeda salsa

2.5 鹽地堿蓬成分指標(biāo)的主成分分析及綜合評(píng)價(jià)

對(duì)來自24份鹽地堿蓬樣本中的成分指標(biāo)含量進(jìn)行主成分分析，得出主成分的特征值、貢獻(xiàn)率、主成分載荷矩陣及特征向量矩陣等。由表6可知，共提取到兩個(gè)特征值大于1的主成分，貢獻(xiàn)率分別為54.427%和25.292%，累積貢獻(xiàn)率為79.719%，表明這兩個(gè)主成分可以基本反映本試驗(yàn)測(cè)定的所有成分指標(biāo)的信息。各主成分特征值、對(duì)原始性狀載荷及對(duì)表型的貢獻(xiàn)率如表6所示。

表6 總方差解釋Table 6 Explanation of total variance

載荷值反映了各品質(zhì)指標(biāo)對(duì)主成分矩陣中的權(quán)重，數(shù)值絕對(duì)值的大小反映了原始變量對(duì)于因子影響的強(qiáng)度，正負(fù)反映了方向[26]。圖1、表7顯示了主成分分析載荷值的分布，水分、抗壞血酸、亞硝酸鹽和黃酮在第一主成分上有較高載荷，說明第一主成分主要反映了這些元素的信息，其中水分、亞硝酸鹽在第一主成分上呈正向分布，載荷值分別為0.821和0.794，抗壞血酸和黃酮呈負(fù)向分布，載荷值分別為?0.781和?0.866?？梢?，負(fù)向作用大于正向作用。且抗壞血酸含量與水分、蛋白質(zhì)，黃酮與亞硝酸鹽均呈顯著相關(guān)性（表5），故選擇抗壞血酸和黃酮為第一主成分的代表指標(biāo)。氨基酸、蛋白質(zhì)在第二主成分上有較高載荷，其中氨基酸呈正向分布，載荷值為0.801，蛋白質(zhì)呈負(fù)向分布，載荷值?0.610?？梢姡蜃饔么笥谪?fù)向作用，故選擇氨基酸作為第二主成分的代表指標(biāo)。最終確定抗壞血酸、黃酮和氨基酸作為不同區(qū)域鹽地堿蓬品質(zhì)評(píng)價(jià)的核心指標(biāo)。

圖1 主成分分析荷載圖Fig.1 load diagram of PCA results

表7 主成分載荷矩陣Table 7 Principal component loading matrix

為了消除不同單位及數(shù)據(jù)量綱的影響，對(duì)品質(zhì)指標(biāo)原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，采用F1和F2兩個(gè)新變量來代替原來的6個(gè)變量，得出以下主成分表達(dá)式：

其中，F(xiàn)1和F2為2個(gè)主成分，X1～X6為6個(gè)品質(zhì)指標(biāo)含量經(jīng)標(biāo)準(zhǔn)化后的數(shù)據(jù)。

以每個(gè)主成分分別對(duì)應(yīng)的方差貢獻(xiàn)率占累積方差貢獻(xiàn)率的比例作為權(quán)重，計(jì)算綜合評(píng)價(jià)得分F綜。

根據(jù)綜合評(píng)價(jià)模型計(jì)算不同區(qū)域鹽地堿蓬的得分并進(jìn)行排序，如表8所示，總得分越高表明該地區(qū)鹽地堿蓬的綜合品質(zhì)越好，8個(gè)區(qū)域鹽地堿蓬中，排名前三的樣本編號(hào)依次為S6、S1、S3，說明南堡區(qū)、一場(chǎng)及七場(chǎng)的鹽地堿蓬品質(zhì)表現(xiàn)好，但南堡區(qū)鹽地堿蓬樣本檢測(cè)出亞硝酸鹽含量超過國(guó)家標(biāo)準(zhǔn)（小于4 mg/kg）[12]，S2即二場(chǎng)鹽地堿蓬主成分綜合得分最低，品質(zhì)表現(xiàn)最差。

表8 主成分綜合得分和排序Table 8 Comprehensive scores and ranking of principal components

2.6 鹽地堿蓬成分指標(biāo)的聚類分析

鹽地堿蓬各成分指標(biāo)數(shù)據(jù)較為離散，不易直觀地進(jìn)行相近程度的比較。聚類分析可同時(shí)對(duì)大量性狀進(jìn)行綜合考察[27]。利用系統(tǒng)聚類法對(duì)不同區(qū)域鹽地堿蓬進(jìn)行聚類分析，構(gòu)建樹系圖（圖2），將樣本劃為不同類群進(jìn)行評(píng)價(jià)。由圖可知，在平方歐式距離15處可將24份樣本分為三類：第一類為S1、S3、S5、S6，第二類為S4、S7、S8，第三類為S2。從總體數(shù)據(jù)來看，不同地區(qū)樣品間具有一定差異，其中S2與其他樣品差異較大，距離為 25 時(shí)就可以歸為不同的類型；在距離為1時(shí)所有樣本完全分離，說明雖然24份堿蓬樣本來自于不同區(qū)域，但仍有較強(qiáng)相似性。結(jié)合主成分分析綜合得分和排名可知，第一類的品質(zhì)為優(yōu)，第二類的品質(zhì)為一般，第三類的品質(zhì)為差，即一場(chǎng)、七場(chǎng)、柳贊鎮(zhèn)和南堡區(qū)的鹽地堿蓬品質(zhì)較好，二場(chǎng)的鹽地堿蓬品質(zhì)最差，聚類分析結(jié)果與主成分分析結(jié)果基本一致，表明聚類分析和主成分分析均可用來分析鹽地堿蓬的成分指標(biāo)，綜合評(píng)價(jià)不同區(qū)域鹽地堿蓬品質(zhì)優(yōu)劣。研究結(jié)果與劉丙花[27]和李躍紅[28]等采用兩種方法分別對(duì)不同品種藍(lán)莓品質(zhì)分析和不同產(chǎn)地紅心獼猴桃品質(zhì)的評(píng)判結(jié)果一致。

圖2 不同地區(qū)鹽地堿蓬的聚類分析圖Fig.2 Cluster analysis diagram of Suaeda salsa from different producing areas

3 結(jié)論

不同區(qū)域鹽地堿蓬品質(zhì)存在明顯差異，亞硝酸鹽含量差異最大，水分含量差異最小。采用主成分分析對(duì)8個(gè)區(qū)域鹽地堿蓬6個(gè)品質(zhì)性狀進(jìn)行主成分分析研究，提取前2個(gè)主成分，累積貢獻(xiàn)率達(dá)到79.72%，篩選出抗壞血酸、黃酮和氨基酸含量3個(gè)核心指標(biāo)。結(jié)合系統(tǒng)聚類結(jié)果可知，南堡區(qū)、一場(chǎng)及七場(chǎng)的鹽地堿蓬品質(zhì)較好，二場(chǎng)的鹽地堿蓬品質(zhì)最差，綜合運(yùn)用主成分分析和聚類分析方法可以簡(jiǎn)化鹽地堿蓬成分評(píng)價(jià)指標(biāo)，有助于優(yōu)良品質(zhì)特性的鹽地堿蓬的評(píng)價(jià)與篩選。結(jié)合土壤鹽分和速效養(yǎng)分結(jié)果分析，表明土壤鹽分、速效鉀含量高，能夠顯著提升鹽地堿蓬蛋白質(zhì)和亞硝酸鹽含量，但抑制了總黃酮含量的積累，建議鹽地堿蓬栽培發(fā)展中著重關(guān)注土壤鹽分和速效養(yǎng)分含量。本試驗(yàn)對(duì)于人工質(zhì)控栽培鹽地堿蓬提供依據(jù)，可為鹽地堿蓬的質(zhì)量評(píng)價(jià)提供可靠的分析方法與參考依據(jù)，為進(jìn)一步的綜合開發(fā)利用和產(chǎn)業(yè)發(fā)展提供數(shù)據(jù)支持。但鹽地堿蓬品質(zhì)差異還可能與多糖、脂質(zhì)、膳食纖維含量或者其他無機(jī)離子等因素有關(guān)，在后續(xù)的研究中需要加入更多的評(píng)價(jià)指標(biāo)，以期更好地為鹽地堿蓬品質(zhì)篩選提供理論指導(dǎo)。