基于种实性状的无患子优良单株选择1)

刁松锋 邵文豪 姜景民 董汝湘 肖可青

(中国林业科学研究院亚热带林业研究所，富阳，311400) (天台县满园春农林开发有限公司)

基于种实性状的无患子优良单株选择1)

刁松锋 邵文豪 姜景民 董汝湘 肖可青

(中国林业科学研究院亚热带林业研究所，富阳，311400) (天台县满园春农林开发有限公司)

应用主成分分析法，从样本相关矩阵出发，以102株8年生无患子(SapindusmukurossiGaertn.)初选优株的12个主要种实性状(果实质量、果指数、果体积、果皮厚、果皮质量、种子质量、种指数、种体积、种仁质量、单株产量、果皮皂苷质量分数、种仁油脂质量分数)进行分析，以性状累积方差贡献率大于80%为标准，确定了4个反映无患子种实主要经济性状的主成分(特征值>1)，即产量因子、油脂因子、产量皂苷因子和种实形态因子，其相应贡献率分别为47.271%、15.285%、10.206%和8.711%。通过对样本重要主成分值的比较分析，按5%入选率分别筛选出不同利用方向的6种类型的优良单株，即高产型、高油脂型、高皂苷型、高产高油脂型、高产高皂苷型和复合型。不同类型优良单株主成分得分值的现实增益为165%～311%，这说明利用主成分分析法选出的不同类型优良单株具有显著的优异性，可在生产中推广应用。

无患子；种实性状；主成分分析；优树选择

Sapindusmukorossi; Fruit and seed traits; Principal component analysis; Superior individual selection

无患子(SapindusmukorossiGaertn.)，又名肥皂树、洗手果等，为无患子科(Sapindaceae)无患子属落叶大乔木，具有广泛的立地适应性，主要生长在我国东部、南部、西南部，树形美观、生长迅速、材质优良，是我国南方城市主要的园林绿化树种[1-2]。无患子果皮富含皂苷达10.76%[3]，其表面活性成分主要为萜类皂苷和倍半萜糖苷，具有良好的起泡性和去污能力，是天然的环境友好型洗涤剂[4]。无患子种仁油脂含量达42.7%，油脂不饱和脂肪酸含量高达86.63%，为重要的生物质能源树种[5-6]。

随着国家对过度使用化工洗涤剂引起水土污染治理力度的加强，以及对生物质能源需求的增加。可以预见，依托无患子皂素和油脂资源发展无患子生物质产业，必将具有良好的生态效益和巨大的市场前景。我国大部分无患子资源处于野生状态，对其开发利用尚处于起步阶段，目前对无患子的研究多集中在园林绿化[7]、组织培养和体细胞发生[8-9]、化学成分及其提取工艺[5，10]和医药应用等方面[11-12]，对其遗传多样性[13]和遗传图谱构建[14]研究也略有涉及，但其作为传统的洗涤用品原料和新型生物柴油原料，以往并无定向培育的意识，完全是对无患子野生资源的利用，对无患子以种实为利用目的的优良单株选择研究还未见报道。

主成分分析是从多个存在一定相关关系的变量中选出几个新的综合变量，而新的综合变量相互独立又能反映原来多个变量所提供的主要信息，从而简化数据结构，寻找变量间线性关系[15]。所以用主成分值作为优种选择指标，可较准确的了解各性状的综合表现。因此，主成分分析法已成功应用在多种农产品资源评价[16-17]、品种选育[18]和优良单株选择[19]。本文以果用无患子优良单株为选择目标，通过对初选优树性状变异及其相关性进行研究，应用主成分分析法对其进行综合评价，初步筛选出不同利用方向的无患子优良单株，为无患子果用林定向培育提供了基础材料。

1 试验地概况

试验材料来自浙江省天台县满园春农林开发有限公司无患子原料林基地，基地面积约100 hm2，林分为2003年利用当地无患子老树种子繁育营建，株行距多为4 m×4 m。海拔75～110 m；年平均气温16.8 ℃，降水量1 320 mm，年无霜期平均232 d；基地土壤以黏性红黄壤为主，土层厚度20～70 cm，pH值5.3～6.0。林分郁闭度较大，平均0.75，林下部分区块生长一些耐荫植物，主要有枸骨(Ilexcornuta)、小果蔷薇(Rosacymosa)、络石(Trachelospermumjasminoides)、山莓(Rubuscorchorifolius)、铁芒萁(Dicranopterislinearis)等。

