供试的麻栎种子于2019年9月至10月分别采自麻栎自然分布区(25°21′～39°49′ N，100°54′～121°59′ E)的7个天然群体。麻栎林中选择优势树种作为母树，采种母树之间相距50 m以上，每一群体选取20～30株母树采集种子作为1份，每份种子不少于300粒，将采集的所有种子进行充分混合，然后将其用作群体的种子样品，为增加实验的准确性，将每个地理种源都使用了GPS定位软件进行定位。种子来源情况详见表1。

对不同群体的种子长、种子宽、种子长宽比、表面积和体积等利用扫描仪Epson/V700 PHOTO测定表型性状，采用winseedle Pro 2013b软件可获得表型性状指标，使用1/100电子天平测定千粒质量，从每组样品分别随机选取50粒种子称质量，重复3次，取其平均值。

采用巢式线性模型对麻栎群体间和群体内的表型变异进行方差分析。用变异系数（CV）表示表型性状的离散程度。计算表型分化系数（Vst）来表示麻栎群体的表型分化程度，即群体间方差分量与方差总量的比值。计算公式分别为：CV = (δ/${{\overline x}}$) × 100%，其中，δ为标准差、${{\overline x}}$为群体平均值，Vst = [δ²_t/s /(δ²_t/s + δ²_s)] × 100%，其中，δ²_t/s为群体间的方差分量，δ²_s为群体内的方差分量^[11]。采用Excel 2010软件对数据进行统计，用SPSS 21软件对各项指标进行分析，其中，采用Pearson相关分析双变量相关性，并根据各性状在群体间的欧氏距离进行聚类分析。

由表2的方差分析结果可知：在群体间和群体内，种子长、宽、长宽比、表面积、体积及千粒质量均差异极显著（P < 0.01），说明麻栎种子表型性状在群体间存在一定程度的变异。

根据表3可以看出：各群体种子长度的平均值为1.63～2.32 cm，其中，最大的是云南祥云县群体，最小的是贵州贵定县群体；各群体种子宽度的平均值为1.33～1.84 cm，贵州贵定县群体最小，云南祥云县群体最大；长宽比最大的是山西中条山群体，最小的是河北联峰山群体；种子体积的平均值为1.48～4.15 cm³ ，其中，最小的是贵州贵定县群体，最大的是云南祥云县群体；种子表面积平均值为5.19～10.41 cm²，最小的是贵州贵定县群体，最大的是云南祥云县群体；种子千粒质量差异最大，平均变幅为1.74～4.82 kg，贵州贵定县群体最小，云南祥云县群体最大。

分析和比较不同表型性状时，性状的离散程度用变异系数来表示^[12]，变异系数越大，性状的变异程度越大；变异系数越小，性状的变异程度越小。由表4可知：在麻栎的6个表型性状中，平均变异系数最大的为种子体积，达28.18%，其次是种子表面积(18.79%)，变异系数最小的为种子长(10.34%)，表明种子长的变异程度最小，而种子体积这一性状的变异程度最大。

由表5可以看出：在表型性状的总方差分量占比中，群体间方差分量占比51.29%，群体内方差分量占比13.64%，机误占比35.07%。从群体间表型性状的变异幅度看，群体间种子长宽比的方差分量百分比最小，为41.67%；而千粒质量的方差分量百分比最大，为72.20%。从群体内表型性状的变异幅度看，种子宽的方差分量百分比最大，其次为种子长，最小的是长宽比。种子表型分化系数(Vst)为73.17%～84.34%，其中，千粒质量的Vst最大，种子宽的Vst最小；表型分化系数的平均值为78.84%，表明麻栎群体78.84%的变异来自群体间。

根据表6可知：各种子表型性状中两两性状具有极显著的相关关系，其中，千粒质量与长宽比的相关系数较小，为0.289；而种子宽和种子体积的相关系数达0.907。说明选择的6个表型性状进行麻栎种子表型性状分析是可行的。

利用欧氏平方距离系数对麻栎群体进行聚类，图1表明：以欧氏距离10为阈值，7个麻栎群体可以划分为2大类，第1类包括山西中条山、辽宁普兰店区、河南桐柏县、云南祥云县群体，第2类包括贵州贵定县、山东黄岛区、河北联峰山群体。

