Study on Phenotypic Variations of Elite Trees in Natural Populations of Hippophae rhamnoides L. subsp. sinensis Rousi in Shanxi

2019年9月至 11月，完成中国沙棘选优、采样工作，共计11个种群，110个优良单株。每个种群内取样植株间距在50 m以上。各取样地点的基本情况见表1。

表型性状指标主要包括果实横径、果实纵径、果形指数、果柄长、百果质量、种子千粒质量、结实密度、当年枝长、刺长、刺密度。

将同一种群采集的果实充分混匀，采用四分法，50个果实1组，随机抽取3组，共150个果实；利用数显游标卡尺（精度0.01 mm）分别测量各组果实纵径、横径、果柄长。根据果实纵径和横径计算果形指数（果形指数 = 果实纵径/果实横径）。从每个种群中随机抽取500粒果实，分成5组，每组100粒，利用电子天平（精度0.001 g）测定每组果实的质量，取其平均值即为百果质量。将同一种群采集的果实中剥取种子，充分混匀，采用四分法，100粒1组，随机抽取3组，共300粒种子，用电子天平（精度0.001 g）测定每组种子的质量，取其平均值乘以10，即为种子千粒质量。从每个种群内单株的东、西、南、北四个方向分别随机抽取10 cm长的结果枝，测定结果和棘刺数量，选3株优良单株进行重复，计算其平均值即为结实密度和刺密度。从每个种群内单株的东、西、南、北四个方向分别随机抽取、测定当年生枝长，选3株优良单株进行重复，计算其平均值即为当年枝长。从每个种群内单株上随机选取10个棘刺，测量其长度，重复3次，计算刺长。

采用SPSS 22.0对数据进行分析。通过单因素方差分析、Duncan新复极差法进行多重比较。计算表型性状的均值、标准差、变异系数（CV），其中，变异系数（CV）= 标准差/均值。

采用相关分析法计算各表型性状间及其与生态因子的Pearson相关系数。对各种群优树的表型性状指标进行主成分分析，鉴于刺长和刺密度越小越好，将刺长、刺密度通过变换倒数方式进行正向转换。采用方差最大化正交旋转方法进行主成分分析^[28]。

中国沙棘天然种群优树表型性状利用Y值进行综合评价，其计算公式如下：

　　式中：λ₁、λ₂、λ_n分别代表第1主成分、第2主成分、第n主成分的特征值；Y₁、Y₂、Y_n分别代表第1主成分、第2主成分、第n主成分。

采用离差平方和法，根据欧氏距离对各种群优树的主成分综合得分值进行聚类分析，以掌握山西省境内中国沙棘天然种群优树表型性状的地理分化特点。

中国沙棘天然种群优树表型性状的方差分析结果表明：种子千粒质量、果柄长在不同种群间的差异呈显著水平（p < 0.05），其余8个表型性状指标均呈极显著水平（p < 0.01）。可见，中国沙棘优树表型性状在不同群体间存在差异。多重比较结果见表2。

由表2 可知：隰县种群优树的果实纵径、果形指数、果柄长、结实密度、刺长均最大，百果质量、种子千粒质量较大；广灵种群优树的果实横径、百果质量最大，果实纵径、种子千粒质量、结实密度均较大，果柄最短；兴县种群优树的种子千粒质量最大，果实横径、果实纵径、百粒质量也较大；文水种群优树的当年枝长最长，刺密度最大，果形指数最小；右玉种群优树的果实横径、果实纵径最小，百果质量、种子千粒质量也较小；岢岚种群优树的百果质量、当年枝长、刺长、刺密度均最小；五寨种群优树的结实密度最小；云州种群优树的种子千粒质量最小。

中国沙棘天然种群优树表型性状的变异情况见表3。研究表明，中国沙棘种群优树的表型性状变异大小因性状而异，10个表型性状的平均变异系数为22.10%，其变化幅度为5.23%～37.42%，以刺密度的变异系数最大，果形指数的变异系数最小。营养生长性状的变异大于生殖生长性状的变异。

由表3还得知：右玉种群优树的果实横径、果实纵径、百果质量、种子千粒质量变异系数均最大，反映出右玉群体优树的果实性状变异较明显，遗传稳定性较低，果实多样性最丰富；广灵种群优树除种子千粒质量变异系数较小外，其余表型性状的变异系数均较大，反映出该群体表型性状变异分化较明显，遗传稳定性较低。隰县种群优树的10个表型性状的变异系数均较小，反映出该种群优树表型性状遗传稳定性较高，多样性相对较单一。