2 材料与方法

2.1 性状测量及方法

在无患子果实成熟期，2010—2012年连续3 a对候选优树进行跟踪观测，选出单株产量高、果实大、果枝粗、生长势好、无病虫害等综合表现优良的102棵为初选单株。选择果实质量、果指数、果体积、果皮厚、果皮质量、种子质量、种指数、种体积、种仁质量、单株产量、果皮皂苷质量分数、种仁油脂质量分数等12个种实性状对无患子优树选择观测分析，同时测量树高、胸径、冠幅和枝下高等4个生长性状。

单株产量的估测方法：在南北两面树冠的上、中、下3个位置各选3个能代表结果平均状况的样枝，统计每个样枝果实质量，计算3个样枝的平均质量，单株产量=结果枝数量×样枝平均结果量。果实成熟后，在树冠南面的中上部采样，分株采收、存放，约50粒·株-1。果实自然风干后，采用GB/T 3543.7—1995的百粒四分法随机抽取30粒果实，用电子游标卡尺测量果实和种子纵径、横径和侧径，测量精度0.02 mm；计算果实大小指数和种子大小指数，指数=纵径/横径；用1/1 000电子天平称量30粒果实质量，称量精度为0.01 g；剥去果皮和种皮称果皮质量、种子质量和种仁质量，其平均值分别为果实质量、果皮质量、种子质量和种仁质量；用游标卡尺测量果皮随机部位的厚度，6次平均值为果皮厚。

无患子皂苷测定采用高效液相色谱分析法[20]。将30个果皮混合均匀后用微型粉碎机将其粉碎成粉末状，再放入电热恒温鼓风干燥箱中70 ℃保持干燥。称取0.3 g果皮粉末，装入20 mL的带密封塞子的试管中，将甲醇用移液管定容至15 mL，立即盖上试管塞，浸泡20 h后，用漏斗过滤，并对滤液进行HPLC检测。HPLC检测以0.4 g/L常春藤皂苷元甲醇溶液为对照[21]。色谱条件：色谱柱为SymmetryTMC18(3.9 mm×150 mm)；柱温为40 ℃；流速为1 mL/min；30 min内流动相体积变化为V(CH3CN)∶V(H2O)=0.1∶0.9～V(CH3CN)∶V(H2O)=0.8∶0.2；检测波长为210 nm。

果皮皂苷质量分数=(纯皂苷量/果皮干质量)×100%。

无患子油脂采用索氏提取法[5]。将30个种仁混合均匀后用微型粉碎机粉碎，在电热鼓风干燥箱中80 ℃烘干至恒质量，称取10 g种仁粉末，置于索氏提取器中，加正己烷200 mL，90 ℃回流提取9h，至回流液无色(取回流液滴于滤纸上检查，待溶剂挥发尽，滤纸上不留油迹为提取终点)。用旋转蒸发仪回收正己烷，经干燥后称质量得到黄色清亮无患子籽油。

种仁油脂质量分数=((提取前种仁质量-提取后残渣质量)/提取前种仁质量)×100%。

2.2 数据统计和分析

应用Excel 2010和SPSS18.0分析软件，从样本相关矩阵出发，对12个种实性状指标进行主成分分析，根据性状累积方差贡献率达到80%以上，且特征值大于1确定主成分的个数[15]。由各性状相关矩阵的特征向量，依据各主成分的向量模型，计算重要主成分值，选出优良单株。

3 结果与分析

3.1 无患子种实性状的变异

无患子12个种实性状中，单株产量是每株的果实总质量，果皮皂苷和种仁油脂质量分数的性状值是每株30个果实混合后测定的值，对另外9个表型性状进行方差分析：果实质量(F=60.364)、果指数(F=2.778)、果体积(F=61.445)、果皮厚(F=23.115)、果皮质量(F=48.726)、种子质量(F=15.892)、种指数(F=2.635)、种体积(F=52.280)、种仁质量(F=10.457)。各性状在102棵初选单株间差异均极显著(p<0.01)，其中果实性状的变异程度大于种子性状的变异程度，果皮质量的变异程度大于果皮厚度的变异程度。

