根据西藏野生光核桃的自然集中分布区，兼顾其生长的海拔、气候和成熟期的差异性，对西藏光核桃的5个地级市18个种群进行实地调查采样，调查采样点覆盖了光核桃的主要分布区。通过实地调查，采用GPS定位记录采样地的经纬度和海拔，根据光核桃褪去绿色，着黄、粉红色和形态判断成熟，并根据开始成熟时间记录各种群光核桃的成熟期。从每个种群中选取13株长势良好的光核桃，各单株间距大于50 m，种群间距离大于3 km^[16]。分别从每个植株的阴面上、阴面下、阳面上和阳面下4方位采集成熟果实，从每个方位采集的成熟果实中随机选择5枚，共20枚，测量其表型性状。18个种群的位置及概况信息见图1和表1，其中年降雨量和年均温数据来自中国气象网及当地的气象局。

图  1  光核桃18个种群的分布图

Figure 1.  Distribution of the 18 Amygdalus mira populations

光核桃果实表型性状中果实质量（FWT）和果核质量（NWT）用电子天平（精度0.01 g）测得，果实的纵径（FL）、横径（FWH）、侧径（FT）及果核的长（NL）、宽（NWH）、厚（NT）用电子游标卡尺（精度0.01 mm）测得。纵径为果实基部到顶部的距离；横径为果实腹缝线到果实边缘最大宽度；侧径为垂直于纵径、横径处的最厚距离，果核的长、宽、厚取值方式同果实纵径、横径和侧径。果实的果肉厚（PT）= FT－NT，出核率（NR）= NWT/FWT × 100%。

通过巢式方差分析、方差分量及种群间表型分化系数（V_ST）的计算^[17]，研究18个种群光核桃果实表型性状种群间与种群内变异；使用均值、变异系数（coefficient of variation，CV）分析光核桃种群间与种群内表型性状的离散程度。

式中：s为种群数目；t为种群内的株数；$ {\sigma }_{\frac{t}{s}}^{2} $为种群间方差分量；$ {\sigma }_{s}^{2} $为种群内方差分量。

式中：S为标准差，X为平均值。

采用Pearson相关系数分析果实性状间的相关性。采用主成分分析（principal component analysis，PCA）果实表型性状与种群的相对关系，根据排序图筛选优良性状的主要代表种群。冗余分析（redundancy analysis，RDA）生态因子对果实表型性状总体变异的影响，并采用多元回归分析生态因子对特征性状的影响规律。PCA、RDA使用R语言vegan软件包分析，其它使用SAS及EXCEL软件分析。

实地调查18个种群光核桃的成熟期，发现种群间有明显差异（表1）。光核桃的成熟期主要从8月中旬到9月下旬，林芝市朗县朗镇（海拔3 118 m）光核桃成熟最早，日喀则亚东县上亚东乡（海拔3 218 m）成熟较晚。即使在同一地区不同的种群间也存在差异，如山南市贡嘎县江塘镇（海拔3 770 m）光核桃的成熟期比山南市其他种群的成熟期晚；昌都市芒康县木许乡（海拔2 280 m）和八宿县林卡乡（海拔2 984 m）与山南市加查县的安绕镇（海拔3 265 m）和冷达乡（海拔3 168 m）属于不同地区，海拔相差也较大，但却有相近的成熟期。

对18个种群的光核桃果实表型性状进行了巢式方差分析、方差分量及V_ST的计算（表2、3），结果显示10个果实表型性状在种群间和种群内均是极显著差异（P < 0.01）；种群间的表型变异占总体变异的60.98%，种群内的占28.09%；表型分化系数为40.48% ~ 87.46%，果核厚表型分化系数最大，果核质量表型分化系数最小。表型分化系数的平均值为67.99%，表明光核桃果实表型变异主要来源于种群间。

对18个种群光核桃的果实表型性状进行相关性分析（表4），结果显示果实质量、纵径、横径、侧径、果肉厚相互之间存在显著性正相关；果实质量、纵径、横径、侧径、果肉厚与出核率呈显著负相关，果实的大小与果核的大小相关性不显著（r = − 0.13 ~ 0.20），表明果核的大小受果实大小的影响较小。