中国沙棘优树表型性状间的相关分析结果（表4）表明：中国沙棘优树的10个表型性状指标间存在不同程度的相关性，其中，果实横径与果实纵径、百果质量、种子千粒质量、当年枝长均呈极显著（p < 0.01）正相关，与果形指数、结实密度呈极显著（p < 0.01）负相关；果形指数与果实纵径、果柄长、结实密度呈极显著（p < 0.01）正相关，与果实横径呈极显著（p < 0.01）负相关；种子千粒质量与果实横径、果实纵径、百果质量呈极显著（p < 0.01）正相关；当年枝长与刺长、刺密度呈极显著（p < 0.01）正相关；刺密度与果实纵径、百果质量、当年枝长、刺长呈极显著（p < 0.01）正相关。

中国沙棘优树表型性状与生态因子间的相关分析结果（表5）表明：中国沙棘优树的表型性状与生态因子间存在不同程度的相关性，其中，果实横径、果实纵径、百果质量、种子千粒质量、当年枝长与海拔因子呈极显著（p < 0.01）正相关；果实纵径、果形指数、果柄长、当年枝长与纬度因子呈极显著（p < 0.01）负相关；果形指数与年降水量呈极显著（p < 0.01）正相关。

对中国沙棘天然种群优树的表型性状进行主成分分析，结果（表6），得到3个特征根大于1的主成分，第1主成分中百果质量、果实纵径、当年枝长载荷度较高，主要表征的是果实大小性状和当年生长量，贡献率为46.15%；第2主成分中果形指数、结实密度载荷度较高，主要表征的是果实的长宽比和结实性状，贡献率为26.13%；第3主成分中刺密度、种子千粒质量、刺长载荷度较高，主要表征的是棘刺的性状和种子大小性状，贡献率为15.38%。上述3个主成分的累积方差贡献率为87.66%，已能反应中国沙棘优树表型性状10个测定指标的绝大多数信息，故选取前3个主成分为中国沙棘优树表型性状的重要主成分。

Y₁、Y₂、Y₃主成分的函数式是根据样本成分相关矩阵的因子得分矩阵得出：

Y₁=0.349X₁+0.112X₂−0.174X₃+0.264X₄+0.349X₅−0.076X₆+0.146X₇ +0.046X₈+0.092X₉+0.120X₁₀

Y₂=−0.196X₁+0.212X₂+0.366X₃+0.001X₄−0.028X₅+0.304X₆+0.237X₇ +0.065X₈−0.020X₉+0.049X₁₀

Y₃=−0.059X₁−0.044X₂−0.094X₃−0.019X₄+0.150X₅−0.020X₆+0.228X₇−0.240X₈+0.385X₉+0.415X₁₀

式中：X₁代表果实横径；X₂代表果实纵径；X₃代表果形指数；X₄代表百果质量；X₅代表种子千粒质量；X₆代表果柄长；X₇代表结实密度；X₈代表当年枝长；X₉代表刺长的倒数；X₁₀代表刺密度的倒数。

通过公式（1），计算得出中国沙棘天然种群优树表型性状的综合评价结果见图1。结果表明：中国沙棘天然种群优树表型性状的综合得分从高到低依次为隰县、岢岚、广灵、原平、兴县、文水、云州、宁武、右玉、和顺、五寨，说明隰县种群优树的表型性状表现最佳，五寨种群优树的综合表型性状最差。为掌握山西省境内中国沙棘优树的表型性状地理分化特点，进一步将主成分综合得分进行聚类分析，结果见图2。

Figure 1.  Comprehensive principal component scores of elite trees in natural populations of H. rhamnoides subsp. sinensis

Figure 2.  Cluster analysis of phenotypic traits of elite trees in natural populations of H. rhamnoides subsp. sinensis

聚类分析结果表明，在距离系数为5时，可将11个中国沙棘天然种群优树分为4类，第1类为隰县种群优树，综合性状表现最好；第2类分别为广灵、原平、岢岚种群优树，综合表现较好；第3类分别为云州、文水、宁武、兴县种群优树，综合表现一般；第4类分别为五寨、右玉、和顺种群优树，综合表现较差。

丰富的优良种质资源是林木育种研究重要的物质基础，其中，表型性状分析对揭示基因型差异以及鉴定种群生态学分析具有重要意义^[23]。生物的变异和多样性是开展育种工作的前提，变异系数可以反映性状值的离散特征，变异系数越大，性状值的离散程度越大，表型多样性越高；相反，变异系数越小，说明该种群的性状离散程度越低，物种多样性越低^[12]。有研究表明，植物群体间的变异程度能反映出植物本身对环境的适应程度，变异系数越大，说明植物适应环境的范围越广^[29]。

