汪 穎,王 尖,馮子珊,汪寶根,吳新義,魯忠富,孫玉燕,董文其,李國景,吳曉花
(浙江省農(nóng)業(yè)科學院 蔬菜研究所,浙江 杭州 310021)
瓠瓜[Lagenariasiceraria(Molina) Standl.](2n=2x=22)別名瓠子、長瓜、夜開花、葫蘆等,隸屬于葫蘆科葫蘆屬,為一年生草本植物,具有重要的經(jīng)濟價值。瓠瓜幼嫩果實可作蔬菜,成熟果實用作器具、飾品等,為西瓜、甜瓜等作物嫁接用砧木[1-2],是我國南方地區(qū)重要特色蔬菜作物之一,近年來在山東等北方省份栽培面積迅速增長,目前全國栽植面積約13萬hm2。瓠瓜栽培效益較高,產(chǎn)值可達22.5萬~45.0萬元·hm-2。
目前,便于運輸?shù)拈L棒形瓠瓜市場占有率高,但追求風味品質(zhì)的消費者則更傾向于購買市場占有率較低的近圓形、牛腿形、葫蘆形或短棒形瓠瓜。隨著人民生活水平的提高,果實品質(zhì)已成為瓠瓜市場競爭力的核心,高品質(zhì)育種已成為瓠瓜育種未來發(fā)展的趨勢,但當前瓠瓜市場缺乏高端產(chǎn)品質(zhì)量標準,瓠瓜品質(zhì)受多種因素影響且各因素影響效果未明,缺乏品質(zhì)代表性指標,本課題組前期研究找到游離谷氨酸含量是影響?zhàn)硝r味的主要因子之一,然而,由于瓠瓜品質(zhì)為復雜性狀,更多的品質(zhì)因子尚未鑒定,缺乏標準的客觀的品質(zhì)評價體系,因此,尋求簡單、準確的瓠瓜品質(zhì)評估方法是瓠瓜產(chǎn)業(yè)提質(zhì)增效的重要任務之一。
在園藝作物品質(zhì)評價指標的選擇上,學者們大多采用主成分分析法,如臍橙品質(zhì)評價因子的選擇研究[3],桃、楊梅、棗等的品質(zhì)綜合分析以及加工型黃瓜、番茄等的品質(zhì)分析和分類[4-8]。因子分析方法,與主成分分析方法類似,也是一種較為常用的進行多指標綜合的多元統(tǒng)計方法[9],其基本原理是通過降維步驟,將多個測量指標轉(zhuǎn)化為少數(shù)的、相互獨立的新指標,再根據(jù)各樣品的因子得分,客觀、合理地綜合評價[10]。目前,在砂梨、制干棗、無花果品質(zhì)指標分析與綜合評價、甜瓜果實品質(zhì)綜合評價和種質(zhì)材料選擇等方面已有應用[11-14]。
本研究以206份瓠瓜自然群體為試材,對可能影響?zhàn)瞎麑嵠焚|(zhì)的瓜長、瓜橫徑、肉厚、單瓜重、含水量、可溶性固形物含量、硬度、蛋白質(zhì)含量以及20種氨基酸,包括天冬氨酸(Asp)、蘇氨酸(Thr)、絲氨酸(Ser)、谷氨酸(Glu)、甘氨酸(Gly)、丙氨酸(Ala)、半胱氨酸(Cys)、纈氨酸(Val)、蛋氨酸(Met)、異亮氨酸(Ile)、亮氨酸(Leu)、酪氨酸(Tyr)、苯丙氨酸(Phe)、γ-氨基丁酸(g-ABA)、賴氨酸(Lys)、游離氨(NH3)、組氨酸(His)、精氨酸(Arg)、脯氨酸(Pro)、總氨基酸含量等28項指標,綜合運用相關(guān)性分析、因子分析、多元線性回歸等多種多元統(tǒng)計分析方法,明確了各指標間的相互關(guān)系,選擇了瓠瓜果實品質(zhì)的代表性指標,構(gòu)建了瓠瓜品質(zhì)量化評分模型并進一步綜合評價了206份瓠瓜種質(zhì)材料的品質(zhì)。
本研究所用的206份瓠瓜種質(zhì)材料[15]均由浙江省農(nóng)業(yè)科學院蔬菜研究所豇豆瓠瓜研究室提供。所有材料均于2021年春季種植于浙江省農(nóng)業(yè)科學院海寧楊渡實驗創(chuàng)新基地(30°N,120°E)。采用人字架栽培,畦寬連溝1.6 m。雙行種植,行距75 cm,株距50 cm,每份材料種植20株。相同栽培管理措施。
參照Wu等[16]的測定方法,在瓠瓜果實商品成熟期(花后7~10 d),采摘瓜型勻稱、大小一致的嫩果,每份材料采摘9個商品瓜,3個果實一組,每組均取果實中間段果肉,混合制成勻漿,每份材料3個重復。檢測時,每份勻漿樣品稱取1.5 g,放入15 mL離心管中,加入3%水楊酸準確定容至8 mL,室溫下振蕩1 h,12 000 r·min-1離心10 min,用一次性注射器抽取上清液,用400目水相過濾器過濾上清液,用日立L-8900氨基酸自動分析儀檢測樣品中各游離氨基酸的含量。
基于質(zhì)構(gòu)儀穿刺檢測瓠瓜果肉硬度,具體操作步驟為:在瓠瓜果實商品成熟期,采摘瓜型勻稱、大小一致的嫩果,每份材料采摘3個商品瓜,每個商品瓜均取果實中間段果肉,并切成3個邊長約2 cm的正方體塊,分別將每個小正方體塊置于探頭下,將TMS-Pro型質(zhì)構(gòu)儀設置為穿刺模式,使用100 N穿刺程序,調(diào)最大的力為96 N,起始力為0.05 N,穿刺距離10 mm,穿刺速度30 mm·min-1,回程距離20 mm,回程速度50 mm·min-1,最終取穿刺果肉的最大力作為衡量果肉硬度的指標。
采用硫酸消解,凱氏定氮法測定粗蛋白含量,具體步驟包括:準確稱取凍干樣品0.2 g(精確值0.1 mg)左右,加入消化管,再加入混合催化劑3 g K2SO4、0.2 g CuSO4,加入10 mL濃硫酸。利用石墨消解爐進行消解,將消化管放在石墨爐上,蓋上排氣罩,連接廢氣吸收系統(tǒng),消化過程采用曲線升溫模式,設定消解參數(shù)。消解樣品到澄清狀態(tài),K1100全自動凱氏定氮儀測定。
本試驗所有數(shù)據(jù)均使用 Excel 軟件整理,利用SPSS 22.0統(tǒng)計軟件進行28項瓠瓜果實品質(zhì)指標間的相關(guān)性分析、因子分析、多元線性回歸分析[17]。瓜長、瓜橫徑、單果重數(shù)據(jù)為6個重復,其他指標均設3個重復。
瓠瓜果實品質(zhì)28項指標的平均值、最大值、最小值、標準差和變異系數(shù)見表1,各指標變異系數(shù)在2.71%~63.60%。瓠瓜果實品質(zhì)各指標中,變異系數(shù)最大的是游離谷氨酸(Glu)含量,為63.60%,變異系數(shù)最小的是含水量,僅為2.71%,其他指標的變異系數(shù)在23.36%~54.67%。從表1可見,瓠瓜果實的游離氨基酸含量和氨基酸總量的變異系數(shù)均較大,在33.75%~63.60%,此外,硬度和瓜長的變異系數(shù)也較大。除瓠瓜果實含水量之外,蛋白質(zhì)含量的變異系數(shù)最小,為23.36%。
對瓠瓜28項不同品質(zhì)性狀相關(guān)指標進行了相關(guān)性分析,從圖1可知,28項果實品質(zhì)指標間表現(xiàn)出不同程度的正相關(guān)性和負相關(guān)性。其中20種游離氨基酸兩兩之間呈極顯著正相關(guān)(P<0.01);瓜橫徑、肉厚、單瓜重兩兩之間呈極顯著正相關(guān)(P<0.01);瓜長、瓜橫徑與肉厚呈極顯著負相關(guān),此外,瓜長還與硬度、Glu含量、Gly含量、Met含量等呈極顯著負相關(guān);可溶性固形物含量與含水量呈極顯著負相關(guān),與硬度呈極顯著正相關(guān);蛋白質(zhì)含量與瓜長、含水量呈極顯著正相關(guān),與Glu含量呈顯著正相關(guān),而與肉厚呈極顯著負相關(guān)。Glu含量與瓜長、瓜橫徑、肉厚、單瓜重存在不同程度的負相關(guān)性,而與硬度和蛋白質(zhì)含量呈不同程度的正相關(guān)性;總游離氨基酸除了與不同游離氨基酸呈極顯著正相關(guān)外,還與瓜橫徑、肉厚和單瓜重呈極顯著負相關(guān),與硬度和蛋白質(zhì)呈極顯著正相關(guān)(圖1)。
1~28分別代表瓜長、瓜橫徑、肉厚、單瓜重、含水量、可溶性固形物含量、硬度、蛋白質(zhì)含量,以及天冬氨酸、蘇氨酸、絲氨酸、谷氨酸、甘氨酸、丙氨酸、半胱氨酸、纈氨酸、蛋氨酸、異亮氨酸、亮氨酸、酪氨酸、苯丙氨酸、γ-氨基丁酸、賴氨酸、游離氨、組氨酸、精氨酸、脯氨酸、總氨基酸含量。1-28 represents fruit length, fruit diameter, flesh thickness, single fruit weight, water content, soluble solids content, hardness, protein content, and the contents of Asp, Thr, Ser, Glu, Gly, Ala, Cys, Val, Met, Ile, Leu, Tyr, Phe, GABA, Lys, NH3, His, Arg, Pro, total amino acids, respectively.
將瓠瓜果實品質(zhì)性狀相關(guān)的28項指標進行因子分析,標準為旋轉(zhuǎn)后的因子載荷矩陣和起始特征值大于1,如表2,得到F1、F2、F3、F4、F5、F6、F7、F8和F9共9個公因子。以旋轉(zhuǎn)元件矩陣中的載荷值的絕對值大于0.8的指標為解釋指標(下同)。F1解釋指標較多,主要集中于游離氨基酸類,分別為Thr、Glu、Gly、Val、Ile、Leu、His和總氨基酸含量;F2解釋指標為瓜長;F3解釋指標為瓜橫徑和肉厚;F5解釋指標為可溶性固形物含量;F7解釋指標為蛋白質(zhì)含量;F8解釋指標為含水量;F9解釋指標為硬度。按照該標準F4和F6無解釋指標。F1~F9方差貢獻率分別為38.46%、5.64%、9.34%、3.46%、4.36%、4.90%、4.31%、5.14%和11.17%,累計方差貢獻率 86.78%(>85%),可以反映瓠瓜果實品質(zhì)性狀相關(guān)指標的足夠信息[17]。因此,依據(jù)各因子方差貢獻率,建立瓠瓜品質(zhì)量化評分模型Y品質(zhì)性狀=0.3846F1+0.0564F2+0.0934F3+0.0346F4+0.0436F5+0.0490F6+0.0431F7+0.0514F8+0.1117F9。進一步計算出各瓠瓜種質(zhì)材料的Y品質(zhì)性狀得分,以瓠瓜種質(zhì)材料的Y品質(zhì)性狀得分為因變量,9個公因子所包含的15個指標為自變量,利用逐步回歸分析,得到3個有效預測模型(R2>0.9,下同),任選其中一個有效預測模型,計算各瓠瓜種質(zhì)材料的Y品質(zhì)性狀-2得分。利用相關(guān)性分析分析各瓠瓜種質(zhì)材料的Y品質(zhì)性狀得分和Y品質(zhì)性狀-2得分的相關(guān)性,結(jié)果表明,Y品質(zhì)性狀得分和Y品質(zhì)性狀-2得分極顯著相關(guān)(0.998 2**),說明瓠瓜品質(zhì)量化評分模型Y具有較高的預測準確性,可以用該模型對瓠瓜206份種質(zhì)材料進行評價,Thr、Glu、Gly、Val、Ile、Leu、His含量,以及總氨基酸含量、瓜長、瓜橫徑、肉厚、可溶性固形物含量、蛋白質(zhì)含量、含水量、硬度這15項指標可以代表瓠瓜果實品質(zhì)性狀及其在不同材料間的差異。
表2 瓠瓜果實品質(zhì)性狀旋轉(zhuǎn)后的因子載荷值
利用瓠瓜果實品質(zhì)性狀相關(guān)的28項指標數(shù)據(jù),采用歐式距離、組間平均鏈接法[18]和主成分分析等結(jié)果將瓠瓜果實各品質(zhì)指標進行聚類分析,利用SPSS22.0軟件計算輸出聚類分析樹形圖,如圖2。從圖中可以看出,聚類分析將瓠瓜果實28個品質(zhì)相關(guān)因子聚成5類,可溶性固形物含量和瓜長2個指標聚成一類,Pro、Lys、Met、g-ABA、Asp、Cys含量和硬度這7個指標聚成一類,Leu、Val、Ile、His、Phe、Tyr、Arg、總氨基酸、Ser、Ala、Gly、Thr、NH3和Glu含量這14個指標聚成一類,蛋白質(zhì)和含水量這2個指標聚成一類,單瓜重、肉厚和瓜橫徑這3個指標聚成一類。不同類的各指標間具有相對獨立性,同一類的各指標間具有較強的正相關(guān)性。
圖2 瓠瓜不同品質(zhì)指標的聚類分析Fig.2 Cluster analysis of different quality traits of bottle gourd
本著指標易用和避免重復的原則,根據(jù)瓠瓜果實品質(zhì)性狀相關(guān)的28項指標數(shù)據(jù)的相關(guān)性分析和聚類分析結(jié)果,結(jié)合前人研究經(jīng)驗,進行了指標簡化??扇苄怨绦挝锖亢凸祥L聚成一類,選擇可溶性固形物含量為代表性指標;Pro、Lys、Met、g-ABA、Asp、Cys含量和硬度這7個指標聚成一類,相關(guān)性分析結(jié)果顯示,除Thr含量之外,硬度與各個游離氨基酸及總氨基酸含量均呈極顯著正相關(guān),所以這類選擇硬度為代表性指標;Leu、Val、Ile、His、Phe、Tyr、Arg、總氨基酸、Ser、Ala、Gly、Thr、NH3和Glu含量這14個指標聚成一類,Glu含量是前期鑒定的瓠瓜鮮味關(guān)鍵因子之一[16],而且其與各個游離氨基酸及總氨基酸含量均呈極顯著正相關(guān),故這類選擇Glu含量為代表性指標;蛋白質(zhì)和含水量這2個指標聚成一類,相關(guān)性分析顯示,二者呈極顯著正相關(guān),選擇蛋白質(zhì)為代表性指標;單瓜重、肉厚和瓜橫徑這3個指標聚成一類,相關(guān)性分析顯示,肉厚與單瓜重和瓜橫徑呈極顯著正相關(guān),選擇肉厚為代表性指標。最終,將Thr、Glu、Gly、Val、Ile、Leu、His、總氨基酸含量,以及瓜長、瓜橫徑、肉厚、可溶性固形物含量、蛋白質(zhì)含量、含水量、硬度這15項指標簡化為肉厚、可溶性固形物含量、硬度、Glu含量、蛋白質(zhì)含量這5項代表性指標,用于綜合評估瓠瓜果實品質(zhì)。
根據(jù)各因子得分及對應方差貢獻率,分別以第1公因子和第2公因子、第1公因子和第3公因子作二維排序圖(圖3)。從圖3-A可見,材料Z-126、Z-127 第1公因子較大,即各種游離氨基酸含量比較高。結(jié)合第2公因子,材料Z-25和Z-179瓜長最長,實際數(shù)據(jù)Z-25和Z-179瓜長確實較長,分別為52.67和56.33 cm。從圖3-B可見,材料Z-170第3公因子最小,即瓜橫徑和肉厚較大,實際數(shù)據(jù)測得材料Z-170瓜橫徑為13.67 cm,肉厚為1.93 cm,為所有材料靠前。結(jié)合瓠瓜品質(zhì)量化評分模型Y品質(zhì)性狀=0.3846F1+0.0564F2+0.0934F3+0.0346F4+0.0436F5+0.0490F6+0.0431F7+0.0514F8+0.1117F9,206份瓠瓜種質(zhì)材料綜合品質(zhì)得分見表3,各因子得分與綜合得分排名與二維排序結(jié)果基本一致,進一步驗證了本品質(zhì)評價模型的準確性。此外,總排名靠前的為品質(zhì)較優(yōu)質(zhì)的種質(zhì)資源,可在育種中選擇應用。
圖3 206份瓠瓜種質(zhì)資源各因子二維排序圖Fig.3 Two dimensional ordered graph of each factor of 206 bottle gourd germplasm collections
表3 206份瓠瓜種質(zhì)資源綜合品質(zhì)得分表
果實品質(zhì)是決定瓠瓜市場競爭力的關(guān)鍵因素,其是由多因子構(gòu)成的綜合、復雜性狀。瓠瓜商品成熟果實生食時有微甜味,清炒或作湯后呈現(xiàn)較強的鮮味。鮮味是瓠瓜最重要的品質(zhì)性狀之一,本課題組前期研究發(fā)現(xiàn)游離谷氨酸含量對瓠瓜鮮味起重要作用[16]。然而,游離谷氨酸含量僅是呈鮮物質(zhì)中的一種,實地分析考察發(fā)現(xiàn)瓠瓜風味受眾多因素影響,包括含水量、硬度、可溶性固形物含量、蛋白質(zhì)含量、其他氨基酸含量等。此外,還與瓜長、瓜橫徑、單瓜重等因素相關(guān),如不同果形的瓠瓜品質(zhì)不同,一般來說,長棒形瓠瓜的品質(zhì)低于近圓形、牛腿形、葫蘆形或短棒形等瓠瓜。篩選瓠瓜品質(zhì)構(gòu)成因子及建立瓠瓜品質(zhì)綜合評價技術(shù)體系將有助于主要瓠瓜種質(zhì)資源的綜合評價和高品質(zhì)優(yōu)良品種的選育。
本研究研究了206份瓠瓜28項果實品質(zhì)相關(guān)指標,變異分析發(fā)現(xiàn)各指標變異系數(shù)在2.71%~63.60%。變異系數(shù)最大的是游離谷氨酸含量,為63.60%,硬度和瓜長的變異系數(shù)也較大。變異系數(shù)最小的是含水量,變異系數(shù)僅為2.71%,除瓠瓜果實含水量之外,變異系數(shù)最小的為蛋白質(zhì)含量,為23.36%。在前人瓜類作物研究中,單果重變異系數(shù)一般最大,如閆洪朗等[19]研究發(fā)現(xiàn)65個甜瓜品種的數(shù)量性狀指標中,變異系數(shù)最大的為單果重,為52.33%。尚建立等[20]評測了厚皮甜瓜和薄皮甜瓜遺傳多樣性,結(jié)果顯示,單果重變異系數(shù)也是最大,達74.3%。本研究中瓠瓜單瓜重的變異系數(shù)為36.91%(表1),較前人同為瓜類的研究偏低,可能與試驗選用的瓠瓜種質(zhì)材料中無野生資源有關(guān),野生資源一般偏小,單瓜重偏低,可增大群體的遺傳多樣性。
研究中206份瓠瓜種質(zhì)資源的28項果實品質(zhì)相關(guān)指標相關(guān)性分析顯示絕大多數(shù)品質(zhì)相關(guān)指標存在顯著或極顯著相關(guān)性,說明這些性狀間既相對獨立,又密切相關(guān)。目前瓠瓜綜合品質(zhì)評價指標過多且仍以主觀判斷為主,缺乏系統(tǒng)、科學的種質(zhì)材料瓠瓜評價技術(shù)體系。因子分析作為最常用的統(tǒng)計分析方法之一,可最大程度保留原有信息,將多個變量綜合為少數(shù)幾個相互獨立的“因子”,從而有效代表數(shù)目較多且彼此相關(guān)的指標[9-10],在園藝作物品質(zhì)綜合評價和種質(zhì)材料選擇等方面已有應用[11-14]。本研究利用因子分析方法,首先建立了瓠瓜品質(zhì)量化評分模型,并篩選出Thr、Glu、Gly、Val、Ile、Leu、His、總氨基酸含量,以及瓜長、瓜橫徑、肉厚、可溶性固形物含量、蛋白質(zhì)含量、含水量、硬度這15項品質(zhì)相關(guān)指標,代表了瓠瓜果實品質(zhì)性狀及其在不同材料間的差異。進一步,我們利用28項品質(zhì)指標數(shù)據(jù)的相關(guān)性分析和聚類分析結(jié)果,將指標簡化為肉厚、可溶性固形物含量、硬度、Glu含量、蛋白質(zhì)含量這5項代表性指標,并綜合評估了206份瓠瓜種質(zhì)材料果實品質(zhì)。后續(xù)在瓠瓜高品質(zhì)優(yōu)良品種的選育中,可根據(jù)不同種質(zhì)材料的不同優(yōu)勢,選育和推廣滿足不同市場需求的瓠瓜新品種。