王珮璇，陳 靜，郭鈺柬，周開兵

（海南大學(xué) 熱帶作物新品種選育教育部工程研究中心，海南 ?？?570228）

油茶是山茶科Theaceae山茶屬Camellia植物中油脂含量較高且具有一定栽培規(guī)模的一類植物的總稱[1]，是中國(guó)和世界四大木本油料作物之一，被稱為“東方的橄欖樹”[2]。茶籽油是優(yōu)質(zhì)保健食用油和高級(jí)天然化妝品原料，其副產(chǎn)品在農(nóng)業(yè)、醫(yī)療、日化等方面均有應(yīng)用價(jià)值[3-7]。越南油茶Camellia vietnamensis是中國(guó)的廣東省、廣西壯族自治區(qū)、海南省和越南、老撾等國(guó)家和地區(qū)特產(chǎn)的油茶樹種[8]，也是我國(guó)海南省種植的主要油茶樹種，其栽培歷史至少有500多年[9]，海南省不乏野生原生種和老林分布。在長(zhǎng)期適應(yīng)南亞熱帶環(huán)境的過(guò)程中，越南油茶已形成了與內(nèi)陸油茶不同的特點(diǎn)，具有抗逆性強(qiáng)、耐熱性強(qiáng)、壽命長(zhǎng)與樹冠高大、結(jié)果早、果大籽多等優(yōu)點(diǎn)，是很有發(fā)展前途的油茶樹種之一；同時(shí)能用于山地造林，具有生態(tài)保護(hù)功能；并且越南油茶籽油在口感、品質(zhì)、生產(chǎn)加工、市場(chǎng)潛力等方面都有獨(dú)特的優(yōu)勢(shì)。但是，長(zhǎng)期以來(lái)，由于地理的隔離及相關(guān)科研、產(chǎn)業(yè)發(fā)展的滯后，海南省的越南油茶產(chǎn)業(yè)并沒有得到應(yīng)有的重視[10-11]。

目前，正逢海南國(guó)際自由貿(mào)易港建設(shè)機(jī)遇，人們對(duì)海南產(chǎn)區(qū)越南油茶油特有食用和保健價(jià)值的認(rèn)識(shí)不斷加深，對(duì)其需求量加大，選育海南當(dāng)?shù)貎?yōu)良的越南油茶品種已成為海南省越南油茶產(chǎn)業(yè)發(fā)展的當(dāng)務(wù)之急[10]。前人已有一些包括越南油茶在內(nèi)的油茶豐產(chǎn)單株的選擇和產(chǎn)區(qū)比較的研究報(bào)道，但其研究局限于同一產(chǎn)區(qū)不同單株間、不同年份間和同一年份不同產(chǎn)區(qū)間的比較分析[12-13]；也有研究者針對(duì)個(gè)別產(chǎn)區(qū)油茶的農(nóng)藝性狀、產(chǎn)量性狀、茶籽油理化性質(zhì)和脂肪酸組成等品質(zhì)特點(diǎn)，探討了優(yōu)株篩選的標(biāo)準(zhǔn)和方法[14-16]。然而，針對(duì)海南省的越南油茶單株及產(chǎn)區(qū)的評(píng)價(jià)研究尚未見諸報(bào)道。

灰色系統(tǒng)理論以其簡(jiǎn)單便捷又不失科學(xué)合理等優(yōu)點(diǎn)，被廣泛應(yīng)用于如水稻、小麥、玉米、花生等[17-18]農(nóng)作物品種的綜合評(píng)價(jià)及其性狀間的相關(guān)分析中。在對(duì)油茶的研究中，有關(guān)研究者采用這種分析方法研究過(guò)其他地區(qū)的普通油茶的各性狀指標(biāo)與產(chǎn)油量之間關(guān)系的密切程度、氣象等因素與茶籽油產(chǎn)量及品質(zhì)間的關(guān)系、茶油中活性成分之間的關(guān)系等問題[19-21]。而采用此法研究越南油茶的報(bào)道極少，僅見有聶浪等[22]對(duì)包括越南油茶在內(nèi)的廣東油茶種質(zhì)資源篩選及評(píng)價(jià)的研究報(bào)道。為給越南油茶種質(zhì)資源育種的開發(fā)利用提供數(shù)據(jù)支持，同時(shí)為海南省高產(chǎn)油茶基地的規(guī)劃建設(shè)提供參考依據(jù)，本研究對(duì)海南省的6個(gè)越南油茶產(chǎn)區(qū)實(shí)生老林中的32個(gè)單株的農(nóng)藝性狀、產(chǎn)量性狀和茶籽油品質(zhì)進(jìn)行了比較分析，并采用灰色關(guān)聯(lián)分析法對(duì)優(yōu)良單株與豐產(chǎn)穩(wěn)產(chǎn)產(chǎn)區(qū)進(jìn)行了篩選和評(píng)價(jià)，現(xiàn)將研究結(jié)果分析報(bào)道如下。

