陳詩強(qiáng)，陳 亮，蘇素嬌，李 娜，朱 思，牛廣俊，梁一池

（福建中醫(yī)藥大學(xué)，福建 福州 350122）

油茶Camellia oleiferaAbel.屬山茶科山茶屬植物，為多年生喬木或灌木。油茶是我國重要的木本油料植物，其良種選育歷史悠久[1]。隨著油茶科研的深入和推廣，我國油茶產(chǎn)業(yè)發(fā)展迅速，國務(wù)院批準(zhǔn)頒布的《全國油茶產(chǎn)業(yè)發(fā)展規(guī)劃（2009～2020）》中提出，2009年全國油茶林面積突破3×106hm2，2020年全國茶油產(chǎn)量將達(dá)到250萬t，按銷售價(jià)4萬元/t計(jì)算，年產(chǎn)值達(dá)到1 000億元。

典范相關(guān)分析是研究2組多元變量之間相關(guān)關(guān)系的一種統(tǒng)計(jì)方法，能同時(shí)處理2組不同屬性的變量，并能找出這2類屬性間的關(guān)系。筆者根據(jù)典范分析原理研究油茶生物性狀與質(zhì)量性狀2組變量之間的關(guān)系，旨在尋找影響油茶質(zhì)量性狀的主要因素，從整體上考察多個(gè)生物學(xué)性狀與質(zhì)量性狀的關(guān)系，通過對(duì)生物學(xué)性狀的綜合選擇，為油茶栽培與加工方面的研究提供參考。

1 材料與方法

1.1 材 料

20個(gè)無性系油茶優(yōu)良單株采自福建省三明市后樓村油茶園豐產(chǎn)示范基地，經(jīng)福建中醫(yī)藥大學(xué)藥學(xué)院梁一池教授鑒定，為山茶科山茶屬普通油茶Camellia oleiferaAbe1。20個(gè)無性系編號(hào)為：YX1～YX20。測(cè)定其9個(gè)方面的生物學(xué)性狀：果形指數(shù)（高徑比）、單果質(zhì)量、單果籽數(shù)、單果籽質(zhì)量、果皮厚度、鮮出籽率、干出仁率、干仁含油率、單株果質(zhì)量，分別用X1～X9表示；6個(gè)方面的質(zhì)量性狀：單株產(chǎn)油量、棕櫚酸含量、亞油酸含量、油酸含量、硬脂酸含量、角鯊烯含量，分別用Y1～Y6表示。

Agilent 6890N氣相色譜法（美國安捷倫公司），F(xiàn)ID檢測(cè)器。對(duì)照品均購自成都瑞芬生物科技有限公司，棕櫚酸甲酯（批號(hào)：Z-060-120712），油酸甲酯（批號(hào)：Y-172-120622），亞油酸甲酯（批號(hào)：Y-173-120929），硬脂酸甲酯（批號(hào)：Y-171-120926）。

1.2 方 法

1.2.1 供試品的制備

依據(jù)GB/T 5009.6-2003《食品中脂肪的測(cè)定》中的索式抽提法進(jìn)行供試品的制備。取烘干至恒質(zhì)量粉碎的油茶籽仁20 g，石油醚為溶劑，料液比1∶20，趁熱減壓抽濾，過濾，合并濾液，60 ℃減壓濃縮至無石油醚，于100 ℃干燥至恒質(zhì)量[2]。

取油脂0.1 g，置于10 mL容量瓶中，加入正己烷溶解至刻度，取5 mL溶液于15 mL具塞試管中，加入2 mL苯搖勻，加入4%氫氧化鈉甲醇溶液5 mL，60 ℃水浴20 min，靜置分層，取上清液加入無水硫酸鈉，22 μm濾膜過濾，即得供試品。

1.2.2 油脂含量的測(cè)定

取一定量供試品，采用氣相色譜法進(jìn)行測(cè)定。色譜條件為：DB-WAX 石英玻璃毛細(xì)管柱（30 m×0.25 mm×0.5 μm）；載氣N2（純度99.999%），空氣400 mL·min-1，氫氣40 mL·min-1；柱流量為0.77 mL·min-1，分流比1∶50，進(jìn)樣量1 μL；進(jìn)樣口溫度240 ℃，檢測(cè)器溫度230 ℃；柱溫采用等度升溫190 ℃，維持20 min，進(jìn)樣口壓力100 kPa。所得結(jié)果采取面積歸一化法進(jìn)行計(jì)算。

1.2.3 角鯊烯含量的測(cè)定[3]

（1）對(duì)照品溶液的制備 精密稱量20 mg角鯊烯標(biāo)準(zhǔn)品放入10 mL容量瓶，加正己烷溶液至刻度值，配成質(zhì)量濃度為2 mg·mL-1的標(biāo)準(zhǔn)溶液。