变异系数的大小可间接反映出群体的表型多样性丰富程度，变异系数大，说明该群体表型多样性丰富。无患子初选优株的12个种实性状存在着不同的变异(表1)，范围为0.05～0.66，由大到小依次为：单株产量、果皮质量、种仁质量、果皮皂苷质量分数、果皮厚、果实质量、果体积、种体积、种子质量、种仁油脂质量分数、果指数和种指数。其中，单株产量、果实质量、果皮质量、种仁质量、果皮皂苷质量分数等主要经济性状的变异系数均相对较大，分别为0.66、0.22、0.29、0.28和0.24。表征果实和种子形态果指数和种指数变异系数最小，均为0.05，说明无患子初选优树的种实形态变异较小，相对稳定。种仁油脂质量分数作为重要的经济性状其变异系数并不大仅为0.08，但其变异幅度却较大，为26.50%～42.20%，可见在不同单株间选择高油脂优树仍具有可行性。综上可知，无患子初选优株种实的主要经济性状变异丰富，主要经济性状具有较高的选择空间。

表1 无患子初选优株种实性状变异

果皮质量变幅/g均值±标准差/g极差/g变异系数/%种子质量变幅/g均值±标准差/g极差/g变异系数/%种体积变幅/mm3均值±标准差/mm3极差/mm3变异系数/%1．32～6．403．15±0．095．08291．17～2．781．88±0．031．61161616．15～4432．063040．01±55．122815．9118

种仁质量变幅/g均值±标准差/g极差/g变异系数/%单株产量变幅/kg均值±标准差/kg极差/kg变异系数/%果皮皂苷质量分数/%变幅均值±标准差极差变异系数0．32～1．801．09±0．031．48281．50～50．0011．05±0．7248．50661．00～10．275．33±0．139．2724

果皮厚变幅/mm均值±标准差/mm极差/mm变异系数/%种指数变幅均值±标准差极差变异系数/%种仁油脂质量分数/%变幅均值±标准差极差变异系数/%1．23～4．212．29±0．052．98230．89～1．100．99±00．21526．50～42．2035．03±0．2815．708

3.2 无患子种实性状相关性

由表2可知，果实质量与果体积、果皮厚、果皮质量、种子质量、种体积等性状间具有极显著的正相关关系(p<0.01)，其中与果皮质量最为密切(r=0.973)，与果体积次之(r=0.938)，而与种仁质量具有显著正相关(p<0.05)，说明在以上几种性状中果皮质量和果体积对果实质量的贡献率最大，种仁质量最小。果实质量与果指数、种指数及单株产量具有极显著的负相关关系(p<0.01)，表明果实和种子越趋向于圆球体果实质量越小，调查中亦发现，一般果实相对扁圆的类型其质量较大。单株产量除了与果实质量具有极显著的负相关外，还与果体积、果皮质量呈现极显著的负相关(p<0.01)，单株产量高的植株其单果质量反而越小是由于果实较大的个体，其单果序结果数一般较少，产量也相对较低。果皮皂苷质量分数与种指数(r=-0.305)和种仁油脂质量分数(r=-0.303)具有极显著的负相关关系(p<0.01)；种仁油脂质量分数除了与果皮皂苷质量分数具有极显著的负相关关系，还和种仁质量(r=0.413)呈极显著的正相关关系(p<0.01)。果皮皂苷质量分数、种仁油脂质量分数与其它多数性状间相关性不显著性，表明两者性状在不同初选优株间变异相对较小。

表2 无患子初选优株种实性状相关分析

注：** 表示极显著(p<0.01)；*表示显著(p<0.05)。

3.3 无患子种实性状的主成分分析