1 材料與方法

1.1 單株的選擇與果實(shí)采樣地的概況

分別從海南省的保梅嶺生態(tài)保護(hù)區(qū)（簡(jiǎn)稱為“保梅嶺”）、霸王嶺長(zhǎng)臂猿保護(hù)區(qū)（簡(jiǎn)稱為“霸王嶺”）、澄邁縣福山鎮(zhèn)長(zhǎng)安村（簡(jiǎn)稱為“澄邁”）、白沙縣福門村（簡(jiǎn)稱為“白沙”）、儋州市王五鎮(zhèn)興文村（簡(jiǎn)稱為“儋州”）、瓊海市陽(yáng)江鎮(zhèn)上科村（簡(jiǎn)稱為“瓊?！保┑仍侥嫌筒璁a(chǎn)區(qū)中選定2017和2018年連續(xù)兩年豐產(chǎn)（株產(chǎn)果量在10 kg以上）、穩(wěn)產(chǎn)單株32株作為試驗(yàn)材料，各產(chǎn)區(qū)中的環(huán)境概況和豐產(chǎn)單株的編號(hào)與樹齡等信息見表1。

表1 海南省6個(gè)越南油茶產(chǎn)區(qū)的環(huán)境概況和豐產(chǎn)單株的編號(hào)與樹齡?Table 1 Environmental profiles of 6 Vietnamese oil tea producing areas in Hainan Province and the number of high yield trees with tree age

1.2 研究方法

在除瓊海產(chǎn)區(qū)（因?yàn)樵摦a(chǎn)區(qū)只篩選出1株高產(chǎn)單株樣樹，所以在產(chǎn)區(qū)評(píng)價(jià)時(shí)，此產(chǎn)區(qū)的相關(guān)數(shù)據(jù)不包含在內(nèi)）外的5個(gè)產(chǎn)區(qū)內(nèi)，選擇越南油茶豐產(chǎn)單株，對(duì)各產(chǎn)區(qū)不同年份間和同一年份不同產(chǎn)區(qū)間其主要農(nóng)藝性狀、產(chǎn)量性狀、茶籽油品質(zhì)性狀進(jìn)行比較，并采用灰色關(guān)聯(lián)分析法，對(duì)所有單株和產(chǎn)區(qū)越南油茶的主要產(chǎn)量性狀、茶籽油的品質(zhì)性狀進(jìn)行了調(diào)查與分析，并對(duì)其中的優(yōu)良單株與豐產(chǎn)產(chǎn)區(qū)進(jìn)行了初步的篩選與評(píng)價(jià)。

1.3 取樣方法與樣品處理

分別于2017和2018年的10月即海南省的越南油茶成熟采收期，在選定的6個(gè)產(chǎn)區(qū)的32個(gè)供試單株樹冠中部外圍的東、西、南、北4個(gè)方位隨機(jī)選取當(dāng)年生春梢，在其中部采集無(wú)病蟲害的成熟葉片，每個(gè)單株各采集樣葉30片；在每株樣樹樹冠中部的東、西、南、北4個(gè)方位隨機(jī)采摘無(wú)病蟲害的果實(shí)30個(gè)。將所采葉片和果實(shí)樣品按照單株分裝帶回實(shí)驗(yàn)室，并將所采果實(shí)置于室溫下后熟至果皮開裂時(shí)取出種子以備用。

1.4 測(cè)定指標(biāo)與測(cè)定方法

每株樣樹各選葉片樣品30片、果實(shí)樣品30個(gè)，采用常規(guī)方法調(diào)查各單株樹體特征與葉片、果實(shí)的大小和形狀等農(nóng)藝性狀及單果質(zhì)量、出籽率、出仁率、含油率和單位樹冠面積果實(shí)產(chǎn)量、單位面積毛油產(chǎn)量等產(chǎn)量性狀表現(xiàn)，以某一產(chǎn)區(qū)中所有樣株的平均值作為該產(chǎn)區(qū)的觀測(cè)值，所有數(shù)值均為（平均值±標(biāo)準(zhǔn)差）。

采用索氏提取法提取茶籽油，參照國(guó)家標(biāo)準(zhǔn)《食品中過(guò)氧化值的測(cè)定》（GB 5009.227—2016）[23]、《動(dòng)植物油脂碘值的測(cè)定》（GB/T 5532—2008）[24]、《食品中酸價(jià)的測(cè)定》（GB 5009.229—2016）[25]、《動(dòng)植物油脂皂化值的測(cè)定》（GB/T 5534—2008）[26]分別測(cè)定茶籽油的過(guò)氧化值、碘值、酸價(jià)、皂化值等理化性質(zhì)指標(biāo)值；按照《食品中脂肪酸的測(cè)定》（GB 5009.168—2016）[27]中說(shuō)明的具體要求制備脂肪酸甲酯，采用氣相色譜法檢測(cè)脂肪酸的種類與含量。每個(gè)單株的待測(cè)樣品各制備3份，取其平均值作為該單株的觀測(cè)值，所有數(shù)值亦均為（平均值±標(biāo)準(zhǔn)差）。