（2）樣品含量的測(cè)定 取一定量供試品，采用氣相色譜法進(jìn)行測(cè)定。色譜條件：HP-5 石英玻璃毛細(xì)管柱（30 m×0.25 mm×0.5 μm）；載氣N2（純度99.999%），空氣 400 mL·min-1，H240 mL·min-1；柱流量為1.34 mL·min-1，分流比為1∶10，進(jìn)樣量1 μL；進(jìn)樣口溫度270 ℃，檢測(cè)器溫度290 ℃；柱溫采用程序升溫，初溫150 ℃，維持1 min，以15 ℃·min-1升溫至270 ℃，維持10 min，進(jìn)樣口壓力204 kPa，時(shí)長19 min。

（3）標(biāo)準(zhǔn)曲線的繪制 分別取出一定量2 mg·mL-1的溶液配成質(zhì)量濃度分別為0.2、0.4、0.6、0.8、1.0 mg·mL-1的角鯊烯標(biāo)準(zhǔn)溶液，按樣品含量測(cè)定方法分別進(jìn)行角鯊烯含量的測(cè)定，并繪制標(biāo)準(zhǔn)曲線。

2 結(jié)果與分析

20個(gè)無性系油茶優(yōu)良單株的生物學(xué)性狀及質(zhì)量性狀指標(biāo)測(cè)定結(jié)果見表1，利用SPSS18.0軟件對(duì)這些生物學(xué)性狀和質(zhì)量性狀指標(biāo)進(jìn)行相關(guān)分析[4]，得到相關(guān)系數(shù)矩陣，如表2所示。由表2可以看出，15個(gè)成分指標(biāo)間共有105對(duì)組合，其中9對(duì)組合極顯著相關(guān)，10對(duì)組合顯著相關(guān)。生物學(xué)性狀和質(zhì)量性狀指標(biāo)間有4對(duì)極顯著相關(guān)，1對(duì)顯著相關(guān)。單果質(zhì)量和單果籽質(zhì)量、果皮厚度，干出仁率和干仁含油率、單株產(chǎn)油量、油酸含量，干仁含油率和油酸含量、角鯊烯含量，亞油酸含量和油酸含量，棕櫚酸含量和亞油酸含量呈極顯著相關(guān)；單果籽數(shù)和單果籽質(zhì)量、油酸含量，單果籽質(zhì)量和鮮出籽率，果皮厚度和鮮出籽率，單株果質(zhì)量和亞油酸含量，單株產(chǎn)油量和油酸含量、角鯊烯含量，棕櫚酸含量和油酸含量、角鯊烯含量，油酸含量和硬脂酸含量、角鯊烯含量呈顯著相關(guān)。

利用DPS軟件對(duì)20個(gè)無性系油茶優(yōu)良單株的生物學(xué)性狀和質(zhì)量性狀進(jìn)行典范相關(guān)分析[5]，典范相關(guān)系數(shù)的檢驗(yàn)結(jié)果見表3，典范變量分析結(jié)果見表4。由表3可以看出，第1對(duì)典范相關(guān)系數(shù)最大為0.949 4，Wilk's統(tǒng)計(jì)量為0.007 3，df＝54，p＝0.047 9，具有統(tǒng)計(jì)學(xué)意義，其它典范相關(guān)系數(shù)均無統(tǒng)計(jì)學(xué)意義。據(jù)此可以認(rèn)為第1對(duì)典范變量U1和V1已經(jīng)把2組變量X和Y之間存在的全部相關(guān)信息進(jìn)行分解并提取出來。

表1 油茶優(yōu)良單株的生物學(xué)性狀及質(zhì)量性狀Table 1 The biological and quality characters of the fine trees in Camellia oleifera

表2 生物學(xué)性狀與質(zhì)量性狀之間的相關(guān)系數(shù)矩陣?Table 2 The correlation coefficient matrix between the biological and quality characters

表3 典范相關(guān)系數(shù)的檢驗(yàn)Table 3 Test of canonical correlation coefficient

由表4可以看出，在生物學(xué)性狀線性組合中，果形指數(shù)、單果籽質(zhì)量、干出仁率、干仁含油率、單株果質(zhì)量呈正效應(yīng)；單果質(zhì)量、單果籽數(shù)、果皮厚度、鮮出籽率[6]呈負(fù)效應(yīng)。根據(jù)其相關(guān)系數(shù)得出因子權(quán)重，由高到低依次為：單果籽質(zhì)量、單果質(zhì)量、單株果質(zhì)量、干出仁率、鮮出籽率、干仁含油率、果形指數(shù)、果皮厚度。與之對(duì)應(yīng)的質(zhì)量性狀線性組合中，油酸含量、硬脂酸含量呈負(fù)相關(guān)；單株產(chǎn)油量、棕櫚酸含量、亞油酸含量、角鯊烯含量呈正相關(guān)。根據(jù)其相關(guān)系數(shù)得出因子權(quán)重由高到低依次為：油酸含量、亞油酸含量、單株產(chǎn)油量、角鯊烯含量、硬脂酸含量、棕櫚酸含量[7]。