特征根和贡献率是选择主成分的依据。由样本相关矩阵出发，对原始数据经标准化处理，计算性状相关矩阵的特征根和特征向量，并根据性状特征根和原变量数，计算各个主成分的贡献率和方差贡献率。由于无患子9个种实表型性状在单株间都具极显著差异，为此结合单株产量、果皮含皂率和种仁油脂质量分数等3个主要经济性状做主成分分析(表3)。由表3可知，第1主成分的特征值为5.727，方差贡献率47.727%，是最重要的主因子，且果实质量、果体积、果皮质量和单株产量在第1主成分中分别表现为显著性(p<0.05)和极显著性(p<0.01)，则把第1主成分定义为构成单株和种实质量的因子，即产量因子；第2主成分的特征值为1.834，方差贡献率15.285%，其中种子质量、种体积在第2主成分的相关系数为显著(p<0.05)，种仁质量和种仁油脂质量分数则为极显著性(p<0.01)，则把第2主成分定义为构成油脂质量分数高低的因子，即油脂因子；第3主成分的特征根为1.225，方差贡献率为10.206%，其中单株产量和果皮皂苷质量分数在第3主成分的相关系数为极显著(p<0.01)，则称其为产量皂苷因子；第4主成分的特征根为1.045，方差贡献率为8.711%，其中果指数和种指数在第3主成分的相关系数极显著，而种仁油脂质量分数为显著，因种仁油脂质量分数在第2主成分表现为极显著，则称第4主成分为种实形态因子。其他因子特征值均小于1。可见，前4个主成分的累积贡献率为81.926%，能较好地反映无患子果实大小、果皮皂苷质量分数及种仁油脂质量分数等主要经济性状特征。

表3 无患子种实性状主成分向量和具有显著意义的因子

注：** 表示极显著(p<0.01);*表示显著(p<0.05)。

3.4 不同目标的优株选择

主成分是原性状的线性组合函数，通过计算样本相关矩阵的特征向量得出主成分函数式，即可计算每棵初选优株的主成分得分值(Yi)和综合得分(Y)。

Yi=∑Fi×Xj，其中i=1～4，j=1～12，F为主成分向量值，X各主成分对应的性状。

Y=∑Yi×yi，其中i=1～4，Yi为主成分的分值，yi为对应的贡献率。

ΔG=[(M-m)/m]×100%，其中ΔG为入选优良单株的现实增益，M为入选优良单株的主成分平均值，m为所有单株主成分平均值。

依据单株产量、果皮皂苷质量分数和种仁油脂质量分数等主要经济性状筛选不同的优良单株进行良种培育，将为无患子产业发展提供基础保障。根据不同的选育目标，本文设计高产型、高油脂型、高皂苷型、高产高油脂型、高产高皂苷型以及复合型共6种优树类型，其中复合型是指单株产量、果皮皂苷质量分数、种仁油脂质量分数等重要经济性状综合表现优良的个体类型。

对每种类型单株选取前5%为优良单株，即得分前6位的单株为优良单株。从表4可以看出，入选不同类型的优良单株的现实增益均比较大，其中最高的为高产型优良单株，ΔG=311%；高油脂型优良单株的增益最低，但ΔG仍高达165%。这说明按照5%入选的不同类优良单株具有显著的优异性，可在生产中推广应用。

表4 不同类型优良单株及其得分值和现实增益

注：ΔG为入选优良单株的现实增益，M为入选优良单株的主成分平均值，m为所有单株主成分平均值。

4 结论与讨论

表型多样性是基因型和环境互作的结果，表型变异越大，可能存在的遗传变异越大，表型多样性研究是进行性状选择的重要基础[21]。植物表型多样性是遗传多样性的一个重要组成部分，植物表型性状包括种实、花、叶及树体等，而对于经济林产果类植物来说主要体现为种实性状[22-23]。

本研究在同一生态环境下，根据102棵无患子单株连续3 a的种实性状调查数据进行研究，结果表明，无患子不同优株间种实性状变异丰富，主要经济性状具有较高的选择空间。对无患子12个种实性状分析表明，果实质量、果皮厚、果皮质量、种子质量等变异程度较大，而果指数、种指数和种体积等变异较小，这说明无患子的质量性状变异较大而形态性状相对稳定，这可能与形状间发育程度不同有关[24]；在重要的经济性状方面，单株产量变异幅度最大，其变异系数高达0.66，而次生代谢物中果皮皂苷质量分数和种仁油脂质量分数的变异幅度差异较大，其变异系数分别为0.24和0.08，这可能和次级代谢物产生的不同部位有关。无患子种实性状丰富的变异为育种工作者选择不同性状优良单株提供丰富空间，为其资源良种化和产业化提供了广阔的前景。