1.5 數(shù)據(jù)處理方法

采用SAS軟件中的GLM過(guò)程作方差分析，采用DUNCAN法對(duì)不同產(chǎn)區(qū)作多重比較分析，用TTEST過(guò)程對(duì)各產(chǎn)區(qū)年份間的差異顯著性進(jìn)行分析。用Excel軟件編輯公式，采用灰色關(guān)聯(lián)分析法[17-21]分別根據(jù)豐產(chǎn)性和茶籽油品質(zhì)等性狀指標(biāo)的觀測(cè)值對(duì)不同單株和產(chǎn)區(qū)進(jìn)行評(píng)價(jià)及排序。

2 結(jié)果與分析

2.1 同一年份不同產(chǎn)區(qū)間和同一產(chǎn)區(qū)不同年份間越南油茶農(nóng)藝性狀的比較

同一年份不同產(chǎn)區(qū)間和同一產(chǎn)區(qū)不同年份間越南油茶農(nóng)藝性狀的比較結(jié)果見表2。2017年，不同產(chǎn)區(qū)間越南油茶的葉片長(zhǎng)度、果實(shí)縱徑、果形指數(shù)、果皮厚度的差異均顯著，而其余農(nóng)藝性狀的差異均不顯著。2018年，不同產(chǎn)區(qū)間越南油茶的果縱徑、果形指數(shù)、果皮厚度均無(wú)顯著差異，而其余農(nóng)藝性狀的差異均顯著。2017年，越南油茶的葉片長(zhǎng)度，霸王嶺產(chǎn)區(qū)的顯著最長(zhǎng)，而儋州產(chǎn)區(qū)的顯著最短；儋州產(chǎn)區(qū)與保梅嶺、澄邁產(chǎn)區(qū)間的差異均不顯著，白沙產(chǎn)區(qū)與其他4個(gè)產(chǎn)區(qū)間的差異均不顯著。2018年，越南油茶的葉片寬度，霸王嶺產(chǎn)區(qū)的顯著最大，而澄邁產(chǎn)區(qū)的顯著最??；澄邁產(chǎn)區(qū)與白沙產(chǎn)區(qū)間的差異不顯著，保梅嶺和儋州產(chǎn)區(qū)的居中，且其與最大、最小的兩組間均無(wú)顯著差異。2018年，越南油茶的葉形指數(shù)，白沙產(chǎn)區(qū)的顯著最大，儋州產(chǎn)區(qū)的顯著最?。毁僦莓a(chǎn)區(qū)與其他3個(gè)產(chǎn)區(qū)間的差異均不顯著。2017年越南油茶的果實(shí)縱徑和2018年的果實(shí)橫徑，白沙產(chǎn)區(qū)的均顯著最大，儋州產(chǎn)區(qū)的均顯著最小。2017年越南油茶的果形指數(shù)和果皮厚度，澄邁產(chǎn)區(qū)的均顯著最小，白沙產(chǎn)區(qū)的均顯著最大；白沙產(chǎn)區(qū)與其他3個(gè)產(chǎn)區(qū)間的差異均不顯著。將不同年份間的單株觀測(cè)值進(jìn)行比較，僅有保梅嶺產(chǎn)區(qū)的葉長(zhǎng)、霸王嶺產(chǎn)區(qū)的果實(shí)橫徑和儋州產(chǎn)區(qū)的果皮厚度均表現(xiàn)出顯著差異。這一觀測(cè)結(jié)果表明，不論同一產(chǎn)區(qū)不同年份還是同一年份不同產(chǎn)區(qū)，越南油茶的農(nóng)藝性狀都會(huì)出現(xiàn)不同的變化情況，說(shuō)明其農(nóng)藝性狀的表現(xiàn)是因環(huán)境因素的變化而變化的，意味著越南油茶新品種的選育研究宜基于多年多點(diǎn)的田間試驗(yàn)結(jié)果進(jìn)行綜合分析，還意味著越南油茶栽培技術(shù)的標(biāo)準(zhǔn)化是十分重要的。

表2 同一年份不同產(chǎn)區(qū)間和同一產(chǎn)區(qū)不同年份間越南油茶農(nóng)藝性狀的比較結(jié)果?Table 2 Comparison results of agronomic traits of C.vietnamensis between different producings of the same year and between different years of the same producing

2.2 同一年份不同產(chǎn)區(qū)間和同一產(chǎn)區(qū)不同年份間越南油茶產(chǎn)量性狀的比較