图  1  麻栎群体表型性状聚类树系图

Figure 1.  Cluster tree of phenotypic traits of Quercus acutissima population

表型变异既是生物遗传变异最直接的表现形式，也是遗传多样性研究的重要内容之一^[13-14]。通常分布范围广泛的树种都具有丰富的表型变异，本研究以麻栎群体种子为材料进行表型性状变异的分析，可以在一定程度上反映其群体遗传多样性水平^[15-16]。

由巢式方差分析结果可知：7个麻栎群体种子表型性状在群体间和群体内均差异极显著，表明麻栎种子遗传变异较丰富，不同群体的地理环境不同使得麻栎种子性状变异差异明显，为下一步进行优良种源和家系的选择奠定了基础。麻栎群体间的变异大于群体内，这与许洋等^[17]、刘志龙等^[18]及张元燕等^[19]的研究结论相反，可能是与遗传因素和环境因素的共同作用有关。此外，麻栎结实有明显的“大小年”现象，而且不同群体间种子的成熟期不一致，可能会造成个体的差异，所以，群体间的变异程度也可能与个体的差异有关。各个不同性状中，种子体积的变异系数最大，为28.18%，种子长的变异系数最小，为10.34%，说明种子长这一性状较稳定，这与常恩福等^[20]、梁德洋等^[21]及厉月桥等^[22]对锥连栎和蒙古栎的研究基本上一致。7个群体中，平均变异系数最大的是青岛黄岛区群体，为20.58%，平均变异系数最小的是贵州贵定县群体，为13.17%，说明山东黄岛区麻栎群体的生境条件相对复杂。麻栎群体平均表型分化系数为78.84%，表明麻栎群体表型变异在群体间的贡献率为78.84%，变异来源以群体间为主，群体间差异反映了地理隔离和生殖隔离的差异，群体间多样性是物种内多样性的重要组成部分^[23]。因此，在麻栎种质资源保护工作中应当注意对群体间的变异进行合理利用。

相关性分析结果表明：麻栎种子长和种子宽、千粒质量均呈极显著正相关关系，这与厉月桥等^[24]对不同种源蒙古栎种子表型性状的研究结果相近，说明种子长增加，其种子宽和千粒质量也极显著的增加。根据聚类分析结果，将7个麻栎群体划分成2个组群，发现麻栎群体并没有完全按照地理距离远近而聚类，与地域没有明显相关性，这可能与不同的地理环境有关，此结果与栎属的其他树种白栎^[25]、铁橡栎^[26]和栓皮栎^[27]的种子表型性状的研究基本吻合。千粒质量可以作为检测种子品质的一个指标，千粒质量越大，种子品质越好，将麻栎群体聚成了2类，第1类代表千粒质量大的，包括山西中条山、辽宁普兰店区、河南桐柏县、云南祥云县群体；第2类代表千粒质量较小的，包括贵州贵定县、山东黄岛区、河北联峰山群体。不同种源麻栎种子千粒质量具有一定的地理变异模式，可能与分布区域的环境因素有关，如光照、温度、经纬度、降雨及海拔等因素相互作用造成了不同群体麻栎种子千粒质量的不同变异模式^[28]。

表型变异既受遗传控制，也受环境因素的影响^[29-31]，表型多样性是遗传和环境影响的结果。然而，表型变异中有多少遗传成分可以被利用，需要进一步深入研究，下一步有必要选择合适的分子标记开展麻栎分子水平的变异分析，为有效保护麻栎种质资源提供可靠的理论依据，提高麻栎的利用价值。

麻栎天然群体种子在群体间和群体内都达到了极显著的水平，其遗传变异较丰富，是遗传因素和环境因素共同作用的结果，且群体间的变异大于群体内。研究结果表明，各表型性状间呈极显著相关关系。聚类分析将7个群体分为2个组群，与各个群体的地理位置分布无明显相关性。麻栎天然群体表型变异的研究，为麻栎多样性种质资源保护和利用提供了理论依据。