无论是以果皮皂苷质量分数为经济目的，还是以种仁油脂质量分数为经济目的，均需建立在单株高产的基础之上。因此，在优良单株选择时，应优先考虑单株产量。相关性分析表明，果实质量和果体积与单株产量均为极显著的负相关关系(p<0.01)，相关系数分别为r=-0.336和r=-0.335，说明单果质量越大、单果体积越大，其单株产量越低。而果指数和种指数与单株产量均为极显著的正相关(p<0.01)，其相关系数分别为r=0.221和r=0.217，这表明果实和种子越趋向于圆球形，单株产量越高。因此，在高产优良单株选择中应以果实小而圆为主要选择性状。同理，由于种指数与皂苷质量分数呈极显著的负相关关系，在高皂苷优株选择时应特别考虑种子形状为椭圆形的单株；种仁质量与油脂质量分数呈极显著的正相关关系，皂苷含量与油脂质量分数则为极显著的负向相关系，所以在选择高油脂优株时应选择种仁质量大而皂苷质量分数低的单株。

经济效益最大化是经济树种选择的目标，而影响单株经济效益的性状众多，且性状之间存在着相关性。如果仅凭少数几个性状的表现型对种质资源进行评价、选择亲本，势必带有主观性。而主成分分析法中的各主成分是一个相对独立的指标体系，它们之间不存在相关，并且数值直观，容易分析[15]。本文利用主成分分析法将无患子12个种实主要经济性状转化为4个主成分，4个主成分提供了原性状的81.926%的信息，并具有明确的生物学意义。

本文按照5%的入选率筛选出高产、高油脂、高皂苷、高产高油脂、高产高皂苷和复合型等6种不同利用类型的无患子优良单株各6株，分别为99、88、81、4、64、48；1、23、91、69、58、2；68、22、78、27、9、76；99、88、81、50、62、91；99、88、81、62、4、27和99、88、81、62、91、84。以99号单株为例，其不仅仅在几种主成分得分值位居前列，而且在单个主要经济性状中的表现也具有显著的优异性，如102株的单果质量、单株产量、果皮皂苷质量分数和种仁油脂质量分数等几种性状的均值分别为：5.02 g、4.82 kg、5.33%和35.03%，而99号单株的几种性状相对应的值为：7.98 g、24.50 kg、9.62%和41.75%，现实增益分别59%、122%、81%和20%，由于102棵初选优株种仁油脂质量分数的变异系数仅为0.08，因此99号单株种仁油脂质量分数20%的现实增益仍然具有显著的优异性。这表明，本文利用主成分分析法选出的不同利用目标的优良单株可以作为扩繁的基础材料。

本文仅研究了浙江省天台县无患子实生林初选优株的种实部分性状，未能对不同种源无患子种实性状分析，也未涉及初选优株的树体生长因子，如树高、胸径、冠幅和枝下高等。这将是今后研究的重点，以便从中选出各性状指标均优良的单株，为选育出适应不同地区生长的新品种奠定基础。

[1] 《中国植物志》编辑委员会.中国植物志：47卷[M].北京:科学出版社,1998.

[2] 《中国树木志》编辑委员会.中国树木志：6卷[M].北京:中国林业出版社,2004.

[3] 邵文豪,姜景民,董汝湘,等.不同产地无患子果皮含皂率的地理变异研究[J].植物研究,2012,32(5):627-631.

[4] Rupesh Kumar, Surek Ha, Balu A, et al. Study of functional properties of Sapindus mukorossi as a potential bio-surfactant[J]. Indian Journal of Science and Technology,2011,4(5):530-533.

[5] 黄素梅,王敬文,杜孟浩,等.无患子籽油脂肪酸成分分析[J].中国油脂,2009,34(12):74-76.

[6] Sun Shangde, Ke Xiaoqiao, Yang Longlong, et al. Enzymatic epoxidation ofSapindusmukorossiseed oil by perstearic acid optimized using response surface methodology[J]. Industrial Crops and products,2011,33(3):676-682.

[7] 黄顺，潘文明.优良造景树无患子[J].园林,2004(9):42.

[8] Kim Hyuntae, Yang Byeunghoon, Park Young Goo, et al. Somatic embryogenesis in leaf tissue culture of Soapberry (SapindusmukorossiGaertn)[J]. Plant Biotechnology,2012,29(3):311-314.

[9] Upma Dobhal, Snehlata Bhandari, Shivani Bisht, et al. Highly efficient plant regeneration through somatic embryogenesis inSapindusmukorossi[J]. Plant Archives,2009,2(6):219-221.

[10] 刘淑明,孙丙寅,贺安乾,等.西部地区文冠果种群种子特征及主要化学成分的地理变化[J].林业科学,2012,48(4):43-48.