同一年份不同產(chǎn)區(qū)間和同一產(chǎn)區(qū)不同年份間越南油茶產(chǎn)量性狀的比較結(jié)果見表3。2017年，不同產(chǎn)區(qū)間越南油茶單果質(zhì)量和干種仁含油率的差異均顯著，儋州產(chǎn)區(qū)的單果質(zhì)量顯著最小，而其干種仁含油率顯著最大，其余4個(gè)產(chǎn)區(qū)間這2個(gè)性狀指標(biāo)的差異均不顯著；而不同產(chǎn)區(qū)間其余產(chǎn)量性狀均無(wú)顯著差異。2018年，不同產(chǎn)區(qū)間越南油茶單果質(zhì)量的差異顯著；保梅嶺產(chǎn)區(qū)的單果質(zhì)量顯著最大，保梅嶺產(chǎn)區(qū)與白沙產(chǎn)區(qū)間單果質(zhì)量的差異不顯著，而儋州產(chǎn)區(qū)的顯著最小，其余2個(gè)產(chǎn)區(qū)間其單果質(zhì)量的差異不顯著。在除瓊海產(chǎn)區(qū)外的5個(gè)產(chǎn)區(qū)中，僅有儋州產(chǎn)區(qū)的越南油茶干種仁含油率在不同年份間的差異顯著，而儋州產(chǎn)區(qū)除干種仁含油率外的其余產(chǎn)量性狀在不同年份間均無(wú)顯著差異。可見，各產(chǎn)區(qū)豐產(chǎn)單株群體的產(chǎn)量性狀表現(xiàn)基本穩(wěn)定。

表3 同一年份不同產(chǎn)區(qū)間和同一產(chǎn)區(qū)不同年份間越南油茶主要產(chǎn)量性狀的比較結(jié)果Table 3 Comparison results of yield traits of C.vietnamensis between different producings of the same year and between different years of the same producing

2.3 同一年份不同產(chǎn)區(qū)間和同一產(chǎn)區(qū)不同年份間越南油茶茶籽油品質(zhì)的比較

2.3.1 茶籽油的理化性質(zhì)

同一年份不同產(chǎn)區(qū)間和同一產(chǎn)區(qū)不同年份間茶籽油理化性質(zhì)的比較結(jié)果見表4。由表4可知，在2017和2018年連續(xù)兩年里，不同產(chǎn)區(qū)的茶籽油其碘值、酸價(jià)都符合GB /T 11765—2018《油茶籽油》[28]中的標(biāo)準(zhǔn)，且不同產(chǎn)區(qū)間均無(wú)顯著差異，說(shuō)明各產(chǎn)區(qū)的茶籽油其耐貯藏能力和新鮮度都較好，其不飽和脂肪酸的不飽和程度、耐酸敗能力均相當(dāng)。2017和2018年不同產(chǎn)區(qū)間越南油茶茶籽油過(guò)氧化值的差異顯著；2017年儋州產(chǎn)區(qū)的茶籽油其過(guò)氧化值顯著最高，且與白沙產(chǎn)區(qū)的差異不顯著；其他3個(gè)產(chǎn)區(qū)間其過(guò)氧化值的差異均不顯著；而2018年霸王嶺產(chǎn)區(qū)的茶籽油其過(guò)氧化值顯著最高，其他4個(gè)產(chǎn)區(qū)間其過(guò)氧化值的差異均不顯著。將同一產(chǎn)區(qū)不同年份間的觀測(cè)值進(jìn)行比較，霸王嶺、白沙和儋州3個(gè)產(chǎn)區(qū)的茶籽油過(guò)氧化值在不同年份間的差異均顯著，卻都遠(yuǎn)低于GB /T 11765—2018《油茶籽油》[28]中規(guī)定的一級(jí)油閾值標(biāo)準(zhǔn)，說(shuō)明各產(chǎn)區(qū)不同年份的茶籽油其抗氧化性都很強(qiáng)，均有較高的保健價(jià)值且都耐貯藏。2017年，白沙產(chǎn)區(qū)的茶籽油其皂化值顯著最高，各產(chǎn)區(qū)的茶籽油其皂化值都符合GB /T 11765—2018《油茶籽油》[28]中規(guī)定的合格級(jí)別的閾值標(biāo)準(zhǔn)?？傊Ｄ喜煌a(chǎn)區(qū)的茶籽油其理化性質(zhì)均符合GB /T 11765—2018《油茶籽油》[28]中規(guī)定的標(biāo)準(zhǔn)，同一年份不同產(chǎn)區(qū)間和同一產(chǎn)區(qū)不同年份間茶籽油的過(guò)氧化值和皂化值均表現(xiàn)出明顯的差異，說(shuō)明這兩項(xiàng)理化性質(zhì)指標(biāo)與油茶生長(zhǎng)的環(huán)境變化均有密切的關(guān)系。

