摘 要：研究以4个不同产地的38份亚麻资源为供试材料，开展9项产量性状指标和6项品质性状指标检测，通过多样性?相关性和主成分分析等研究不同亚麻资源单株产量性状及籽粒品质性状的差异?结果发现：不同亚麻资源15项经济性状的变异幅度为7.94%～47.57%，棕榈酸的多样性指数达到最大；相关性分析得出，株果数和单株粒重呈极显著正相关，相关系数达到最大，为0.887；主成分分析得出5个主成分累计方差贡献率达80.01%，其中可以把单株粒重?株果数?千粒重?单株粒数?株高?工艺长度?硬脂酸?分枝数?棕榈酸和油酸作为38份亚麻资源经济性状的综合评价指标；采用隶属函数法计算出每个亚麻资源的D值，通过D值分析得出经济性状最佳的亚麻资源为陇亚5号?
关键词：亚麻；遗传多样性；产量性状；品质性状；综合评价
亚麻(Linum usitatissimum L.)，别名胡麻，是亚麻科(Linaceae)亚麻属(Linum)一年生草本纤维经济植物，具有悠久的栽培历史，可分为油用亚麻?纤维用亚麻和油纤兼用亚麻[1-2]?其中亚麻籽中富含α-亚麻酸?棕榈酸和膳食性纤维等多种优质营养成分，具有排毒减肥?降压减脂?防癌抗癌等功效，在食用油和保健食品领域中备受消费者青睐[3-5]?目前，油用亚麻资源的经济性状鉴定通常集中在部分有代表性的产量性状和品质性状指标上，产量性状指标包括株高?分支数?千粒重等，品质性状指标包括籽粒含油量和各类脂肪酸组分等?与国外相比，我国拥有的亚麻种质资源数量小?类型不足?研究不够深入，并且亚麻产地较多，各产地间自然条件差异较大，因此深入分析不同地区?不同亚麻种质资源的产量性状和品质性状指标具有重要的意义[5-6]?
不同亚麻资源经济性状鉴定需要将产量性状和品质性状等指标相结合来进行综合评价，因此，简单有效的鉴定步骤和评价方法是筛选优质资源的关键?目前，众多学者从不同角度研究并提出了品种鉴定选育的指标选择和评价方法，针对不同作物?不同生育期?不同地区等种质资源提出了更加全面的分析及评价方法，目前已在大豆[7]?玉米[8]?西瓜[9]?水稻[10]等作物中应用，如主成分分析法?隶属函数法?聚类分析法等，将多指标多方法相结合能更好地制定评价方案，评价结果更加真实可靠?
本研究选取伊犁综合试验站提供的来自不同地区的38份亚麻资源为试验材料，对其单株产量性状和品质性状进行比较分析和综合评价，以期通过研究筛选特异种质，为亚麻种质资源创新和品质育种提供参考?
1材料与方法
1.1试验材料
本研究选取来自河北?甘肃和新疆3个产地的38份亚麻资源，具体资源信息见表1?
表1 38份亚麻资源信息表
1.2测定项目与方法
参考《亚麻种质资源描述规范和数据标准》，在亚麻植株收获前取样，测定各小区亚麻植株株高?工艺长度?分茎数?分枝数?株果数；用天平?种子分析仪等测定亚麻植株单株粒数?单株粒重?每果粒数?千粒重等产量性状指标；亚麻籽粒含油率参照CB/T14488.1-2008中8.3.6进行处理，根据步骤9进行测定；亚麻籽粒各类脂肪酸含量测定参考孟桂元等[11]的方法?
1.3数据统计与分析
隶属函数值按式(1)计算?
式中：Xj—第j个综合指标，j=1，2，3，…，n；Xmin?Xmax—第j个综合指标的最值；U(Xj)—第j个综合指标的隶属函数值?
权重按式(2)计算?
式中：Wj—第j个综合指标的权重；Rj—第j个综合指标的贡献率?
综合评价值按式(3)计算?
式中：D—不同产地亚麻经济性状指标的综合评价值?
采用Excel2017软件进行数据统计整理；Shannon-Weave多样性指数(H’)计算参考关虎等[12]方法；用IBMSPSS Statistics25软件进行相关性分析?聚类分析和主成分分析(PCA)，聚类分析采用组间连接方法，遗传距离为平方欧式距离?
2结果与分析
2.1 38份亚麻资源经济性状多样性分析
参试不同亚麻资源的经济性状测定值统计见下表2?38份亚麻籽粒的亚麻酸含量平均值为44.13%，亚麻酸含量最高的是陇西白胡麻，为54.85%，亚麻酸含量最低的是陇亚5号，为33.82%，变异系数仅为7.95%，说明亚麻酸含量在亚麻资源间比较稳定?从表2可以看出，变异系数最大的是分茎数，变异系数达到47.57%，其平均值为1.53，标准差为0.73；棕榈酸变异系数最小，只有7.94%，其平均值和标准差分别为5.64和0.45?其他13个经济性状变异系数由高到低排序依次为：单株粒数>单株粒重>株果数>分枝数>千粒重>工艺长度>每果粒数>油酸>含油率>硬脂酸>亚油酸>株高>亚麻酸?其中分茎数?单株粒数?单株粒重?株果数?分枝数这5个性状的变异系数均超过了30%?综上所述，供试亚麻资源间经济性状差别很大，具有丰富的遗传背景，为亚麻育种提供优异的种质基础?
研究表明，15个经济性状遗传多样性指数变幅在5.06～5.24，平均为5.21，总体来看，15个经济性状的多样性指数都比较高，其中棕榈酸和亚麻酸的遗传多样性指数最高，均为5.24，分茎数的遗传多样性指数最低，为5.06?多样性指数从大到小排列依次为：棕榈酸>亚麻酸>株高>亚油酸>硬脂酸>含油率>油酸>每果粒数>工艺长度>千粒重>分枝数>株果数>单株粒重>单株粒数>分茎数?
2.2 38份亚麻资源经济性状的相关性分析
由表3可知，15对品质性状组合中有1对性状组合呈极显著负相关，3对性状呈显著正相关，1对性状呈显著负相关(表3)?亚麻酸含量与棕榈酸?亚油酸含量呈显著正相关，相关系数分别为0.336和0.375，与硬脂酸含量呈显著负相关，与油酸含量呈极显著负相关，表明随着亚麻酸含量的升高，棕榈酸和亚油酸含量升高，油酸和硬脂酸含量降低?含油率与棕榈酸?亚油酸和亚麻酸含量呈负相关，与油酸和硬脂酸含量呈正相关?亚油酸含量与油酸?硬脂酸含量呈负相关，与棕榈酸含量呈显著正相关，相关系数为0.383?硬脂酸含量与棕榈酸呈含量负相关，与油酸含量呈正相关，说明随着硬脂酸含量的提高，棕榈酸含量降低，油酸含量增加?油酸含量和棕榈酸含量呈负相关?
表2 38份亚麻单株产量性状和籽粒品质性状的统计结果
表3 38份亚麻籽粒品质性状间相关性分析
注：“∗∗”代表在p<0.01水平极显著相关；“∗”代表在p<0.05水平显著相关?下同
36对产量性状组合中有9对性状组合呈极显著正相关，有1对呈极显著负相关(表4)?株高与每果粒数和千粒重呈负相关，与分茎数?株果数和工艺长度呈(极)显著正相关，相关系数分别为0.411?0.375和0.825?工艺长度与分茎数?分枝数和株果数呈正相关，与单株粒数?单株粒重?每果粒数和千粒重呈负相关，均未达到显著水平?分茎数与株果数?单株粒数和单株粒重呈极显著正相关，相关系数分别为0.772?0.700和0.655，与分枝数呈正相关，与每果粒数和千粒重呈负相关?分枝数与株果数和单株粒重呈极显著正相关，相关系数分别为0.563和0.591，与单株粒数呈正相关，与千粒重呈显著正相关，与每果粒数呈显著负相关?株果数与单株粒数和单株粒重呈极显著正相关，与每果粒数呈负相关，与千粒重呈正相关?单株粒数与单株粒重呈极显著正相关，相关系数为0.744，与每果粒数呈显著正相关，与千粒重呈显著负相关，表明随着单株粒数的增加，单株粒重和每果粒数增加，千粒重降低?单株粒重与千粒重呈正比，与每果粒数呈负相关?每果粒数与千粒重呈极显著负相关，相关系数为-0.730?
从相关性来看，株果数与株高?分茎数?分枝数?单株粒数?单株粒重关系更加密切，除了与每果粒数呈较弱的负相关外，与其他产量性状都呈现正相关，且大部分呈极显著正相关?因此，用株果数可以表征亚麻产量性状，但所测的6个品质性状里，单个品质性状无法衡量亚麻总体品质性状?
表4 38份亚麻单株产量性状间相关性分析
2.3 38份亚麻资源经济性状的主成分分析
2.3.1主成分提取
以38份亚麻资源的15种经济性状指标构成38×15的矩阵，对其进行主成分分析，结果见表5?由表5可知，前5个主成分对应特征值均大于1，其方差贡献率分别为29.181%?21.687%?13.973%?7.938%和7.233%，累计方差贡献率达80.013%，说明可以用前5个主成分代替15个经济性状指标对38份亚麻资源进行综合评价?由表6可知，5个主成分对应的特征向量为：
由表5?6可知，PC1的方差贡献率为29.181%，其中单株粒重和株果数的特征向量值绝对值较大，分别为0.419和0.396，可知PC1主要由株粒重和株果数决定；PC2的方差贡献率为21.687%，其中千粒重和单株粒数的特征向量值绝对值较大，分别为0.412和0.397，可知PC2主要由千粒重和单株粒数决定；PC3的方差贡献率为13.973%，其中株高和工艺长度的特征向量值绝对值较大，分别为0.459和0.577，可知PC3主要由株高和工艺长度决定；PC4的方差贡献率为7.938%，其中硬脂酸和分枝数的特征向量值绝对值较大，分别为0.733和0.325，可知PC4主要由硬脂酸含量和分枝数决定；PC5的方差贡献率为7.233%，其中棕榈酸和油酸的特征向量值绝对值较大，分别为0.619和0.486，可知PC5主要由棕榈酸和油酸含量决定?综上所述，株粒重?株果数?千粒重?单株粒数?株高?工艺长度?硬脂酸?分枝数?棕榈酸和油酸可以作为38份亚麻资源经济性状的综合评价指标?
表5 38份亚麻经济性状指标主成分分析的特征值?贡献率和累积贡献率
表6 38份亚麻经济性状指标性状主成分的特征向量
2.3.2综合评价
(1)隶属函数分析
通过IBMSPSS Statistics 25软件将数据进行标准化处理?根据式(1)可计算出每个亚麻资源综合指标的隶属函数值?5个综合指标的隶属函数值不同，PC1对应的亚麻资源869和陇西白胡麻的隶属函数值分别为1.000和0，说明869在PC1上表现出的品质最优，陇西白胡麻在PC1上表现出的品质最差?同理可得，P.I177451和庆阳胡麻在PC2上表现分别为品质最优和最差；静宁红胡麻和BGOLD在PC3上表现分别为品质最优和最差；陇亚5号和外引6在PC4上表现分别为品质最优和最差；庆阳老和0216-16-8在PC5上表现分别为品质最优和最差?
(2)各综合指标权重的确定
根据综合指标贡献率的大小，第1综合指标贡献率为29.181%，第2综合指标贡献率为21.687%，第3综合指标贡献率为13.973%，第4综合指标贡献率为7.938%，第5综合指标贡献率为7.233%?可根据式(2)求出权重?经计算，5个综合指标权重分别为0.292?0.217?0.140?0.079和0.072(表7)?
(3)经济指标数据的综合评价
根据式(3)计算出D值大小，通过对D值排序，得出38份亚麻资源优劣顺序为：陇亚5号>869>安西红胡麻>P.I177451>861>天水市老胡麻>88(88)-1-10>西礼红胡麻>民勤胡麻>陇西红胡麻>陇亚4号>武威红胡麻>礼县底脚>88(125-1)>灵台五星>P.I1181058>线胡麻>定西17号>陇西白胡麻>灵台转那>定亚18号>伊亚2号>No.547>榆中红胡麻>78-18-1-34-1-2>FIELDNo.17>静宁红胡麻>BGOLD>0216-16-8>LINAGRDSSES>8513>定亚21>DEHISCTNT>庆阳老>BLANC>庆阳胡麻>AC.Emerson>外引6?
2.3.3聚类分析
采用组间连接法对表7中的D值进行聚类分析，建立了聚类树状图(图1)，从图中可以看出，当平方欧式距离为12时，聚类分析将其划分为3大类，其综合评价值依次降低?第I类包括陇亚5号?869?安西红胡麻等16个资源；第II类包括线胡麻?定西17号?陇西白胡麻等15个资源；第Ⅲ类包括定亚21?庆阳老?庆阳胡麻等7个资源?
表7 38份亚麻经济性状指标评价综合指标值?隶属函数值?D值及排序
3讨论
性状指标的合理选择是进行作物综合评价的关键，大量学者从不同角度提出了综合评价方案，涉及的主要性状类别不同，如陈婷婷等[13]研究发现，影响花生经济性状的主要指标包括蛋白和总氨基酸含量等，而徐澜等[14]认为影响小麦经济性状的主要指标包括千粒重和成穗数等?针对不同作物及作物利用方向的不同，所选评价指标也大不一样，完善油用亚麻经济性状综合评价体系在亚麻品种选育上至关重要?采用隶属函数的方法可对不同隶属关系的性状指标进行综合评价，通过提取主成分来提高评价结果准确性，在品种选育上被广大学者认为是一种可靠的评价方法[15]?
图1 聚类树状图
本研究多样性分析显示，38份供试油用亚麻资源的15个经济性状的变异性较大，变异幅度达到7.94%～47.57%，仅棕榈酸?亚油酸?亚麻酸和株高4个指标变异系数小于10，其他均大于10，说明这些亚麻性状变异丰富，具有很大研究潜力?相关性分析显示，绝大多数性状之间都有显著性或者极显著性差异，如油酸和亚麻酸之间存在极显著负相关，这与谢冬微等[16]的研究结果一致，表明这些性状之间互相影响，存在复杂的相关性，在关注目标性状的同时也要兼顾其他性状?
通过主成分分析法，将供试的38份油用亚麻15个主要经济性状指标简化为5个具有代表性的主成分，累计贡献率达80%以上，能更科学地反映供试材料经济性状的综合特征，同时大大减少了数据处理量，数据更加直观和容易处理?根据各个资源的D值大小进行聚类分析，划分等级，可将亚麻经济性状表现较好的归为一类，供后续试验?最终D值排名最高的是陇亚5号，其亚麻酸含量只有33.82%，但变异系数小，说明亚麻酸含量在亚麻资源间比较稳定，根据D值可进一步进行陇亚5号?869和安西红胡麻的单株选择，以期更科学?全面?客观地遴选出优异的亚麻资源?
本研究仅对38份油用亚麻资源经济性状指标进行了分析和评价，但在亚麻育种研究中，必须增加植株抗性?田间适应性等关键指标进行综合评价，才能更好适应生产需要?
4结论
38份亚麻资源变异系数最大的指标是分茎数，为47.57%，最小的指标是棕榈酸，为7.94%，其中分茎数?分枝数?株果数?单株粒数和单株粒重的变异系数均超过30%，说明供试亚麻资源间经济性状差别很大，具有丰富的遗传背景?而15个经济性状指标的遗传多样性指数变化幅度不大，分茎数和棕榈酸的多样性指数为两个最值，分别为最小值5.062和最大值5.243?相关性分析发现，在15个经济性状指标内，株果数和单株粒重呈极显著正相关，相关系数达到最大，为0.887?主成分分析提取出5个主成分，累计方差贡献率达80.013%，其中株粒重?株果数?千粒重?单株粒数?株高?工艺长度?硬脂酸?分枝数?棕榈酸和油酸可以作为38份亚麻资源经济性状的综合评价指标?通过隶属函数法计算出每个亚麻资源的D值，通过D值进行亚麻资源优劣排序为：陇亚5号>869>安西红胡麻>P.I177451>861>天水市老胡麻>88(88)-1-10>西礼红胡麻>民勤胡麻>陇西红胡麻>陇亚4号>武威红胡麻>礼县底脚>88(125-1)>灵台五星>P.I1181058>线胡麻>定西17号>陇西白胡麻>灵台转那>定亚18号>伊亚2号>No.547>榆中红胡麻>78-18-1-34-1-2>FIELDNo.17>静宁红胡麻>BGOLD>0216-16-8>LINAGRDSSES>8513>定亚21>DEHISCTNT>庆阳老>BLANC>庆阳胡麻>AC.Emerson>外引6，得出经济性状最佳的亚麻资源为陇亚5号，最差亚麻资源为外引6，此评价结果与各亚麻资源间产量性状指标数据变化趋势基本一致?将各亚麻资源D值进行聚类分析，把D值前16位分到第Ⅰ类群，D值后7位分到第Ⅲ类群，其他分在第Ⅱ类群?
参考文献
[1]张丽丽，刘晶晶，乔海明，等.从俄罗斯引进亚麻种质资源的农艺性状评价[J].中国油料作物学报，2017，39(5)：698-703.
[2]崔翠，周清元，王利鹃，等.亚麻种质主要农艺性状主成分分析与综合评价[J].西南大学学报(自然科学版)，2016，38(12)：10-18.
[3]邓乾春，禹晓，黄庆德，等.亚麻籽油的营养特性研究进展[J].天然产物研究与开发，2010，22(4)：715-721.
[4]郭永利，范丽娟.亚麻籽的保健功效和药用价值[J].中国麻业科学，2007(3)：147-149.
[5]兴丽，赵凤敏，曹有福，等.不同产地亚麻籽矿物质元素及脂肪酸组成的主成分分析[J].光谱学与光谱分析，2014，34(9)：2538-2543.
[6]李秋芝，姜颖，鲁振家，等.300份亚麻种质资源主要农艺性状的鉴定及评价[J].中国麻业科学，2017，39(4)：172-179.
[7]宋江峰，刘春泉，姜晓青，等.基于主成分与聚类分析的菜用大豆品质综合评价[J].食品科学，2015，36(13)：12-17.
[8]于康珂，刘源，李亚明，等.玉米花期耐高温品种的筛选与综合评价[J].玉米科学，2016，24(2)：62-71.
[9]侯东颖，苏东涛，郝科星，等.基于主成分和聚类分析的无籽西瓜果实性状的综合评价[J].中国瓜菜，2022，35(5)：37-41.
[10]陶维旭，程生海，冀俊超，等.水稻品种资源耐盐性综合评价及耐盐指标筛选[J].江苏农业科学，2022，50(18)：180-187.
[11]孟桂元，孙方，周静，等.亚麻种质脂肪酸成分差异及其相关性研究[J].分子植物育种，2016，14(9)：2502-2508.
[12]关虎，王振华，曹禹，等.亚麻品种主要农艺性状遗传多样性分析[J].新疆农业科学，2011，48(11)：2035-2040.
[13]陈婷婷，王苗苗，黄杨，等.花生种质农艺?产量和品质性状的综合评价[J].花生学报，2020，49(4)：38-46.
[14]徐澜，秦秀娟，安伟，等.基于主成分分析的南引小麦品种综合评价[J].甘肃农业大学学报，2022，57(3)：32-41.
[15]JIN D， XU Y， GUI H， et al. Evaluation of cotton (Gossypium hirsutum L.) leaf abscission sensitivity triggered by thidiazuron through membership function value[J]. Plants (Basel)， 2020， 10(1)： 49.
[16]谢冬微，赵德宝，杨学，等.不同类型亚麻种质资源脂肪酸组分相关及变异特点分析[J].中国麻业科学，2015，37(4)：173-177.
文章摘自：陈宇，守合热提·牙地卡尔，祖勒胡玛尔·乌斯满江等.38份油用亚麻资源经济性状多样性分析与综合评价[J].中国麻业科学，2023，45(02)：49-57.