[11] Kanchana T, Kasai R, Yamasaki K, et al. Acetylated triterpene saponins from the thai medicinal plant, sapindus emarginatus[J]. Chemical and Pharmaceutical Bulletin,2001,49(9):1195-1197.

[12] 郭英,谢建平,柳爱华,等.无患子药理作用概述[J].中国病原生物学杂志,2011,6(1):873-874.

[13] Pooja Asthana, Jaiswal V S, Jaiswal Uma. Micropropagation ofSapindustrifoliatusL and assessment of genetic fidelity of micropropagated plants using RAPD analysis[J]. Acta Physiol Plant,2011,33:1821-1829.

[14] Mahar K S, Rana T S, Ranade S A, et al. Estimation of genetic variability and population structure inSapindustrifoliatusL, using DNA fingerprinting methods[J]. Trees,2013,27(1):85-96.

[15] 裴鑫德.多元统计分析及其应用[M].北京:中国农业大学出版社,1991：196-247.

[16] 苗锦山,刘彩霞,戴振建.葱种质资源数量性状的聚类分析、相关性和主成分分析[J].中国农业大学学报,2010,15(3):41-49.

[17] 李清华,黄金堂,陈海玲,等.27份花生种质资源的主成分分析及遗传距离测定[J].植物遗传资源学报,2011,12(4):519-524.

[18] 郭宝林,杨俊霞,李永慈,等.主成分分析法在仁用杏品种主要经济性状选中上的应用[J].林业科学,2000,36(6):53-56.

[19] 敖妍.因子分析法在文冠果优良单株选择中的应用[J].华南农业大学学报：自然科学版,2009,30(4):70-73.

[20] 刘志,阮长春,刘天志,等.HPLC法同时测定林下参、鲜人参、生晒参和红参中14种人参皂苷[J].中草药,2012,43(12):2431-2434.

[21] Andrew Lowe, Stephen H, Paul Ashton. Ecological genetics: design, analysis and application[M]. London: Blackwell Publishing,1981.

[22] 黄勇.小果油茶遗传多样性分析及杂交渐渗研究[D].北京：中国林业科学研究院,2009.

[23] Wheeler N C, Guries R P. Population structure, genic diversity, and morphological variation inPinuscontortaDougl[J]. Canadian Journal of Forest Research,1982，12:595-606.

[24] 康永祥,赵宝鑫,玉洁,等.毛梾天然群体种实表型多样性研究[J].西北农林科技大学学报：自然科学版,2011,39(9):107-117.

刁松锋，男， 1989年2月生，中国林业科学研究院亚热带林业研究所，硕士研究生。

姜景民，中国林业科学研究院亚热带林业研究所，研究员。E-mail:jmjiang6001@126.com。

2013年6月20日。

S722.1+3； S722.5

Superior Individual Selection ofSapindusmukorossiBased on Fruit and Seed Traits/Diao Songfeng, Shao Wenhao, Jiang Jingmin, Dong Ruxiang(Research Institute of Subtropical Forestry, Chinese Academy of Forestry, Fuyang 311400, P. R. China); Xiao Keqing(Tiantai Manyuanchun Agroforestry Development Company)//Journal of Northeast Forestry University.-2014,42(4).-6～10，45

1) 浙江省重大科技专项重点农业项目(2011C12015)；林业公益性行业科研专项项目(200804032)。

责任编辑：潘 华。

Based on a specimen correlation matrix, with 102 eight-year-oldSapindusmukorossiGaertn. individuals, the experiment was conducted to determine 12 major economic characteristics by the principal component analysis, including the weight per fruit, fruit finger, the volume of fruit, the thick of pericarp, the weight of pericarp, the weight per seed, seed finger, the volume of seed, the weight of kernel, the weight per kernel, the yield per plant, the saponin contents in the peep and the oil contents, in the kernel. Four factors (eigenvalues>1) were extracted from the converted data matrix with their cumulative contribution of 81.929%, including 47.727% of yield factor, 15.285% of oil factor yield, 10.206% of saponin factor and 8.711% of form factor. The superior individuals of six types are selected, in top 5%, including high-yield type, high-oil type, high-saponin type, high-yield and high-oil type, high-yield and high-saponin type, and compound type. The realized gains of superior individuals in each type are in 165%-311%. The principal component analysis can be used to select superior individual based on the main fruit traits and seed traits ofS.mukorossi. The superior individuals can be applied in production.