表4 同一年份不同產(chǎn)區(qū)間和同一產(chǎn)區(qū)不同年份間茶籽油理化性質(zhì)的比較結(jié)果Table 4 Comparison results of physical and chemical indicators of C.vietnamensis seed oil between different producings of the same year and between different years of the same producing

2.3.2 茶籽油的脂肪酸組成

不同產(chǎn)區(qū)間和不同年份間茶籽油脂肪酸組成的比較結(jié)果見表5。由表5可知，從不同產(chǎn)區(qū)的茶籽油中共檢測(cè)到8種主要脂肪酸，具體包括棕櫚酸、硬脂酸、花生酸、肉豆蔻酸等飽和脂肪酸和油酸、亞油酸、亞麻酸、花生—烯酸等不飽和脂肪酸，其中油酸含量位居8種脂肪酸之首，所有產(chǎn)區(qū)茶籽油的油酸含量均高于75%，且所有產(chǎn)區(qū)的茶籽油其主要脂肪酸組分含量均符合GB /T 11765—2018《油茶籽油》[28]中規(guī)定的標(biāo)準(zhǔn)。

表5 同一年份不同產(chǎn)區(qū)間和同一產(chǎn)區(qū)不同年份間茶籽油脂肪酸組分的比較結(jié)果Table 5 Comparison results of fatty acid composition of C.vietnamensis seed oil between different producings of the same year and between different years of the same producing %

在2017年產(chǎn)的茶籽油中，硬脂酸、油酸、亞油酸、亞麻酸、花生—烯酸等脂肪酸組分的含量在不同產(chǎn)區(qū)間的差異均不顯著，而肉豆蔻酸、花生酸、棕櫚酸等脂肪酸組分的含量在不同產(chǎn)區(qū)間的差異均顯著。肉豆蔻酸僅在儋州產(chǎn)區(qū)和保梅嶺產(chǎn)區(qū)的茶籽油中被發(fā)現(xiàn)，其在儋州產(chǎn)區(qū)的茶籽油中的含量顯著較高；花生酸含量，儋州產(chǎn)區(qū)的顯著最高，而保梅嶺產(chǎn)區(qū)的顯著最低，霸王嶺、澄邁、白沙3個(gè)產(chǎn)區(qū)的居中，且此3組與最高、最低的2組之間均無(wú)顯著差異；棕櫚酸含量，儋州產(chǎn)區(qū)的顯著最低，而霸王嶺、澄邁、白沙、保梅嶺這4個(gè)產(chǎn)區(qū)間的差異均不顯著。在2018年產(chǎn)的茶籽油中，棕櫚酸僅在保梅嶺、霸王嶺、儋州產(chǎn)區(qū)的茶籽油中被發(fā)現(xiàn)，保梅嶺產(chǎn)區(qū)的顯著最高，而其他2個(gè)產(chǎn)區(qū)間的差異不顯著；其余各種脂肪酸組分的含量在各產(chǎn)區(qū)間均無(wú)顯著差異。在保梅嶺產(chǎn)區(qū)的茶籽油中，肉豆蔻酸和亞麻酸的含量在不同年份間存在顯著差異；在澄邁和白沙產(chǎn)區(qū)的茶籽油中，棕櫚酸的含量在不同年份間也存在顯著差異；在霸王嶺產(chǎn)區(qū)的茶籽油中，亞麻酸和棕櫚酸等含量在不同年份間也存在顯著差異；而在其余產(chǎn)區(qū)的茶籽油中，其他脂肪酸組分的含量在不同年份間均無(wú)顯著差異。這一結(jié)果說(shuō)明，除亞麻酸以外的不飽和脂肪酸組分的含量都有較強(qiáng)的遺傳穩(wěn)定性，而多數(shù)飽和脂肪酸組分的含量則相反，這意味著采用不同的栽培技術(shù)可能會(huì)調(diào)節(jié)飽和脂肪酸和亞麻酸等脂肪酸組分的含量。因此，在新品種的選育研究中，只有在相同的環(huán)境條件下才能對(duì)茶籽油的品質(zhì)作出準(zhǔn)確的選擇，不飽和脂肪酸含量的穩(wěn)定性在越南油茶種質(zhì)資源鑒定和新品種選育中或可在其表現(xiàn)型標(biāo)記中用作參考。

2.4 灰色關(guān)聯(lián)分析結(jié)果

2.4.1 各個(gè)單株的灰色關(guān)聯(lián)度及其排序