变异系数（CV）反映了果实表型性状种群内与种群间的离散程度^[18]，结果种群内各果实表型性状变异系数的平均值DB最大，为18.53%；BM最小，为8.43%。果实表型性状种群间平均变异系数为7.86% ~ 25.25%，果实质量最大，果实纵径最小（表5）。

PCA分析通过降维形式直观的表现光核桃果实表型性状与种群的相对关系^[19]，结果显示PC₁、PC₂共解释了90.81%观测到的果实表型性状变异，PC₁和PC₂分别解释58.54%和32.27%（表6）。PC₁主要解释果实质量、纵径、横径、侧径、果肉厚和出核率；PC₂主要解释果核长、果核宽和果核厚。根据主成分PC₁和PC₂绘制线性排序图，结合利用部位果肉和果核，可以看出果实大、果核小、果肉厚、出核率低的代表种群为JT、QN、LZ，其表型性状果实质量24.43 ~ 28.48 g，纵径33.14 ~ 36.56 mm，横径34.55 ~ 36.76 mm，侧径33.34 ~ 36.17 mm，果核质量1.77 ~ 1.88 g，果核长14.94 ~ 20.93 mm，果核宽11.19 ~ 17.05 mm，果核厚5.24 ~ 11.07 mm，果肉厚25.10 ~ 28.10 mm，出核率6.87% ~ 7.56%；果实小、果核大、果肉薄、高出核率代表种群有GY、SYD、QD，其表型性状果实质量8.37 ~ 9.30 g，纵径27.39 ~ 28.24 mm，横径23.25 ~ 25.54 mm，侧径22.84 ~ 23.99 mm，果核质量1.31 ~ 1.57 g，果核长19.99 ~ 22.11 mm，果核宽15.68 ~ 15.93 mm，果核厚11.10 ~ 11.43 mm，果肉厚11.41 ~ 12.89 mm，出核率15.60% ~ 19.33%（图2，表5）。

图  2  18个种群光核桃果实表型性状PCA分析

Figure 2.  Segregation of the 18 populations according to fruit phenotypic traits determined by PCA analysis

RDA分析结果显示，5个生态因子共解释56.00%的果实表型性状变异，RDA₁、RDA₂分别解释44.82%和9.20%，共54.02%（表7，图3)。生态因子与主成分的相关性分析结果显示只有年均温对10个表型性状的总体变异有显著影响，年均温对果实性状的影响大小顺序为NR > NL > NWT > NWH > NT > FL > FWT > FT > FWH > PT（表8，图3）。

图  3  5个生态因子与18个种群光核桃果实表型性状RDA分析

Figure 3.  RDA analysis of 18 populations phenotypic traits with 5 ecological factors

根据光核桃果实质量、纵径、横径、侧径与果肉厚存在显著性正相关，果实的大小与果核的大小相关性不显著，并结合利用部位果肉和果核，筛选出果肉厚和出核率两个特征性状。果肉厚和出核率与生态因子进行多元回归分析（表9），结果显示两方程调整后判定系数分别为0.599 3，0.798 1，F检验P值均小于0.01，表明两个方程均回归显著。果肉厚与经度、海拔、年均降雨量呈负相关，与纬度和年均温呈正相关；出核率与经度、海拔、年均降雨量和年均温呈正相关，与纬度呈负相关。两个回归方程中自变量的决定系数通过t检验显示5个生态因子中经度、纬度、海拔对果肉厚和出核率的影响极显著。

不同种群光核桃的成熟期存在明显差异，主要从8月中旬到9月下旬。植物果实成熟期受环境条件和遗传特性的影响，经过长期的自然适应每个种群形成相对稳定的成熟期。Dicenta等人研究表明成熟期是可以遗传的^[20]。因此，对于成熟期的调查为早熟或晚熟品种选育提供数据支持，成熟期的差异有助于满足不同时期桃子的市场需求。