本研究发现，山西省中国沙棘天然种群优树间存在丰富的表型变异，平均变异系数为22.10%，其中，刺密度的变异系数最大（37.42%），果形指数的变异系数最小（5.23%），营养生长性状的变异大于生殖生长性状的变异。黄铨等^[24]研究发现，山西省关帝山的中国沙棘每个性状都有很大幅度的变异，其中，营养器官的变异多大于生殖器官的变异，而在生殖器官中，果柄长度和鲜果百粒质量的变异又大于其它性状。吴琼等^[25]研究发现，山西省中国沙棘天然种群各表型性状的平均变异系数介于3.1%～26.5%，其中，营养性状变异幅度最大，种子性状最小。本研究结果与其一致，可能是因为营养器官主要通过光合作用为植物提供养分，容易受到外界环境因素的影响，变异性较强；而果实属于生殖器官，主要受自身遗传因素控制，性状相对稳定，这也保证了物种遗传上的相似性和稳定性^[6]。

另外，本研究还发现，中国沙棘优树的同一性状在不同种群的变异幅度存在差异性。余天蓝等^[23]在对新疆6个居群野生沙棘果实性状多样性研究也得出类似结论。这在一定程度上也表明，生长环境的异质性可能是导致中国沙棘种群优树表型变异的主要原因。

不同植物对环境因子的适宜机制和敏感程度的差异导致其呈现不同的地理变异规律^[8]。中国沙棘优树表型性状相关分析表明，种子千粒质量与果实横径、果实纵径、百果质量呈极显著（p < 0.01）正相关，表明沙棘种子大小与果实大小存在关系，果实的纵、横径越大，种子也越大，该研究结果与吴琼^[30]对中国沙棘表型性状相关性分析结论一致。中国沙棘优树表型性状与生态因子间的相关分析显示，果形指数与纬度因子呈极显著负相关，表明随纬度降低，果实纵径越大，果形指数增大；果实横径、果实纵径、百果质量、种子千粒质量与海拔因子呈极显著正相关，表明随海拔升高，沙棘果实、种子增大，该研究结果也与吴琼^[30]对中国沙棘表型性状相关性分析结论一致。通常情况下，植物表型性状的地理变异比较复杂，而且因物种而异，如长柄扁桃种群表型变异主要受年平均气温、纬度和无霜期3个主要地理气候因子影响^[11]，无患子表型性状的变异与地理纬度和年平均气温2个主要因子密切相关^[28]，川西云杉种实性状存在经度和纬度并存的变异模式^[29]。本研究结果表明，中国沙棘天然种群优树的表型性状变异与地理纬度和海拔2个主要因子密切相关。

由于分布和环境条件等因子的综合作用，植物会形成连续变异、不连续变异以及随机变异等多种地理变异模式^[31]。主成分分析和聚类分析结果发现，遗传距离大致在5时，分界线将山西省11个中国沙棘种群优树分为4类，各种群并未严格按照地理距离而聚类，说明山西省中国沙棘天然种群优树表型性状未形成连续变异，具有随机变异的特点。吴琼等^[25]研究也发现，山西省境内中国沙棘没有表现出明显的地理梯度变异式样，表型性状分化有多态性的特点。可能是由于中国沙棘主要分布在深山沟壑中，种群间大多受吕梁山、太行山等高大山脉的隔离，种子和花粉难以在种群间传播，即使地理距离较近的种群间难以进行基因交流，增加了各种群独立分化的机会。

沙棘是集水土保持、医药保健、工业日化等为一体的多功能生态树种，在我国生态建设和乡村振兴中发挥着巨大的作用。中国沙棘在山西省境内分布跨度大，生长环境条件复杂，通过长期的地理隔离、自然选择，使得山西省中国沙棘种群优树间存在丰富的变异，为中国沙棘遗传改良奠定了重要物质基础。为了充分保护和利用现有中国沙棘天然资源，有效地进行遗传改良，培育优质、高产及功能性新品种，建议如下：（1）中国沙棘表型变异主要与纬度和海拔相关，尽可能丰富这些因子的采样范围，进行随机采样，同时增加采样数和叶片、种子的评测指标，使研究更系统、全面，才能更有效地开展中国沙棘的保护和利用工作；（2）开展生态遗传学研究，将形态学、分子标记和地理生态结合分析种群遗传多样性，掌握各种群的亲缘关系；（3）加强对现有中国沙棘天然种群和优良单株的保护力度，通过收集种质资源，建立基因库，有效避免采种对现有遗传资源的破坏。

山西省中国沙棘天然种群优树间存在丰富的表型变异，平均变异系数为22.10%，营养生长性状的变异大于生殖生长性状的变异，其中，隰县中国沙棘种群优树的果实纵径、结实密度最大，表型性状的变异系数较小，综合性状表现最好，可参考作为中国沙棘优良种质资源收集的重点区域。中国沙棘天然种群优树的表型性状变异与地理纬度和海拔2个主要因子密切相关。山西省中国沙棘天然种群优树表型性状未形成连续变异，具有随机变异的特点。以上研究为中国沙棘优良种质资源收集和品种创制奠定了重要物质基础，也为其栽培区划提供了技术参考。