在灰色關(guān)聯(lián)度分析中，首先需要結(jié)合研究目標(biāo)構(gòu)建理想個(gè)體，具備可以達(dá)到研究目的的最優(yōu)參考指標(biāo)。本研究以篩選海南越南油茶豐產(chǎn)、優(yōu)質(zhì)單株為目的，首先需要確定理想單株產(chǎn)量性狀的參考指標(biāo)，即兩年調(diào)查觀察到的產(chǎn)量性狀平均數(shù)的最大值，依據(jù)單果質(zhì)量和樹冠單位面積產(chǎn)果量、果實(shí)出籽率、種子出仁率、種仁含油率、樹冠單位面積產(chǎn)油量等產(chǎn)量指標(biāo)值越高越好的原則，即其最大值依次為56.722 g、2.0 45 kg·m-2、43.664%、65.422%、54.016%、1 001.108 kg·hm-2；茶油品質(zhì)性狀參考指標(biāo)的選擇有所不同，并不全是越大越好，需要按照國(guó)家質(zhì)量標(biāo)準(zhǔn)GB /T 11765—2018《油茶籽油》[28]來(lái)選擇，用測(cè)得的最優(yōu)值作為理想個(gè)體的參考指標(biāo)，得到茶籽油的過(guò)氧化值、酸價(jià)和棕櫚酸、硬脂酸、油酸、亞油酸、亞麻酸、花生酸的含量依次為0.010%、0.500 mg·g-1、9.200%、2.550%、77.500%、8.900%、0.160%、0.360%?；诟餍誀钪笜?biāo)兩年觀測(cè)值的平均值對(duì)32個(gè)單株作灰色關(guān)聯(lián)分析，結(jié)果見表6。由表6可知，產(chǎn)量性狀排名前10位的單株編號(hào)依次為8、14、2、19、9、21、16、13、1、20；比較其灰色關(guān)聯(lián)度的大小可知，2和19、16和13號(hào)單株相差并不明顯，即其產(chǎn)量性狀的綜合表現(xiàn)相似；結(jié)合樹冠單位面積產(chǎn)油量的排名情況進(jìn)行綜合分析，結(jié)果表明，8、14、19、9、21號(hào)等單株產(chǎn)量性狀的綜合表現(xiàn)均佳。從表6中還可看出，品質(zhì)性狀排名前10位的單株編號(hào)依次為24、6、16、20、9、8、11、15、18、23；比較其灰色關(guān)聯(lián)度的大小可知，20、9、8、11號(hào)單株之間，15、18、23號(hào)單株之間，其相差均不明顯，即其品質(zhì)性狀的綜合表現(xiàn)相似，表明這10個(gè)單株茶籽油的綜合品質(zhì)都較好。綜上所述，產(chǎn)量和品質(zhì)性狀綜合表現(xiàn)均進(jìn)入前10位的單株有8、9號(hào)單株，說(shuō)明這兩個(gè)單株品種有望上升為豐產(chǎn)、優(yōu)質(zhì)的品種。

表6 不同單株產(chǎn)量和品質(zhì)性狀的灰色關(guān)聯(lián)分析結(jié)果Table 6 The result of grey relational degree analysis of different single tree’s yield and quality traits

2.4.2 各個(gè)產(chǎn)區(qū)的灰色關(guān)聯(lián)度及其排序

與各個(gè)單株的灰色關(guān)聯(lián)度分析一樣，首先要以篩選優(yōu)勢(shì)產(chǎn)區(qū)為研究目標(biāo)構(gòu)建產(chǎn)區(qū)的理想個(gè)體，以某一產(chǎn)區(qū)內(nèi)全部單株各性狀指標(biāo)的平均值作為該產(chǎn)區(qū)的各性狀指標(biāo)值。采用相同的辦法確定理想產(chǎn)區(qū)的產(chǎn)量性狀和茶籽油品質(zhì)性狀指標(biāo)值，即理想產(chǎn)區(qū)的單果質(zhì)量、樹冠單位面積產(chǎn)果量、果實(shí)出籽率、種子出仁率、種仁含油率、樹冠單位面積產(chǎn)油量等產(chǎn)量性狀指標(biāo)值分別為40.914 g、1.01 kg·m-2、41.148%、48.949%、52.271%、477.099 kg·hm-2，理想產(chǎn)區(qū)的茶籽油的過(guò)氧化值、酸價(jià)和棕櫚酸、硬脂酸、油酸、亞油酸、亞麻酸、花生酸含量等茶油品質(zhì)性狀指標(biāo)值分別為0.02%、1 mg·g-1、9.2%、2.55%、77.5%、8.9%、0.1%、0.4%?；诟鳟a(chǎn)區(qū)全部單株的各性狀兩年平均觀測(cè)值對(duì)5個(gè)產(chǎn)區(qū)進(jìn)行灰色關(guān)聯(lián)分析，結(jié)果見表7。由表7可知，不同產(chǎn)區(qū)的產(chǎn)量性狀的灰色關(guān)聯(lián)位序依次為霸王嶺產(chǎn)區(qū)、保梅嶺產(chǎn)區(qū)、澄邁產(chǎn)區(qū)、白沙產(chǎn)區(qū)和儋州產(chǎn)區(qū)；比較其灰色關(guān)聯(lián)度的大小可知，霸王嶺產(chǎn)區(qū)和保梅嶺產(chǎn)區(qū)、白沙產(chǎn)區(qū)和儋州產(chǎn)區(qū)的差別并不明顯，即其產(chǎn)量性狀的綜合表現(xiàn)相似。這個(gè)排名與其樹冠單位面積產(chǎn)油量的排名一致，表明霸王嶺產(chǎn)區(qū)和保梅嶺產(chǎn)區(qū)的產(chǎn)量性狀綜合表現(xiàn)均佳。不同產(chǎn)區(qū)的品質(zhì)性狀排名依次為澄邁產(chǎn)區(qū)、保梅嶺產(chǎn)區(qū)、白沙產(chǎn)區(qū)、霸王嶺產(chǎn)區(qū)、儋州產(chǎn)區(qū)；比較其灰色關(guān)聯(lián)度的大小可知，5個(gè)產(chǎn)區(qū)相差并不明顯，說(shuō)明其品質(zhì)性狀的綜合表現(xiàn)相似。綜上所述，產(chǎn)量和品質(zhì)性狀綜合排名靠前的依次是霸王嶺產(chǎn)區(qū)和保梅嶺產(chǎn)區(qū)，居末位的是儋州產(chǎn)區(qū)，說(shuō)明霸王嶺產(chǎn)區(qū)和保梅嶺產(chǎn)區(qū)可能均為海南省越南油茶的優(yōu)勢(shì)產(chǎn)區(qū)。