表4 生物學(xué)性狀與質(zhì)量性狀的典范變量分析Table 4 Analysis of canonical correlation between the biological and quality characters

由于X與Y之間有極顯著正相關(guān)關(guān)系，從油茶當(dāng)前產(chǎn)油量及其品質(zhì)需要方面考慮，應(yīng)以單果籽質(zhì)量、單果質(zhì)量、單株果質(zhì)量為主要選擇因子，同時(shí)由于單果籽質(zhì)量與單果質(zhì)量呈極顯著正相關(guān)（見表2），因此應(yīng)主要考慮單株果質(zhì)量與單果質(zhì)量，能獲得較好的效果。綜合以上分析，可以構(gòu)造1個(gè)綜合選擇指數(shù)方程：

式中為標(biāo)準(zhǔn)化基因型值，I表示綜合選擇指數(shù)，X1～X9分別表示果形指數(shù)（高徑比）、單果質(zhì)量、單果籽數(shù)、單果籽質(zhì)量、果皮厚度、鮮出籽率、干出仁率、干仁含油率、單株果質(zhì)量。由該公式可算出20個(gè)無性系的油茶優(yōu)良單株綜合指數(shù)[8]。用聚類分析方法，根據(jù)各優(yōu)樹得分對(duì)其進(jìn)行分類，將20個(gè)無性系分為3類，結(jié)果見表5。

表5 油茶樹綜合指數(shù)聚類Table 5 The cluster of the fine trees in Camellia oleifera based on comprehensive index

3 結(jié)論與討論

油茶生物學(xué)性狀與質(zhì)量性狀的關(guān)系復(fù)雜，運(yùn)用典范相關(guān)方法分析2組性狀之間的關(guān)系，可以使它們的關(guān)系簡單化。典范相關(guān)分析結(jié)果有5對(duì)變量，彼此獨(dú)立，第1對(duì)變量達(dá)極顯著水平，幾乎占了全部信息量的90.43%。第1對(duì)典范變量中的主導(dǎo)因子是單株果質(zhì)量和單果質(zhì)量，它們與質(zhì)量性狀因子油酸含量、亞油酸含量和單株產(chǎn)油量之間關(guān)系緊密，說明在油茶產(chǎn)量和品質(zhì)選擇中要優(yōu)先考慮果質(zhì)量[8]。

根據(jù)典范變量系數(shù)所構(gòu)造的綜合選擇指數(shù)方程I＝－0.112 0X1＋0.611 2X2＋0.309 9X3－0.952 9X4＋0.049 7X5＋0.338 2X6－0.456 1X7－0.314 0X8－0.495 2X9，得出20株優(yōu)樹的生物學(xué)綜合得分值，依據(jù)綜合得分值，進(jìn)行聚類分析，將它們分成3類。其中，I類的生物學(xué)性狀與質(zhì)量性狀均很好，是選育所需要的優(yōu)良植株，II類性狀較好，III類最差[9]。

參考文獻(xiàn)：

[1] 張 可,錢 和,張 添.油茶籽的綜合開發(fā)利用[J].食品科技,2003,3(4):85－86．

[2] 陳毓靜,王琳琳,陳小鵬,等.廣西河池、百色和梧州油茶種子含油率和脂肪酸組成測(cè)定[J].食品科學(xué),2011,5(8):172－176.

[3] 張國武,闕龍善,賴筍芽,等.6個(gè)油茶優(yōu)良無性系種子脂肪酸組分與含量分析[J].江西科學(xué),2007,6(1):35－57.

[4] 梁新華,鄭彩霞,張風(fēng)俠,等．甘草角鯊烯的提取及高效液相色譜分析[J].北京林業(yè)大學(xué)學(xué)報(bào),2010,8(2):123.

[5] 王文浩,蘇淑釵,陳 鳳.遼榛7號(hào)結(jié)果母枝各性狀因子間的相關(guān)分析[J].經(jīng)濟(jì)林研究,2013,31(2):119－153.

[6] 陳 鳳,蘇淑釵,白 倩.平歐雜交榛枝條主要性狀因子間的相關(guān)分析[J].經(jīng)濟(jì)林研究,2013,31(2):21－25.

[7] 王 瑞,陳永忠,陳隆升,等.油茶葉片SPAD值與葉綠素含量的相關(guān)分析 [J].中南林業(yè)科技大學(xué),2013,33(2):76－80.

[8] 曹永慶,姚小華,任華東,等.油茶果實(shí)礦質(zhì)元素含量和油脂積累的相關(guān)性 [J].中南林業(yè)科技大學(xué),2013,33(10):38－41.

[9] 吳志莊,梁一池,劉淑英.增廣隨機(jī)區(qū)組設(shè)計(jì)在肉桂苗期測(cè)定中的應(yīng)用[J].經(jīng)濟(jì)林研究,1999,27(2):35－37.