光核桃表型性状在种群间和种群内存在显著差异，表明光核桃表型多样性丰富。表型分化系数平均为67.99%，种群间的变异大于种群内的变异，光核桃的表型变异主要来源于种群间，这与包文泉在林芝、山南范围内分析光核桃表型变异的研究结论相一致^[14]。光核桃在西藏分布广泛，生境多样，各种群间有山地阻隔或距离较远，限制了花粉的传播和基因交流，也反映了地理和生殖隔离对表型变异的影响^[17,21]；另一方面，由于种群生境的差异造成了局部适应，致使不同生境下种群间有效基因流的降低，从而形成了适应性隔离^[22-24]。种群变异包括表型变异与遗传变异，表型变异又与遗传变异密切相关，植物表型变异大表明其具有遗传多样性潜力^[25]。因此，光核桃果实表型变异大，可以选育具有特征性状的光核桃。

表型变异系数越小离散程度越小，性状稳定性越好^[26]，种群内各果实表型性状变异系数的平均值DB最大为18.53%，BM最小为8.43%。种群内部环境差异大造成表型变异大，DB处在海拔变化较大的山谷，光照、水分等差异明显。种群间10个果实表型性状变异系数平均值为7.86% ~ 25.25%，果实质量最大，果实纵径最小，说明果实纵径较为稳定。

结合光核桃利用部位果肉和果核，筛选出果实大、果核小、果肉厚、出核率低的果用型和果实小、果核大、果肉薄、高出核率的核用型两种光核桃种质资源。果用型光核桃可用于加工果酒、果醋、果脯等，或作为水果食用^[1,9]，代表性种群有JT、QN、LZ；核用型光核桃桃核可做活性炭、生物吸附剂外；桃仁有炒制食用和浸泡榨油食用的历史，其种子可加工或做坚果^[12]，代表种群有GY、SYD、QD。

表型差异是通过基因与环境的相互作用实现的，生态因子对植物表型的塑造起着重要的作用^[27]。结合RDA分析，经度、纬度、海拔、年均降雨量、年均温5个生态因子解释56.00%的表型性状变异，其他44.00%的表型变异可能来自于生物遗传特性^[28]、土壤养分^[29]、光^[30]、昼夜温差^[31]等条件。5个因子中只有年均温对光核桃的10个表型性状影响显著，这与李洪果研究的温度因子在气候因子中对种群表型性状的影响起主导作用的结论相一致^[32]，也同时说明了局部环境的温度受地形、海拔、降雨和光照等因素的影响，是环境变化的集中表现。

特征性状果肉厚与经度、海拔、年均降雨量呈负相关，可能与西藏降雨规律有关。光核桃分布区随着经度的增加，降雨量有所增加，且降雨集中在5—10月^[33]，通过影响温度、光照等进一步影响果实有机物质的积累与生长^[34-35]。此外，随着海拔的增加，温度降低，植物生长缓慢，从而影响果实的生长。果肉厚与纬度和年均温呈正相关，随着纬度的升高，降雨量相对较少^[33]，白天阳光充足，夜晚温度低，有利于有机物质的积累^[34]；年均温的增加也有利于果实有机物质的积累^[35]。出核率与经度、海拔、年均降雨量和年均温呈正相关，与纬度呈负相关。随着经度、海拔、年均降雨量的增加，果实果肉厚的减小，出核率增加；同时年均温的增加，果核也增大，与出核率呈现出正相关；随着纬度的升高，果肉厚增加，出核率降低。生态因子对特征性状的影响规律对光核桃优良种质资源推广种植具有指导意义。

根据西藏野生光核桃的集中分布区，兼顾其生长的海拔、气候和成熟期的差异性，对5个地级市18个种群进行实地调查采样，研究光核桃果实表型性状种群变异规律，根据利用部位对光核桃种质资源进行筛选，得出如下结论：光核桃果实种群间表型性状多样，存在丰富的遗传变异潜力。依据光核桃果实利用部位及种群间果实表型性状筛选出果用型和核用型两类种质资源，山南贡嘎县江塘镇、林芝米林县羌纳乡、林芝朗县朗镇为果用型光核桃的代表种群，林芝察隅县古玉乡、日喀则亚东县上亚东乡、昌都芒康县曲登乡为核用型光核桃的代表种群。经度、纬度、海拔、年均降雨量、年均温5个生态因子中年均温对10个果实表型性状的总体影响显著，其中经度、纬度、海拔对特征性状果肉厚和出核率的影响极显著。研究光核桃果实表型种群变异，筛选优良种质资源，为光核桃果实的开发利用提供思路，也为光核桃的良种选育和推广提供参考。