表7 不同產(chǎn)區(qū)油茶產(chǎn)量和品質(zhì)性狀的灰色關(guān)聯(lián)分析結(jié)果Table 7 The result of grey relational degree analysis of the different producing regions’ yield and quality traits

3 討 論

本研究結(jié)果表明，不同產(chǎn)區(qū)間及不同年份間越南油茶的農(nóng)藝、產(chǎn)量和品質(zhì)性狀的差異普遍存在，但茶籽油的品質(zhì)性狀都能達(dá)到GB /T11765—2018《油茶籽油》[28]中的要求，因此海南省越南油茶產(chǎn)區(qū)的茶籽油無(wú)品質(zhì)問題，其主要的育種目標(biāo)仍然是豐產(chǎn)性、穩(wěn)產(chǎn)性的遺傳改良，這與本課題組前期多年的研究結(jié)論[13]一致。同時(shí)，不同單株茶籽油品質(zhì)差異的存在為選育具有地方特色的高產(chǎn)品種提供了良好的遺傳物質(zhì)基礎(chǔ)，這與前人對(duì)越南油茶的研究結(jié)論[22]一致。

本研究結(jié)果還表明，在不同年份或不同產(chǎn)區(qū)間越南油茶不同單株茶籽油的多個(gè)產(chǎn)量和品質(zhì)性狀均存在顯著差異。這一研究結(jié)果對(duì)其相關(guān)研究和栽培實(shí)踐具有以下兩個(gè)方面的指導(dǎo)意義：一方面，越南油茶實(shí)生選種的研究需要開展多年連續(xù)、多點(diǎn)同時(shí)進(jìn)行的田間試驗(yàn)；另一方面，將性狀與其較大影響因素結(jié)合起來(lái)進(jìn)行研究，有利于指導(dǎo)栽培實(shí)踐，因?yàn)檫@些性狀主要受環(huán)境條件和不同栽培技術(shù)的影響，意味著人為地改善油茶的栽培環(huán)境，如改善土質(zhì)、加強(qiáng)栽培技術(shù)指導(dǎo)、精準(zhǔn)施肥，或者建立優(yōu)質(zhì)商品生產(chǎn)基地，都利于提高越南油茶的果實(shí)產(chǎn)量和茶籽油的品質(zhì)，這與前人的研究結(jié)果[29-32]基本一致。

霸王嶺產(chǎn)區(qū)和保梅嶺產(chǎn)區(qū)的產(chǎn)量性狀的關(guān)聯(lián)度相差不明顯，分別排在第一與第二名，而其品質(zhì)性狀的關(guān)聯(lián)度分別排在第四與第二名，但是，由于5個(gè)產(chǎn)區(qū)的品質(zhì)性狀相差并不明顯，且均符合GB/T 11765—2018《油茶籽油》[28]中的要求，所以霸王嶺產(chǎn)區(qū)和保梅嶺產(chǎn)區(qū)可能均為越南油茶的優(yōu)勢(shì)產(chǎn)區(qū)。霸王嶺產(chǎn)區(qū)和保梅嶺產(chǎn)區(qū)的產(chǎn)量性狀相似，且其油茶產(chǎn)量均高于其余產(chǎn)區(qū)的，這得益于其相似的氣候條件、土壤類型、樹齡和管理程度，與其余3個(gè)產(chǎn)區(qū)相比，兩者所處山地?zé)釒в炅值娜站鶞馗摺囟饶贻^差大、雨量多、濕度大、日照短，這樣的氣候條件可能更適合越南油茶的栽培，這與前人篩選出的普通油茶最優(yōu)產(chǎn)區(qū)的環(huán)境條件有明顯差異[31]，可能因?yàn)樵侥嫌筒枧c普通油茶對(duì)環(huán)境條件的要求不同，這是不同物種遺傳差異的重要表現(xiàn)。

研究中發(fā)現(xiàn)，產(chǎn)量性狀和品質(zhì)性狀的綜合表現(xiàn)均進(jìn)入前10位的8、9號(hào)單株都出自保梅嶺產(chǎn)區(qū)，而均非出自產(chǎn)量性狀排名第一的霸王嶺產(chǎn)區(qū)，其原因可能是，保梅嶺產(chǎn)區(qū)的油茶實(shí)生老林殘存面積達(dá)60 hm2以上，是我國(guó)面積最大的油茶老林之一，其中的單株數(shù)量大，因而更有可能具有豐富的遺傳多樣性，也更有可能選育出豐產(chǎn)性和優(yōu)質(zhì)性兼?zhèn)涞膬?yōu)良單株。

供試越南油茶豐產(chǎn)單株的取樣地雖然限于海南部分產(chǎn)區(qū)，但對(duì)越南油茶單株和產(chǎn)區(qū)作關(guān)聯(lián)分析時(shí)所涉及的觀測(cè)指標(biāo)均是國(guó)家標(biāo)準(zhǔn)《油茶良種選育技術(shù)規(guī)程》（GB/T 28991—2020）[33]、林業(yè)行業(yè)標(biāo)準(zhǔn)《油茶第1部分：優(yōu)樹選擇和優(yōu)良無(wú)性系選育技術(shù)規(guī)程》（LY/T 1730.1—2008）[34]中列出的關(guān)鍵指標(biāo)，其茶籽油的全部品質(zhì)性狀指標(biāo)基本都達(dá)到了《油茶籽油》（GB /T 11765—2018）[28]中的要求，因此，本研究對(duì)其產(chǎn)量和品質(zhì)等性狀指標(biāo)的等權(quán)關(guān)聯(lián)度的排序基本上是合理的。隨著越南油茶種質(zhì)資源研究的不斷深入，相關(guān)研究所需考察的農(nóng)藝和經(jīng)濟(jì)性狀指標(biāo)將會(huì)越來(lái)越多，這樣就能更加全面地了解不同性狀對(duì)茶籽油產(chǎn)量和品質(zhì)性狀的影響的相對(duì)重要性及其差異程度，針對(duì)不同重要性的性狀給予不同的權(quán)重，然后對(duì)種質(zhì)資源作出更加全面的比較分析，研究結(jié)果也將更加可靠更加準(zhǔn)確，同時(shí)，灰色關(guān)聯(lián)分析法在今后有關(guān)越南油茶品種的綜合評(píng)價(jià)、各因子對(duì)產(chǎn)量的影響度的排名、主導(dǎo)影響因子的確定等方面的研究中也將有望得以推廣應(yīng)用。

4 結(jié) 論

雖然海南產(chǎn)區(qū)不同越南油茶豐產(chǎn)單株的某些農(nóng)藝、產(chǎn)量和茶籽油品質(zhì)等性狀在不同產(chǎn)區(qū)間與不同年份間普遍存在明顯差異，但其茶籽油品質(zhì)性狀基本上都達(dá)到了《油茶籽油》（GB/T 11765—2018）[28]中的要求，因此，豐產(chǎn)性仍是越南油茶的主要育種目標(biāo)。越南油茶實(shí)生選種需要開展連續(xù)多年、多點(diǎn)同時(shí)進(jìn)行的試驗(yàn)研究。研究選取的8和9號(hào)單株均為豐產(chǎn)、優(yōu)質(zhì)單株，可以推薦其進(jìn)入優(yōu)株篩查的復(fù)選階段；霸王嶺產(chǎn)區(qū)、保梅嶺產(chǎn)區(qū)可能均為優(yōu)勢(shì)產(chǎn)區(qū)，即海南省西南部熱帶雨林油茶產(chǎn)區(qū)很可能是越南油茶的優(yōu)勢(shì)產(chǎn)區(qū)，這一產(chǎn)區(qū)蘊(yùn)含著豐產(chǎn)、優(yōu)質(zhì)單株。采用灰色關(guān)聯(lián)分析法對(duì)越南油茶單株和產(chǎn)區(qū)作綜合評(píng)價(jià)研究，在一定程度上能提高油茶研究的效率和準(zhǔn)確性，在今后有關(guān)越南油茶種質(zhì)資源的評(píng)價(jià)研究中可以推廣應(yīng)用。